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 2.19 by greg, Tue Aug 16 18:09:53 2011 UTC vs.
Revision 2.20 by greg, Fri Sep 14 13:34:02 2012 UTC

# Line 21 | Line 21 | static const char      RCSid[] = "$Id$";
21   #define  MINSAMPLES     16              /* minimum number of pretest samples */
22   #define  STESTMAX       32              /* maximum seeks per sample */
23  
24 + #define FEQ(a,b)        ((a)-(b)+FTINY >= 0 && (b)-(a)+FTINY >= 0)
25  
26 +
27   static OBJECT  *vobject;                /* virtual source objects */
28   static int  nvobjects = 0;              /* number of virtual source objects */
29  
30  
31 < extern void
31 > static int
32 > isident4(MAT4 m)
33 > {
34 >        int     i, j;
35 >
36 >        for (i = 4; i--; )
37 >                for (j = 4; j--; )
38 >                        if (!FEQ(m[i][j], i==j))
39 >                                return(0);
40 >        return(1);
41 > }
42 >
43 >
44 > void
45   markvirtuals(void)                      /* find and mark virtual sources */
46   {
47 <        register OBJREC  *o;
48 <        register int  i;
47 >        OBJREC  *o;
48 >        int  i;
49                                          /* check number of direct relays */
50          if (directrelay <= 0)
51                  return;
# Line 69 | Line 84 | markvirtuals(void)                     /* find and mark virtual sources
84   }
85  
86  
87 < extern void
87 > void
88   addvirtuals(            /* add virtuals associated with source */
89          int  sn,
90          int  nr
91   )
92   {
93 <        register int  i;
93 >        int  i;
94                                  /* check relay limit first */
95          if (nr <= 0)
96                  return;
# Line 88 | Line 103 | addvirtuals(           /* add virtuals associated with source *
103   }
104  
105  
106 < extern void
106 > void
107   vproject(               /* create projected source(s) if they exist */
108          OBJREC  *o,
109          int  sn,
110          int  n
111   )
112   {
113 <        register int  i;
114 <        register VSMATERIAL  *vsmat;
113 >        int  i;
114 >        VSMATERIAL  *vsmat;
115          MAT4  proj;
116          int  ns;
117  
# Line 117 | Line 132 | vproject(              /* create projected source(s) if they exist
132   }
133  
134  
135 < extern OBJREC *
135 > OBJREC *
136   vsmaterial(                     /* get virtual source material pointer */
137          OBJREC  *o
138   )
139   {
140 <        register int  i;
141 <        register OBJREC  *m;
140 >        int  i;
141 >        OBJREC  *m;
142  
143          i = o->omod;
144          m = findmaterial(objptr(i));
# Line 137 | Line 152 | vsmaterial(                    /* get virtual source material pointer *
152   }
153  
154  
155 < extern int
155 > int
156   makevsrc(               /* make virtual source if reasonable */
157          OBJREC  *op,
158 <        register int  sn,
158 >        int  sn,
159          MAT4  pm
160   )
161   {
# Line 148 | Line 163 | makevsrc(              /* make virtual source if reasonable */
163          double  maxrad2, d;
164          int  nsflags;
165          SPOT  theirspot, ourspot;
166 <        register int  i;
167 <
166 >        int  i;
167 >                                        /* check for no-op */
168 >        if (isident4(pm))
169 >                return(0);
170          nsflags = source[sn].sflags | (SVIRTUAL|SSPOT|SFOLLOW);
171                                          /* get object center and max. radius */
172          maxrad2 = getdisk(ocent, op, sn);
# Line 263 | Line 280 | memerr:
280   }
281  
282  
283 < extern double
283 > double
284   getdisk(                /* get visible object disk */
285          FVECT  oc,
286          OBJREC  *op,
287 <        register int  sn
287 >        int  sn
288   )
289   {
290          double  rad2, roffs, offs, d, rd, rdoto;
# Line 296 | Line 313 | getdisk(               /* get visible object disk */
313   }
314  
315  
316 < extern int
316 > int
317   vstestvis(              /* pretest source visibility */
318          int  f,                 /* virtual source flags */
319          OBJREC  *o,             /* relay object */
320          FVECT  oc,              /* relay object center */
321          double  or2,            /* relay object radius squared */
322 <        register int  sn        /* target source number */
322 >        int  sn /* target source number */
323   )
324   {
325          RAY  sr;
# Line 312 | Line 329 | vstestvis(             /* pretest source visibility */
329          double  or, d, d1;
330          int  stestlim, ssn;
331          int  nhit, nok;
332 <        register int  i, n;
332 >        int  i, n;
333                                  /* return if pretesting disabled */
334          if (vspretest <= 0)
335                  return(f);
# Line 410 | Line 427 | vstestvis(             /* pretest source visibility */
427          
428  
429   #ifdef DEBUG
430 < extern void
430 > void
431   virtverb(       /* print verbose description of virtual source */
432 <        register int  sn,
432 >        int  sn,
433          FILE  *fp
434   )
435   {

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines