ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/rt/virtuals.c
(Generate patch)

Comparing ray/src/rt/virtuals.c (file contents):
Revision 1.4 by greg, Thu Jun 20 16:36:48 1991 UTC vs.
Revision 1.8 by greg, Tue Jun 25 08:52:04 1991 UTC

# Line 11 | Line 11 | static char SCCSid[] = "$SunId$ LBL";
11  
12   #include  "ray.h"
13  
14 + #include  "octree.h"
15 +
16   #include  "otypes.h"
17  
18   #include  "source.h"
19  
20 + #include  "random.h"
21  
19 double  intercircle();
22  
23 + double  getdisk();
24 +
25   static OBJECT  *vobject;                /* virtual source objects */
26   static int  nvobjects = 0;              /* number of virtual source objects */
27  
# Line 55 | Line 59 | markvirtuals()                 /* find and mark virtual sources */
59   #endif
60                                          /* append virtual sources */
61          for (i = nsources; i-- > 0; )
62 <                if (!(source[i].sflags & SSKIP))
59 <                        addvirtuals(i, directrelay);
62 >                addvirtuals(i, directrelay);
63                                          /* done with our object list */
64          free((char *)vobject);
65          nvobjects = 0;
# Line 71 | Line 74 | int  nr;
74                                  /* check relay limit first */
75          if (nr <= 0)
76                  return;
77 +        if (source[sn].sflags & SSKIP)
78 +                return;
79                                  /* check each virtual object for projection */
80          for (i = 0; i < nvobjects; i++)
81                                          /* vproject() calls us recursively */
# Line 97 | Line 102 | int  n;
102                  if ((*vsmat->vproj)(proj, o, &source[sn], i))
103                          if ((ns = makevsrc(o, sn, proj)) >= 0) {
104   #ifdef DEBUG
105 <                                virtverb(&source[sn], stderr);
105 >                                virtverb(ns, stderr);
106   #endif
107                                  addvirtuals(ns, n);
108                          }
# Line 110 | Line 115 | OBJREC  *op;
115   register int  sn;
116   MAT4  pm;
117   {
118 <        register int  nsn;
114 <        FVECT  nsloc, ocent, nsnorm;
115 <        int  nsflags;
118 >        FVECT  nsloc, nsnorm, ocent;
119          double  maxrad2;
120 +        int  nsflags;
121          double  d1;
122          SPOT  theirspot, ourspot;
123          register int  i;
124  
125 <        nsflags = (source[sn].sflags|(SVIRTUAL|SFOLLOW)) & ~SSPOT;
125 >        nsflags = source[sn].sflags | (SVIRTUAL|SSPOT|SFOLLOW);
126                                          /* get object center and max. radius */
127 <        if (sfun[op->otype].of->getdisk != NULL) {
128 <                maxrad2 = (*sfun[op->otype].of->getdisk)(ocent, op);
129 <                if (maxrad2 <= FTINY)                   /* too small? */
126 <                        return(NULL);
127 <                nsflags |= SSPOT;
128 <        }
127 >        maxrad2 = getdisk(ocent, op, sn);
128 >        if (maxrad2 <= FTINY)                   /* too small? */
129 >                return(-1);
130                                          /* get location and spot */
131          if (source[sn].sflags & SDISTANT) {             /* distant source */
132                  if (source[sn].sflags & SPROX)
133 <                        return(NULL);           /* should never get here! */
133 >                        return(-1);             /* should never get here! */
134                  multv3(nsloc, source[sn].sloc, pm);
135 <                if (nsflags & SSPOT) {
136 <                        VCOPY(ourspot.aim, ocent);
137 <                        ourspot.siz = PI*maxrad2;
137 <                        ourspot.flen = 0.;
138 <                }
135 >                VCOPY(ourspot.aim, ocent);
136 >                ourspot.siz = PI*maxrad2;
137 >                ourspot.flen = 0.;
138                  if (source[sn].sflags & SSPOT) {
139                          copystruct(&theirspot, source[sn].sl.s);
140                          multp3(theirspot.aim, source[sn].sl.s->aim, pm);
141 <                        if (nsflags & SSPOT &&
142 <                                !commonbeam(&ourspot, &theirspot, nsloc))
144 <                                return(NULL);           /* no overlap */
141 >                        if (!commonbeam(&ourspot, &theirspot, nsloc))
142 >                                return(-1);             /* no overlap */
143                  }
144          } else {                                /* local source */
145                  multp3(nsloc, source[sn].sloc, pm);
146 <                if (nsflags & SSPOT) {
147 <                        for (i = 0; i < 3; i++)
148 <                                ourspot.aim[i] = ocent[i] - nsloc[i];
149 <                        if ((d1 = normalize(ourspot.aim)) == 0.)
150 <                                return(NULL);           /* at source!! */
151 <                        if (source[sn].sflags & SPROX &&
152 <                                        d1 > source[sn].sl.prox)
153 <                                return(NULL);           /* too far away */
156 <                        ourspot.siz = 2.*PI*(1. - d1/sqrt(d1*d1+maxrad2));
157 <                        ourspot.flen = 0.;
158 <                } else if (source[sn].sflags & SPROX) {
159 <                        FVECT  norm;
160 <                        double  offs;
161 <                                                /* use distance from plane */
162 <                        offs = (*sfun[op->otype].of->getpleq)(norm, op);
163 <                        d1 = DOT(norm, nsloc) - offs;
164 <                        if (d1 < 0.) d1 = -d1;
165 <                        if (d1 > source[sn].sl.prox)
166 <                                return(NULL);           /* too far away */
167 <                }
146 >                for (i = 0; i < 3; i++)
147 >                        ourspot.aim[i] = ocent[i] - nsloc[i];
148 >                if ((d1 = normalize(ourspot.aim)) == 0.)
149 >                        return(-1);             /* at source!! */
150 >                if (source[sn].sflags & SPROX && d1 > source[sn].sl.prox)
151 >                        return(-1);             /* too far away */
152 >                ourspot.siz = 2.*PI*(1. - d1/sqrt(d1*d1+maxrad2));
153 >                ourspot.flen = 0.;
154                  if (source[sn].sflags & SSPOT) {
155                          copystruct(&theirspot, source[sn].sl.s);
156                          multv3(theirspot.aim, source[sn].sl.s->aim, pm);
157 <                        if (nsflags & SSPOT) {
158 <                                if (!commonspot(&ourspot, &theirspot, nsloc))
159 <                                        return(NULL);   /* no overlap */
174 <                                ourspot.flen = theirspot.flen;
175 <                        }
157 >                        if (!commonspot(&ourspot, &theirspot, nsloc))
158 >                                return(-1);     /* no overlap */
159 >                        ourspot.flen = theirspot.flen;
160                  }
161                  if (source[sn].sflags & SFLAT) {        /* behind source? */
162                          multv3(nsnorm, source[sn].snorm, pm);
163 <                        if (nsflags & SSPOT && checkspot(&ourspot, nsnorm) < 0)
164 <                                return(NULL);
163 >                        if (checkspot(&ourspot, nsnorm) < 0)
164 >                                return(-1);
165                  }
166          }
167 <                                        /* everything is OK, make source */
168 <        if ((nsn = newsource()) < 0)
167 >                                        /* pretest visibility */
168 >        nsflags = vstestvis(nsflags, op, ocent, maxrad2, sn);
169 >        if (nsflags & SSKIP)
170 >                return(-1);     /* obstructed */
171 >                                        /* it all checks out, so make it */
172 >        if ((i = newsource()) < 0)
173                  goto memerr;
174 <        source[nsn].sflags = nsflags;
175 <        VCOPY(source[nsn].sloc, nsloc);
174 >        source[i].sflags = nsflags;
175 >        VCOPY(source[i].sloc, nsloc);
176          if (nsflags & SFLAT)
177 <                VCOPY(source[nsn].snorm, nsnorm);
178 <        source[nsn].ss = source[sn].ss; source[nsn].ss2 = source[sn].ss2;
179 <        if ((nsflags | source[sn].sflags) & SSPOT) {
180 <                if ((source[nsn].sl.s = (SPOT *)malloc(sizeof(SPOT))) == NULL)
181 <                        goto memerr;
194 <                if (nsflags & SSPOT)
195 <                        copystruct(source[nsn].sl.s, &ourspot);
196 <                else
197 <                        copystruct(source[nsn].sl.s, &theirspot);
198 <                source[nsn].sflags |= SSPOT;
199 <        }
177 >                VCOPY(source[i].snorm, nsnorm);
178 >        source[i].ss = source[sn].ss; source[i].ss2 = source[sn].ss2;
179 >        if ((source[i].sl.s = (SPOT *)malloc(sizeof(SPOT))) == NULL)
180 >                goto memerr;
181 >        copystruct(source[i].sl.s, &ourspot);
182          if (nsflags & SPROX)
183 <                source[nsn].sl.prox = source[sn].sl.prox;
184 <        source[nsn].sa.svnext = sn;
185 <        source[nsn].so = op;
186 <        return(nsn);
183 >                source[i].sl.prox = source[sn].sl.prox;
184 >        source[i].sa.svnext = sn;
185 >        source[i].so = op;
186 >        return(i);
187   memerr:
188          error(SYSTEM, "out of memory in makevsrc");
189   }
190  
191  
192 < commonspot(sp1, sp2, org)       /* set sp1 to intersection of sp1 and sp2 */
193 < register SPOT  *sp1, *sp2;
194 < FVECT  org;
192 > double
193 > getdisk(oc, op, sn)             /* get visible object disk */
194 > FVECT  oc;
195 > OBJREC  *op;
196 > register int  sn;
197   {
198 <        FVECT  cent;
199 <        double  rad2, cos1, cos2;
200 <
201 <        cos1 = 1. - sp1->siz/(2.*PI);
202 <        cos2 = 1. - sp2->siz/(2.*PI);
203 <        if (sp2->siz >= 2.*PI-FTINY)            /* BIG, just check overlap */
204 <                return(DOT(sp1->aim,sp2->aim) >= cos1*cos2 -
205 <                                        sqrt((1.-cos1*cos1)*(1.-cos2*cos2)));
206 <                                /* compute and check disks */
207 <        rad2 = intercircle(cent, sp1->aim, sp2->aim,
208 <                        1./(cos1*cos1) - 1.,  1./(cos2*cos2) - 1.);
209 <        if (rad2 <= FTINY || normalize(cent) == 0.)
210 <                return(0);
211 <        VCOPY(sp1->aim, cent);
212 <        sp1->siz = 2.*PI*(1. - 1./sqrt(1.+rad2));
213 <        return(1);
198 >        double  rad2, roffs, offs, d, rd, rdoto;
199 >        FVECT  rnrm, nrm;
200 >                                /* first, use object getdisk function */
201 >        rad2 = (*sfun[op->otype].of->getdisk)(oc, op);
202 >        if (!(source[sn].sflags & SVIRTUAL))
203 >                return(rad2);           /* all done for normal source */
204 >                                /* check for correct side of relay surface */
205 >        roffs = (*sfun[source[sn].so->otype].of->getpleq)(rnrm, source[sn].so);
206 >        rd = DOT(rnrm, source[sn].sloc);        /* source projection */
207 >        if (!(source[sn].sflags & SDISTANT))
208 >                rd -= roffs;
209 >        d = DOT(rnrm, oc) - roffs;      /* disk distance to relay plane */
210 >        if ((d > 0.) ^ (rd > 0.))
211 >                return(rad2);           /* OK if opposite sides */
212 >        if (d*d >= rad2)
213 >                return(.0);             /* no relay is possible */
214 >                                /* we need a closer look */
215 >        offs = (*sfun[op->otype].of->getpleq)(nrm, op);
216 >        rdoto = DOT(rnrm, nrm);
217 >        if (d*d >= rad2*(1.-rdoto*rdoto))
218 >                return(0.);             /* disk entirely on projection side */
219 >                                /* should shrink disk but I'm lazy */
220 >        return(rad2);
221   }
222  
223  
224 < commonbeam(sp1, sp2, dir)       /* set sp1 to intersection of sp1 and sp2 */
225 < register SPOT  *sp1, *sp2;
226 < FVECT  dir;
224 > int
225 > vstestvis(f, o, oc, or2, sn)            /* pretest source visibility */
226 > int  f;                 /* virtual source flags */
227 > OBJREC  *o;             /* relay object */
228 > FVECT  oc;              /* relay object center */
229 > double  or2;            /* relay object radius squared */
230 > register int  sn;       /* target source number */
231   {
232 <        FVECT  cent, c1, c2;
233 <        double  rad2, d;
234 <        register int  i;
235 <                                        /* move centers to common plane */
236 <        d = DOT(sp1->aim, dir);
237 <        for (i = 0; i < 3; i++)
238 <                c1[i] = sp1->aim[i] - d*dir[i];
239 <        d = DOT(sp2->aim, dir);
240 <        for (i = 0; i < 3; i++)
241 <                c2[i] = sp2->aim[i] - d*dir[i];
242 <                                        /* compute overlap */
243 <        rad2 = intercircle(cent, c1, c2, sp1->siz/PI, sp2->siz/PI);
244 <        if (rad2 <= FTINY)
245 <                return(0);
246 <        VCOPY(sp1->aim, cent);
247 <        sp1->siz = PI*rad2;
248 <        return(1);
249 < }
255 <
256 <
257 < checkspot(sp, nrm)              /* check spotlight for behind source */
258 < register SPOT  *sp;
259 < FVECT  nrm;
260 < {
261 <        double  d, d1;
262 <
263 <        d = DOT(sp->aim, nrm);
264 <        if (d > FTINY)                  /* center in front? */
265 <                return(0);
266 <                                        /* else check horizon */
267 <        d1 = 1. - sp->siz/(2.*PI);
268 <        return(1.-FTINY-d*d > d1*d1);
269 < }
270 <
271 <
272 < double
273 < intercircle(cc, c1, c2, r1s, r2s)       /* intersect two circles */
274 < FVECT  cc;                      /* midpoint (return value) */
275 < FVECT  c1, c2;                  /* circle centers */
276 < double  r1s, r2s;               /* radii squared */
277 < {
278 <        double  a2, d2, l;
279 <        FVECT  disp;
280 <        register int  i;
281 <
282 <        for (i = 0; i < 3; i++)
283 <                disp[i] = c2[i] - c1[i];
284 <        d2 = DOT(disp,disp);
285 <                                        /* circle within overlap? */
286 <        if (r1s < r2s) {
287 <                if (r2s >= r1s + d2) {
288 <                        VCOPY(cc, c1);
289 <                        return(r1s);
290 <                }
232 >        RAY  sr;
233 >        FVECT  onorm;
234 >        FVECT  offsdir;
235 >        double  or, d;
236 >        int  infront;
237 >        int  ssn;
238 >        int  nok, nhit;
239 >        register int  i, n;
240 >                                /* return if pretesting disabled */
241 >        if (vspretest <= 0)
242 >                return(f);
243 >                                /* get surface normal */
244 >        (*sfun[o->otype].of->getpleq)(onorm, o);
245 >                                /* set number of rays to sample */
246 >        if (source[sn].sflags & SDISTANT) {
247 >                n = (2./3.*PI*PI)*or2/(thescene.cusize*thescene.cusize)*
248 >                                vspretest + .5;
249 >                infront = DOT(onorm, source[sn].sloc) > 0.;
250          } else {
251 <                if (r1s >= r2s + d2) {
252 <                        VCOPY(cc, c2);
253 <                        return(r2s);
254 <                }
251 >                n = or2/dist2(oc,source[sn].sloc)*vspretest + .5;
252 >                for (i = 0; i < 3; i++)
253 >                        offsdir[i] = source[sn].sloc[i] - oc[i];
254 >                infront = DOT(onorm, offsdir) > 0.;
255          }
256 <        a2 = .25*(2.*(r1s+r2s) - d2 - (r2s-r1s)*(r2s-r1s)/d2);
257 <                                        /* no overlap? */
258 <        if (a2 <= 0.)
259 <                return(0.);
260 <                                        /* overlap, compute center */
261 <        l = sqrt((r1s - a2)/d2);
262 <        for (i = 0; i < 3; i++)
263 <                cc[i] = c1[i] + l*disp[i];
264 <        return(a2);
256 >        if (n < 1) n = 1;
257 >                                /* sample */
258 >        or = sqrt(or2);
259 >        ssn = 7*n;
260 >        nhit = nok = 0;
261 >        while (n-- > 0) {
262 >                                        /* get sample point */
263 >                do {
264 >                        if (--ssn < 0)
265 >                                return(f);      /* too small a target! */
266 >                        for (i = 0; i < 3; i++)
267 >                                offsdir[i] = or*(1. -
268 >                                                2.*urand(931*i+5827+ssn));
269 >                        for (i = 0; i < 3; i++)
270 >                                sr.rorg[i] = oc[i] + offsdir[i];
271 >                        d = DOT(offsdir,onorm);
272 >                        if (infront)
273 >                                for (i = 0; i < 3; i++) {
274 >                                        sr.rorg[i] -= (d-.0001)*onorm[i];
275 >                                        sr.rdir[i] = -onorm[i];
276 >                                }
277 >                        else
278 >                                for (i = 0; i < 3; i++) {
279 >                                        sr.rorg[i] -= (d+.0001)*onorm[i];
280 >                                        sr.rdir[i] = onorm[i];
281 >                                }
282 >                        rayorigin(&sr, NULL, PRIMARY, 1.0);
283 >                } while (!(*ofun[o->otype].funp)(o, &sr));
284 >                                        /* check against source */
285 >                samplendx++;
286 >                if (srcray(&sr, NULL, sn) == 0.0)
287 >                        continue;
288 >                sr.revf = srcvalue;
289 >                rayvalue(&sr);
290 >                if (bright(sr.rcol) <= FTINY)
291 >                        continue;
292 >                nok++;
293 >                                        /* check against obstructions */
294 >                srcray(&sr, NULL, sn);
295 >                rayvalue(&sr);
296 >                if (bright(sr.rcol) <= FTINY)
297 >                        continue;
298 >                nhit++;
299 >        }
300 >                                /* interpret results */
301 >        if (nhit == 0)
302 >                return(f | SSKIP);      /* 0% hit rate:  totally occluded */
303 >        if (nhit == nok)
304 >                return(f & ~SFOLLOW);   /* 100% hit rate:  no occlusion */
305 >        return(f);              /* no comment */
306   }
307 +        
308  
308
309   #ifdef DEBUG
310 < virtverb(vs, fp)        /* print verbose description of virtual source */
311 < register SRCREC  *vs;
310 > virtverb(sn, fp)        /* print verbose description of virtual source */
311 > register int  sn;
312   FILE  *fp;
313   {
314          register int  i;
315  
316          fprintf(fp, "%s virtual source %d in %s %s\n",
317 <                        vs->sflags & SDISTANT ? "distant" : "local",
318 <                        vs-source, ofun[vs->so->otype].funame, vs->so->oname);
317 >                        source[sn].sflags & SDISTANT ? "distant" : "local",
318 >                        sn, ofun[source[sn].so->otype].funame,
319 >                        source[sn].so->oname);
320          fprintf(fp, "\tat (%f,%f,%f)\n",
321 <                        vs->sloc[0], vs->sloc[1], vs->sloc[2]);
321 >                source[sn].sloc[0], source[sn].sloc[1], source[sn].sloc[2]);
322          fprintf(fp, "\tlinked to source %d (%s)\n",
323 <                        vs->sa.svnext, source[vs->sa.svnext].so->oname);
324 <        if (vs->sflags & SFOLLOW)
323 >                source[sn].sa.svnext, source[source[sn].sa.svnext].so->oname);
324 >        if (source[sn].sflags & SFOLLOW)
325                  fprintf(fp, "\talways followed\n");
326          else
327                  fprintf(fp, "\tnever followed\n");
328 <        if (!(vs->sflags & SSPOT))
328 >        if (!(source[sn].sflags & SSPOT))
329                  return;
330          fprintf(fp, "\twith spot aim (%f,%f,%f) and size %f\n",
331 <                        vs->sl.s->aim[0], vs->sl.s->aim[1], vs->sl.s->aim[2],
332 <                        vs->sl.s->siz);
331 >                        source[sn].sl.s->aim[0], source[sn].sl.s->aim[1],
332 >                        source[sn].sl.s->aim[2], source[sn].sl.s->siz);
333   }
334   #endif

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines