ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/rt/m_brdf.c
Revision: 2.11
Committed: Wed Dec 21 09:51:46 1994 UTC (29 years, 4 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.10: +11 -3 lines
Log Message:
added -bv option for back face visibility (normally on)

File Contents

# User Rev Content
1 greg 1.7 /* Copyright (c) 1991 Regents of the University of California */
2 greg 1.1
3     #ifndef lint
4     static char SCCSid[] = "$SunId$ LBL";
5     #endif
6    
7     /*
8     * Shading for materials with arbitrary BRDF's
9     */
10    
11     #include "ray.h"
12    
13     #include "data.h"
14    
15     #include "otypes.h"
16    
17 greg 2.2 #include "func.h"
18    
19 greg 2.11 extern int backvis; /* back faces visible? */
20    
21 greg 1.1 /*
22     * Arguments to this material include the color and specularity.
23     * String arguments include the reflection function and files.
24     * The BRDF is currently used just for the specular component to light
25     * sources. Reflectance values or data coordinates are functions
26 greg 2.7 * of the direction to the light source. (Data modification functions
27     * are passed the source direction as args 2-4.)
28 greg 1.1 * We orient the surface towards the incoming ray, so a single
29     * surface can be used to represent an infinitely thin object.
30     *
31     * Arguments for MAT_PFUNC and MAT_MFUNC are:
32 greg 1.4 * 2+ func funcfile transform
33 greg 1.1 * 0
34 greg 1.4 * 4+ red grn blu specularity A5 ..
35 greg 1.1 *
36     * Arguments for MAT_PDATA and MAT_MDATA are:
37 greg 1.4 * 4+ func datafile funcfile v0 .. transform
38 greg 1.1 * 0
39 greg 1.4 * 4+ red grn blu specularity A5 ..
40 greg 1.5 *
41     * Arguments for MAT_TFUNC are:
42     * 2+ func funcfile transform
43     * 0
44     * 4+ red grn blu rspec trans tspec A7 ..
45     *
46     * Arguments for MAT_TDATA are:
47     * 4+ func datafile funcfile v0 .. transform
48     * 0
49     * 4+ red grn blu rspec trans tspec A7 ..
50     *
51     * Arguments for the more general MAT_BRTDF are:
52     * 10+ rrefl grefl brefl
53     * rtrns gtrns btrns
54     * rbrtd gbrtd bbrtd
55     * funcfile transform
56     * 0
57 greg 2.6 * 9+ rdf gdf bdf
58     * rdb gdb bdb
59     * rdt gdt bdt A10 ..
60 greg 1.5 *
61     * In addition to the normal variables available to functions,
62     * we define the following:
63     * NxP, NyP, NzP - perturbed surface normal
64     * RdotP - perturbed ray dot product
65 greg 2.6 * CrP, CgP, CbP - perturbed material color (or pattern)
66 greg 1.1 */
67    
68     typedef struct {
69     OBJREC *mp; /* material pointer */
70     RAY *pr; /* intersected ray */
71 greg 1.5 DATARRAY *dp; /* data array for PDATA, MDATA or TDATA */
72 greg 2.6 COLOR mcolor; /* material (or pattern) color */
73     COLOR rdiff; /* diffuse reflection */
74     COLOR tdiff; /* diffuse transmission */
75     double rspec; /* specular reflectance (1 for BRDTF) */
76     double trans; /* transmissivity (.5 for BRDTF) */
77     double tspec; /* specular transmittance (1 for BRDTF) */
78 greg 1.1 FVECT pnorm; /* perturbed surface normal */
79     double pdot; /* perturbed dot product */
80     } BRDFDAT; /* BRDF material data */
81    
82    
83     dirbrdf(cval, np, ldir, omega) /* compute source contribution */
84     COLOR cval; /* returned coefficient */
85     register BRDFDAT *np; /* material data */
86     FVECT ldir; /* light source direction */
87     double omega; /* light source size */
88     {
89     double ldot;
90     double dtmp;
91     COLOR ctmp;
92 greg 1.4 FVECT ldx;
93 greg 2.3 double lddx[3], pt[MAXDIM];
94 greg 2.7 double vldx[4];
95 greg 1.5 register char **sa;
96 greg 1.1 register int i;
97    
98     setcolor(cval, 0.0, 0.0, 0.0);
99    
100     ldot = DOT(np->pnorm, ldir);
101    
102 greg 1.5 if (ldot <= FTINY && ldot >= -FTINY)
103     return; /* too close to grazing */
104 greg 2.6
105 greg 1.5 if (ldot < 0.0 ? np->trans <= FTINY : np->trans >= 1.0-FTINY)
106 greg 1.1 return; /* wrong side */
107    
108 greg 2.6 if (ldot > 0.0) {
109 greg 1.1 /*
110     * Compute and add diffuse reflected component to returned
111     * color. The diffuse reflected component will always be
112     * modified by the color of the material.
113     */
114 greg 2.6 copycolor(ctmp, np->rdiff);
115     dtmp = ldot * omega / PI;
116 greg 1.1 scalecolor(ctmp, dtmp);
117     addcolor(cval, ctmp);
118 greg 2.6 } else {
119 greg 1.1 /*
120 greg 1.5 * Diffuse transmitted component.
121 greg 1.1 */
122 greg 2.6 copycolor(ctmp, np->tdiff);
123     dtmp = -ldot * omega / PI;
124 greg 1.5 scalecolor(ctmp, dtmp);
125     addcolor(cval, ctmp);
126 greg 1.1 }
127 greg 1.5 if (ldot > 0.0 ? np->rspec <= FTINY : np->tspec <= FTINY)
128     return; /* no specular component */
129     /* set up function */
130 greg 1.10 setbrdfunc(np);
131 greg 1.5 sa = np->mp->oargs.sarg;
132     errno = 0;
133     /* transform light vector */
134     multv3(ldx, ldir, funcxf.xfm);
135     for (i = 0; i < 3; i++)
136 greg 2.3 lddx[i] = ldx[i]/funcxf.sca;
137 greg 1.5 /* compute BRTDF */
138     if (np->mp->otype == MAT_BRTDF) {
139 greg 2.6 if (sa[6][0] == '0') /* special case */
140     colval(ctmp,RED) = 0.0;
141     else
142     colval(ctmp,RED) = funvalue(sa[6], 3, lddx);
143 greg 1.7 if (!strcmp(sa[7],sa[6]))
144 greg 1.5 colval(ctmp,GRN) = colval(ctmp,RED);
145     else
146 greg 2.3 colval(ctmp,GRN) = funvalue(sa[7], 3, lddx);
147 greg 1.7 if (!strcmp(sa[8],sa[6]))
148 greg 1.5 colval(ctmp,BLU) = colval(ctmp,RED);
149 greg 1.7 else if (!strcmp(sa[8],sa[7]))
150 greg 1.5 colval(ctmp,BLU) = colval(ctmp,GRN);
151     else
152 greg 2.3 colval(ctmp,BLU) = funvalue(sa[8], 3, lddx);
153 greg 1.5 dtmp = bright(ctmp);
154     } else if (np->dp == NULL) {
155 greg 2.3 dtmp = funvalue(sa[0], 3, lddx);
156 greg 1.5 setcolor(ctmp, dtmp, dtmp, dtmp);
157     } else {
158     for (i = 0; i < np->dp->nd; i++)
159 greg 2.3 pt[i] = funvalue(sa[3+i], 3, lddx);
160 greg 2.7 vldx[0] = datavalue(np->dp, pt);
161     vldx[1] = lddx[0]; vldx[2] = lddx[1]; vldx[3] = lddx[2];
162     dtmp = funvalue(sa[0], 4, vldx);
163 greg 1.5 setcolor(ctmp, dtmp, dtmp, dtmp);
164     }
165 greg 2.2 if (errno) {
166     objerror(np->mp, WARNING, "compute error");
167     return;
168     }
169 greg 1.5 if (dtmp <= FTINY)
170     return;
171     if (ldot > 0.0) {
172     /*
173     * Compute reflected non-diffuse component.
174     */
175 greg 2.6 if (np->mp->otype == MAT_MFUNC | np->mp->otype == MAT_MDATA)
176 greg 1.6 multcolor(ctmp, np->mcolor);
177     dtmp = ldot * omega * np->rspec;
178 greg 1.5 scalecolor(ctmp, dtmp);
179     addcolor(cval, ctmp);
180     } else {
181     /*
182     * Compute transmitted non-diffuse component.
183     */
184 greg 2.6 if (np->mp->otype == MAT_TFUNC | np->mp->otype == MAT_TDATA)
185 greg 1.6 multcolor(ctmp, np->mcolor);
186 greg 1.5 dtmp = -ldot * omega * np->tspec;
187     scalecolor(ctmp, dtmp);
188     addcolor(cval, ctmp);
189     }
190 greg 1.1 }
191    
192    
193 greg 2.6 m_brdf(m, r) /* color a ray which hit a BRDTF material */
194 greg 1.1 register OBJREC *m;
195     register RAY *r;
196     {
197     BRDFDAT nd;
198 greg 2.6 RAY sr;
199 greg 1.7 double transtest, transdist;
200 greg 2.6 int hasrefl, hastrans;
201 greg 1.1 COLOR ctmp;
202 greg 2.6 register MFUNC *mf;
203 greg 1.1 register int i;
204 greg 1.5 /* check arguments */
205 greg 2.6 if (m->oargs.nsargs < 10 | m->oargs.nfargs < 9)
206     objerror(m, USER, "bad # arguments");
207     nd.mp = m;
208     nd.pr = r;
209     /* dummy values */
210     nd.rspec = nd.tspec = 1.0;
211     nd.trans = 0.5;
212     /* diffuse reflectance */
213     if (r->rod > 0.0)
214     setcolor(nd.rdiff, m->oargs.farg[0],
215     m->oargs.farg[1],
216     m->oargs.farg[2]);
217     else
218     setcolor(nd.rdiff, m->oargs.farg[3],
219     m->oargs.farg[4],
220     m->oargs.farg[5]);
221     /* diffuse transmittance */
222     setcolor(nd.tdiff, m->oargs.farg[6],
223     m->oargs.farg[7],
224     m->oargs.farg[8]);
225     /* get modifiers */
226     raytexture(r, m->omod);
227     nd.pdot = raynormal(nd.pnorm, r); /* perturb normal */
228     if (r->rod < 0.0) { /* orient perturbed values */
229     nd.pdot = -nd.pdot;
230     for (i = 0; i < 3; i++) {
231     nd.pnorm[i] = -nd.pnorm[i];
232     r->pert[i] = -r->pert[i];
233     }
234 greg 1.5 }
235 greg 2.6 copycolor(nd.mcolor, r->pcol); /* get pattern color */
236     multcolor(nd.rdiff, nd.mcolor); /* modify diffuse values */
237     multcolor(nd.tdiff, nd.mcolor);
238     hasrefl = bright(nd.rdiff) > FTINY;
239     hastrans = bright(nd.tdiff) > FTINY;
240     /* load cal file */
241     nd.dp = NULL;
242     mf = getfunc(m, 9, 0x3f, 0);
243     /* compute transmitted ray */
244     setbrdfunc(&nd);
245     transtest = 0;
246 greg 2.8 transdist = r->rot;
247 greg 2.6 errno = 0;
248     setcolor(ctmp, evalue(mf->ep[3]),
249     evalue(mf->ep[4]),
250     evalue(mf->ep[5]));
251     if (errno)
252     objerror(m, WARNING, "compute error");
253     else if (rayorigin(&sr, r, TRANS, bright(ctmp)) == 0) {
254     if (!(r->crtype & SHADOW) &&
255     DOT(r->pert,r->pert) > FTINY*FTINY) {
256     for (i = 0; i < 3; i++) /* perturb direction */
257     sr.rdir[i] = r->rdir[i] - .75*r->pert[i];
258     if (normalize(sr.rdir) == 0.0) {
259     objerror(m, WARNING, "illegal perturbation");
260     VCOPY(sr.rdir, r->rdir);
261     }
262     } else {
263     VCOPY(sr.rdir, r->rdir);
264     transtest = 2;
265     }
266     rayvalue(&sr);
267     multcolor(sr.rcol, ctmp);
268     addcolor(r->rcol, sr.rcol);
269     transtest *= bright(sr.rcol);
270     transdist = r->rot + sr.rt;
271     }
272     if (r->crtype & SHADOW) /* the rest is shadow */
273 greg 2.10 return(1);
274 greg 2.6 /* compute reflected ray */
275     setbrdfunc(&nd);
276     errno = 0;
277     setcolor(ctmp, evalue(mf->ep[0]),
278     evalue(mf->ep[1]),
279     evalue(mf->ep[2]));
280     if (errno)
281     objerror(m, WARNING, "compute error");
282     else if (rayorigin(&sr, r, REFLECTED, bright(ctmp)) == 0) {
283     for (i = 0; i < 3; i++)
284     sr.rdir[i] = r->rdir[i] + 2.0*nd.pdot*nd.pnorm[i];
285     rayvalue(&sr);
286     multcolor(sr.rcol, ctmp);
287     addcolor(r->rcol, sr.rcol);
288     }
289     /* compute ambient */
290     if (hasrefl) {
291     if (nd.pdot < 0.0)
292     flipsurface(r);
293     ambient(ctmp, r);
294     multcolor(ctmp, nd.rdiff);
295     addcolor(r->rcol, ctmp); /* add to returned color */
296     if (nd.pdot < 0.0)
297     flipsurface(r);
298     }
299     if (hastrans) { /* from other side */
300     if (nd.pdot > 0.0)
301     flipsurface(r);
302     ambient(ctmp, r);
303     multcolor(ctmp, nd.tdiff);
304     addcolor(r->rcol, ctmp);
305     if (nd.pdot > 0.0)
306     flipsurface(r);
307     }
308     if (hasrefl | hastrans || m->oargs.sarg[6][0] != '0')
309     direct(r, dirbrdf, &nd); /* add direct component */
310     /* check distance */
311     if (transtest > bright(r->rcol))
312     r->rt = transdist;
313 greg 2.10
314     return(1);
315 greg 2.6 }
316    
317    
318    
319     m_brdf2(m, r) /* color a ray which hit a BRDF material */
320     register OBJREC *m;
321     register RAY *r;
322     {
323     BRDFDAT nd;
324     COLOR ctmp;
325     double dtmp;
326     /* always a shadow */
327     if (r->crtype & SHADOW)
328 greg 2.10 return(1);
329 greg 2.6 /* check arguments */
330     if (m->oargs.nsargs < (hasdata(m->otype)?4:2) | m->oargs.nfargs <
331     (m->otype==MAT_TFUNC|m->otype==MAT_TDATA?6:4))
332 greg 1.1 objerror(m, USER, "bad # arguments");
333     nd.mp = m;
334     nd.pr = r;
335 greg 2.6 /* get material color */
336     setcolor(nd.mcolor, m->oargs.farg[0],
337     m->oargs.farg[1],
338     m->oargs.farg[2]);
339 greg 1.5 /* get specular component */
340     nd.rspec = m->oargs.farg[3];
341 greg 2.6 /* compute transmittance */
342     if (m->otype == MAT_TFUNC | m->otype == MAT_TDATA) {
343 greg 1.5 nd.trans = m->oargs.farg[4]*(1.0 - nd.rspec);
344     nd.tspec = nd.trans * m->oargs.farg[5];
345 greg 2.6 dtmp = nd.trans - nd.tspec;
346     setcolor(nd.tdiff, dtmp, dtmp, dtmp);
347     } else {
348     nd.tspec = nd.trans = 0.0;
349     setcolor(nd.tdiff, 0.0, 0.0, 0.0);
350     }
351     /* compute reflectance */
352     dtmp = 1.0 - nd.trans - nd.rspec;
353     setcolor(nd.rdiff, dtmp, dtmp, dtmp);
354 greg 2.11 /* check for back side */
355     if (r->rod < 0.0) {
356     if (!backvis && m->otype != MAT_TFUNC
357     && m->otype != MAT_TDATA) {
358     raytrans(r);
359     return(1);
360     }
361     flipsurface(r); /* reorient if backvis */
362     }
363 greg 1.5 /* get modifiers */
364     raytexture(r, m->omod);
365     nd.pdot = raynormal(nd.pnorm, r); /* perturb normal */
366     multcolor(nd.mcolor, r->pcol); /* modify material color */
367 greg 2.6 multcolor(nd.rdiff, nd.mcolor);
368     multcolor(nd.tdiff, nd.mcolor);
369 greg 1.1 /* load auxiliary files */
370 greg 2.2 if (hasdata(m->otype)) {
371 greg 1.1 nd.dp = getdata(m->oargs.sarg[1]);
372 greg 2.6 getfunc(m, 2, 0, 0);
373 greg 1.1 } else {
374     nd.dp = NULL;
375 greg 2.6 getfunc(m, 1, 0, 0);
376 greg 1.1 }
377     /* compute ambient */
378 greg 2.6 if (nd.trans < 1.0-FTINY) {
379 greg 1.1 ambient(ctmp, r);
380 greg 2.6 scalecolor(ctmp, 1.0-nd.trans);
381 greg 1.1 multcolor(ctmp, nd.mcolor); /* modified by material color */
382     addcolor(r->rcol, ctmp); /* add to returned color */
383 greg 1.5 }
384 greg 2.6 if (nd.trans > FTINY) { /* from other side */
385 greg 1.5 flipsurface(r);
386     ambient(ctmp, r);
387 greg 2.6 scalecolor(ctmp, nd.trans);
388 greg 1.5 multcolor(ctmp, nd.mcolor);
389     addcolor(r->rcol, ctmp);
390     flipsurface(r);
391 greg 1.1 }
392     /* add direct component */
393     direct(r, dirbrdf, &nd);
394 greg 2.10
395     return(1);
396 greg 1.10 }
397    
398    
399     setbrdfunc(np) /* set up brdf function and variables */
400     register BRDFDAT *np;
401     {
402     FVECT vec;
403    
404     if (setfunc(np->mp, np->pr) == 0)
405     return(0); /* it's OK, setfunc says we're done */
406     /* else (re)assign special variables */
407     multv3(vec, np->pnorm, funcxf.xfm);
408     varset("NxP", '=', vec[0]/funcxf.sca);
409     varset("NyP", '=', vec[1]/funcxf.sca);
410     varset("NzP", '=', vec[2]/funcxf.sca);
411 greg 1.11 varset("RdotP", '=', np->pdot <= -1.0 ? -1.0 :
412     np->pdot >= 1.0 ? 1.0 : np->pdot);
413 greg 1.10 varset("CrP", '=', colval(np->mcolor,RED));
414     varset("CgP", '=', colval(np->mcolor,GRN));
415     varset("CbP", '=', colval(np->mcolor,BLU));
416     return(1);
417 greg 1.1 }