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

Comparing ray/src/hd/rholo.c (file contents):
Revision 3.45 by gwlarson, Tue Jan 5 11:34:52 1999 UTC vs.
Revision 3.86 by greg, Thu Jun 5 19:10:31 2025 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1998 Silicon Graphics, Inc. */
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ SGI";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6 –
4   /*
5   * Radiance holodeck generation controller
6   */
7  
8 < #include "rholo.h"
9 < #include "random.h"
8 > #include <stdio.h>
9 > #include <time.h>
10   #include <signal.h>
14 – #include <sys/types.h>
11   #include <sys/stat.h>
12 + #include <string.h>
13  
14 + #include "platform.h"
15 + #include "rterror.h"
16 + #include "resolu.h"
17 + #include "rholo.h"
18 + #include "random.h"
19 +
20   #ifndef FRAGWARN
21   #define FRAGWARN        20              /* fragmentation for warning (%) */
22   #endif
23   #ifndef MAXQTIME
24 < #define MAXQTIME        5               /* target maximum seconds in queue */
24 > #define MAXQTIME        3               /* target maximum seconds in queue */
25   #endif
26 +                                        /* manual cache flushing frequency */
27 + #ifndef RTFLUSH
28 + #if MAXQTIME
29 + #define RTFLUSH         (300/MAXQTIME*totqlen)  /* <= 5 minutes */
30 + #else
31 + #define RTFLUSH         (50*totqlen)            /* just guess */
32 + #endif
33 + #endif
34                          /* the following must be consistent with rholo.h */
35   int     NVARS = NRHVARS;                /* total number of variables */
36  
37   VARIABLE        vv[] = RHVINIT;         /* variable-value pairs */
38  
28 – char    *progname;              /* our program name */
39   char    *hdkfile;               /* holodeck file name */
40   char    froot[256];             /* root file name */
41  
# Line 35 | Line 45 | char   *outdev = NULL;         /* output device name */
45  
46   int     readinp = 0;            /* read commands from stdin */
47  
48 < int     force = 0;              /* allow overwrite of holodeck */
48 > int     force = 0;              /* allow overwrite of holodeck (-1 == read-only) */
49  
50   time_t  starttime;              /* time we got started */
51   time_t  endtime;                /* time we should end by */
52   time_t  reporttime;             /* time for next report */
53  
54 < long    maxdisk;                /* maximum file space (bytes) */
54 > off_t   maxdisk;                /* maximum file space (bytes) */
55  
56   int     rtargc = 1;             /* rtrace command */
57 < char    *rtargv[128] = {"rtrace", NULL};
57 > char    *rtargv[MAXRTARGC] = {"rtrace", NULL};
58  
59   int     orig_mode = -1;         /* original file mode (-1 if unchanged) */
60  
# Line 52 | Line 62 | long   nraysdone = 0L;         /* number of rays done */
62   long    npacksdone = 0L;        /* number of packets done */
63  
64   PACKET  *freepacks;             /* available packets */
65 < int     totqlen;                /* maximum queue length when full */
65 > int     totqlen;                /* maximum queue length (number of packets) */
66  
67   char  *sigerr[NSIG];            /* signal error messages */
68  
69   extern int      nowarn;         /* turn warnings off? */
70  
71 < extern time_t   time();
71 > static void onsig(int  signo);
72 > static void sigdie(int  signo, char  *msg);
73 > static int resfmode(int fd, int mod);
74 > static void initrholo(void);
75 > static int rholo(void);
76 > static void setdefaults(HDGRID  *gp);
77 > static void creatholo(HDGRID    *gp);
78 > static gethfunc headline;
79 > static void loadholo(void);
80 > static void rootname(char       *rn, char       *fn);
81 > static void badvalue(int        vc);
82  
83  
84 < main(argc, argv)
85 < int     argc;
86 < char    *argv[];
84 > int
85 > main(
86 >        int     argc,
87 >        char    *argv[]
88 > )
89   {
90          int     i;
91  
92 <        initurand(16384);                       /* initialize urand */
71 <        progname = argv[0];                     /* get arguments */
92 >        fixargv0(argv[0]);                      /* get arguments */
93          for (i = 1; i < argc && argv[i][0] == '-'; i++)
94                  switch (argv[i][1]) {
95                  case 'w':                       /* turn off warnings */
96                          nowarn++;
97                          break;
98                  case 'f':                       /* force overwrite */
99 <                        force++;
99 >                        force = 1;
100                          break;
101 +                case 'r':                       /* read-only mode */
102 +                        force = -1;
103 +                        break;
104                  case 'i':                       /* read input from stdin */
105                          readinp++;
106                          break;
# Line 119 | Line 143 | char   *argv[];
143                  HDGRID  hdg[HDMAX];
144                                                          /* set defaults */
145                  setdefaults(hdg);
146 +                                                        /* check read-only */
147 +                if (force < 0)
148 +                        error(USER, "cannot create read-only holodeck");
149                                                          /* holodeck exists? */
150                  if (!force && access(hdkfile, R_OK|W_OK) == 0)
151                          error(USER,
# Line 136 | Line 163 | char   *argv[];
163          quit(0);
164   userr:
165          fprintf(stderr,
166 < "Usage: %s [-n nprocs][-o disp][-w][-f] output.hdk [control.hif|+|- [VAR=val ..]]\n",
166 > "Usage: %s [-n nprocs][-o disp][-i][-w][-r|-f] output.hdk [control.hif|+|- [VAR=val ..]]\n",
167                          progname);
168          quit(1);
169 +        return 1; /* pro forma return */
170   }
171  
172  
173 < onsig(signo)                            /* fatal signal */
174 < int  signo;
173 > static void
174 > onsig(                          /* fatal signal */
175 >        int  signo
176 > )
177   {
178          static int  gotsig = 0;
179  
180 <        if (gotsig++)                   /* two signals and we're gone! */
180 >        if (gotsig > 1)                 /* we're going as fast as we can! */
181 >                return;
182 >        if (gotsig++) {                 /* two signals and we split */
183 >                hdsync(NULL, 0);        /* don't leave w/o saying goodbye */
184                  _exit(signo);
185 <
186 <        alarm(120);                     /* allow 2 minutes to clean up */
154 <        signal(SIGALRM, SIG_DFL);       /* make certain we do die */
185 >        }
186 >        alarm(300);                     /* allow 10 minutes to clean up */
187          eputs("signal - ");
188          eputs(sigerr[signo]);
189          eputs("\n");
# Line 159 | Line 191 | int  signo;
191   }
192  
193  
194 < sigdie(signo, msg)                      /* set fatal signal */
195 < int  signo;
196 < char  *msg;
194 > static void
195 > sigdie(                 /* set fatal signal */
196 >        int  signo,
197 >        char  *msg
198 > )
199   {
200          if (signal(signo, onsig) == SIG_IGN)
201                  signal(signo, SIG_IGN);
# Line 169 | Line 203 | char  *msg;
203   }
204  
205  
206 < int
207 < resfmode(fd, mod)               /* restrict open file access mode */
208 < int     fd, mod;
206 > static int
207 > resfmode(               /* restrict open file access mode */
208 >        int     fd,
209 >        int     mod
210 > )
211   {
212          struct stat     stbuf;
213                                          /* get original mode */
# Line 189 | Line 225 | int    fd, mod;
225   }
226  
227  
228 < initrholo()                     /* get our holodeck running */
228 > static void
229 > initrholo(void)                 /* get our holodeck running */
230   {
231 <        extern int      global_packet();
195 <        register int    i;
231 >        int     i;
232                                                  /* close holodeck on exec() */
233          fcntl(hdlist[0]->fd, F_SETFD, FD_CLOEXEC);
234  
# Line 202 | Line 238 | initrholo()                    /* get our holodeck running */
238                  init_global();
239                                                  /* record disk space limit */
240          if (!vdef(DISKSPACE))
241 <                maxdisk = 0;
241 >                maxdisk = ((off_t)1<<(sizeof(off_t)*8-2)) - 1024;
242          else
243                  maxdisk = 1024.*1024.*vflt(DISKSPACE);
244                                                  /* set up memory cache */
245          if (outdev == NULL)
246 <                hdcachesize = 0;        /* manual flushing */
246 >                hdcachesize = 0;                /* manual flushing */
247          else if (vdef(CACHE))
248                  hdcachesize = 1024.*1024.*vflt(CACHE);
249                                                  /* open report file */
250          if (vdef(REPORT)) {
251 <                register char   *s = sskip2(vval(REPORT), 1);
251 >                char    *s = sskip2(vval(REPORT), 1);
252                  if (*s && freopen(s, "a", stderr) == NULL)
253                          quit(2);
254          }
# Line 254 | Line 290 | initrholo()                    /* get our holodeck running */
290          }
291                                          /* set up signal handling */
292          sigdie(SIGINT, "Interrupt");
257 –        sigdie(SIGHUP, "Hangup");
293          sigdie(SIGTERM, "Terminate");
294 + #ifdef SIGHUP
295 +        sigdie(SIGHUP, "Hangup");
296 + #endif
297 + #ifdef SIGPIPE
298          sigdie(SIGPIPE, "Broken pipe");
299 + #endif
300 + #ifdef SIGALRM
301          sigdie(SIGALRM, "Alarm clock");
302 + #endif
303   #ifdef  SIGXCPU
304          sigdie(SIGXCPU, "CPU limit exceeded");
305 + #endif
306 + #ifdef  SIGXFSZ
307          sigdie(SIGXFSZ, "File size exceeded");
308   #endif
309 <                                                /* protect holodeck file */
310 <        orig_mode = resfmode(hdlist[0]->fd, ncprocs>0 ? 0 : 0444);
309 >                                        /* protect holodeck file */
310 >        orig_mode = resfmode(hdlist[0]->fd, (ncprocs>0) & (force>=0) ? 0 : 0444);
311          return;
312   memerr:
313          error(SYSTEM, "out of memory in initrholo");
314   }
315  
316  
317 < rholo()                         /* holodeck main loop */
317 > static int
318 > rholo(void)                             /* holodeck main loop */
319   {
320 <        static long     nextfragwarn = 100*(1L<<20);
320 >        static off_t    nextfragwarn = 100L<<20;
321          static int      idle = 0;
322          PACKET  *pl = NULL, *plend;
323 <        long    fsiz;
323 >        off_t   fsiz;
324          int     pksiz;
325 <        register PACKET *p;
325 >        PACKET  *p;
326          time_t  t;
327                                          /* check display */
328          if (nprocs <= 0)
# Line 291 | Line 336 | rholo()                                /* holodeck main loop */
336                  return(0);              /* all done */
337          fsiz = hdfilen(hdlist[0]->fd);  /* check file size */
338          if (maxdisk > 0 && fsiz >= maxdisk) {
339 <                error(WARNING, "file limit exceeded");
339 >                error(USER, "file limit exceeded");
340                  done_rtrace();
341                  return(1);      /* comes back */
342          }
343   #if FRAGWARN
344 <        if (fsiz >= nextfragwarn &&
345 <                (fsiz-hdfiluse(hdlist[0]->fd,0))/(fsiz/100) > FRAGWARN) {
346 <                sprintf(errmsg, "holodeck file fragmentation is %.0f%%",
347 <                                100.*(fsiz-hdfiluse(hdlist[0]->fd,1))/fsiz);
348 <                error(WARNING, errmsg);
349 <                nextfragwarn = fsiz + (fsiz>>2);        /* decent interval */
344 >        if (fsiz >= nextfragwarn) {
345 >                double  pctfrag = 100.*(fsiz-hdfiluse(hdlist[0]->fd))/fsiz;
346 >                if (pctfrag >= (double)FRAGWARN) {
347 >                        sprintf(errmsg, "holodeck file fragmentation is %.0f%%",
348 >                                        pctfrag);
349 >                        error(WARNING, errmsg);
350 >                        nextfragwarn = fsiz + (fsiz>>2);
351 >                } else
352 >                        nextfragwarn = fsiz + (10L<<20);
353          }
354   #endif
355          t = time(NULL);                 /* check time */
356          if (endtime > 0 && t >= endtime) {
357 <                error(WARNING, "time limit exceeded");
357 >                error(USER, "time limit exceeded");
358                  done_rtrace();
359                  return(1);      /* comes back */
360          }
361          if (reporttime > 0 && t >= reporttime)
362                  report(t);
363                                          /* figure out good packet size */
364 +        pksiz = RPACKSIZ;
365   #if MAXQTIME
366 <        pksiz = nraysdone*MAXQTIME/(totqlen*(t - starttime + 1L));
367 <        if (pksiz < 1)
368 <                pksiz = 1;
369 <        else if (pksiz > RPACKSIZ)
366 >        if (!chunkycmp) {
367 >                pksiz = nraysdone*MAXQTIME/(totqlen*(t - starttime + 1L));
368 >                if (pksiz < 1) pksiz = 1;
369 >                else if (pksiz > RPACKSIZ) pksiz = RPACKSIZ;
370 >        }
371   #endif
322 –                pksiz = RPACKSIZ;
372          idle = 0;                       /* get packets to process */
373          while (freepacks != NULL) {
374                  p = freepacks; freepacks = p->next; p->next = NULL;
# Line 338 | Line 387 | rholo()                                /* holodeck main loop */
387   }
388  
389  
390 < setdefaults(gp)                 /* set default values */
391 < register HDGRID *gp;
390 > static void
391 > setdefaults(                    /* set default values */
392 >        HDGRID  *gp
393 > )
394   {
395 <        extern char     *atos();
345 <        register int    i;
395 >        int     i;
396          int     n;
397          double  len[3], d;
348 –        char    buf[64];
398  
399          if (!vdef(SECTION)) {
400                  sprintf(errmsg, "%s must be defined", vnam(SECTION));
# Line 363 | Line 412 | register HDGRID        *gp;
412          }
413                                  /* append rendering options */
414          if (vdef(RENDER))
415 <                rtargc += wordstring(rtargv+rtargc, vval(RENDER));
415 >                rtargc += wordstring(rtargv+rtargc, MAXRTARGC-rtargc, vval(RENDER));
416          
417          if (gp == NULL)         /* already initialized? */
418                  return;
# Line 395 | Line 444 | register HDGRID        *gp;
444   }
445  
446  
447 < creatholo(gp)                   /* create a holodeck output file */
448 < HDGRID  *gp;
447 > static void
448 > creatholo(                      /* create a holodeck output file */
449 >        HDGRID  *gp
450 > )
451   {
452          extern char     VersionID[];
453 <        int4    lastloc, nextloc;
453 >        off_t   lastloc, nextloc;
454          int     n;
455          int     fd;
456          FILE    *fp;
457                                          /* open & truncate file */
458 <        if ((fp = fopen(hdkfile, "w+")) == NULL) {
458 >        if ((fp = fopen(hdkfile, "wb+")) == NULL) {
459                  sprintf(errmsg, "cannot open \"%s\" for writing", hdkfile);
460                  error(SYSTEM, errmsg);
461          }
462                                          /* write information header */
463          newheader("RADIANCE", fp);
464          fprintf(fp, "SOFTWARE= %s\n", VersionID);
465 +        fputnow(fp);
466          printvars(fp);
467 +        fputendian(fp);
468          fputformat(HOLOFMT, fp);
469          fputc('\n', fp);
470          putw(HOLOMAGIC, fp);            /* put magic number */
471          fd = dup(fileno(fp));
472          fclose(fp);                     /* flush and close stdio stream */
473 <        lastloc = lseek(fd, 0L, 2);
473 >        lastloc = lseek(fd, (off_t)0, SEEK_END);
474          for (n = vdef(SECTION); n--; gp++) {    /* initialize each section */
475                  nextloc = 0L;
476                  write(fd, (char *)&nextloc, sizeof(nextloc));
# Line 425 | Line 478 | HDGRID *gp;
478                  if (!n)
479                          break;
480                  nextloc = hdfilen(fd);          /* write section pointer */
481 <                if (lseek(fd, (long)lastloc, 0) < 0)
481 >                if (lseek(fd, lastloc, SEEK_SET) < 0)
482                          error(SYSTEM,
483                                  "cannot seek on holodeck file in creatholo");
484                  write(fd, (char *)&nextloc, sizeof(nextloc));
485 <                lseek(fd, (long)(lastloc=nextloc), 0);
485 >                lseek(fd, (lastloc=nextloc), SEEK_SET);
486          }
487   }
488  
489  
490 < int
491 < headline(s)                     /* process information header line */
492 < char    *s;
490 > static int
491 > headline(                       /* process information header line */
492 >        char    *s,
493 >        void    *p
494 > )
495   {
496          extern char     FMTSTR[];
497 <        register char   *cp;
498 <        char    fmt[32];
497 >        char    *cp;
498 >        char    fmt[MAXFMTLEN];
499  
500          if (formatval(fmt, s)) {
501                  if (strcmp(fmt, HOLOFMT)) {
# Line 460 | Line 515 | char   *s;
515   }
516  
517  
518 < loadholo()                      /* start loading a holodeck from fname */
518 > static void
519 > loadholo(void)                  /* start loading a holodeck from fname */
520   {
465 –        extern long     ftell();
521          FILE    *fp;
522          int     fd;
523          int     n;
524 <        int4    nextloc;
525 <                                        /* open holodeck file */
526 <        if ((fp = fopen(hdkfile, ncprocs>0 ? "r+" : "r")) == NULL) {
527 <                sprintf(errmsg, "cannot %s \"%s\"",
528 <                                ncprocs>0 ? "append" : "read", hdkfile);
529 <                error(SYSTEM, errmsg);
524 >        off_t   nextloc;
525 >        
526 >        if ((ncprocs > 0) & (force >= 0))
527 >                fp = fopen(hdkfile, "rb+");
528 >        else
529 >                fp = NULL;
530 >        if (fp == NULL) {
531 >                if ((fp = fopen(hdkfile, "rb")) == NULL) {
532 >                        sprintf(errmsg, "cannot open \"%s\"", hdkfile);
533 >                        error(SYSTEM, errmsg);
534 >                }
535 >                if (ncprocs > 0) {
536 >                        sprintf(errmsg, "\"%s\" is read-only", hdkfile);
537 >                        if (outdev == NULL)
538 >                                error(USER, errmsg);
539 >                        strcat(errmsg, "; new rays will be discarded");
540 >                        error(WARNING, errmsg);
541 >                        force = -1;
542 >                }
543          }
544                                          /* load variables from header */
545          getheader(fp, headline, NULL);
# Line 485 | Line 553 | loadholo()                     /* start loading a holodeck from fname */
553          fd = dup(fileno(fp));
554          fclose(fp);                             /* done with stdio */
555          for (n = 0; nextloc > 0L; n++) {        /* initialize each section */
556 <                lseek(fd, (long)nextloc, 0);
556 >                lseek(fd, nextloc, SEEK_SET);
557                  read(fd, (char *)&nextloc, sizeof(nextloc));
558                  hdinit(fd, NULL);
559          }
# Line 497 | Line 565 | loadholo()                     /* start loading a holodeck from fname */
565   }
566  
567  
568 < done_packets(pl)                /* handle finished packets */
569 < PACKET  *pl;
568 > void
569 > done_packets(           /* handle finished packets */
570 >        PACKET  *pl
571 > )
572   {
573          static int      n2flush = 0;
574 <        register PACKET *p;
574 >        PACKET  *p;
575  
576          while (pl != NULL) {
577                  p = pl; pl = p->next; p->next = NULL;
578                  if (p->nr > 0) {                /* add to holodeck */
579 <                        bcopy((char *)p->ra,
580 <                                (char *)hdnewrays(hdlist[p->hd],p->bi,p->nr),
579 >                        memcpy( (void *)hdnewrays(hdlist[p->hd],p->bi,p->nr),
580 >                                (void *)p->ra,
581                                  p->nr*sizeof(RAYVAL));
582                          if (outdev != NULL)     /* display it */
583                                  disp_packet((PACKHEAD *)p);
584 <                        if (hdcachesize <= 0)   /* manual flushing */
584 >                        if (hdcachesize <= 0)
585                                  n2flush++;
586                          nraysdone += p->nr;
587                          npacksdone++;
# Line 520 | Line 590 | PACKET *pl;
590                  p->next = freepacks;            /* push onto free list */
591                  freepacks = p;
592          }
593 < #if MAXQTIME
594 <        if (n2flush > 300/MAXQTIME*totqlen) {
595 < #else
526 <        if (n2flush > 50*totqlen) {
527 < #endif
528 <                if (outdev == NULL)
529 <                        hdflush(NULL);          /* flush holodeck buffers */
593 >        if (n2flush >= RTFLUSH) {
594 >                if (outdev != NULL)
595 >                        hdsync(NULL, 1);
596                  else
597 <                        hdsync(NULL, 1);        /* sync holodeck file */
597 >                        hdflush(NULL);
598                  n2flush = 0;
599          }
600   }
601  
602  
603 < rootname(rn, fn)                /* remove tail from end of fn */
604 < register char   *rn, *fn;
603 > static void
604 > rootname(               /* remove tail from end of fn */
605 >        char    *rn,
606 >        char    *fn
607 > )
608   {
609          char    *tp, *dp;
610  
611 <        for (tp = NULL, dp = rn; *rn = *fn++; rn++)
611 >        for (tp = NULL, dp = rn; (*rn = *fn++); rn++) {
612                  if (*rn == '/')
613                          dp = rn;
614                  else if (*rn == '.')
615                          tp = rn;
616 +        }
617          if (tp != NULL && tp > dp)
618                  *tp = '\0';
619   }
620  
621  
622 < badvalue(vc)                    /* report bad variable value and exit */
623 < int     vc;
622 > static void
623 > badvalue(                       /* report bad variable value and exit */
624 >        int     vc
625 > )
626   {
627          sprintf(errmsg, "bad value for variable '%s'", vnam(vc));
628          error(USER, errmsg);
629   }
630  
631  
632 < eputs(s)                        /* put error message to stderr */
633 < register char  *s;
632 > void
633 > eputs(const char *s)            /* put error message to stderr */
634   {
635          static int  midline = 0;
636  
# Line 576 | Line 648 | register char  *s;
648   }
649  
650  
651 < quit(ec)                        /* exit program gracefully */
652 < int     ec;
651 > void
652 > quit(int ec)                    /* exit program gracefully */
653   {
654          int     status = 0;
655  
656          if (hdlist[0] != NULL) {        /* close holodeck */
657                  if (nprocs > 0)
658                          status = done_rtrace();         /* calls hdsync() */
659 <                if (ncprocs > 0 && vdef(REPORT)) {
660 <                        long    fsiz, fuse;
659 >                if ((ncprocs > 0) & (force >= 0) && vdef(REPORT)) {
660 >                        off_t   fsiz, fuse;
661                          fsiz = hdfilen(hdlist[0]->fd);
662 <                        fuse = hdfiluse(hdlist[0]->fd, 1);
662 >                        fuse = hdfiluse(hdlist[0]->fd);
663                          fprintf(stderr,
664                          "%s: %.1f Mbyte holodeck file, %.1f%% fragmentation\n",
665                                          hdkfile, fsiz/(1024.*1024.),

Diff Legend

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