ViewVC Help

View File | Revision Log | Show Annotations | Download File | Root Listing

root/Development/ray/src/rt/normal.c

Comparing ray/src/rt/normal.c (file contents):
Revision 2.37 by gwlarson, Wed Dec 16 18:14:58 1998 UTC vs.
Revision 2.41 by greg, Wed Mar 12 04:59:05 2003 UTC

#	Line 1 | Line 1
1	–	/* Copyright (c) 1998 Silicon Graphics, Inc. */
2	–
1		#ifndef lint
2	<	static char SCCSid[] = "$SunId$ SGI";
2	>	static const char RCSid[] = "$Id$";
3		#endif
6	–
4		/*
5		*  normal.c - shading function for normal materials.
6		*
#	Line 14 | Line 11 | static char SCCSid[] = "$SunId$ SGI";
11		*     Later changes described in delta comments.
12		*/
13
14	+	#include "copyright.h"
15	+
16		#include  "ray.h"
17
18		#include  "otypes.h"
19
20		#include  "random.h"
21
23	–	extern double  specthresh;              /* specular sampling threshold */
24	–	extern double  specjitter;              /* specular sampling jitter */
25	–
26	–	extern int  backvis;                    /* back faces visible? */
27	–
22		#ifndef  MAXITER
23		#define  MAXITER        10              /* maximum # specular ray attempts */
24		#endif
25	+	/* estimate of Fresnel function */
26	+	#define  FRESNE(ci)     (exp(-6.0*(ci)) - 0.00247875217)
27
28	<	static  gaussamp();
28	>	static void  gaussamp();
29
30		/*
31		*      This routine implements the isotropic Gaussian
#	Line 69 | Line 65 | typedef struct {
65		}  NORMDAT;             /* normal material data */
66
67
68	+	static void
69		dirnorm(cval, np, ldir, omega)          /* compute source contribution */
70		COLOR  cval;                    /* returned coefficient */
71		register NORMDAT  *np;          /* material data */
#	Line 76 | Line 73 | FVECT  ldir;                   /* light source direction */
73		double  omega;                  /* light source size */
74		{
75		double  ldot;
76	+	double  ldiff;
77		double  dtmp, d2;
78		FVECT  vtmp;
79		COLOR  ctmp;
#	Line 87 | Line 85 | double  omega;                 /* light source size */
85		if (ldot < 0.0 ? np->trans <= FTINY : np->trans >= 1.0-FTINY)
86		return;         /* wrong side */
87
88	<	if (ldot > FTINY && np->rdiff > FTINY) {
88	>	/* Fresnel estimate */
89	>	ldiff = np->rdiff;
90	>	if (np->specfl & SP_PURE && (np->rspec > FTINY & ldiff > FTINY))
91	>	ldiff *= 1. - FRESNE(fabs(ldot));
92	>
93	>	if (ldot > FTINY && ldiff > FTINY) {
94		/*
95		*  Compute and add diffuse reflected component to returned
96		*  color.  The diffuse reflected component will always be
97		*  modified by the color of the material.
98		*/
99		copycolor(ctmp, np->mcolor);
100	<	dtmp = ldot * omega * np->rdiff / PI;
100	>	dtmp = ldot * omega * ldiff / PI;
101		scalecolor(ctmp, dtmp);
102		addcolor(cval, ctmp);
103		}
#	Line 154 | Line 157 | double  omega;                 /* light source size */
157		}
158
159
160	+	int
161		m_normal(m, r)                  /* color a ray that hit something normal */
162		register OBJREC  *m;
163		register RAY  *r;
164		{
165		NORMDAT  nd;
166	+	double  fest;
167		double  transtest, transdist;
168		double  mirtest, mirdist;
169		int     hastexture;
#	Line 177 | Line 182 | register RAY  *r;
182		raytrans(r);
183		return(1);
184		}
185	+	raytexture(r, m->omod);
186		flipsurface(r);                 /* reorient if backvis */
187	<	}
187	>	} else
188	>	raytexture(r, m->omod);
189		nd.mp = m;
190		nd.rp = r;
191		/* get material color */
#	Line 190 | Line 197 | register RAY  *r;
197		nd.alpha2 = m->oargs.farg[4];
198		if ((nd.alpha2 *= nd.alpha2) <= FTINY)
199		nd.specfl |= SP_PURE;
200	<	if (r->ro != NULL && isflat(r->ro->otype))
201	<	nd.specfl |= SP_FLAT;
195	<	/* get modifiers */
196	<	raytexture(r, m->omod);
197	<	if (hastexture = DOT(r->pert,r->pert) > FTINY*FTINY)
200	>
201	>	if (hastexture = (DOT(r->pert,r->pert) > FTINY*FTINY)) {
202		nd.pdot = raynormal(nd.pnorm, r);       /* perturb normal */
203	<	else {
203	>	} else {
204		VCOPY(nd.pnorm, r->ron);
205		nd.pdot = r->rod;
206	+	if (r->ro != NULL && isflat(r->ro->otype))
207	+	nd.specfl |= SP_FLAT;
208		}
209		if (nd.pdot < .001)
210		nd.pdot = .001;                 /* non-zero for dirnorm() */
#	Line 206 | Line 212 | register RAY  *r;
212		mirtest = transtest = 0;
213		mirdist = transdist = r->rot;
214		nd.rspec = m->oargs.farg[3];
215	+	/* compute Fresnel approx. */
216	+	if (nd.specfl & SP_PURE && nd.rspec > FTINY) {
217	+	fest = FRESNE(r->rod);
218	+	nd.rspec += fest*(1. - nd.rspec);
219	+	} else
220	+	fest = 0.;
221		/* compute transmission */
222		if (m->otype == MAT_TRANS) {
223		nd.trans = m->oargs.farg[5]*(1.0 - nd.rspec);
#	Line 254 | Line 266 | register RAY  *r;
266		if (nd.rspec > FTINY) {
267		nd.specfl |= SP_REFL;
268		/* compute specular color */
269	<	if (m->otype == MAT_METAL)
269	>	if (m->otype != MAT_METAL) {
270	>	setcolor(nd.scolor, nd.rspec, nd.rspec, nd.rspec);
271	>	} else if (fest > FTINY) {
272	>	d = nd.rspec*(1. - fest);
273	>	for (i = 0; i < 3; i++)
274	>	nd.scolor[i] = fest + nd.mcolor[i]*d;
275	>	} else {
276		copycolor(nd.scolor, nd.mcolor);
277	<	else
278	<	setcolor(nd.scolor, 1.0, 1.0, 1.0);
261	<	scalecolor(nd.scolor, nd.rspec);
277	>	scalecolor(nd.scolor, nd.rspec);
278	>	}
279		/* check threshold */
280		if (!(nd.specfl & SP_PURE) && specthresh >= nd.rspec-FTINY)
281		nd.specfl |= SP_RBLT;
#	Line 333 | Line 350 | register RAY  *r;
350		}
351
352
353	<	static
353	>	static void
354		gaussamp(r, np)                 /* sample gaussian specular */
355		RAY  *r;
356		register NORMDAT  *np;

Diff Legend

–	Removed lines
+	Added lines
<	Changed lines (old)
>	Changed lines (new)