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

Comparing ray/src/rt/srcdraw.c (file contents):
Revision 2.3 by greg, Sat Jun 8 21:16:04 1996 UTC vs.
Revision 2.9 by greg, Wed Oct 1 22:07:19 2003 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1996 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   * Draw small sources into image in case we missed them.
6 + *
7 + *  External symbols declared in ray.h
8   */
9  
10 + #include "copyright.h"
11 +
12   #include  "ray.h"
13  
14   #include  "view.h"
# Line 26 | Line 27 | typedef struct splist {
27          struct splist   *next;                  /* next source in list */
28          int     sn;                             /* source number */
29          short   nv;                             /* number of vertices */
30 <        FLOAT   vl[3][2];                       /* vertex array (last) */
30 >        RREAL   vl[3][2];                       /* vertex array (last) */
31   } SPLIST;                               /* source polygon list */
32  
33   extern VIEW     ourview;                /* our view parameters */
# Line 36 | Line 37 | static SPLIST  *sphead = NULL;         /* our list of source p
37  
38   static int
39   inregion(p, cv, crit)                   /* check if vertex is in region */
40 < FLOAT   p[2];
40 > RREAL   p[2];
41   double  cv;
42   int     crit;
43   {
# Line 56 | Line 57 | int    crit;
57  
58   static
59   clipregion(a, b, cv, crit, r)           /* find intersection with boundary */
60 < register FLOAT  a[2], b[2];
60 > register RREAL  a[2], b[2];
61   double  cv;
62   int     crit;
63 < FLOAT   r[2];   /* return value */
63 > RREAL   r[2];   /* return value */
64   {
65          switch (crit) {
66          case CLIP_ABOVE:
# Line 78 | Line 79 | FLOAT  r[2];   /* return value */
79  
80   static int
81   hp_clip_poly(vl, nv, cv, crit, vlo)     /* clip polygon to half-plane */
82 < FLOAT   vl[][2];
82 > RREAL   vl[][2];
83   int     nv;
84   double  cv;
85   int     crit;
86 < FLOAT   vlo[][2];       /* return value */
86 > RREAL   vlo[][2];       /* return value */
87   {
88 <        FLOAT   *s, *p;
88 >        RREAL   *s, *p;
89          register int    j, nvo;
90  
91          s = vl[nv-1];
# Line 105 | Line 106 | FLOAT  vlo[][2];       /* return value */
106  
107   static int
108   box_clip_poly(vl, nv, xl, xr, yb, ya, vlo)      /* clip polygon to box */
109 < FLOAT   vl[MAXVERT][2];
109 > RREAL   vl[MAXVERT][2];
110   int     nv;
111   double  xl, xr, yb, ya;
112 < FLOAT   vlo[MAXVERT][2];        /* return value */
112 > RREAL   vlo[MAXVERT][2];        /* return value */
113   {
114 <        FLOAT   vlt[MAXVERT][2];
114 >        RREAL   vlt[MAXVERT][2];
115          int     nvt, nvo;
116  
117          nvt = hp_clip_poly(vl, nv, yb, CLIP_BELOW, vlt);
# Line 124 | Line 125 | FLOAT  vlo[MAXVERT][2];        /* return value */
125  
126   static double
127   minw2(vl, nv, ar2)                      /* compute square of minimum width */
128 < FLOAT   vl[][2];
128 > RREAL   vl[][2];
129   int     nv;
130   double  ar2;
131   {
132          double  d2, w2, w2min, w2max;
133 <        register FLOAT  *p0, *p1, *p2;
133 >        register RREAL  *p0, *p1, *p2;
134          int     i, j;
135                                  /* find minimum for all widths */
136          w2min = FHUGE;
# Line 157 | Line 158 | double ar2;
158  
159   static
160   convex_center(vl, nv, cv)               /* compute center of convex polygon */
161 < register FLOAT  vl[][2];
161 > register RREAL  vl[][2];
162   int     nv;
163 < FLOAT   cv[2];          /* return value */
163 > RREAL   cv[2];          /* return value */
164   {
165          register int    i;
166                                          /* simple average (suboptimal) */
# Line 175 | Line 176 | FLOAT  cv[2];          /* return value */
176  
177   static double
178   poly_area(vl, nv)                       /* compute area of polygon */
179 < register FLOAT  vl[][2];
179 > register RREAL  vl[][2];
180   int     nv;
181   {
182          double  a;
183 <        FLOAT   v0[2], v1[2];
183 >        RREAL   v0[2], v1[2];
184          register int    i;
185  
186          a = 0.;
# Line 197 | Line 198 | int    nv;
198  
199   static int
200   convex_hull(vl, nv, vlo)                /* compute polygon's convex hull */
201 < FLOAT   vl[][2];
201 > RREAL   vl[][2];
202   int     nv;
203 < FLOAT   vlo[][2];       /* return value */
203 > RREAL   vlo[][2];       /* return value */
204   {
205          int     nvo, nvt;
206 <        FLOAT   vlt[MAXVERT][2];
206 >        RREAL   vlt[MAXVERT][2];
207          double  voa, vta;
208          register int    i, j;
209                                          /* start with original polygon */
# Line 234 | Line 235 | FLOAT  vlo[][2];       /* return value */
235   static
236   spinsert(sn, vl, nv)                    /* insert new source polygon */
237   int     sn;
238 < FLOAT   vl[][2];
238 > RREAL   vl[][2];
239   int     nv;
240   {
241          register SPLIST *spn;
# Line 243 | Line 244 | int    nv;
244          if (nv < 3)
245                  return;
246          if (nv > 3)
247 <                spn = (SPLIST *)malloc(sizeof(SPLIST)+sizeof(FLOAT)*2*(nv-3));
247 >                spn = (SPLIST *)malloc(sizeof(SPLIST)+sizeof(RREAL)*2*(nv-3));
248          else
249                  spn = (SPLIST *)malloc(sizeof(SPLIST));
250          if (spn == NULL)
# Line 260 | Line 261 | int    nv;
261   int
262   sourcepoly(sn, sp)                      /* compute image polygon for source */
263   int     sn;
264 < FLOAT   sp[MAXVERT][2];
264 > RREAL   sp[MAXVERT][2];
265   {
266 <        static char     cubeord[8][6] = {{1,3,2,6,4,5},{0,4,5,7,3,2},
266 >        static short    cubeord[8][6] = {{1,3,2,6,4,5},{0,4,5,7,3,2},
267                                           {0,1,3,7,6,4},{0,1,5,7,6,2},
268                                           {0,2,6,7,5,1},{0,4,6,7,3,1},
269                                           {0,2,3,7,5,4},{1,5,4,6,2,3}};
270          register SRCREC *s = source + sn;
271          FVECT   ap, ip;
272 <        FLOAT   pt[6][2];
272 >        RREAL   pt[6][2];
273          int     dir;
274          register int    i, j;
275  
# Line 284 | Line 285 | FLOAT  sp[MAXVERT][2];
285                  for (j = 0; j < 4; j++) {       /* four corners */
286                          for (i = 0; i < 3; i++) {
287                                  ap[i] = s->sloc[i];
288 <                                if (j==1|j==2) ap[i] += s->ss[SU][i];
288 >                                if ((j==1)|(j==2)) ap[i] += s->ss[SU][i];
289                                  else ap[i] -= s->ss[SU][i];
290 <                                if (j==2|j==3) ap[i] += s->ss[SV][i];
290 >                                if ((j==2)|(j==3)) ap[i] += s->ss[SV][i];
291                                  else ap[i] -= s->ss[SV][i];
292                                  if (s->sflags & SDISTANT) {
293                                          ap[i] *= 1. + ourview.vfore;
# Line 330 | Line 331 | FLOAT  sp[MAXVERT][2];
331   init_drawsources(rad)
332   int     rad;                            /* source sample size */
333   {
334 <        FLOAT   spoly[MAXVERT][2];
334 >        RREAL   spoly[MAXVERT][2];
335          int     nsv;
336          register SPLIST *sp;
337          register int    i;
338                                          /* free old source list if one */
339          for (sp = sphead; sp != NULL; sp = sphead) {
340                  sphead = sp->next;
341 <                free((char *)sp);
341 >                free((void *)sp);
342          }
343                                          /* loop through all sources */
344          for (i = nsources; i--; ) {
# Line 356 | Line 357 | int    rad;                            /* source sample size */
357          }
358   }
359  
360 <                        /* add sources smaller than rad to computed subimage */
360 > void                    /* add sources smaller than rad to computed subimage */
361   drawsources(pic, zbf, x0, xsiz, y0, ysiz)
362   COLOR   *pic[];                         /* subimage pixel value array */
363   float   *zbf[];                         /* subimage distance array (opt.) */
364   int     x0, xsiz, y0, ysiz;             /* origin and size of subimage */
365   {
366 <        FLOAT   spoly[MAXVERT][2], ppoly[MAXVERT][2];
366 >        RREAL   spoly[MAXVERT][2], ppoly[MAXVERT][2];
367          int     nsv, npv;
368          int     xmin, xmax, ymin, ymax, x, y;
369 <        FLOAT   cxy[2];
370 <        double  pa;
369 >        RREAL   cxy[2];
370 >        double  w;
371          RAY     sr;
372          register SPLIST *sp;
373          register int    i;
# Line 408 | Line 409 | int    x0, xsiz, y0, ysiz;             /* origin and size of subimag
409                                  if (source[sp->sn].sflags & SSPOT &&
410                                                  spotout(&sr, source[sp->sn].sl.s))
411                                          continue;       /* outside spot */
412 <                                rayorigin(&sr, NULL, SHADOW, 1.0);
413 <                                sr.rsrc = sp->sn;
412 >                                w = poly_area(ppoly, npv) * hres * vres;
413 >                                if (w < .95) {          /* subpixel source */
414 >                                        rayorigin(&sr, NULL, SHADOW, 1.0);
415 >                                        sr.rsrc = sp->sn;
416 >                                } else
417 >                                        rayorigin(&sr, NULL, PRIMARY, 1.0);
418                                  rayvalue(&sr);          /* compute value */
419                                  if (bright(sr.rcol) <= FTINY)
420                                          continue;       /* missed/blocked */
421                                                          /* modify pixel */
422                                  if (zbf[y-y0] != NULL &&
423 <                                                sr.rt < zbf[y-y0][x-x0])
423 >                                                sr.rt < 0.999*zbf[y-y0][x-x0])
424                                          zbf[y-y0][x-x0] = sr.rt;
425 <                                pa = poly_area(ppoly, npv);
426 <                                scalecolor(sr.rcol, pa*hres*vres);
427 <                                scalecolor(pic[y-y0][x-x0], (1.-pa*hres*vres));
425 >                                else if (!bigdiff(sr.rcol, pic[y-y0][x-x0],
426 >                                                0.01))  /* source sample */
427 >                                        setcolor(pic[y-y0][x-x0], 0., 0., 0.);
428 >                                scalecolor(sr.rcol, w);
429 >                                scalecolor(pic[y-y0][x-x0], 1.-w);
430                                  addcolor(pic[y-y0][x-x0], sr.rcol);
431                          }
432          }

Diff Legend

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