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

Comparing ray/src/rt/source.c (file contents):
Revision 2.50 by greg, Tue Apr 19 01:15:07 2005 UTC vs.
Revision 2.59 by greg, Thu Oct 28 09:38:46 2010 UTC

# Line 55 | Line 55 | findmaterial(register OBJREC *o)
55                          ao = objptr(aobj);
56                          if (ismaterial(ao->otype))
57                                  return(ao);
58 +                        if (ao->otype == MOD_ALIAS) {
59 +                                o = ao;
60 +                                continue;
61 +                        }
62                  }
63                  if (o->omod == OVOID)
64                          return(NULL);
# Line 82 | Line 86 | marksources(void)                      /* find and mark source objects */
86                          continue;
87                                          /* find material */
88                  m = findmaterial(objptr(o->omod));
89 <                if (m == NULL || !islight(m->otype))
89 >                if (m == NULL)
90 >                        continue;
91 >                if (m->otype == MAT_CLIP) {
92 >                        markclip(m);    /* special case for antimatter */
93 >                        continue;
94 >                }
95 >                if (!islight(m->otype))
96                          continue;       /* not source modifier */
97          
98                  if (m->oargs.nfargs != (m->otype == MAT_GLOW ? 4 :
# Line 174 | Line 184 | srcray(                                /* send a ray to a source, return domega */
184          SRCINDEX  *si                   /* source sample index */
185   )
186   {
187 <    double  d;                          /* distance to source */
188 <    register SRCREC  *srcp;
187 >        double  d;                              /* distance to source */
188 >        register SRCREC  *srcp;
189  
190 <    rayorigin(sr, SHADOW, r, NULL);             /* ignore limits */
190 >        rayorigin(sr, SHADOW, r, NULL);         /* ignore limits */
191  
192 <    while ((d = nextssamp(sr, si)) != 0.0) {
193 <        sr->rsrc = si->sn;                      /* remember source */
194 <        srcp = source + si->sn;
195 <        if (srcp->sflags & SDISTANT) {
196 <                if (srcp->sflags & SSPOT && spotout(sr, srcp->sl.s))
197 <                        continue;
198 <                return(1);              /* sample OK */
199 <        }
192 >        if (r == NULL)
193 >                sr->rmax = 0.0;
194 >
195 >        while ((d = nextssamp(sr, si)) != 0.0) {
196 >                sr->rsrc = si->sn;                      /* remember source */
197 >                srcp = source + si->sn;
198 >                if (srcp->sflags & SDISTANT) {
199 >                        if (srcp->sflags & SSPOT && spotout(sr, srcp->sl.s))
200 >                                continue;
201 >                        return(1);              /* sample OK */
202 >                }
203                                  /* local source */
204                                                  /* check proximity */
205 <        if (srcp->sflags & SPROX && d > srcp->sl.prox)
193 <                continue;
194 <                                                /* check angle */
195 <        if (srcp->sflags & SSPOT) {
196 <                if (spotout(sr, srcp->sl.s))
205 >                if (srcp->sflags & SPROX && d > srcp->sl.prox)
206                          continue;
207 +                                                /* check angle */
208 +                if (srcp->sflags & SSPOT) {
209 +                        if (spotout(sr, srcp->sl.s))
210 +                                continue;
211                                          /* adjust solid angle */
212 <                si->dom *= d*d;
213 <                d += srcp->sl.s->flen;
214 <                si->dom /= d*d;
212 >                        si->dom *= d*d;
213 >                        d += srcp->sl.s->flen;
214 >                        si->dom /= d*d;
215 >                }
216 >                return(1);                      /* sample OK */
217          }
218 <        return(1);                      /* sample OK */
204 <    }
205 <    return(0);                  /* no more samples */
218 >        return(0);                      /* no more samples */
219   }
220  
221  
# Line 296 | Line 309 | sourcehit(                     /* check to see if ray hit distant source
309                   * If it's a glow or transparent illum, just remember it.
310                   */
311                  if (source[i].sflags & SSKIP) {
312 <                        glowsrc = i;
312 >                        if (glowsrc < 0)
313 >                                glowsrc = i;
314                          continue;
315                  }
316                  if (transillum(source[i].so->omod)) {
317 <                        transrc = i;
317 >                        if (transrc < 0)
318 >                                transrc = i;
319                          continue;
320                  }
321                  r->ro = source[i].so;   /* otherwise, use first hit */
# Line 382 | Line 397 | direct(                                        /* add direct component */
397                  scp->sno = sr.rsrc;
398                                                  /* compute coefficient */
399                  (*f)(scp->coef, p, sr.rdir, si.dom);
400 <                cntord[sn].brt = bright(scp->coef);
400 >                cntord[sn].brt = intens(scp->coef);
401                  if (cntord[sn].brt <= 0.0)
402                          continue;
403   #if SHADCACHE
# Line 395 | Line 410 | direct(                                        /* add direct component */
410                  VCOPY(scp->dir, sr.rdir);
411                  copycolor(sr.rcoef, scp->coef);
412                                                  /* compute potential */
413 <                srcvalue(&sr);
413 >                sr.revf = srcvalue;
414 >                rayvalue(&sr);
415                  multcolor(sr.rcol, sr.rcoef);
416                  copycolor(scp->val, sr.rcol);
417 <                cntord[sn].brt = bright(scp->val);
417 >                cntord[sn].brt = bright(sr.rcol);
418          }
419                                                  /* sort contributions */
420          qsort(cntord, sn, sizeof(CNTPTR), cntcmp);
# Line 437 | Line 453 | direct(                                        /* add direct component */
453                  scp = srccnt + cntord[sn].sndx;
454                                                  /* test for hit */
455                  rayorigin(&sr, SHADOW, r, NULL);
456 +                copycolor(sr.rcoef, scp->coef);
457                  VCOPY(sr.rdir, scp->dir);
458                  sr.rsrc = scp->sno;
459                                                  /* keep statistics */
# Line 451 | Line 468 | direct(                                        /* add direct component */
468                          raycont(&sr);
469                          if (trace != NULL)
470                                  (*trace)(&sr);  /* trace execution */
454 –                        rayparticipate(&sr);
471                          if (bright(sr.rcol) <= FTINY) {
472   #if SHADCACHE
473                                  if ((scp <= srccnt || scp[-1].sno != scp->sno)
# Line 461 | Line 477 | direct(                                        /* add direct component */
477   #endif
478                                  continue;       /* missed! */
479                          }
480 +                        rayparticipate(&sr);
481 +                        multcolor(sr.rcol, sr.rcoef);
482                          copycolor(scp->val, sr.rcol);
483 <                        multcolor(scp->val, scp->coef);
483 >                } else if (trace != NULL &&
484 >                        (source[scp->sno].sflags & (SDISTANT|SVIRTUAL|SFOLLOW))
485 >                                                == (SDISTANT|SFOLLOW) &&
486 >                                sourcehit(&sr) && rayshade(&sr, sr.ro->omod)) {
487 >                        (*trace)(&sr);          /* trace execution */
488 >                        /* skip call to rayparticipate() & scp->val update */
489                  }
490                                                  /* add contribution if hit */
491                  addcolor(r->rcol, scp->val);
# Line 484 | Line 507 | direct(                                        /* add direct component */
507                  scp = srccnt + cntord[sn].sndx;
508                  prob = hwt * (double)source[scp->sno].nhits /
509                                  (double)source[scp->sno].ntests;
510 <                if (prob > 1.0)
511 <                        prob = 1.0;
489 <                scalecolor(scp->val, prob);
510 >                if (prob < 1.0)
511 >                        scalecolor(scp->val, prob);
512                  addcolor(r->rcol, scp->val);
513          }
514   }
# Line 533 | Line 555 | srcscatter(                    /* compute source scattering into ray */
555                          sr.rorg[0] = r->rorg[0] + r->rdir[0]*t;
556                          sr.rorg[1] = r->rorg[1] + r->rdir[1]*t;
557                          sr.rorg[2] = r->rorg[2] + r->rdir[2]*t;
536 –                        sr.rmax = 0.;
558                          initsrcindex(&si);      /* sample ray to this source */
559                          si.sn = r->slights[i];
560                          nopart(&si, &sr);

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines