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

Comparing ray/src/rt/pmapsrc.c (file contents):
Revision 2.1 by greg, Tue Feb 24 19:39:27 2015 UTC vs.
Revision 2.8 by greg, Tue Sep 1 16:27:53 2015 UTC

# Line 1 | Line 1
1 + #ifndef lint
2 + static const char RCSid[] = "$Id$";
3 + #endif
4   /*
5     ==================================================================
6     Photon map support routines for emission from light sources
7  
8     Roland Schregle (roland.schregle@{hslu.ch, gmail.com})
9     (c) Fraunhofer Institute for Solar Energy Systems,
10 <       Lucerne University of Applied Sciences & Arts  
10 >   (c) Lucerne University of Applied Sciences and Arts,
11 >   supported by the Swiss National Science Foundation (SNSF, #147053)
12     ==================================================================
13 <  
10 <   $Id$
13 >
14   */
15  
16  
# Line 46 | Line 49 | static int flatPhotonPartition2 (EmissionMap* emap, un
49        if (!emap -> partitions)
50           error(USER, "can't allocate source partitions");
51  
52 <      bzero(emap -> partitions + (emap -> maxPartitions >> 2),
53 <            emap -> maxPartitions >> 2);
52 >      memset(emap -> partitions + (emap -> maxPartitions >> 2), 0,
53 >             emap -> maxPartitions >> 2);
54     }
55    
56     if (du2 * dv2 <= 1) {                /* hit limit? */
# Line 100 | Line 103 | static void flatPhotonPartition (EmissionMap* emap)
103     RREAL    *vp;
104     double   du2, dv2;
105  
106 <   bzero(emap -> partitions, emap -> maxPartitions >> 1);
106 >   memset(emap -> partitions, 0, emap -> maxPartitions >> 1);
107     emap -> partArea = srcsizerat * thescene.cusize;
108     emap -> partArea *= emap -> partArea;
109     vp = emap -> src -> ss [SU];
# Line 132 | Line 135 | static void sourcePhotonPartition (EmissionMap* emap)
135    
136     else {
137        /* No photon ports defined, so partition scene cube faces */
138 <      bzero(emap -> partitions, emap -> maxPartitions >> 1);
138 >      memset(emap -> partitions, 0, emap -> maxPartitions >> 1);
139        setpart(emap -> partitions, 0, S0);
140        emap -> partitionCnt = 0;
141        emap -> numPartitions = 1 / srcsizerat;
# Line 150 | Line 153 | static void spherePhotonPartition (EmissionMap* emap)
153   {
154     unsigned numTheta, numPhi;
155    
156 <   bzero(emap -> partitions, emap -> maxPartitions >> 1);
156 >   memset(emap -> partitions, 0, emap -> maxPartitions >> 1);
157     setpart(emap -> partitions, 0, S0);
158     emap -> partArea = 4 * PI * sqr(emap -> src -> srad);
159     emap -> numPartitions = emap -> partArea /
# Line 180 | Line 183 | static int cylPhotonPartition2 (EmissionMap* emap, uns
183                                                     emap -> maxPartitions >> 1);
184        if (!emap -> partitions)
185           error(USER, "can't allocate source partitions");
186 <      bzero(emap -> partitions + (emap -> maxPartitions >> 2),
186 >        
187 >      memset(emap -> partitions + (emap -> maxPartitions >> 2), 0,
188              emap -> maxPartitions >> 2);
189     }
190    
# Line 222 | Line 226 | static void cylPhotonPartition (EmissionMap* emap)
226   {
227     double d2;
228  
229 <   bzero(emap -> partitions, emap -> maxPartitions >> 1);
229 >   memset(emap -> partitions, 0, emap -> maxPartitions >> 1);
230     d2 = srcsizerat * thescene.cusize;
231     d2 = PI * emap -> src -> ss2 / (2 * emap -> src -> srad * sqr(d2));
232     d2 *= d2 * DOT(emap -> src -> ss [SU], emap -> src -> ss [SU]);
# Line 590 | Line 594 | void initPhotonEmission (EmissionMap *emap, float numP
594     unsigned i, t, p;
595     double phi, cosTheta, sinTheta, du, dv, dOmega, thetaScale;
596     EmissionSample* sample;
597 <   const OBJREC* mod =  objptr(emap -> src -> so -> omod);
597 >   const OBJREC* mod =  findmaterial(emap -> src -> so);
598     static RAY r;
599   #if 0  
600     static double lastCosNorm = FHUGE;
601     static SRCREC *lastSrc = NULL, *lastPort = NULL;
602   #endif  
603  
604 +   setcolor(emap -> partFlux, 0, 0, 0);
605 +
606     photonOrigin [emap -> src -> so -> otype] (emap);
607     cosTheta = DOT(emap -> ws, emap -> wh);
608  
# Line 619 | Line 625 | void initPhotonEmission (EmissionMap *emap, float numP
625   #endif
626          
627     /* Need to recompute flux & PDF */
622 –   setcolor(emap -> partFlux, 0, 0, 0);
628     emap -> cdf = 0;
629     emap -> numSamples = 0;
630    
# Line 668 | Line 673 | void initPhotonEmission (EmissionMap *emap, float numP
673          
674        VCOPY(r.rorg, emap -> photonOrg);
675        VCOPY(r.rop, emap -> photonOrg);
676 <      r.rmax = FHUGE;
676 >      r.rmax = 0;
677        
678        for (t = 0; t < emap -> numTheta; t++) {
679           for (p = 0; p < emap -> numPhi; p++) {
# Line 738 | Line 743 | void emitPhoton (const EmissionMap* emap, RAY* ray)
743     unsigned long i, lo, hi;
744     const EmissionSample* sample = emap -> samples;
745     RREAL du, dv, cosTheta, cosThetaSqr, sinTheta, phi;  
746 <   const OBJREC* mod = objptr(emap -> src -> so -> omod);
746 >   const OBJREC* mod = findmaterial(emap -> src -> so);
747    
748     /* Choose a new origin within current partition for every
749        emitted photon to break up clustering artifacts */
750     photonOrigin [emap -> src -> so -> otype] ((EmissionMap*)emap);
751 +   /* If we have a local glow source with a maximum radius, then
752 +      restrict our photon to the specified distance (otherwise no limit) */
753 +   if (mod -> otype == MAT_GLOW && !(emap -> src -> sflags & SDISTANT)
754 +                && mod -> oargs.farg[3] > FTINY)
755 +      ray -> rmax = mod -> oargs.farg[3];
756 +   else
757 +      ray -> rmax = 0;
758     rayorigin(ray, PRIMARY, NULL, NULL);
747 –   ray -> rmax = FHUGE;
759    
760     if (!emap -> numSamples) {
761        /* Source is unmodified and has no port, and either local with

Diff Legend

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