--- ray/src/rt/srcsupp.c 1992/01/04 19:54:20 2.2 +++ ray/src/rt/srcsupp.c 2003/02/22 02:07:29 2.9 @@ -1,13 +1,69 @@ -/* Copyright (c) 1992 Regents of the University of California */ - #ifndef lint -static char SCCSid[] = "$SunId$ LBL"; +static const char RCSid[] = "$Id: srcsupp.c,v 2.9 2003/02/22 02:07:29 greg Exp $"; #endif - /* * Support routines for source objects and materials + * + * External symbols declared in source.h */ +/* ==================================================================== + * The Radiance Software License, Version 1.0 + * + * Copyright (c) 1990 - 2002 The Regents of the University of California, + * through Lawrence Berkeley National Laboratory. All rights reserved. + * + * Redistribution and use in source and binary forms, with or without + * modification, are permitted provided that the following conditions + * are met: + * + * 1. Redistributions of source code must retain the above copyright + * notice, this list of conditions and the following disclaimer. + * + * 2. Redistributions in binary form must reproduce the above copyright + * notice, this list of conditions and the following disclaimer in + * the documentation and/or other materials provided with the + * distribution. + * + * 3. The end-user documentation included with the redistribution, + * if any, must include the following acknowledgment: + * "This product includes Radiance software + * (http://radsite.lbl.gov/) + * developed by the Lawrence Berkeley National Laboratory + * (http://www.lbl.gov/)." + * Alternately, this acknowledgment may appear in the software itself, + * if and wherever such third-party acknowledgments normally appear. + * + * 4. The names "Radiance," "Lawrence Berkeley National Laboratory" + * and "The Regents of the University of California" must + * not be used to endorse or promote products derived from this + * software without prior written permission. For written + * permission, please contact radiance@radsite.lbl.gov. + * + * 5. Products derived from this software may not be called "Radiance", + * nor may "Radiance" appear in their name, without prior written + * permission of Lawrence Berkeley National Laboratory. + * + * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESSED OR IMPLIED + * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES + * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE + * DISCLAIMED. IN NO EVENT SHALL Lawrence Berkeley National Laboratory OR + * ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, + * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT + * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF + * USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND + * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, + * OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT + * OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF + * SUCH DAMAGE. + * ==================================================================== + * + * This software consists of voluntary contributions made by many + * individuals on behalf of Lawrence Berkeley National Laboratory. For more + * information on Lawrence Berkeley National Laboratory, please see + * . + */ + #include "ray.h" #include "otypes.h" @@ -26,13 +82,10 @@ int nsources = 0; /* the number of sources */ SRCFUNC sfun[NUMOTYPE]; /* source dispatch table */ +void initstypes() /* initialize source dispatch table */ { extern VSMATERIAL mirror_vs, direct1_vs, direct2_vs; - extern int fsetsrc(), ssetsrc(), sphsetsrc(), cylsetsrc(), rsetsrc(); - extern int nopart(), flatpart(), cylpart(); - extern double fgetplaneq(), rgetplaneq(); - extern double fgetmaxdisk(), rgetmaxdisk(); static SOBJECT fsobj = {fsetsrc, flatpart, fgetplaneq, fgetmaxdisk}; static SOBJECT ssobj = {ssetsrc, nopart}; static SOBJECT sphsobj = {sphsetsrc, nopart}; @@ -67,6 +120,7 @@ newsource() /* allocate new source in our array */ } +void setflatss(src) /* set sampling for a flat source */ register SRCREC *src; { @@ -86,6 +140,7 @@ register SRCREC *src; } +void fsetsrc(src, so) /* set a face as a source */ register SRCREC *src; OBJREC *so; @@ -119,8 +174,7 @@ OBJREC *so; } src->srad = sqrt(src->srad); /* compute size vectors */ - if (f->nv == 4 || (f->nv == 5 && /* parallelogram case */ - dist2(VERTEX(f,0),VERTEX(f,4)) <= FTINY*FTINY)) + if (f->nv == 4) /* parallelogram case */ for (j = 0; j < 3; j++) { src->ss[SU][j] = .5*(VERTEX(f,1)[j]-VERTEX(f,0)[j]); src->ss[SV][j] = .5*(VERTEX(f,3)[j]-VERTEX(f,0)[j]); @@ -130,6 +184,7 @@ OBJREC *so; } +void ssetsrc(src, so) /* set a source as a source */ register SRCREC *src; register OBJREC *so; @@ -156,6 +211,7 @@ register OBJREC *so; } +void sphsetsrc(src, so) /* set a sphere as a source */ register SRCREC *src; register OBJREC *so; @@ -178,6 +234,7 @@ register OBJREC *so; } +void rsetsrc(src, so) /* set a ring (disk) as a source */ register SRCREC *src; OBJREC *so; @@ -199,6 +256,7 @@ OBJREC *so; } +void cylsetsrc(src, so) /* set a cylinder as a source */ register SRCREC *src; OBJREC *so; @@ -239,16 +297,49 @@ register OBJREC *m; { register SPOT *ns; + if ((ns = (SPOT *)m->os) != NULL) + return(ns); if ((ns = (SPOT *)malloc(sizeof(SPOT))) == NULL) return(NULL); ns->siz = 2.0*PI * (1.0 - cos(PI/180.0/2.0 * m->oargs.farg[3])); VCOPY(ns->aim, m->oargs.farg+4); if ((ns->flen = normalize(ns->aim)) == 0.0) objerror(m, USER, "zero focus vector"); + m->os = (char *)ns; return(ns); } +int +spotout(r, s) /* check if we're outside spot region */ +register RAY *r; +register SPOT *s; +{ + double d; + FVECT vd; + + if (s == NULL) + return(0); + if (s->flen < -FTINY) { /* distant source */ + vd[0] = s->aim[0] - r->rorg[0]; + vd[1] = s->aim[1] - r->rorg[1]; + vd[2] = s->aim[2] - r->rorg[2]; + d = DOT(r->rdir,vd); + /* wrong side? + if (d <= FTINY) + return(1); */ + d = DOT(vd,vd) - d*d; + if (PI*d > s->siz) + return(1); /* out */ + return(0); /* OK */ + } + /* local source */ + if (s->siz < 2.0*PI * (1.0 + DOT(s->aim,r->rdir))) + return(1); /* out */ + return(0); /* OK */ +} + + double fgetmaxdisk(ocent, op) /* get center and squared radius of face */ FVECT ocent; @@ -320,6 +411,7 @@ OBJREC *op; } +int commonspot(sp1, sp2, org) /* set sp1 to intersection of sp1 and sp2 */ register SPOT *sp1, *sp2; FVECT org; @@ -343,6 +435,7 @@ FVECT org; } +int commonbeam(sp1, sp2, dir) /* set sp1 to intersection of sp1 and sp2 */ register SPOT *sp1, *sp2; FVECT dir; @@ -367,6 +460,7 @@ FVECT dir; } +int checkspot(sp, nrm) /* check spotlight for behind source */ register SPOT *sp; /* spotlight */ FVECT nrm; /* source surface normal */ @@ -462,154 +556,4 @@ double r1s, r2s; /* radii squared */ for (i = 0; i < 3; i++) cc[i] = c1[i] + l*disp[i]; return(a2); -} - - -sourcehit(r) /* check to see if ray hit distant source */ -register RAY *r; -{ - int first, last; - register int i; - - if (r->rsrc >= 0) { /* check only one if aimed */ - first = last = r->rsrc; - } else { /* otherwise check all */ - first = 0; last = nsources-1; - } - for (i = first; i <= last; i++) - if ((source[i].sflags & (SDISTANT|SVIRTUAL)) == SDISTANT) - /* - * Check to see if ray is within - * solid angle of source. - */ - if (2.0*PI * (1.0 - DOT(source[i].sloc,r->rdir)) - <= source[i].ss2) { - r->ro = source[i].so; - if (!(source[i].sflags & SSKIP)) - break; - } - - if (r->ro != NULL) { - for (i = 0; i < 3; i++) - r->ron[i] = -r->rdir[i]; - r->rod = 1.0; - r->rox = NULL; - return(1); - } - return(0); -} - - -/**************************************************************** - * The following macros were separated from the m_light() routine - * because they are very nasty and difficult to understand. - */ - -/* wrongillum * - * - * We cannot allow an illum to pass to another illum, because that - * would almost certainly constitute overcounting. - * However, we do allow an illum to pass to another illum - * that is actually going to relay to a virtual light source. - */ - -#define wrongillum(m, r) (!(source[r->rsrc].sflags&SVIRTUAL) && \ - objptr(source[r->rsrc].so->omod)->otype==MAT_ILLUM) - -/* wrongsource * - * - * This source is the wrong source (ie. overcounted) if we are - * aimed to a different source than the one we hit and the one - * we hit is not an illum which should be passed. - */ - -#define wrongsource(m, r) (r->rsrc>=0 && source[r->rsrc].so!=r->ro && \ - (m->otype!=MAT_ILLUM || wrongillum(m,r))) - -/* badspecular * - * - * Any undirected specular ray that hits a light source - * should be discarded. This is because the source contribution to - * specular components is calculated separately to reduce variance. - */ - -#define badspecular(m, r) (r->rsrc<0 && r->crtype&SPECULAR) - -/* distglow * - * - * A distant glow is an object that sometimes acts as a light source, - * but is too far away from the test point to be one in this case. - */ - -#define distglow(m, r) (m->otype==MAT_GLOW && \ - r->rot > m->oargs.farg[3]) - -/* badambient * - * - * We must avoid including counting light sources in the ambient calculation, - * since the direct component is handled separately. Therefore, any - * ambient ray which hits an active light source must be discarded. - */ - -#define badambient(m, r) ((r->crtype&(AMBIENT|SHADOW))==AMBIENT && \ - !distglow(m, r)) - -/* overcount * - * - * All overcounting possibilities are contained here. - */ - -#define overcount(m, r) (badspecular(m,r) || wrongsource(m,r) || \ - badambient(m,r)) - -/* passillum * - * - * An illum passes to another material type when we didn't hit it - * on purpose (as part of a direct calculation), or it is relaying - * a virtual light source. - */ - -#define passillum(m, r) (m->otype==MAT_ILLUM && \ - (r->rsrc<0 || source[r->rsrc].so!=r->ro || \ - source[r->rsrc].sflags&SVIRTUAL)) - -/* srcignore * - * - * The -di flag renders light sources invisible, and here is the test. - */ - -#define srcignore(m, r) (directinvis && !(r->crtype&SHADOW) && \ - !distglow(m, r)) - - -m_light(m, r) /* ray hit a light source */ -register OBJREC *m; -register RAY *r; -{ - /* check for over-counting */ - if (overcount(m, r)) - return; - /* check for passed illum */ - if (passillum(m, r)) { - if (m->oargs.nsargs < 1 || !strcmp(m->oargs.sarg[0], VOIDID)) - raytrans(r); - else - rayshade(r, modifier(m->oargs.sarg[0])); - return; - } - /* otherwise treat as source */ - /* check for behind */ - if (r->rod < 0.0) - return; - /* check for invisibility */ - if (srcignore(m, r)) - return; - /* get distribution pattern */ - raytexture(r, m->omod); - /* get source color */ - setcolor(r->rcol, m->oargs.farg[0], - m->oargs.farg[1], - m->oargs.farg[2]); - /* modify value */ - multcolor(r->rcol, r->pcol); }