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

Comparing ray/src/rt/rtrace.c (file contents):
Revision 2.14 by greg, Tue May 4 13:26:52 1993 UTC vs.
Revision 2.73 by greg, Tue Nov 13 19:58:33 2018 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   *  rtrace.c - program and variables for individual ray tracing.
9 – *
10 – *     6/11/86
6   */
7  
8 + #include "copyright.h"
9 +
10   /*
11   *  Input is in the form:
12   *
# Line 24 | Line 21 | static char SCCSid[] = "$SunId$ LBL";
21   *  irradiance values are desired.
22   */
23  
24 < #include  "ray.h"
24 > #include  <time.h>
25  
26 < #include  "octree.h"
27 <
26 > #include  "platform.h"
27 > #include  "ray.h"
28 > #include  "ambient.h"
29 > #include  "source.h"
30   #include  "otypes.h"
31 <
31 > #include  "otspecial.h"
32   #include  "resolu.h"
33 + #include  "random.h"
34  
35 < int  dimlist[MAXDIM];                   /* sampling dimensions */
36 < int  ndims = 0;                         /* number of sampling dimensions */
37 < int  samplendx = 0;                     /* index for this sample */
35 > extern int  inform;                     /* input format */
36 > extern int  outform;                    /* output format */
37 > extern char  *outvals;                  /* output values */
38  
39 < int  imm_irrad = 0;                     /* compute immediate irradiance? */
39 > extern int  imm_irrad;                  /* compute immediate irradiance? */
40 > extern int  lim_dist;                   /* limit distance? */
41  
42 < int  inform = 'a';                      /* input format */
43 < int  outform = 'a';                     /* output format */
43 < char  *outvals = "v";                   /* output specification */
42 > extern char  *tralist[];                /* list of modifers to trace (or no) */
43 > extern int  traincl;                    /* include == 1, exclude == 0 */
44  
45 < int  hresolu = 0;                       /* horizontal (scan) size */
46 < int  vresolu = 0;                       /* vertical resolution */
45 > extern int  hresolu;                    /* horizontal resolution */
46 > extern int  vresolu;                    /* vertical resolution */
47  
48 < double  dstrsrc = 0.0;                  /* square source distribution */
49 < double  shadthresh = .05;               /* shadow threshold */
50 < double  shadcert = .5;                  /* shadow certainty */
51 < int  directrelay = 1;                   /* number of source relays */
52 < int  vspretest = 512;                   /* virtual source pretest density */
53 < int  directvis = 1;                     /* sources visible? */
54 < double  srcsizerat = .25;               /* maximum ratio source size/dist. */
48 > static int  castonly = 0;
49  
50 < double  specthresh = .15;               /* specular sampling threshold */
51 < double  specjitter = 1.;                /* specular sampling jitter */
50 > #ifndef  MAXTSET
51 > #define  MAXTSET        8191            /* maximum number in trace set */
52 > #endif
53 > OBJECT  traset[MAXTSET+1]={0};          /* trace include/exclude set */
54  
59 – int  maxdepth = 6;                      /* maximum recursion depth */
60 – double  minweight = 4e-3;               /* minimum ray weight */
61 –
62 – COLOR  ambval = BLKCOLOR;               /* ambient value */
63 – double  ambacc = 0.2;                   /* ambient accuracy */
64 – int  ambres = 32;                       /* ambient resolution */
65 – int  ambdiv = 128;                      /* ambient divisions */
66 – int  ambssamp = 0;                      /* ambient super-samples */
67 – int  ambounce = 0;                      /* ambient bounces */
68 – char  *amblist[128];                    /* ambient include/exclude list */
69 – int  ambincl = -1;                      /* include == 1, exclude == 0 */
70 –
71 – extern OBJREC  Lamb;                    /* a Lambertian surface */
72 –
55   static RAY  thisray;                    /* for our convenience */
56  
57 < static int  oputo(), oputd(), oputv(), oputl(), oputL(),
58 <                oputp(), oputn(), oputN(), oputs(), oputw(), oputm();
57 > typedef void putf_t(RREAL *v, int n);
58 > static putf_t puta, putd, putf;
59  
60 < static int  ourtrace(), tabin();
61 < static int  (*ray_out[16])(), (*every_out[16])();
62 < static int  castonly = 0;
60 > typedef void oputf_t(RAY *r);
61 > static oputf_t  oputo, oputd, oputv, oputV, oputl, oputL, oputc, oputp,
62 >                oputr, oputR, oputx, oputX, oputn, oputN, oputs,
63 >                oputw, oputW, oputm, oputM, oputtilde;
64  
65 < static int  puta(), putf(), putd();
65 > static void setoutput(char *vs);
66 > extern void tranotify(OBJECT obj);
67 > static void bogusray(void);
68 > static void raycast(RAY *r);
69 > static void rayirrad(RAY *r);
70 > static void rtcompute(FVECT org, FVECT dir, double dmax);
71 > static int printvals(RAY *r);
72 > static int getvec(FVECT vec, int fmt, FILE *fp);
73 > static void tabin(RAY *r);
74 > static void ourtrace(RAY *r);
75  
76 < static int  (*putreal)();
76 > static oputf_t *ray_out[32], *every_out[32];
77 > static putf_t *putreal;
78  
79  
80 < quit(code)                      /* quit program */
81 < int  code;
80 > void
81 > quit(                   /* quit program */
82 >        int  code
83 > )
84   {
85 < #ifndef  NIX
86 <        headclean();            /* delete header file */
87 <        pfclean();              /* clean up persist files */
85 >        if (ray_pnprocs > 0)    /* close children if any */
86 >                ray_pclose(0);          
87 > #ifndef  NON_POSIX
88 >        else if (!ray_pnprocs) {
89 >                headclean();    /* delete header file */
90 >                pfclean();      /* clean up persist files */
91 >        }
92   #endif
93          exit(code);
94   }
95  
96  
97   char *
98 < formstr(f)                              /* return format identifier */
99 < int  f;
98 > formstr(                                /* return format identifier */
99 >        int  f
100 > )
101   {
102          switch (f) {
103          case 'a': return("ascii");
# Line 109 | Line 109 | int  f;
109   }
110  
111  
112 < rtrace(fname)                           /* trace rays from file */
113 < char  *fname;
112 > extern void
113 > rtrace(                         /* trace rays from file */
114 >        char  *fname,
115 >        int  nproc
116 > )
117   {
118 <        long  vcount = hresolu>1 ? hresolu*vresolu : vresolu;
119 <        long  nextflush = hresolu;
118 >        unsigned long  vcount = (hresolu > 1) ? (unsigned long)hresolu*vresolu
119 >                                              : (unsigned long)vresolu;
120 >        long  nextflush = (vresolu > 0) & (hresolu > 1) ? 0 : hresolu;
121          FILE  *fp;
122 +        double  d;
123          FVECT  orig, direc;
124                                          /* set up input */
125          if (fname == NULL)
# Line 123 | Line 128 | char  *fname;
128                  sprintf(errmsg, "cannot open input file \"%s\"", fname);
129                  error(SYSTEM, errmsg);
130          }
126 – #ifdef MSDOS
131          if (inform != 'a')
132 <                setmode(fileno(fp), O_BINARY);
129 < #endif
132 >                SET_FILE_BINARY(fp);
133                                          /* set up output */
134 +        setoutput(outvals);
135          if (imm_irrad)
136 <                outvals = "v";
137 <        else
138 <                setoutput(outvals);
136 >                castonly = 0;
137 >        else if (castonly)
138 >                nproc = 1;              /* don't bother multiprocessing */
139          switch (outform) {
140          case 'a': putreal = puta; break;
141          case 'f': putreal = putf; break;
# Line 143 | Line 147 | char  *fname;
147          default:
148                  error(CONSISTENCY, "botched output format");
149          }
150 +        if (nproc > 1) {                /* start multiprocessing */
151 +                ray_popen(nproc);
152 +                ray_fifo_out = printvals;
153 +        }
154          if (hresolu > 0) {
155                  if (vresolu > 0)
156                          fprtresolu(hresolu, vresolu, stdout);
# Line 152 | Line 160 | char  *fname;
160          while (getvec(orig, inform, fp) == 0 &&
161                          getvec(direc, inform, fp) == 0) {
162  
163 <                if (normalize(direc) == 0.0) {          /* zero ==> flush */
164 <                        fflush(stdout);
165 <                        continue;
166 <                }
167 <                samplendx++;
168 <                                                        /* compute and print */
169 <                if (imm_irrad)
170 <                        irrad(orig, direc);
171 <                else
172 <                        traceray(orig, direc);
163 >                d = normalize(direc);
164 >                if (d == 0.0) {                         /* zero ==> flush */
165 >                        if (--nextflush <= 0 || !vcount) {
166 >                                if (nproc > 1 && ray_fifo_flush() < 0)
167 >                                        error(USER, "child(ren) died");
168 >                                bogusray();
169 >                                fflush(stdout);
170 >                                nextflush = (vresolu > 0) & (hresolu > 1) ? 0 :
171 >                                                                hresolu;
172 >                        } else
173 >                                bogusray();
174 >                } else {                                /* compute and print */
175 >                        rtcompute(orig, direc, lim_dist ? d : 0.0);
176                                                          /* flush if time */
177 <                if (--nextflush == 0) {
178 <                        fflush(stdout);
179 <                        nextflush = hresolu;
177 >                        if (!--nextflush) {
178 >                                if (nproc > 1 && ray_fifo_flush() < 0)
179 >                                        error(USER, "child(ren) died");
180 >                                fflush(stdout);
181 >                                nextflush = hresolu;
182 >                        }
183                  }
184                  if (ferror(stdout))
185                          error(SYSTEM, "write error");
186 <                if (--vcount == 0)                      /* check for end */
186 >                if (vcount && !--vcount)                /* check for end */
187                          break;
188          }
189 <        fflush(stdout);
190 <        if (vcount > 0)
191 <                error(USER, "read error");
189 >        if (nproc > 1) {                                /* clean up children */
190 >                if (ray_fifo_flush() < 0)
191 >                        error(USER, "unable to complete processing");
192 >                ray_pclose(0);
193 >        }
194 >        if (fflush(stdout) < 0)
195 >                error(SYSTEM, "write error");
196 >        if (vcount)
197 >                error(USER, "unexpected EOF on input");
198          if (fname != NULL)
199                  fclose(fp);
200   }
201  
202  
203 < setoutput(vs)                           /* set up output tables */
204 < register char  *vs;
203 > static void
204 > trace_sources(void)                     /* trace rays to light sources, also */
205   {
206 <        extern int  (*trace)();
207 <        register int (**table)() = ray_out;
206 >        int     sn;
207 >        
208 >        for (sn = 0; sn < nsources; sn++)
209 >                source[sn].sflags |= SFOLLOW;
210 > }
211  
212 +
213 + static void
214 + setoutput(                              /* set up output tables */
215 +        char  *vs
216 + )
217 + {
218 +        oputf_t **table = ray_out;
219 +
220          castonly = 1;
221          while (*vs)
222                  switch (*vs++) {
223 +                case 'T':                               /* trace sources */
224 +                        if (!*vs) break;
225 +                        trace_sources();
226 +                        /* fall through */
227                  case 't':                               /* trace */
228 +                        if (!*vs) break;
229                          *table = NULL;
230                          table = every_out;
231                          trace = ourtrace;
# Line 201 | Line 237 | register char  *vs;
237                  case 'd':                               /* direction */
238                          *table++ = oputd;
239                          break;
240 +                case 'r':                               /* reflected contrib. */
241 +                        *table++ = oputr;
242 +                        castonly = 0;
243 +                        break;
244 +                case 'R':                               /* reflected distance */
245 +                        *table++ = oputR;
246 +                        castonly = 0;
247 +                        break;
248 +                case 'x':                               /* xmit contrib. */
249 +                        *table++ = oputx;
250 +                        castonly = 0;
251 +                        break;
252 +                case 'X':                               /* xmit distance */
253 +                        *table++ = oputX;
254 +                        castonly = 0;
255 +                        break;
256                  case 'v':                               /* value */
257                          *table++ = oputv;
258                          castonly = 0;
259                          break;
260 +                case 'V':                               /* contribution */
261 +                        *table++ = oputV;
262 +                        castonly = 0;
263 +                        if (ambounce > 0 && (ambacc > FTINY || ambssamp > 0))
264 +                                error(WARNING,
265 +                                        "-otV accuracy depends on -aa 0 -as 0");
266 +                        break;
267                  case 'l':                               /* effective distance */
268                          *table++ = oputl;
269                          castonly = 0;
270                          break;
271 +                case 'c':                               /* local coordinates */
272 +                        *table++ = oputc;
273 +                        break;
274                  case 'L':                               /* single ray length */
275                          *table++ = oputL;
276                          break;
# Line 228 | Line 290 | register char  *vs;
290                  case 'w':                               /* weight */
291                          *table++ = oputw;
292                          break;
293 +                case 'W':                               /* coefficient */
294 +                        *table++ = oputW;
295 +                        castonly = 0;
296 +                        if (ambounce > 0 && (ambacc > FTINY || ambssamp > 0))
297 +                                error(WARNING,
298 +                                        "-otW accuracy depends on -aa 0 -as 0");
299 +                        break;
300                  case 'm':                               /* modifier */
301                          *table++ = oputm;
302                          break;
303 +                case 'M':                               /* material */
304 +                        *table++ = oputM;
305 +                        break;
306 +                case '~':                               /* tilde */
307 +                        *table++ = oputtilde;
308 +                        break;
309                  }
310          *table = NULL;
311   }
312  
313  
314 < traceray(org, dir)              /* compute and print ray value(s) */
315 < FVECT  org, dir;
314 > static void
315 > bogusray(void)                  /* print out empty record */
316   {
317 <        register int  (**tp)();
317 >        thisray.rorg[0] = thisray.rorg[1] = thisray.rorg[2] =
318 >        thisray.rdir[0] = thisray.rdir[1] = thisray.rdir[2] = 0.0;
319 >        thisray.rmax = 0.0;
320 >        rayorigin(&thisray, PRIMARY, NULL, NULL);
321 >        printvals(&thisray);
322 > }
323  
244 –        VCOPY(thisray.rorg, org);
245 –        VCOPY(thisray.rdir, dir);
246 –        rayorigin(&thisray, NULL, PRIMARY, 1.0);
247 –        if (castonly)
248 –                localhit(&thisray, &thescene) || sourcehit(&thisray);
249 –        else
250 –                rayvalue(&thisray);
324  
325 <        if (ray_out[0] == NULL)
326 <                return;
327 <        for (tp = ray_out; *tp != NULL; tp++)
328 <                (**tp)(&thisray);
329 <        if (outform == 'a')
330 <                putchar('\n');
325 > static void
326 > raycast(                        /* compute first ray intersection only */
327 >        RAY *r
328 > )
329 > {
330 >        if (!localhit(r, &thescene)) {
331 >                if (r->ro == &Aftplane) {       /* clipped */
332 >                        r->ro = NULL;
333 >                        r->rot = FHUGE;
334 >                } else
335 >                        sourcehit(r);
336 >        }
337   }
338  
339  
340 < irrad(org, dir)                 /* compute immediate irradiance value */
341 < FVECT  org, dir;
340 > static void
341 > rayirrad(                       /* compute irradiance rather than radiance */
342 >        RAY *r
343 > )
344   {
345 <        register int  i;
345 >        void    (*old_revf)(RAY *) = r->revf;
346  
347 <        for (i = 0; i < 3; i++) {
348 <                thisray.rorg[i] = org[i] + dir[i];
349 <                thisray.rdir[i] = -dir[i];
347 >        r->rot = 1e-5;                  /* pretend we hit surface */
348 >        VSUM(r->rop, r->rorg, r->rdir, r->rot);
349 >        r->ron[0] = -r->rdir[0];
350 >        r->ron[1] = -r->rdir[1];
351 >        r->ron[2] = -r->rdir[2];
352 >        r->rod = 1.0;
353 >                                        /* compute result */
354 >        r->revf = raytrace;
355 >        (*ofun[Lamb.otype].funp)(&Lamb, r);
356 >        r->revf = old_revf;
357 > }
358 >
359 >
360 > static void
361 > rtcompute(                      /* compute and print ray value(s) */
362 >        FVECT  org,
363 >        FVECT  dir,
364 >        double  dmax
365 > )
366 > {
367 >                                        /* set up ray */
368 >        rayorigin(&thisray, PRIMARY, NULL, NULL);
369 >        if (imm_irrad) {
370 >                VSUM(thisray.rorg, org, dir, 1.1e-4);
371 >                thisray.rdir[0] = -dir[0];
372 >                thisray.rdir[1] = -dir[1];
373 >                thisray.rdir[2] = -dir[2];
374 >                thisray.rmax = 0.0;
375 >                thisray.revf = rayirrad;
376 >        } else {
377 >                VCOPY(thisray.rorg, org);
378 >                VCOPY(thisray.rdir, dir);
379 >                thisray.rmax = dmax;
380 >                if (castonly)
381 >                        thisray.revf = raycast;
382          }
383 <        rayorigin(&thisray, NULL, PRIMARY, 1.0);
384 <                                        /* pretend we hit surface */
385 <        thisray.rot = 1.0;
386 <        thisray.rod = 1.0;
387 <        VCOPY(thisray.ron, dir);
388 <        for (i = 0; i < 3; i++)         /* fudge factor */
389 <                thisray.rop[i] = org[i] + 1e-4*dir[i];
390 <                                        /* compute and print */
391 <        (*ofun[Lamb.otype].funp)(&Lamb, &thisray);
392 <        oputv(&thisray);
383 >        if (ray_pnprocs > 1) {          /* multiprocessing FIFO? */
384 >                if (ray_fifo_in(&thisray) < 0)
385 >                        error(USER, "lost children");
386 >                return;
387 >        }
388 >        samplendx++;                    /* else do it ourselves */
389 >        rayvalue(&thisray);
390 >        printvals(&thisray);
391 > }
392 >
393 >
394 > static int
395 > printvals(                      /* print requested ray values */
396 >        RAY  *r
397 > )
398 > {
399 >        oputf_t **tp;
400 >
401 >        if (ray_out[0] == NULL)
402 >                return(0);
403 >        for (tp = ray_out; *tp != NULL; tp++)
404 >                (**tp)(r);
405          if (outform == 'a')
406                  putchar('\n');
407 +        return(1);
408   }
409  
410  
411 < getvec(vec, fmt, fp)            /* get a vector from fp */
412 < register FVECT  vec;
413 < int  fmt;
414 < FILE  *fp;
411 > static int
412 > getvec(         /* get a vector from fp */
413 >        FVECT  vec,
414 >        int  fmt,
415 >        FILE  *fp
416 > )
417   {
290 –        extern char  *fgetword();
418          static float  vf[3];
419          static double  vd[3];
420          char  buf[32];
421 <        register int  i;
421 >        int  i;
422  
423          switch (fmt) {
424          case 'a':                                       /* ascii */
# Line 303 | Line 430 | FILE  *fp;
430                  }
431                  break;
432          case 'f':                                       /* binary float */
433 <                if (fread((char *)vf, sizeof(float), 3, fp) != 3)
433 >                if (getbinary(vf, sizeof(float), 3, fp) != 3)
434                          return(-1);
435 <                vec[0] = vf[0]; vec[1] = vf[1]; vec[2] = vf[2];
435 >                VCOPY(vec, vf);
436                  break;
437          case 'd':                                       /* binary double */
438 <                if (fread((char *)vd, sizeof(double), 3, fp) != 3)
438 >                if (getbinary(vd, sizeof(double), 3, fp) != 3)
439                          return(-1);
440 <                vec[0] = vd[0]; vec[1] = vd[1]; vec[2] = vd[2];
440 >                VCOPY(vec, vd);
441                  break;
442          default:
443                  error(CONSISTENCY, "botched input format");
# Line 319 | Line 446 | FILE  *fp;
446   }
447  
448  
449 < static
450 < ourtrace(r)                             /* print ray values */
451 < RAY  *r;
449 > void
450 > tranotify(                      /* record new modifier */
451 >        OBJECT  obj
452 > )
453   {
454 <        register int  (**tp)();
454 >        static int  hitlimit = 0;
455 >        OBJREC   *o = objptr(obj);
456 >        char  **tralp;
457  
458 +        if (obj == OVOID) {             /* starting over */
459 +                traset[0] = 0;
460 +                hitlimit = 0;
461 +                return;
462 +        }
463 +        if (hitlimit || !ismodifier(o->otype))
464 +                return;
465 +        for (tralp = tralist; *tralp != NULL; tralp++)
466 +                if (!strcmp(o->oname, *tralp)) {
467 +                        if (traset[0] >= MAXTSET) {
468 +                                error(WARNING, "too many modifiers in trace list");
469 +                                hitlimit++;
470 +                                return;         /* should this be fatal? */
471 +                        }
472 +                        insertelem(traset, obj);
473 +                        return;
474 +                }
475 + }
476 +
477 +
478 + static void
479 + ourtrace(                               /* print ray values */
480 +        RAY  *r
481 + )
482 + {
483 +        oputf_t **tp;
484 +
485          if (every_out[0] == NULL)
486                  return;
487 +        if (r->ro == NULL) {
488 +                if (traincl == 1)
489 +                        return;
490 +        } else if (traincl != -1 && traincl != inset(traset, r->ro->omod))
491 +                return;
492          tabin(r);
493          for (tp = every_out; *tp != NULL; tp++)
494                  (**tp)(r);
495 <        putchar('\n');
495 >        if (outform == 'a')
496 >                putchar('\n');
497   }
498  
499  
500 < static
501 < tabin(r)                                /* tab in appropriate amount */
502 < RAY  *r;
500 > static void
501 > tabin(                          /* tab in appropriate amount */
502 >        RAY  *r
503 > )
504   {
505 <        register RAY  *rp;
505 >        const RAY  *rp;
506  
507          for (rp = r->parent; rp != NULL; rp = rp->parent)
508                  putchar('\t');
509   }
510  
511  
512 < static
513 < oputo(r)                                /* print origin */
514 < register RAY  *r;
512 > static void
513 > oputo(                          /* print origin */
514 >        RAY  *r
515 > )
516   {
517 <        (*putreal)(r->rorg[0]);
353 <        (*putreal)(r->rorg[1]);
354 <        (*putreal)(r->rorg[2]);
517 >        (*putreal)(r->rorg, 3);
518   }
519  
520  
521 < static
522 < oputd(r)                                /* print direction */
523 < register RAY  *r;
521 > static void
522 > oputd(                          /* print direction */
523 >        RAY  *r
524 > )
525   {
526 <        (*putreal)(r->rdir[0]);
363 <        (*putreal)(r->rdir[1]);
364 <        (*putreal)(r->rdir[2]);
526 >        (*putreal)(r->rdir, 3);
527   }
528  
529  
530 < static
531 < oputv(r)                                /* print value */
532 < register RAY  *r;
530 > static void
531 > oputr(                          /* print mirrored contribution */
532 >        RAY  *r
533 > )
534   {
535 <        COLR  cout;
536 <        
535 >        RREAL   cval[3];
536 >
537 >        cval[0] = colval(r->mcol,RED);
538 >        cval[1] = colval(r->mcol,GRN);
539 >        cval[2] = colval(r->mcol,BLU);
540 >        (*putreal)(cval, 3);
541 > }
542 >
543 >
544 >
545 > static void
546 > oputR(                          /* print mirrored distance */
547 >        RAY  *r
548 > )
549 > {
550 >        (*putreal)(&r->rmt, 1);
551 > }
552 >
553 >
554 > static void
555 > oputx(                          /* print unmirrored contribution */
556 >        RAY  *r
557 > )
558 > {
559 >        RREAL   cval[3];
560 >
561 >        cval[0] = colval(r->rcol,RED) - colval(r->mcol,RED);
562 >        cval[1] = colval(r->rcol,GRN) - colval(r->mcol,GRN);
563 >        cval[2] = colval(r->rcol,BLU) - colval(r->mcol,BLU);
564 >        (*putreal)(cval, 3);
565 > }
566 >
567 >
568 > static void
569 > oputX(                          /* print unmirrored distance */
570 >        RAY  *r
571 > )
572 > {
573 >        (*putreal)(&r->rxt, 1);
574 > }
575 >
576 >
577 > static void
578 > oputv(                          /* print value */
579 >        RAY  *r
580 > )
581 > {
582 >        RREAL   cval[3];
583 >
584          if (outform == 'c') {
585 +                COLR  cout;
586                  setcolr(cout,   colval(r->rcol,RED),
587                                  colval(r->rcol,GRN),
588                                  colval(r->rcol,BLU));
589 <                fwrite((char *)cout, sizeof(cout), 1, stdout);
589 >                putbinary(cout, sizeof(cout), 1, stdout);
590                  return;
591          }
592 <        (*putreal)(colval(r->rcol,RED));
593 <        (*putreal)(colval(r->rcol,GRN));
594 <        (*putreal)(colval(r->rcol,BLU));
592 >        cval[0] = colval(r->rcol,RED);
593 >        cval[1] = colval(r->rcol,GRN);
594 >        cval[2] = colval(r->rcol,BLU);
595 >        (*putreal)(cval, 3);
596   }
597  
598  
599 < static
600 < oputl(r)                                /* print effective distance */
601 < register RAY  *r;
599 > static void
600 > oputV(                          /* print value contribution */
601 >        RAY *r
602 > )
603   {
604 <        (*putreal)(r->rt);
604 >        RREAL   contr[3];
605 >
606 >        raycontrib(contr, r, PRIMARY);
607 >        multcolor(contr, r->rcol);
608 >        (*putreal)(contr, 3);
609   }
610  
611  
612 < static
613 < oputL(r)                                /* print single ray length */
614 < register RAY  *r;
612 > static void
613 > oputl(                          /* print effective distance */
614 >        RAY  *r
615 > )
616   {
617 <        (*putreal)(r->rot);
617 >        RREAL   d = raydistance(r);
618 >
619 >        (*putreal)(&d, 1);
620   }
621  
622  
623 < static
624 < oputp(r)                                /* print point */
625 < register RAY  *r;
623 > static void
624 > oputL(                          /* print single ray length */
625 >        RAY  *r
626 > )
627   {
628 <        if (r->rot < FHUGE) {
408 <                (*putreal)(r->rop[0]);
409 <                (*putreal)(r->rop[1]);
410 <                (*putreal)(r->rop[2]);
411 <        } else {
412 <                (*putreal)(0.0);
413 <                (*putreal)(0.0);
414 <                (*putreal)(0.0);
415 <        }
628 >        (*putreal)(&r->rot, 1);
629   }
630  
631  
632 < static
633 < oputN(r)                                /* print unperturbed normal */
634 < register RAY  *r;
632 > static void
633 > oputc(                          /* print local coordinates */
634 >        RAY  *r
635 > )
636   {
637 <        if (r->rot < FHUGE) {
424 <                (*putreal)(r->ron[0]);
425 <                (*putreal)(r->ron[1]);
426 <                (*putreal)(r->ron[2]);
427 <        } else {
428 <                (*putreal)(0.0);
429 <                (*putreal)(0.0);
430 <                (*putreal)(0.0);
431 <        }
637 >        (*putreal)(r->uv, 2);
638   }
639  
640  
641 < static
642 < oputn(r)                                /* print perturbed normal */
643 < RAY  *r;
641 > static RREAL    vdummy[3] = {0.0, 0.0, 0.0};
642 >
643 >
644 > static void
645 > oputp(                          /* print point */
646 >        RAY  *r
647 > )
648   {
649 +        if (r->rot < FHUGE)
650 +                (*putreal)(r->rop, 3);
651 +        else
652 +                (*putreal)(vdummy, 3);
653 + }
654 +
655 +
656 + static void
657 + oputN(                          /* print unperturbed normal */
658 +        RAY  *r
659 + )
660 + {
661 +        if (r->rot < FHUGE)
662 +                (*putreal)(r->ron, 3);
663 +        else
664 +                (*putreal)(vdummy, 3);
665 + }
666 +
667 +
668 + static void
669 + oputn(                          /* print perturbed normal */
670 +        RAY  *r
671 + )
672 + {
673          FVECT  pnorm;
674  
675          if (r->rot >= FHUGE) {
676 <                (*putreal)(0.0);
443 <                (*putreal)(0.0);
444 <                (*putreal)(0.0);
676 >                (*putreal)(vdummy, 3);
677                  return;
678          }
679          raynormal(pnorm, r);
680 <        (*putreal)(pnorm[0]);
449 <        (*putreal)(pnorm[1]);
450 <        (*putreal)(pnorm[2]);
680 >        (*putreal)(pnorm, 3);
681   }
682  
683  
684 < static
685 < oputs(r)                                /* print name */
686 < register RAY  *r;
684 > static void
685 > oputs(                          /* print name */
686 >        RAY  *r
687 > )
688   {
689          if (r->ro != NULL)
690                  fputs(r->ro->oname, stdout);
# Line 463 | Line 694 | register RAY  *r;
694   }
695  
696  
697 < static
698 < oputw(r)                                /* print weight */
699 < register RAY  *r;
697 > static void
698 > oputw(                          /* print weight */
699 >        RAY  *r
700 > )
701   {
702 <        (*putreal)(r->rweight);
702 >        RREAL   rwt = r->rweight;
703 >        
704 >        (*putreal)(&rwt, 1);
705   }
706  
707  
708 < static
709 < oputm(r)                                /* print modifier */
710 < register RAY  *r;
708 > static void
709 > oputW(                          /* print coefficient */
710 >        RAY  *r
711 > )
712   {
713 +        RREAL   contr[3];
714 +                                /* shadow ray not on source? */
715 +        if (r->rsrc >= 0 && source[r->rsrc].so != r->ro)
716 +                setcolor(contr, 0.0, 0.0, 0.0);
717 +        else
718 +                raycontrib(contr, r, PRIMARY);
719 +
720 +        (*putreal)(contr, 3);
721 + }
722 +
723 +
724 + static void
725 + oputm(                          /* print modifier */
726 +        RAY  *r
727 + )
728 + {
729          if (r->ro != NULL)
730 <                fputs(objptr(r->ro->omod)->oname, stdout);
730 >                if (r->ro->omod != OVOID)
731 >                        fputs(objptr(r->ro->omod)->oname, stdout);
732 >                else
733 >                        fputs(VOIDID, stdout);
734          else
735                  putchar('*');
736          putchar('\t');
737   }
738  
739  
740 < static
741 < puta(v)                         /* print ascii value */
742 < double  v;
740 > static void
741 > oputM(                          /* print material */
742 >        RAY  *r
743 > )
744   {
745 <        printf("%e\t", v);
745 >        OBJREC  *mat;
746 >
747 >        if (r->ro != NULL) {
748 >                if ((mat = findmaterial(r->ro)) != NULL)
749 >                        fputs(mat->oname, stdout);
750 >                else
751 >                        fputs(VOIDID, stdout);
752 >        } else
753 >                putchar('*');
754 >        putchar('\t');
755   }
756  
757  
758 < static
759 < putd(v)                         /* print binary double */
760 < double  v;
758 > static void
759 > oputtilde(                      /* output tilde (spacer) */
760 >        RAY  *r
761 > )
762   {
763 <        fwrite((char *)&v, sizeof(v), 1, stdout);
763 >        fputs("~\t", stdout);
764   }
765  
766  
767 < static
768 < putf(v)                         /* print binary float */
769 < double  v;
767 > static void
768 > puta(                           /* print ascii value(s) */
769 >        RREAL *v, int n
770 > )
771   {
772 <        float f = v;
772 >        if (n == 3) {
773 >                printf("%e\t%e\t%e\t", v[0], v[1], v[2]);
774 >                return;
775 >        }
776 >        while (n--)
777 >                printf("%e\t", *v++);
778 > }
779  
780 <        fwrite((char *)&f, sizeof(f), 1, stdout);
780 >
781 > static void
782 > putd(RREAL *v, int n)           /* print binary double(s) */
783 > {
784 > #ifdef  SMLFLT
785 >        double  da[3];
786 >        int     i;
787 >
788 >        if (n > 3)
789 >                error(INTERNAL, "code error in putd()");
790 >        for (i = n; i--; )
791 >                da[i] = v[i];
792 >        putbinary(da, sizeof(double), n, stdout);
793 > #else
794 >        putbinary(v, sizeof(RREAL), n, stdout);
795 > #endif
796 > }
797 >
798 >
799 > static void
800 > putf(RREAL *v, int n)           /* print binary float(s) */
801 > {
802 > #ifndef SMLFLT
803 >        float   fa[3];
804 >        int     i;
805 >
806 >        if (n > 3)
807 >                error(INTERNAL, "code error in putf()");
808 >        for (i = n; i--; )
809 >                fa[i] = v[i];
810 >        putbinary(fa, sizeof(float), n, stdout);
811 > #else
812 >        putbinary(v, sizeof(RREAL), n, stdout);
813 > #endif
814   }

Diff Legend

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