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

Comparing ray/src/rt/ambient.c (file contents):
Revision 2.7 by greg, Thu Jul 16 13:37:12 1992 UTC vs.
Revision 2.12 by greg, Mon Sep 21 12:07:31 1992 UTC

# Line 1 | Line 1
1 < /* Copyright (c) 1991 Regents of the University of California */
1 > /* Copyright (c) 1992 Regents of the University of California */
2  
3   #ifndef lint
4   static char SCCSid[] = "$SunId$ LBL";
# Line 7 | Line 7 | static char SCCSid[] = "$SunId$ LBL";
7   /*
8   *  ambient.c - routines dealing with ambient (inter-reflected) component.
9   *
10 – *  The macro AMBFLUSH (if defined) is the number of ambient values
11 – *      to wait before flushing to the ambient file.
12 – *
10   *     5/9/86
11   */
12  
# Line 23 | Line 20 | static char SCCSid[] = "$SunId$ LBL";
20  
21   #include  "random.h"
22  
23 < #define  OCTSCALE       0.5     /* ceil((valid rad.)/(cube size)) */
23 > #define  OCTSCALE       0.5     /* ceil((valid rad.)/(cube size)) */
24  
25   typedef struct ambtree {
26 <        AMBVAL  *alist;         /* ambient value list */
27 <        struct ambtree  *kid;   /* 8 child nodes */
26 >        AMBVAL  *alist;         /* ambient value list */
27 >        struct ambtree  *kid;   /* 8 child nodes */
28   }  AMBTREE;                     /* ambient octree */
29  
30   extern CUBE  thescene;          /* contains space boundaries */
31  
32 < #define  MAXASET        511     /* maximum number of elements in ambient set */
33 < OBJECT  ambset[MAXASET+1]={0};  /* ambient include/exclude set */
32 > #define  MAXASET        511     /* maximum number of elements in ambient set */
33 > OBJECT  ambset[MAXASET+1]={0};  /* ambient include/exclude set */
34  
35 < double  maxarad;                /* maximum ambient radius */
36 < double  minarad;                /* minimum ambient radius */
35 > double  maxarad;                /* maximum ambient radius */
36 > double  minarad;                /* minimum ambient radius */
37  
38 < static AMBTREE  atrunk;         /* our ambient trunk node */
38 > static AMBTREE  atrunk;         /* our ambient trunk node */
39  
40   static FILE  *ambfp = NULL;     /* ambient file pointer */
44 – static char  *afname;           /* ambient file name */
45 – static long  ambheadlen;        /* length of ambient file header */
41  
42 < #define  AMBFLUSH       (BUFSIZ/AMBVALSIZ)
42 > #define  AMBFLUSH       (BUFSIZ/AMBVALSIZ)
43  
44 < #define  newambval()    (AMBVAL *)bmalloc(sizeof(AMBVAL))
44 > #define  newambval()    (AMBVAL *)bmalloc(sizeof(AMBVAL))
45  
46 < #define  newambtree()   (AMBTREE *)calloc(8, sizeof(AMBTREE))
46 > #define  newambtree()   (AMBTREE *)calloc(8, sizeof(AMBTREE))
47  
48   extern long  ftell(), lseek();
49 + static int  initambfile(), avsave(), avinsert(), ambsync();
50  
51  
52   setambres(ar)                           /* set ambient resolution */
# Line 74 | Line 70 | int  ar;
70   setambient(afile)                       /* initialize calculation */
71   char  *afile;
72   {
73 <        AMBVAL  amb;
73 >        long  headlen;
74 >        AMBVAL  amb;
75                                                  /* init ambient limits */
76          setambres(ambres);
77                                                  /* open ambient file */
78 <        if ((afname = afile) != NULL)
78 >        if (afile != NULL)
79                  if ((ambfp = fopen(afile, "r+")) != NULL) {
80                          initambfile(0);
81 +                        headlen = ftell(ambfp);
82                          while (readambval(&amb, ambfp))
83                                  avinsert(&amb, &atrunk, thescene.cuorg,
84                                                  thescene.cusize);
85                                                          /* align */
86 <                        fseek(ambfp, -((ftell(ambfp)-ambheadlen)%AMBVALSIZ), 1);
87 <                } else if ((ambfp = fopen(afile, "w")) != NULL)
86 >                        fseek(ambfp, -((ftell(ambfp)-headlen)%AMBVALSIZ), 1);
87 >                } else if ((ambfp = fopen(afile, "w+")) != NULL)
88                          initambfile(1);
89                  else {
90                          sprintf(errmsg, "cannot open ambient file \"%s\"",
# Line 96 | Line 94 | char  *afile;
94   }
95  
96  
99 – initambfile(creat)              /* initialize ambient file */
100 – int  creat;
101 – {
102 –        extern char  *progname, *octname;
103 –
104 –        setbuf(ambfp, bmalloc(BUFSIZ));
105 –        if (creat) {                    /* new file */
106 –                fprintf(ambfp, "%s -av %g %g %g -ab %d -aa %g ",
107 –                                progname, colval(ambval,RED),
108 –                                colval(ambval,GRN), colval(ambval,BLU),
109 –                                ambounce, ambacc);
110 –                fprintf(ambfp, "-ad %d -as %d -ar %d %s\n",
111 –                                ambdiv, ambssamp, ambres,
112 –                                octname==NULL ? "" : octname);
113 –                fputformat(AMBFMT, ambfp);
114 –                putc('\n', ambfp);
115 –                putambmagic(ambfp);
116 –                fflush(ambfp);
117 –        } else if (checkheader(ambfp, AMBFMT, NULL) < 0
118 –                        || !hasambmagic(ambfp)) {
119 –                sprintf(errmsg, "\"%s\" is not an ambient file", afname);
120 –                error(USER, afname);
121 –        }
122 –        ambheadlen = ftell(ambfp);
123 – }
124 –
125 –
97   ambnotify(obj)                  /* record new modifier */
98 < OBJECT  obj;
98 > OBJECT  obj;
99   {
100          static int  hitlimit = 0;
101 <        register OBJREC  *o = objptr(obj);
101 >        register OBJREC  *o = objptr(obj);
102          register char  **amblp;
103  
104          if (hitlimit || !ismodifier(o->otype))
# Line 150 | Line 121 | COLOR  acol;
121   register RAY  *r;
122   {
123          static int  rdepth = 0;                 /* ambient recursion */
124 <        double  d;
124 >        double  d;
125  
126          if (ambdiv <= 0)                        /* no ambient calculation */
127                  goto dumbamb;
# Line 192 | Line 163 | sumambient(acol, r, al, at, c0, s)     /* get interpolated
163   COLOR  acol;
164   register RAY  *r;
165   int  al;
166 < AMBTREE  *at;
166 > AMBTREE  *at;
167   FVECT  c0;
168 < double  s;
168 > double  s;
169   {
170          extern double  sqrt();
171 <        double  d, e1, e2, wt, wsum;
171 >        double  d, e1, e2, wt, wsum;
172          COLOR  ct;
173          FVECT  ck0;
174          int  i;
175          register int  j;
176 <        register AMBVAL  *av;
176 >        register AMBVAL  *av;
177                                          /* do this node */
178          wsum = 0.0;
179          for (av = at->alist; av != NULL; av = av->next) {
# Line 281 | Line 252 | COLOR  acol;
252   register RAY  *r;
253   int  al;
254   {
255 <        AMBVAL  amb;
255 >        AMBVAL  amb;
256          FVECT   gp, gd;
257                                                  /* compute weight */
258          amb.weight = pow(AVGREFL, (double)al);
# Line 306 | Line 277 | int  al;
277  
278   extambient(cr, ap, pv, nv)              /* extrapolate value at pv, nv */
279   COLOR  cr;
280 < register AMBVAL  *ap;
280 > register AMBVAL  *ap;
281   FVECT  pv, nv;
282   {
283          FVECT  v1, v2;
284          register int  i;
285 <        double  d;
285 >        double  d;
286  
287          d = 1.0;                        /* zeroeth order */
288                                          /* gradient due to translation */
# Line 331 | Line 302 | FVECT  pv, nv;
302  
303  
304   static
305 + initambfile(creat)              /* initialize ambient file */
306 + int  creat;
307 + {
308 +        extern char  *progname, *octname, VersionID[];
309 +
310 + #ifdef MSDOS
311 +        setmode(fileno(ambfp), O_BINARY);
312 + #endif
313 +        setbuf(ambfp, bmalloc(BUFSIZ));
314 +        if (creat) {                    /* new file */
315 +                fprintf(ambfp, "%s -av %g %g %g -ab %d -aa %g ",
316 +                                progname, colval(ambval,RED),
317 +                                colval(ambval,GRN), colval(ambval,BLU),
318 +                                ambounce, ambacc);
319 +                fprintf(ambfp, "-ad %d -as %d -ar %d %s\n",
320 +                                ambdiv, ambssamp, ambres,
321 +                                octname==NULL ? "" : octname);
322 +                fprintf(ambfp, "SOFTWARE= %s\n", VersionID);
323 +                fputformat(AMBFMT, ambfp);
324 +                putc('\n', ambfp);
325 +                putambmagic(ambfp);
326 +                fflush(ambfp);
327 + #ifndef  NIX
328 +                sync();                 /* protect against NFS buffering */
329 + #endif
330 +        } else if (checkheader(ambfp, AMBFMT, NULL) < 0 || !hasambmagic(ambfp))
331 +                error(USER, "bad ambient file");
332 + }
333 +
334 +
335 + static
336   avsave(av)                              /* insert and save an ambient value */
337 < AMBVAL  *av;
337 > AMBVAL  *av;
338   {
339          static int  nunflshed = 0;
340  
# Line 354 | Line 356 | writerr:
356  
357   static
358   avinsert(aval, at, c0, s)               /* insert ambient value in a tree */
359 < AMBVAL  *aval;
359 > AMBVAL  *aval;
360   register AMBTREE  *at;
361   FVECT  c0;
362 < double  s;
362 > double  s;
363   {
364          FVECT  ck0;
365          int  branch;
366 <        register AMBVAL  *av;
366 >        register AMBVAL  *av;
367          register int  i;
368  
369          if ((av = newambval()) == NULL)
# Line 389 | Line 391 | memerr:
391   }
392  
393  
394 < #include  <fcntl.h>
394 > #ifdef  NIX
395  
396 + static
397 + ambsync()                       /* flush ambient file */
398 + {
399 +        return(fflush(ambfp));
400 + }
401  
402 + #else
403 +
404 + #include  <fcntl.h>
405 + #include  <sys/types.h>
406 + #include  <sys/stat.h>
407 +
408   static
409   ambsync()                       /* synchronize ambient file */
410   {
411          static FILE  *ambinp = NULL;
412          struct flock  fls;
413 <        AMBVAL  avs;
413 >        struct stat  sts;
414 >        AMBVAL  avs;
415          long  lastpos, flen;
416          register int  n;
417                                  /* gain exclusive access */
# Line 408 | Line 422 | ambsync()                      /* synchronize ambient file */
422          if (fcntl(fileno(ambfp), F_SETLKW, &fls) < 0)
423                  error(SYSTEM, "cannot lock ambient file");
424                                  /* see if file has grown */
425 <        lastpos = lseek(fileno(ambfp), 0L, 2);  /* may move pointer */
426 <        flen = lseek(fileno(ambfp), 0L, 1);     /* new(?) file length */
427 <        if (n = (flen - lastpos)%AMBVALSIZ)     /* assure alignment */
428 <                lseek(fileno(ambfp), flen -= n, 0);
425 >        lastpos = lseek(fileno(ambfp), 0L, 1);  /* get previous position */
426 >        if (fstat(fileno(ambfp), &sts) < 0)     /* get current length */
427 >                error(SYSTEM, "cannot stat ambient file");
428 >        flen = sts.st_size;
429          if (n = (flen - lastpos)/AMBVALSIZ) {   /* file has grown */
430 <                if (ambinp == NULL && (ambinp = fopen(afname, "r")) == NULL)
431 <                        error(SYSTEM, "cannot reopen ambient file");
418 <                fseek(ambinp, lastpos, 0);      /* go to previous position */
430 >                if (ambinp == NULL)             /* use duplicate file */
431 >                        ambinp = fdopen(dup(fileno(ambfp)), "r");
432                  while (n--) {                   /* load contributed values */
433                          readambval(&avs, ambinp);
434                          avinsert(&avs,&atrunk,thescene.cuorg,thescene.cusize);
435 <                }
435 >                }                               /* moves shared file pointer */
436 >                if (n = (flen - lastpos)%AMBVALSIZ)     /* alignment */
437 >                        lseek(fileno(ambfp), flen-n, 0);
438          }
439          n = fflush(ambfp);                      /* calls write() at last */
440          fls.l_type = F_UNLCK;                   /* release file */
441          fcntl(fileno(ambfp), F_SETLKW, &fls);
442          return(n);
443   }
444 +
445 + #endif

Diff Legend

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