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

Comparing ray/src/rt/srcsupp.c (file contents):
Revision 2.5 by greg, Fri Aug 28 14:07:45 1992 UTC vs.
Revision 2.11 by greg, Wed Apr 23 00:52:34 2003 UTC

# Line 1 | Line 1
1 /* Copyright (c) 1992 Regents of the University of California */
2
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ LBL";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6
4   /*
5   *  Support routines for source objects and materials
6 + *
7 + *  External symbols declared in source.h
8   */
9  
10 + #include "copyright.h"
11 +
12   #include  "ray.h"
13  
14   #include  "otypes.h"
# Line 26 | Line 27 | int  nsources = 0;                     /* the number of sources */
27   SRCFUNC  sfun[NUMOTYPE];                /* source dispatch table */
28  
29  
30 + void
31   initstypes()                    /* initialize source dispatch table */
32   {
33          extern VSMATERIAL  mirror_vs, direct1_vs, direct2_vs;
32        extern int  fsetsrc(), ssetsrc(), sphsetsrc(), cylsetsrc(), rsetsrc();
33        extern int  nopart(), flatpart(), cylpart();
34        extern double  fgetplaneq(), rgetplaneq();
35        extern double  fgetmaxdisk(), rgetmaxdisk();
34          static SOBJECT  fsobj = {fsetsrc, flatpart, fgetplaneq, fgetmaxdisk};
35          static SOBJECT  ssobj = {ssetsrc, nopart};
36          static SOBJECT  sphsobj = {sphsetsrc, nopart};
# Line 56 | Line 54 | newsource()                    /* allocate new source in our array */
54          if (nsources == 0)
55                  source = (SRCREC *)malloc(SRCINC*sizeof(SRCREC));
56          else if (nsources%SRCINC == 0)
57 <                source = (SRCREC *)realloc((char *)source,
57 >                source = (SRCREC *)realloc((void *)source,
58                                  (unsigned)(nsources+SRCINC)*sizeof(SRCREC));
59          if (source == NULL)
60                  return(-1);
# Line 67 | Line 65 | newsource()                    /* allocate new source in our array */
65   }
66  
67  
68 + void
69   setflatss(src)                          /* set sampling for a flat source */
70   register SRCREC  *src;
71   {
# Line 86 | Line 85 | register SRCREC  *src;
85   }
86  
87  
88 + void
89   fsetsrc(src, so)                        /* set a face as a source */
90   register SRCREC  *src;
91   OBJREC  *so;
# Line 119 | Line 119 | OBJREC  *so;
119          }
120          src->srad = sqrt(src->srad);
121                                                  /* compute size vectors */
122 <        if (f->nv == 4 || (f->nv == 5 &&        /* parallelogram case */
123 <                        dist2(VERTEX(f,0),VERTEX(f,4)) <= FTINY*FTINY))
122 >        if (f->nv == 4)                         /* parallelogram case */
123                  for (j = 0; j < 3; j++) {
124                          src->ss[SU][j] = .5*(VERTEX(f,1)[j]-VERTEX(f,0)[j]);
125                          src->ss[SV][j] = .5*(VERTEX(f,3)[j]-VERTEX(f,0)[j]);
# Line 130 | Line 129 | OBJREC  *so;
129   }
130  
131  
132 + void
133   ssetsrc(src, so)                        /* set a source as a source */
134   register SRCREC  *src;
135   register OBJREC  *so;
# Line 156 | Line 156 | register OBJREC  *so;
156   }
157  
158  
159 + void
160   sphsetsrc(src, so)                      /* set a sphere as a source */
161   register SRCREC  *src;
162   register OBJREC  *so;
# Line 178 | Line 179 | register OBJREC  *so;
179   }
180  
181  
182 + void
183   rsetsrc(src, so)                        /* set a ring (disk) as a source */
184   register SRCREC  *src;
185   OBJREC  *so;
# Line 199 | Line 201 | OBJREC  *so;
201   }
202  
203  
204 + void
205   cylsetsrc(src, so)                      /* set a cylinder as a source */
206   register SRCREC  *src;
207   OBJREC  *so;
# Line 252 | Line 255 | register OBJREC  *m;
255   }
256  
257  
258 < spotout(r, s, dist)             /* check if we're outside spot region */
258 > int
259 > spotout(r, s)                   /* check if we're outside spot region */
260   register RAY  *r;
261   register SPOT  *s;
258 int  dist;
262   {
263          double  d;
264          FVECT  vd;
265          
266          if (s == NULL)
267                  return(0);
268 <        if (dist) {                     /* distant source */
268 >        if (s->flen < -FTINY) {         /* distant source */
269                  vd[0] = s->aim[0] - r->rorg[0];
270                  vd[1] = s->aim[1] - r->rorg[1];
271                  vd[2] = s->aim[2] - r->rorg[2];
# Line 353 | Line 356 | OBJREC  *op;
356   }
357  
358  
359 + int
360   commonspot(sp1, sp2, org)       /* set sp1 to intersection of sp1 and sp2 */
361   register SPOT  *sp1, *sp2;
362   FVECT  org;
# Line 376 | Line 380 | FVECT  org;
380   }
381  
382  
383 + int
384   commonbeam(sp1, sp2, dir)       /* set sp1 to intersection of sp1 and sp2 */
385   register SPOT  *sp1, *sp2;
386   FVECT  dir;
# Line 400 | Line 405 | FVECT  dir;
405   }
406  
407  
408 + int
409   checkspot(sp, nrm)              /* check spotlight for behind source */
410   register SPOT  *sp;     /* spotlight */
411   FVECT  nrm;             /* source surface normal */
# Line 495 | Line 501 | double  r1s, r2s;              /* radii squared */
501          for (i = 0; i < 3; i++)
502                  cc[i] = c1[i] + l*disp[i];
503          return(a2);
498 }
499
500
501 sourcehit(r)                    /* check to see if ray hit distant source */
502 register RAY  *r;
503 {
504        int  first, last;
505        register int  i;
506
507        if (r->rsrc >= 0) {             /* check only one if aimed */
508                first = last = r->rsrc;
509        } else {                        /* otherwise check all */
510                first = 0; last = nsources-1;
511        }
512        for (i = first; i <= last; i++)
513                if ((source[i].sflags & (SDISTANT|SVIRTUAL)) == SDISTANT)
514                        /*
515                         * Check to see if ray is within
516                         * solid angle of source.
517                         */
518                        if (2.0*PI * (1.0 - DOT(source[i].sloc,r->rdir))
519                                        <= source[i].ss2) {
520                                r->ro = source[i].so;
521                                if (!(source[i].sflags & SSKIP))
522                                        break;
523                        }
524
525        if (r->ro != NULL) {
526                for (i = 0; i < 3; i++)
527                        r->ron[i] = -r->rdir[i];
528                r->rod = 1.0;
529                r->rox = NULL;
530                return(1);
531        }
532        return(0);
504   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines