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.42 by greg, Thu May 26 06:55:22 2005 UTC vs.
Revision 2.72 by greg, Thu Nov 8 00:54:07 2018 UTC

# Line 28 | Line 28 | static const char      RCSid[] = "$Id$";
28   #include  "ambient.h"
29   #include  "source.h"
30   #include  "otypes.h"
31 + #include  "otspecial.h"
32   #include  "resolu.h"
33 + #include  "random.h"
34  
35 < CUBE  thescene;                         /* our scene */
36 < OBJECT  nsceneobjs;                     /* number of objects in our scene */
35 > extern int  inform;                     /* input format */
36 > extern int  outform;                    /* output format */
37 > extern char  *outvals;                  /* output values */
38  
39 < int  dimlist[MAXDIM];                   /* sampling dimensions */
40 < int  ndims = 0;                         /* number of sampling dimensions */
38 < int  samplendx = 0;                     /* index for this sample */
39 > extern int  imm_irrad;                  /* compute immediate irradiance? */
40 > extern int  lim_dist;                   /* limit distance? */
41  
42 < int  imm_irrad = 0;                     /* compute immediate irradiance? */
43 < int  lim_dist = 0;                      /* limit distance? */
42 > extern char  *tralist[];                /* list of modifers to trace (or no) */
43 > extern int  traincl;                    /* include == 1, exclude == 0 */
44  
45 < int  inform = 'a';                      /* input format */
46 < int  outform = 'a';                     /* output format */
45 < char  *outvals = "v";                   /* output specification */
45 > extern int  hresolu;                    /* horizontal resolution */
46 > extern int  vresolu;                    /* vertical resolution */
47  
48 < int  do_irrad = 0;                      /* compute irradiance? */
48 > static int  castonly = 0;
49  
50 < void  (*trace)() = NULL;                /* trace call */
51 <
52 < char  *tralist[128];                    /* list of modifers to trace (or no) */
52 < int  traincl = -1;                      /* include == 1, exclude == 0 */
53 < #define  MAXTSET        511             /* maximum number in trace set */
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  
56 – int  hresolu = 0;                       /* horizontal (scan) size */
57 – int  vresolu = 0;                       /* vertical resolution */
58 –
59 – double  dstrsrc = 0.0;                  /* square source distribution */
60 – double  shadthresh = .03;               /* shadow threshold */
61 – double  shadcert = .75;                 /* shadow certainty */
62 – int  directrelay = 2;                   /* number of source relays */
63 – int  vspretest = 512;                   /* virtual source pretest density */
64 – int  directvis = 1;                     /* sources visible? */
65 – double  srcsizerat = .2;                /* maximum ratio source size/dist. */
66 –
67 – COLOR  cextinction = BLKCOLOR;          /* global extinction coefficient */
68 – COLOR  salbedo = BLKCOLOR;              /* global scattering albedo */
69 – double  seccg = 0.;                     /* global scattering eccentricity */
70 – double  ssampdist = 0.;                 /* scatter sampling distance */
71 –
72 – double  specthresh = .15;               /* specular sampling threshold */
73 – double  specjitter = 1.;                /* specular sampling jitter */
74 –
75 – int  backvis = 1;                       /* back face visibility */
76 –
77 – int  maxdepth = 8;                      /* maximum recursion depth */
78 – double  minweight = 2e-3;               /* minimum ray weight */
79 –
80 – char  *ambfile = NULL;                  /* ambient file name */
81 – COLOR  ambval = BLKCOLOR;               /* ambient value */
82 – int  ambvwt = 0;                        /* initial weight for ambient value */
83 – double  ambacc = 0.15;                  /* ambient accuracy */
84 – int  ambres = 256;                      /* ambient resolution */
85 – int  ambdiv = 1024;                     /* ambient divisions */
86 – int  ambssamp = 512;                    /* ambient super-samples */
87 – int  ambounce = 0;                      /* ambient bounces */
88 – char  *amblist[AMBLLEN];                /* ambient include/exclude list */
89 – int  ambincl = -1;                      /* include == 1, exclude == 0 */
90 –
91 – static int  castonly = 0;
92 –
55   static RAY  thisray;                    /* for our convenience */
56  
57 < typedef void putf_t(double v);
57 > typedef void putf_t(RREAL *v, int n);
58   static putf_t puta, putd, putf;
59  
60   typedef void oputf_t(RAY *r);
61 < static oputf_t  oputo, oputd, oputv, oputl, oputL, oputc, oputp,
61 > static oputf_t  oputo, oputd, oputv, oputV, oputl, oputL, oputc, oputp,
62                  oputn, oputN, oputs, oputw, oputW, oputm, oputM, oputtilde;
63  
64   static void setoutput(char *vs);
65 < static void tranotify(OBJECT obj);
65 > extern void tranotify(OBJECT obj);
66   static void bogusray(void);
67 < static void rad(FVECT  org, FVECT  dir, double  dmax);
68 < static void irrad(FVECT  org, FVECT  dir);
69 < static void printvals(RAY  *r);
70 < static int getvec(FVECT  vec, int  fmt, FILE  *fp);
71 < static void tabin(RAY  *r);
72 < static void ourtrace(RAY  *r);
67 > static void raycast(RAY *r);
68 > static void rayirrad(RAY *r);
69 > static void rtcompute(FVECT org, FVECT dir, double dmax);
70 > static int printvals(RAY *r);
71 > static int getvec(FVECT vec, int fmt, FILE *fp);
72 > static void tabin(RAY *r);
73 > static void ourtrace(RAY *r);
74  
75   static oputf_t *ray_out[16], *every_out[16];
76   static putf_t *putreal;
77  
115 – void  (*addobjnotify[])() = {ambnotify, tranotify, NULL};
78  
117 –
79   void
80   quit(                   /* quit program */
81          int  code
82   )
83   {
84 < #ifndef  NON_POSIX /* XXX we don't clean up elsewhere? */
85 <        headclean();            /* delete header file */
86 <        pfclean();              /* clean up persist files */
84 >        if (ray_pnprocs > 0)    /* close children if any */
85 >                ray_pclose(0);          
86 > #ifndef  NON_POSIX
87 >        else if (!ray_pnprocs) {
88 >                headclean();    /* delete header file */
89 >                pfclean();      /* clean up persist files */
90 >        }
91   #endif
92          exit(code);
93   }
94  
95  
96 < extern char *
96 > char *
97   formstr(                                /* return format identifier */
98          int  f
99   )
# Line 145 | Line 110 | formstr(                               /* return format identifier */
110  
111   extern void
112   rtrace(                         /* trace rays from file */
113 <        char  *fname
113 >        char  *fname,
114 >        int  nproc
115   )
116   {
117 <        long  vcount = hresolu>1 ? hresolu*vresolu : vresolu;
118 <        long  nextflush = hresolu;
117 >        unsigned long  vcount = (hresolu > 1) ? (unsigned long)hresolu*vresolu
118 >                                              : (unsigned long)vresolu;
119 >        long  nextflush = (vresolu > 0) & (hresolu > 1) ? 0 : hresolu;
120          FILE  *fp;
121          double  d;
122          FVECT  orig, direc;
# Line 160 | Line 127 | rtrace(                                /* trace rays from file */
127                  sprintf(errmsg, "cannot open input file \"%s\"", fname);
128                  error(SYSTEM, errmsg);
129          }
163 – #ifdef _WIN32
130          if (inform != 'a')
131                  SET_FILE_BINARY(fp);
166 – #endif
132                                          /* set up output */
133          setoutput(outvals);
134 +        if (imm_irrad)
135 +                castonly = 0;
136 +        else if (castonly)
137 +                nproc = 1;              /* don't bother multiprocessing */
138          switch (outform) {
139          case 'a': putreal = puta; break;
140          case 'f': putreal = putf; break;
# Line 177 | Line 146 | rtrace(                                /* trace rays from file */
146          default:
147                  error(CONSISTENCY, "botched output format");
148          }
149 +        if (nproc > 1) {                /* start multiprocessing */
150 +                ray_popen(nproc);
151 +                ray_fifo_out = printvals;
152 +        }
153          if (hresolu > 0) {
154                  if (vresolu > 0)
155                          fprtresolu(hresolu, vresolu, stdout);
# Line 188 | Line 161 | rtrace(                                /* trace rays from file */
161  
162                  d = normalize(direc);
163                  if (d == 0.0) {                         /* zero ==> flush */
164 <                        bogusray();
165 <                        if (--nextflush <= 0 || vcount <= 0) {
164 >                        if (--nextflush <= 0 || !vcount) {
165 >                                if (nproc > 1 && ray_fifo_flush() < 0)
166 >                                        error(USER, "child(ren) died");
167 >                                bogusray();
168                                  fflush(stdout);
169 <                                nextflush = hresolu;
170 <                        }
171 <                } else {
172 <                        samplendx++;
173 <                                                        /* compute and print */
174 <                        if (imm_irrad)
200 <                                irrad(orig, direc);
201 <                        else
202 <                                rad(orig, direc, lim_dist ? d : 0.0);
169 >                                nextflush = (vresolu > 0) & (hresolu > 1) ? 0 :
170 >                                                                hresolu;
171 >                        } else
172 >                                bogusray();
173 >                } else {                                /* compute and print */
174 >                        rtcompute(orig, direc, lim_dist ? d : 0.0);
175                                                          /* flush if time */
176 <                        if (--nextflush == 0) {
176 >                        if (!--nextflush) {
177 >                                if (nproc > 1 && ray_fifo_flush() < 0)
178 >                                        error(USER, "child(ren) died");
179                                  fflush(stdout);
180                                  nextflush = hresolu;
181                          }
182                  }
183                  if (ferror(stdout))
184                          error(SYSTEM, "write error");
185 <                if (--vcount == 0)                      /* check for end */
185 >                if (vcount && !--vcount)                /* check for end */
186                          break;
187          }
188 <        fflush(stdout);
189 <        if (vcount > 0)
188 >        if (nproc > 1) {                                /* clean up children */
189 >                if (ray_fifo_flush() < 0)
190 >                        error(USER, "unable to complete processing");
191 >                ray_pclose(0);
192 >        }
193 >        if (fflush(stdout) < 0)
194 >                error(SYSTEM, "write error");
195 >        if (vcount)
196                  error(USER, "unexpected EOF on input");
197          if (fname != NULL)
198                  fclose(fp);
# Line 231 | Line 211 | trace_sources(void)                    /* trace rays to light sources,
211  
212   static void
213   setoutput(                              /* set up output tables */
214 <        register char  *vs
214 >        char  *vs
215   )
216   {
217 <        register oputf_t **table = ray_out;
217 >        oputf_t **table = ray_out;
218  
219          castonly = 1;
220          while (*vs)
# Line 260 | Line 240 | setoutput(                             /* set up output tables */
240                          *table++ = oputv;
241                          castonly = 0;
242                          break;
243 +                case 'V':                               /* contribution */
244 +                        *table++ = oputV;
245 +                        if (ambounce > 0 && (ambacc > FTINY || ambssamp > 0))
246 +                                error(WARNING,
247 +                                        "-otV accuracy depends on -aa 0 -as 0");
248 +                        break;
249                  case 'l':                               /* effective distance */
250                          *table++ = oputl;
251                          castonly = 0;
# Line 318 | Line 304 | bogusray(void)                 /* print out empty record */
304  
305  
306   static void
307 < rad(            /* compute and print ray value(s) */
308 <        FVECT  org,
323 <        FVECT  dir,
324 <        double  dmax
307 > raycast(                        /* compute first ray intersection only */
308 >        RAY *r
309   )
310   {
311 <        VCOPY(thisray.rorg, org);
312 <        VCOPY(thisray.rdir, dir);
313 <        thisray.rmax = dmax;
314 <        rayorigin(&thisray, PRIMARY, NULL, NULL);
315 <        if (castonly) {
316 <                if (!localhit(&thisray, &thescene)) {
317 <                        if (thisray.ro == &Aftplane) {  /* clipped */
334 <                                thisray.ro = NULL;
335 <                                thisray.rot = FHUGE;
336 <                        } else
337 <                                sourcehit(&thisray);
338 <                }
339 <        } else
340 <                rayvalue(&thisray);
341 <        printvals(&thisray);
311 >        if (!localhit(r, &thescene)) {
312 >                if (r->ro == &Aftplane) {       /* clipped */
313 >                        r->ro = NULL;
314 >                        r->rot = FHUGE;
315 >                } else
316 >                        sourcehit(r);
317 >        }
318   }
319  
320  
321   static void
322 < irrad(                  /* compute immediate irradiance value */
323 <        FVECT  org,
348 <        FVECT  dir
322 > rayirrad(                       /* compute irradiance rather than radiance */
323 >        RAY *r
324   )
325   {
326 <        register int  i;
326 >        void    (*old_revf)(RAY *) = r->revf;
327  
328 <        for (i = 0; i < 3; i++) {
329 <                thisray.rorg[i] = org[i] + dir[i];
330 <                thisray.rdir[i] = -dir[i];
331 <        }
332 <        thisray.rmax = 0.0;
328 >        r->rot = 1e-5;                  /* pretend we hit surface */
329 >        VSUM(r->rop, r->rorg, r->rdir, r->rot);
330 >        r->ron[0] = -r->rdir[0];
331 >        r->ron[1] = -r->rdir[1];
332 >        r->ron[2] = -r->rdir[2];
333 >        r->rod = 1.0;
334 >                                        /* compute result */
335 >        r->revf = raytrace;
336 >        (*ofun[Lamb.otype].funp)(&Lamb, r);
337 >        r->revf = old_revf;
338 > }
339 >
340 >
341 > static void
342 > rtcompute(                      /* compute and print ray value(s) */
343 >        FVECT  org,
344 >        FVECT  dir,
345 >        double  dmax
346 > )
347 > {
348 >                                        /* set up ray */
349          rayorigin(&thisray, PRIMARY, NULL, NULL);
350 <                                        /* pretend we hit surface */
351 <        thisray.rot = 1.0-1e-4;
352 <        thisray.rod = 1.0;
353 <        VCOPY(thisray.ron, dir);
354 <        for (i = 0; i < 3; i++)         /* fudge factor */
355 <                thisray.rop[i] = org[i] + 1e-4*dir[i];
356 <                                        /* compute and print */
357 <        (*ofun[Lamb.otype].funp)(&Lamb, &thisray);
350 >        if (imm_irrad) {
351 >                VSUM(thisray.rorg, org, dir, 1.1e-4);
352 >                thisray.rdir[0] = -dir[0];
353 >                thisray.rdir[1] = -dir[1];
354 >                thisray.rdir[2] = -dir[2];
355 >                thisray.rmax = 0.0;
356 >                thisray.revf = rayirrad;
357 >        } else {
358 >                VCOPY(thisray.rorg, org);
359 >                VCOPY(thisray.rdir, dir);
360 >                thisray.rmax = dmax;
361 >                if (castonly)
362 >                        thisray.revf = raycast;
363 >        }
364 >        if (ray_pnprocs > 1) {          /* multiprocessing FIFO? */
365 >                if (ray_fifo_in(&thisray) < 0)
366 >                        error(USER, "lost children");
367 >                return;
368 >        }
369 >        samplendx++;                    /* else do it ourselves */
370 >        rayvalue(&thisray);
371          printvals(&thisray);
372   }
373  
374  
375 < static void
375 > static int
376   printvals(                      /* print requested ray values */
377          RAY  *r
378   )
379   {
380 <        register oputf_t **tp;
380 >        oputf_t **tp;
381  
382          if (ray_out[0] == NULL)
383 <                return;
383 >                return(0);
384          for (tp = ray_out; *tp != NULL; tp++)
385                  (**tp)(r);
386          if (outform == 'a')
387                  putchar('\n');
388 +        return(1);
389   }
390  
391  
392   static int
393   getvec(         /* get a vector from fp */
394 <        register FVECT  vec,
394 >        FVECT  vec,
395          int  fmt,
396          FILE  *fp
397   )
# Line 394 | Line 399 | getvec(                /* get a vector from fp */
399          static float  vf[3];
400          static double  vd[3];
401          char  buf[32];
402 <        register int  i;
402 >        int  i;
403  
404          switch (fmt) {
405          case 'a':                                       /* ascii */
# Line 406 | Line 411 | getvec(                /* get a vector from fp */
411                  }
412                  break;
413          case 'f':                                       /* binary float */
414 <                if (fread((char *)vf, sizeof(float), 3, fp) != 3)
414 >                if (getbinary(vf, sizeof(float), 3, fp) != 3)
415                          return(-1);
416 <                vec[0] = vf[0]; vec[1] = vf[1]; vec[2] = vf[2];
416 >                VCOPY(vec, vf);
417                  break;
418          case 'd':                                       /* binary double */
419 <                if (fread((char *)vd, sizeof(double), 3, fp) != 3)
419 >                if (getbinary(vd, sizeof(double), 3, fp) != 3)
420                          return(-1);
421 <                vec[0] = vd[0]; vec[1] = vd[1]; vec[2] = vd[2];
421 >                VCOPY(vec, vd);
422                  break;
423          default:
424                  error(CONSISTENCY, "botched input format");
# Line 422 | Line 427 | getvec(                /* get a vector from fp */
427   }
428  
429  
430 < static void
430 > void
431   tranotify(                      /* record new modifier */
432          OBJECT  obj
433   )
434   {
435          static int  hitlimit = 0;
436 <        register OBJREC  *o = objptr(obj);
437 <        register char  **tralp;
436 >        OBJREC   *o = objptr(obj);
437 >        char  **tralp;
438  
439          if (obj == OVOID) {             /* starting over */
440                  traset[0] = 0;
# Line 456 | Line 461 | ourtrace(                              /* print ray values */
461          RAY  *r
462   )
463   {
464 <        register oputf_t **tp;
464 >        oputf_t **tp;
465  
466          if (every_out[0] == NULL)
467                  return;
# Line 490 | Line 495 | oputo(                         /* print origin */
495          RAY  *r
496   )
497   {
498 <        (*putreal)(r->rorg[0]);
494 <        (*putreal)(r->rorg[1]);
495 <        (*putreal)(r->rorg[2]);
498 >        (*putreal)(r->rorg, 3);
499   }
500  
501  
# Line 501 | Line 504 | oputd(                         /* print direction */
504          RAY  *r
505   )
506   {
507 <        (*putreal)(r->rdir[0]);
505 <        (*putreal)(r->rdir[1]);
506 <        (*putreal)(r->rdir[2]);
507 >        (*putreal)(r->rdir, 3);
508   }
509  
510  
# Line 512 | Line 513 | oputv(                         /* print value */
513          RAY  *r
514   )
515   {
516 +        RREAL   cval[3];
517 +
518          if (outform == 'c') {
519                  COLR  cout;
520                  setcolr(cout,   colval(r->rcol,RED),
521                                  colval(r->rcol,GRN),
522                                  colval(r->rcol,BLU));
523 <                fwrite((char *)cout, sizeof(cout), 1, stdout);
523 >                putbinary(cout, sizeof(cout), 1, stdout);
524                  return;
525          }
526 <        (*putreal)(colval(r->rcol,RED));
527 <        (*putreal)(colval(r->rcol,GRN));
528 <        (*putreal)(colval(r->rcol,BLU));
526 >        cval[0] = colval(r->rcol,RED);
527 >        cval[1] = colval(r->rcol,GRN);
528 >        cval[2] = colval(r->rcol,BLU);
529 >        (*putreal)(cval, 3);
530   }
531  
532  
533   static void
534 + oputV(                          /* print value contribution */
535 +        RAY *r
536 + )
537 + {
538 +        RREAL   contr[3];
539 +
540 +        raycontrib(contr, r, PRIMARY);
541 +        multcolor(contr, r->rcol);
542 +        (*putreal)(contr, 3);
543 + }
544 +
545 +
546 + static void
547   oputl(                          /* print effective distance */
548          RAY  *r
549   )
550   {
551 <        (*putreal)(r->rt);
551 >        (*putreal)(&r->rt, 1);
552   }
553  
554  
# Line 540 | Line 557 | oputL(                         /* print single ray length */
557          RAY  *r
558   )
559   {
560 <        (*putreal)(r->rot);
560 >        (*putreal)(&r->rot, 1);
561   }
562  
563  
# Line 549 | Line 566 | oputc(                         /* print local coordinates */
566          RAY  *r
567   )
568   {
569 <        (*putreal)(r->uv[0]);
553 <        (*putreal)(r->uv[1]);
569 >        (*putreal)(r->uv, 2);
570   }
571  
572  
573 + static RREAL    vdummy[3] = {0.0, 0.0, 0.0};
574 +
575 +
576   static void
577   oputp(                          /* print point */
578          RAY  *r
579   )
580   {
581 <        if (r->rot < FHUGE) {
582 <                (*putreal)(r->rop[0]);
583 <                (*putreal)(r->rop[1]);
584 <                (*putreal)(r->rop[2]);
566 <        } else {
567 <                (*putreal)(0.0);
568 <                (*putreal)(0.0);
569 <                (*putreal)(0.0);
570 <        }
581 >        if (r->rot < FHUGE)
582 >                (*putreal)(r->rop, 3);
583 >        else
584 >                (*putreal)(vdummy, 3);
585   }
586  
587  
# Line 576 | Line 590 | oputN(                         /* print unperturbed normal */
590          RAY  *r
591   )
592   {
593 <        if (r->rot < FHUGE) {
594 <                (*putreal)(r->ron[0]);
595 <                (*putreal)(r->ron[1]);
596 <                (*putreal)(r->ron[2]);
583 <        } else {
584 <                (*putreal)(0.0);
585 <                (*putreal)(0.0);
586 <                (*putreal)(0.0);
587 <        }
593 >        if (r->rot < FHUGE)
594 >                (*putreal)(r->ron, 3);
595 >        else
596 >                (*putreal)(vdummy, 3);
597   }
598  
599  
# Line 596 | Line 605 | oputn(                         /* print perturbed normal */
605          FVECT  pnorm;
606  
607          if (r->rot >= FHUGE) {
608 <                (*putreal)(0.0);
600 <                (*putreal)(0.0);
601 <                (*putreal)(0.0);
608 >                (*putreal)(vdummy, 3);
609                  return;
610          }
611          raynormal(pnorm, r);
612 <        (*putreal)(pnorm[0]);
606 <        (*putreal)(pnorm[1]);
607 <        (*putreal)(pnorm[2]);
612 >        (*putreal)(pnorm, 3);
613   }
614  
615  
# Line 626 | Line 631 | oputw(                         /* print weight */
631          RAY  *r
632   )
633   {
634 <        (*putreal)(r->rweight);
634 >        RREAL   rwt = r->rweight;
635 >        
636 >        (*putreal)(&rwt, 1);
637   }
638  
639  
640   static void
641 < oputW(                          /* print contribution */
641 > oputW(                          /* print coefficient */
642          RAY  *r
643   )
644   {
645 <        COLOR   contr;
645 >        RREAL   contr[3];
646 >                                /* shadow ray not on source? */
647 >        if (r->rsrc >= 0 && source[r->rsrc].so != r->ro)
648 >                setcolor(contr, 0.0, 0.0, 0.0);
649 >        else
650 >                raycontrib(contr, r, PRIMARY);
651  
652 <        raycontrib(contr, r, PRIMARY);
641 <        (*putreal)(colval(contr,RED));
642 <        (*putreal)(colval(contr,GRN));
643 <        (*putreal)(colval(contr,BLU));
652 >        (*putreal)(contr, 3);
653   }
654  
655  
# Line 688 | Line 697 | oputtilde(                     /* output tilde (spacer) */
697  
698  
699   static void
700 < puta(                           /* print ascii value */
701 <        double  v
700 > puta(                           /* print ascii value(s) */
701 >        RREAL *v, int n
702   )
703   {
704 <        printf("%e\t", v);
704 >        if (n == 3) {
705 >                printf("%e\t%e\t%e\t", v[0], v[1], v[2]);
706 >                return;
707 >        }
708 >        while (n--)
709 >                printf("%e\t", *v++);
710   }
711  
712  
713   static void
714 < putd(v)                         /* print binary double */
701 < double  v;
714 > putd(RREAL *v, int n)           /* print binary double(s) */
715   {
716 <        fwrite((char *)&v, sizeof(v), 1, stdout);
716 > #ifdef  SMLFLT
717 >        double  da[3];
718 >        int     i;
719 >
720 >        if (n > 3)
721 >                error(INTERNAL, "code error in putd()");
722 >        for (i = n; i--; )
723 >                da[i] = v[i];
724 >        putbinary(da, sizeof(double), n, stdout);
725 > #else
726 >        putbinary(v, sizeof(RREAL), n, stdout);
727 > #endif
728   }
729  
730  
731   static void
732 < putf(v)                         /* print binary float */
709 < double  v;
732 > putf(RREAL *v, int n)           /* print binary float(s) */
733   {
734 <        float f = v;
734 > #ifndef SMLFLT
735 >        float   fa[3];
736 >        int     i;
737  
738 <        fwrite((char *)&f, sizeof(f), 1, stdout);
738 >        if (n > 3)
739 >                error(INTERNAL, "code error in putf()");
740 >        for (i = n; i--; )
741 >                fa[i] = v[i];
742 >        putbinary(fa, sizeof(float), n, stdout);
743 > #else
744 >        putbinary(v, sizeof(RREAL), n, stdout);
745 > #endif
746   }

Diff Legend

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