#ifndef lint
static const char	RCSid[] = "$Id: m_mirror.c,v 2.12 2009/06/06 02:11:43 greg Exp $";
#endif
/*
 * Routines for mirror material supporting virtual light sources
 */

#include "copyright.h"

#include  "ray.h"
#include  "otypes.h"
#include  "rtotypes.h"
#include  "source.h"

/*
 * The real arguments for MAT_MIRROR are simply:
 *
 *	3 rrefl grefl brefl
 *
 * Additionally, the user may specify a single string argument
 * which is interpreted as the name of the material to use
 * instead of the mirror if the ray being considered is not
 * part of the direct calculation.
 */

static int mir_proj(MAT4  pm, OBJREC  *o, SRCREC  *s, int  n);
static void mirrorproj(MAT4  m, FVECT  nv, double  offs);

VSMATERIAL  mirror_vs = {mir_proj, 1};


extern int
m_mirror(			/* shade mirrored ray */
	register OBJREC  *m,
	register RAY  *r
)
{
	COLOR  mcolor;
	RAY  nr;
	int  rpure = 1;
	register int  i;
					/* check arguments */
	if (m->oargs.nfargs != 3 || m->oargs.nsargs > 1)
		objerror(m, USER, "bad number of arguments");
					/* check for substitute material */
	if (m->oargs.nsargs > 0 &&
			(r->rsrc < 0 || source[r->rsrc].so != r->ro)) {
		if (!strcmp(m->oargs.sarg[0], VOIDID)) {
			raytrans(r);
			return(1);
		}
		return(rayshade(r, lastmod(objndx(m), m->oargs.sarg[0])));
	}
					/* check for bad source ray */
	if (r->rsrc >= 0 && source[r->rsrc].so != r->ro)
		return(1);

	if (r->rod < 0.)		/* back is black */
		return(1);
					/* get modifiers */
	raytexture(r, m->omod);
					/* assign material color */
	setcolor(mcolor, m->oargs.farg[0],
			m->oargs.farg[1],
			m->oargs.farg[2]);
	multcolor(mcolor, r->pcol);
					/* compute reflected ray */
	if (r->rsrc >= 0) {			/* relayed light source */
		rayorigin(&nr, REFLECTED, r, mcolor);
					/* ignore textures */
		for (i = 0; i < 3; i++)
			nr.rdir[i] = r->rdir[i] + 2.*r->rod*r->ron[i];
					/* source we're aiming for next */
		nr.rsrc = source[r->rsrc].sa.sv.sn;
	} else {				/* ordinary reflection */
		FVECT  pnorm;
		double  pdot;

		if (rayorigin(&nr, REFLECTED, r, mcolor) < 0)
			return(1);
		if (DOT(r->pert,r->pert) > FTINY*FTINY) {
			pdot = raynormal(pnorm, r);	/* use textures */
			for (i = 0; i < 3; i++)
				nr.rdir[i] = r->rdir[i] + 2.*pdot*pnorm[i];
			rpure = 0;
		}
						/* check for penetration */
		if (rpure || DOT(nr.rdir, r->ron) <= FTINY)
			for (i = 0; i < 3; i++)
				nr.rdir[i] = r->rdir[i] + 2.*r->rod*r->ron[i];
	}
	rayvalue(&nr);
	multcolor(nr.rcol, nr.rcoef);
	addcolor(r->rcol, nr.rcol);
	if (rpure && r->ro != NULL && isflat(r->ro->otype))
		r->rt = r->rot + nr.rt;
	return(1);
}


static int
mir_proj(		/* compute a mirror's projection */
	MAT4  pm,
	register OBJREC  *o,
	SRCREC  *s,
	int  n
)
{
	double	corr = 1.;
	FVECT  nv, sc;
	double  od, offs;
	int  i;
				/* get surface normal and offset */
	offs = od = getplaneq(nv, o);
	if (s->sflags & SDISTANT)
		offs = 0.;
				/* check for extreme point behind */
	if (s->sflags & SCIR) {
		if (s->sflags & (SFLAT|SDISTANT))
			corr = 1.12837917;	/* correct setflatss() */
		else
			corr = 1.0/0.7236;	/* correct sphsetsrc() */
	}
	VCOPY(sc, s->sloc);
	for (i = s->sflags & SFLAT ? SV : SW; i >= 0; i--)
		if (DOT(nv, s->ss[i]) > offs)
			VSUM(sc, sc, s->ss[i], corr);
		else
			VSUM(sc, sc, s->ss[i], -corr);
	if (DOT(sc, nv) <= offs+FTINY)
		return(0);
				/* everything OK -- compute projection */
	mirrorproj(pm, nv, od);
	return(1);
}


static void
mirrorproj(		/* get mirror projection for surface */
	register MAT4  m,
	FVECT  nv,
	double  offs
)
{
	register int  i, j;
					/* assign matrix */
	setident4(m);
	for (j = 0; j < 3; j++) {
		for (i = 0; i < 3; i++)
			m[i][j] -= 2.*nv[i]*nv[j];
		m[3][j] = 2.*offs*nv[j];
	}
}