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

Comparing src/hd/rholo.c (file contents):
Revision 3.14 by gregl, Mon Dec 1 16:34:21 1997 UTC vs.
Revision 3.75 by greg, Thu Sep 30 15:43:30 2010 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1997 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"
13 < #include "paths.h"
8 > #include <stdio.h>
9 > #include <time.h>
10   #include <signal.h>
15 – #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        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  
# Line 22 | Line 38 | VARIABLE       vv[] = RHVINIT;         /* variable-value pairs */
38  
39   char    *progname;              /* our program name */
40   char    *hdkfile;               /* holodeck file name */
41 < char    froot[MAXPATH];         /* root file name */
41 > char    froot[256];             /* root file name */
42  
27 – int     nowarn = 0;             /* turn warnings off? */
28 –
29 – double  expval = 1.;            /* global exposure value */
30 –
43   int     ncprocs = 0;            /* desired number of compute processes */
44  
45   char    *outdev = NULL;         /* output device name */
46  
47 + int     readinp = 0;            /* read commands from stdin */
48 +
49 + int     force = 0;              /* allow overwrite of holodeck (-1 == read-only) */
50 +
51   time_t  starttime;              /* time we got started */
52   time_t  endtime;                /* time we should end by */
53   time_t  reporttime;             /* time for next report */
54  
55 < long    maxdisk;                /* maximum file space (bytes) */
55 > off_t   maxdisk;                /* maximum file space (bytes) */
56  
57   int     rtargc = 1;             /* rtrace command */
58   char    *rtargv[128] = {"rtrace", NULL};
# Line 47 | Line 63 | long   nraysdone = 0L;         /* number of rays done */
63   long    npacksdone = 0L;        /* number of packets done */
64  
65   PACKET  *freepacks;             /* available packets */
66 + int     totqlen;                /* maximum queue length (number of packets) */
67  
68   char  *sigerr[NSIG];            /* signal error messages */
69  
70 < extern time_t   time();
70 > extern int      nowarn;         /* turn warnings off? */
71  
72 + static void onsig(int  signo);
73 + static void sigdie(int  signo, char  *msg);
74 + static int resfmode(int fd, int mod);
75 + static void initrholo(void);
76 + static int rholo(void);
77 + static void setdefaults(HDGRID  *gp);
78 + static void creatholo(HDGRID    *gp);
79 + static gethfunc headline;
80 + static void loadholo(void);
81 + static void rootname(char       *rn, char       *fn);
82 + static void badvalue(int        vc);
83  
84 < main(argc, argv)
85 < int     argc;
86 < char    *argv[];
84 >
85 > int
86 > main(
87 >        int     argc,
88 >        char    *argv[]
89 > )
90   {
60 –        HDGRID  hdg;
91          int     i;
92 <        int     force = 0;
63 <                                                /* mark start time */
64 <        starttime = time(NULL);
65 <        initurand(10240);                       /* initialize urand */
92 >
93          progname = argv[0];                     /* get arguments */
94          for (i = 1; i < argc && argv[i][0] == '-'; i++)
95                  switch (argv[i][1]) {
# Line 70 | Line 97 | char   *argv[];
97                          nowarn++;
98                          break;
99                  case 'f':                       /* force overwrite */
100 <                        force++;
100 >                        force = 1;
101                          break;
102 +                case 'r':                       /* read-only mode */
103 +                        force = -1;
104 +                        break;
105 +                case 'i':                       /* read input from stdin */
106 +                        readinp++;
107 +                        break;
108                  case 'n':                       /* compute processes */
109                          if (i >= argc-2)
110                                  goto userr;
# Line 85 | Line 118 | char   *argv[];
118                  default:
119                          goto userr;
120                  }
88 –                                                /* do we have a job? */
89 –        if (outdev == NULL && ncprocs <= 0)
90 –                goto userr;
121                                                  /* get root file name */
122 +        if (i >= argc)
123 +                goto userr;
124          rootname(froot, hdkfile=argv[i++]);
125 <                                                /* load... */
126 <        if (i < argc) {                         /* variables */
127 <                loadvars(argv[i++]);
128 <                                                        /* cmdline settings */
129 <                for ( ; i < argc; i++)
130 <                        if (setvariable(argv[i], matchvar) < 0) {
131 <                                sprintf(errmsg, "unknown variable: %s",
132 <                                                argv[i]);
133 <                                error(USER, errmsg);
134 <                        }
135 <                                                        /* check settings */
136 <                checkvalues();
137 <                                                        /* load RIF if any */
138 <                if (vdef(RIF))
139 <                        getradfile(vval(RIF));
125 >                                                /* load variables? */
126 >        if (i < argc)
127 >                if (argv[i][0] != '-' && argv[i][0] != '+')
128 >                        loadvars(argv[i]);      /* load variables from file */
129 >
130 >        if (i >= argc || argv[i][0] == '+')
131 >                loadholo();                     /* load existing holodeck */
132 >
133 >        while (++i < argc)                      /* get command line settings */
134 >                if (setvariable(argv[i], matchvar) < 0) {
135 >                        sprintf(errmsg, "unknown variable: %s", argv[i]);
136 >                        error(USER, errmsg);
137 >                }
138 >                                                /* check settings */
139 >        checkvalues();
140 >                                                /* load rad input file */
141 >        getradfile();
142 >
143 >        if (hdlist[0] == NULL) {                /* create new holodeck */
144 >                HDGRID  hdg[HDMAX];
145                                                          /* set defaults */
146 <                setdefaults(&hdg);
146 >                setdefaults(hdg);
147 >                                                        /* check read-only */
148 >                if (force < 0)
149 >                        error(USER, "cannot create read-only holodeck");
150                                                          /* holodeck exists? */
151                  if (!force && access(hdkfile, R_OK|W_OK) == 0)
152                          error(USER,
153                                  "holodeck file exists -- use -f to overwrite");
154                                                          /* create holodeck */
155 <                creatholo(&hdg);
156 <        } else {                                /* else load holodeck */
157 <                loadholo();
118 <                if (vdef(RIF))                          /* load RIF if any */
119 <                        getradfile(vval(RIF));
120 <        }
155 >                creatholo(hdg);
156 >        } else                                  /* else just set defaults */
157 >                setdefaults(NULL);
158                                                  /* initialize */
159          initrholo();
160 <                                                /* run */
160 >                                                /* main loop */
161          while (rholo())
162                  ;
163                                                  /* done */
164          quit(0);
165   userr:
166          fprintf(stderr,
167 < "Usage: %s {-n nprocs|-o disp} [-w][-f] output.hdk [control.hif [VAR=val ..]]\n",
167 > "Usage: %s [-n nprocs][-o disp][-i][-w][-r|-f] output.hdk [control.hif|+|- [VAR=val ..]]\n",
168                          progname);
169          quit(1);
170 +        return 1; /* pro forma return */
171   }
172  
173  
174 < onsig(signo)                            /* fatal signal */
175 < int  signo;
174 > static void
175 > onsig(                          /* fatal signal */
176 >        int  signo
177 > )
178   {
179          static int  gotsig = 0;
180  
181 <        if (gotsig++)                   /* two signals and we're gone! */
181 >        if (gotsig > 1)                 /* we're going as fast as we can! */
182 >                return;
183 >        if (gotsig++) {                 /* two signals and we split */
184 >                hdsync(NULL, 0);        /* don't leave w/o saying goodbye */
185                  _exit(signo);
186 <
187 <        alarm(30);                      /* allow 30 seconds to clean up */
145 <        signal(SIGALRM, SIG_DFL);       /* make certain we do die */
186 >        }
187 >        alarm(300);                     /* allow 5 minutes to clean up */
188          eputs("signal - ");
189          eputs(sigerr[signo]);
190          eputs("\n");
# Line 150 | Line 192 | int  signo;
192   }
193  
194  
195 < sigdie(signo, msg)                      /* set fatal signal */
196 < int  signo;
197 < char  *msg;
195 > static void
196 > sigdie(                 /* set fatal signal */
197 >        int  signo,
198 >        char  *msg
199 > )
200   {
201          if (signal(signo, onsig) == SIG_IGN)
202                  signal(signo, SIG_IGN);
# Line 160 | Line 204 | char  *msg;
204   }
205  
206  
207 < int
208 < resfmode(fd, mod)               /* restrict open file access mode */
209 < int     fd, mod;
207 > static int
208 > resfmode(               /* restrict open file access mode */
209 >        int     fd,
210 >        int     mod
211 > )
212   {
213          struct stat     stbuf;
214                                          /* get original mode */
215          if (fstat(fd, &stbuf) < 0)
216                  error(SYSTEM, "cannot stat open holodeck file");
217          mod &= stbuf.st_mode;           /* always more restrictive */
218 <        if (mod == stbuf.st_mode)
218 >        if (mod == (stbuf.st_mode & 0777))
219                  return(-1);             /* already set */
220                                          /* else change it */
221          if (fchmod(fd, mod) < 0) {
# Line 180 | Line 226 | int    fd, mod;
226   }
227  
228  
229 < initrholo()                     /* get our holodeck running */
229 > static void
230 > initrholo(void)                 /* get our holodeck running */
231   {
232          extern int      global_packet();
233          register int    i;
234 +                                                /* close holodeck on exec() */
235 +        fcntl(hdlist[0]->fd, F_SETFD, FD_CLOEXEC);
236  
237          if (outdev != NULL)                     /* open output device */
238                  disp_open(outdev);
# Line 191 | Line 240 | initrholo()                    /* get our holodeck running */
240                  init_global();
241                                                  /* record disk space limit */
242          if (!vdef(DISKSPACE))
243 <                maxdisk = 0;
243 >                maxdisk = ((off_t)1<<(sizeof(off_t)*8-2)) - 1024;
244          else
245                  maxdisk = 1024.*1024.*vflt(DISKSPACE);
197 –                                                /* record end time */
198 –        if (!vdef(TIME) || vflt(TIME) <= FTINY)
199 –                endtime = 0;
200 –        else
201 –                endtime = starttime + vflt(TIME)*3600.;
246                                                  /* set up memory cache */
247          if (outdev == NULL)
248 <                hdcachesize = 0;        /* manual flushing */
248 >                hdcachesize = 0;                /* manual flushing */
249          else if (vdef(CACHE))
250                  hdcachesize = 1024.*1024.*vflt(CACHE);
251                                                  /* open report file */
# Line 210 | Line 254 | initrholo()                    /* get our holodeck running */
254                  if (*s && freopen(s, "a", stderr) == NULL)
255                          quit(2);
256          }
257 +                                                /* mark the starting time */
258 +        starttime = time(NULL);
259 +                                                /* compute end time */
260 +        if (!vdef(TIME) || vflt(TIME) <= FTINY)
261 +                endtime = 0;
262 +        else
263 +                endtime = starttime + vflt(TIME)*3600. + .5;
264                                                  /* start rtrace */
265          if (ncprocs > 0) {
266 <                i = start_rtrace();
266 >                totqlen = i = start_rtrace();
267                  if (i < 1)
268 <                        error(USER, "cannot start rtrace process");
268 >                        error(USER, "cannot start rtrace process(es)");
269                  if (vdef(REPORT)) {             /* make first report */
270                          printargs(rtargc, rtargv, stderr);
271                          report(0);
# Line 225 | Line 276 | initrholo()                    /* get our holodeck running */
276                          goto memerr;
277                  freepacks[--i].nr = 0;
278                  freepacks[i].next = NULL;
279 <                if (!vbool(OBSTRUCTIONS)) {
279 >                if (!vdef(OBSTRUCTIONS) || !vbool(OBSTRUCTIONS)) {
280                          freepacks[i].offset = (float *)bmalloc(
281                                          RPACKSIZ*sizeof(float)*(i+1) );
282                          if (freepacks[i].offset == NULL)
# Line 241 | Line 292 | initrholo()                    /* get our holodeck running */
292          }
293                                          /* set up signal handling */
294          sigdie(SIGINT, "Interrupt");
244 –        sigdie(SIGHUP, "Hangup");
295          sigdie(SIGTERM, "Terminate");
296 + #ifdef SIGHUP
297 +        sigdie(SIGHUP, "Hangup");
298 + #endif
299 + #ifdef SIGPIPE
300          sigdie(SIGPIPE, "Broken pipe");
301 + #endif
302 + #ifdef SIGALRM
303          sigdie(SIGALRM, "Alarm clock");
304 + #endif
305   #ifdef  SIGXCPU
306          sigdie(SIGXCPU, "CPU limit exceeded");
307 + #endif
308 + #ifdef  SIGXFSZ
309          sigdie(SIGXFSZ, "File size exceeded");
310   #endif
311 <                                                /* protect holodeck file */
312 <        orig_mode = resfmode(hdlist[0]->fd, ncprocs>0 ? 0 : 0444);
311 >                                        /* protect holodeck file */
312 >        orig_mode = resfmode(hdlist[0]->fd, (ncprocs>0) & (force>=0) ? 0 : 0444);
313          return;
314   memerr:
315          error(SYSTEM, "out of memory in initrholo");
316   }
317  
318  
319 < rholo()                         /* holodeck main loop */
319 > static int
320 > rholo(void)                             /* holodeck main loop */
321   {
322 <        static int      idle = 1;
322 >        static off_t    nextfragwarn = 100L<<20;
323 >        static int      idle = 0;
324          PACKET  *pl = NULL, *plend;
325 +        off_t   fsiz;
326 +        int     pksiz;
327          register PACKET *p;
328          time_t  t;
329 <        long    l;
330 <
331 <        if (outdev != NULL)             /* check display */
329 >                                        /* check display */
330 >        if (nprocs <= 0)
331 >                idle = 1;
332 >        if (outdev != NULL) {
333                  if (!disp_check(idle))
334 <                        return(0);
335 <                                        /* display only? */
336 <        if (ncprocs <= 0)
337 <                return(1);
338 <                                        /* check file size */
339 <        if (maxdisk > 0 && hdfilen(hdlist[0]->fd) >= maxdisk) {
340 <                error(WARNING, "file limit exceeded");
341 <                return(0);
334 >                        return(0);      /* quit request */
335 >                if (nprocs <= 0)
336 >                        return(1);
337 >        } else if (idle)
338 >                return(0);              /* all done */
339 >        fsiz = hdfilen(hdlist[0]->fd);  /* check file size */
340 >        if (maxdisk > 0 && fsiz >= maxdisk) {
341 >                error(USER, "file limit exceeded");
342 >                done_rtrace();
343 >                return(1);      /* comes back */
344          }
345 <                                        /* check time */
346 <        if (endtime > 0 || reporttime > 0)
347 <                t = time(NULL);
345 > #if FRAGWARN
346 >        if (fsiz >= nextfragwarn) {
347 >                double  pctfrag = 100.*(fsiz-hdfiluse(hdlist[0]->fd))/fsiz;
348 >                if (pctfrag >= (double)FRAGWARN) {
349 >                        sprintf(errmsg, "holodeck file fragmentation is %.0f%%",
350 >                                        pctfrag);
351 >                        error(WARNING, errmsg);
352 >                        nextfragwarn = fsiz + (fsiz>>2);
353 >                } else
354 >                        nextfragwarn = fsiz + (10L<<20);
355 >        }
356 > #endif
357 >        t = time(NULL);                 /* check time */
358          if (endtime > 0 && t >= endtime) {
359 <                error(WARNING, "time limit exceeded");
360 <                return(0);
359 >                error(USER, "time limit exceeded");
360 >                done_rtrace();
361 >                return(1);      /* comes back */
362          }
363          if (reporttime > 0 && t >= reporttime)
364                  report(t);
365 <                                        /* get packets to process */
365 >                                        /* figure out good packet size */
366 >        pksiz = RPACKSIZ;
367 > #if MAXQTIME
368 >        if (!chunkycmp) {
369 >                pksiz = nraysdone*MAXQTIME/(totqlen*(t - starttime + 1L));
370 >                if (pksiz < 1) pksiz = 1;
371 >                else if (pksiz > RPACKSIZ) pksiz = RPACKSIZ;
372 >        }
373 > #endif
374 >        idle = 0;                       /* get packets to process */
375          while (freepacks != NULL) {
376                  p = freepacks; freepacks = p->next; p->next = NULL;
377 <                if (!next_packet(p)) {
377 >                if (!next_packet(p, pksiz)) {
378                          p->next = freepacks; freepacks = p;
379 +                        idle = 1;
380                          break;
381                  }
382                  if (pl == NULL) pl = p;
383                  else plend->next = p;
384                  plend = p;
385          }
386 <        idle = pl == NULL && freepacks != NULL;
300 <                                        /* are we out of stuff to do? */
301 <        if (idle && outdev == NULL)
302 <                return(0);
303 <                                        /* else process packets */
386 >                                        /* process packets */
387          done_packets(do_packets(pl));
388          return(1);                      /* and continue */
389   }
390  
391  
392 < report(t)                       /* report progress so far */
393 < time_t  t;
392 > static void
393 > setdefaults(                    /* set default values */
394 >        register HDGRID *gp
395 > )
396   {
312 –        if (t == 0)
313 –                t = time(NULL);
314 –        fprintf(stderr, "%s: %ld packets (%ld rays) done after %.2f hours\n",
315 –                        progname, npacksdone, nraysdone, (t-starttime)/3600.);
316 –        fflush(stderr);
317 –        if (vdef(REPORT))
318 –                reporttime = t + (time_t)(vflt(REPORT)*60.+.5);
319 – }
320 –
321 –
322 – setdefaults(gp)                 /* set default values */
323 – register HDGRID *gp;
324 – {
397          extern char     *atos();
398          register int    i;
399 <        double  len[3], maxlen, d;
400 <        char    buf[64];
399 >        int     n;
400 >        double  len[3], d;
401  
402          if (!vdef(SECTION)) {
403                  sprintf(errmsg, "%s must be defined", vnam(SECTION));
404                  error(USER, errmsg);
405          }
334 –        if (vdef(SECTION) > 1) {
335 –                sprintf(errmsg, "ignoring all but first %s", vnam(SECTION));
336 –                error(WARNING, errmsg);
337 –        }
338 –        if (sscanf(vval(SECTION),
339 –                        "%lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf",
340 –                        &gp->orig[0], &gp->orig[1], &gp->orig[2],
341 –                        &gp->xv[0][0], &gp->xv[0][1], &gp->xv[0][2],
342 –                        &gp->xv[1][0], &gp->xv[1][1], &gp->xv[1][2],
343 –                        &gp->xv[2][0], &gp->xv[2][1], &gp->xv[2][2]) != 12)
344 –                badvalue(SECTION);
345 –        maxlen = 0.;
346 –        for (i = 0; i < 3; i++)
347 –                if ((len[i] = VLEN(gp->xv[i])) > maxlen)
348 –                        maxlen = len[i];
349 –        if (!vdef(GRID)) {
350 –                sprintf(buf, "%.4f", maxlen/8.);
351 –                vval(GRID) = savqstr(buf);
352 –                vdef(GRID)++;
353 –        }
354 –        if ((d = vflt(GRID)) <= FTINY)
355 –                badvalue(GRID);
356 –        for (i = 0; i < 3; i++)
357 –                gp->grid[i] = len[i]/d + (1.-FTINY);
358 –        if (!vdef(EXPOSURE)) {
359 –                sprintf(errmsg, "%s must be defined", vnam(EXPOSURE));
360 –                error(USER, errmsg);
361 –        }
362 –        expval = vval(EXPOSURE)[0] == '-' || vval(EXPOSURE)[0] == '+' ?
363 –                        pow(2., vflt(EXPOSURE)) : vflt(EXPOSURE);
406          if (!vdef(OCTREE)) {
407                  if ((vval(OCTREE) = bmalloc(strlen(froot)+5)) == NULL)
408                          error(SYSTEM, "out of memory");
409                  sprintf(vval(OCTREE), "%s.oct", froot);
410                  vdef(OCTREE)++;
411          }
370 –        if (!vdef(OBSTRUCTIONS)) {
371 –                vval(OBSTRUCTIONS) = "T";
372 –                vdef(OBSTRUCTIONS)++;
373 –        }
412          if (!vdef(VDIST)) {
413                  vval(VDIST) = "F";
414                  vdef(VDIST)++;
415          }
378 –        if (!vdef(OCCUPANCY)) {
379 –                vval(OCCUPANCY) = "U";
380 –                vdef(OCCUPANCY)++;
381 –        }
416                                  /* append rendering options */
417          if (vdef(RENDER))
418                  rtargc += wordstring(rtargv+rtargc, vval(RENDER));
419 +        
420 +        if (gp == NULL)         /* already initialized? */
421 +                return;
422 +                                /* set grid parameters */
423 +        for (n = 0; n < vdef(SECTION); n++, gp++) {
424 +                gp->grid[0] = gp->grid[1] = gp->grid[2] = 0;
425 +                if (sscanf(nvalue(SECTION, n),
426 +                "%lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %hd %hd %hd",
427 +                                &gp->orig[0], &gp->orig[1], &gp->orig[2],
428 +                                &gp->xv[0][0], &gp->xv[0][1], &gp->xv[0][2],
429 +                                &gp->xv[1][0], &gp->xv[1][1], &gp->xv[1][2],
430 +                                &gp->xv[2][0], &gp->xv[2][1], &gp->xv[2][2],
431 +                                &gp->grid[0], &gp->grid[1], &gp->grid[2]) < 12)
432 +                        badvalue(SECTION);
433 +                for (i = 0; i < 3; i++)
434 +                        len[i] = VLEN(gp->xv[i]);
435 +                if (!vdef(GRID)) {
436 +                        d = 2/2e5*( len[0]*len[0]*(len[1]*len[1] +
437 +                                        len[2]*len[2] + 4*len[1]*len[2])
438 +                                + len[1]*len[1]*len[2]*(len[2] + 4*len[0])
439 +                                + 4*len[0]*len[1]*len[2]*len[2] );
440 +                        d = sqrt(sqrt(d));
441 +                } else if ((d = vflt(GRID)) <= FTINY)
442 +                        badvalue(GRID);
443 +                for (i = 0; i < 3; i++)
444 +                        if (gp->grid[i] <= 0)
445 +                                gp->grid[i] = len[i]/d + (1.-FTINY);
446 +        }
447   }
448  
449  
450 < creatholo(gp)                   /* create a holodeck output file */
451 < HDGRID  *gp;
450 > static void
451 > creatholo(                      /* create a holodeck output file */
452 >        HDGRID  *gp
453 > )
454   {
455          extern char     VersionID[];
456 <        long    endloc = 0;
456 >        off_t   lastloc, nextloc;
457 >        int     n;
458          int     fd;
459          FILE    *fp;
460                                          /* open & truncate file */
# Line 403 | Line 468 | HDGRID *gp;
468          printvars(fp);
469          fputformat(HOLOFMT, fp);
470          fputc('\n', fp);
471 <        putw(HOLOMAGIC, fp);            /* put magic number & terminus */
407 <        fwrite(&endloc, sizeof(long), 1, fp);
471 >        putw(HOLOMAGIC, fp);            /* put magic number */
472          fd = dup(fileno(fp));
473          fclose(fp);                     /* flush and close stdio stream */
474 <        hdinit(fd, gp);                 /* allocate and initialize index */
474 >        lastloc = lseek(fd, (off_t)0, SEEK_END);
475 >        for (n = vdef(SECTION); n--; gp++) {    /* initialize each section */
476 >                nextloc = 0L;
477 >                write(fd, (char *)&nextloc, sizeof(nextloc));
478 >                hdinit(fd, gp);                 /* writes beam index */
479 >                if (!n)
480 >                        break;
481 >                nextloc = hdfilen(fd);          /* write section pointer */
482 >                if (lseek(fd, lastloc, SEEK_SET) < 0)
483 >                        error(SYSTEM,
484 >                                "cannot seek on holodeck file in creatholo");
485 >                write(fd, (char *)&nextloc, sizeof(nextloc));
486 >                lseek(fd, (lastloc=nextloc), SEEK_SET);
487 >        }
488   }
489  
490  
491 < headline(s)                     /* process information header line */
492 < char    *s;
491 > static int
492 > headline(                       /* process information header line */
493 >        char    *s,
494 >        void    *p
495 > )
496   {
497          extern char     FMTSTR[];
498          register char   *cp;
# Line 424 | Line 504 | char   *s;
504                                          HOLOFMT, hdkfile, FMTSTR, fmt);
505                          error(USER, errmsg);
506                  }
507 <                return;
507 >                return(0);
508          }
509          for (cp = s; *cp; cp++)         /* take off any comments */
510                  if (*cp == '#') {
# Line 432 | Line 512 | char   *s;
512                          break;
513                  }
514          setvariable(s, matchvar);       /* don't flag errors */
515 +        return(0);
516   }
517  
518  
519 < loadholo()                      /* start loading a holodeck from fname */
519 > static void
520 > loadholo(void)                  /* start loading a holodeck from fname */
521   {
440 –        extern long     ftell();
522          FILE    *fp;
523          int     fd;
524 <        long    fpos;
525 <                                        /* open holodeck file */
526 <        if ((fp = fopen(hdkfile, ncprocs>0 ? "r+" : "r")) == NULL) {
527 <                sprintf(errmsg, "cannot open \"%s\" for %s", hdkfile,
528 <                                ncprocs>0 ? "appending" : "reading");
529 <                error(SYSTEM, errmsg);
524 >        int     n;
525 >        off_t   nextloc;
526 >        
527 >        if ((ncprocs > 0) & (force >= 0))
528 >                fp = fopen(hdkfile, "r+");
529 >        else
530 >                fp = NULL;
531 >        if (fp == NULL) {
532 >                if ((fp = fopen(hdkfile, "r")) == NULL) {
533 >                        sprintf(errmsg, "cannot open \"%s\"", hdkfile);
534 >                        error(SYSTEM, errmsg);
535 >                }
536 >                if (ncprocs > 0) {
537 >                        sprintf(errmsg, "\"%s\" is read-only", hdkfile);
538 >                        if (outdev == NULL)
539 >                                error(USER, errmsg);
540 >                        strcat(errmsg, "; new rays will be discarded");
541 >                        error(WARNING, errmsg);
542 >                        force = -1;
543 >                }
544          }
545                                          /* load variables from header */
546          getheader(fp, headline, NULL);
# Line 455 | Line 550 | loadholo()                     /* start loading a holodeck from fname */
550                                  hdkfile);
551                  error(USER, errmsg);
552          }
553 <        fread(&fpos, sizeof(long), 1, fp);
459 <        if (fpos != 0)
460 <                error(WARNING, "ignoring multiple sections in holodeck file");
461 <        fpos = ftell(fp);                       /* get stdio position */
553 >        nextloc = ftell(fp);                    /* get stdio position */
554          fd = dup(fileno(fp));
555          fclose(fp);                             /* done with stdio */
556 <        lseek(fd, fpos, 0);                     /* align system file pointer */
557 <        hdinit(fd, NULL);                       /* allocate and load index */
556 >        for (n = 0; nextloc > 0L; n++) {        /* initialize each section */
557 >                lseek(fd, nextloc, SEEK_SET);
558 >                read(fd, (char *)&nextloc, sizeof(nextloc));
559 >                hdinit(fd, NULL);
560 >        }
561 >        if (n != vdef(SECTION)) {
562 >                sprintf(errmsg, "number of sections does not match %s setting",
563 >                                vnam(SECTION));
564 >                error(WARNING, errmsg);
565 >        }
566   }
567  
568  
569 < done_packets(pl)                /* handle finished packets */
570 < PACKET  *pl;
569 > extern void
570 > done_packets(           /* handle finished packets */
571 >        PACKET  *pl
572 > )
573   {
574          static int      n2flush = 0;
575          register PACKET *p;
# Line 475 | Line 577 | PACKET *pl;
577          while (pl != NULL) {
578                  p = pl; pl = p->next; p->next = NULL;
579                  if (p->nr > 0) {                /* add to holodeck */
580 <                        bcopy((char *)p->ra,
581 <                                (char *)hdnewrays(hdlist[p->hd],p->bi,p->nr),
580 >                        memcpy( (void *)hdnewrays(hdlist[p->hd],p->bi,p->nr),
581 >                                (void *)p->ra,
582                                  p->nr*sizeof(RAYVAL));
583                          if (outdev != NULL)     /* display it */
584                                  disp_packet((PACKHEAD *)p);
585 <                        if (hdcachesize <= 0)   /* manual flushing */
586 <                                n2flush += p->nr;
585 >                        if (hdcachesize <= 0)
586 >                                n2flush++;
587                          nraysdone += p->nr;
588                          npacksdone++;
589 +                        p->nr = 0;
590                  }
591 <                p->nr = 0;                      /* push onto free list */
489 <                p->next = freepacks;
591 >                p->next = freepacks;            /* push onto free list */
592                  freepacks = p;
593          }
594 <        if (n2flush > 512*RPACKSIZ*ncprocs) {
595 <                hdflush(NULL);                  /* flush holodeck buffers */
594 >        if (n2flush >= RTFLUSH) {
595 >                if (outdev != NULL)
596 >                        hdsync(NULL, 1);
597 >                else
598 >                        hdflush(NULL);
599                  n2flush = 0;
600          }
601   }
602  
603  
604 < getradfile(rfargs)              /* run rad and get needed variables */
605 < char    *rfargs;
604 > static void
605 > rootname(               /* remove tail from end of fn */
606 >        register char   *rn,
607 >        register char   *fn
608 > )
609   {
502 –        static short    mvar[] = {OCTREE,EXPOSURE,-1};
503 –        static char     tf1[] = TEMPLATE;
504 –        char    tf2[64];
505 –        char    combuf[256];
506 –        char    *pippt;
507 –        register int    i;
508 –        register char   *cp;
509 –                                        /* create rad command */
510 –        mktemp(tf1);
511 –        sprintf(tf2, "%s.rif", tf1);
512 –        sprintf(combuf,
513 –                "rad -v 0 -s -e -w %s OPTFILE=%s | egrep '^[ \t]*(NOMATCH",
514 –                        rfargs, tf1);
515 –        cp = combuf;
516 –        while (*cp){
517 –                if (*cp == '|') pippt = cp;
518 –                cp++;
519 –        }                               /* match unset variables */
520 –        for (i = 0; mvar[i] >= 0; i++)
521 –                if (!vdef(mvar[i])) {
522 –                        *cp++ = '|';
523 –                        strcpy(cp, vnam(mvar[i]));
524 –                        while (*cp) cp++;
525 –                        pippt = NULL;
526 –                }
527 –        if (pippt != NULL)
528 –                strcpy(pippt, "> /dev/null");   /* nothing to match */
529 –        else
530 –                sprintf(cp, ")[ \t]*=' > %s", tf2);
531 –        if (system(combuf)) {
532 –                error(SYSTEM, "cannot execute rad command");
533 –                unlink(tf2);                    /* clean up */
534 –                unlink(tf1);
535 –                quit(1);
536 –        }
537 –        if (pippt == NULL) {
538 –                loadvars(tf2);                  /* load variables */
539 –                unlink(tf2);
540 –        }
541 –        rtargc += wordfile(rtargv+rtargc, tf1); /* get rtrace options */
542 –        unlink(tf1);                    /* clean up */
543 – }
544 –
545 –
546 – rootname(rn, fn)                /* remove tail from end of fn */
547 – register char   *rn, *fn;
548 – {
610          char    *tp, *dp;
611  
612 <        for (tp = NULL, dp = rn; *rn = *fn++; rn++)
613 <                if (ISDIRSEP(*rn))
612 >        for (tp = NULL, dp = rn; (*rn = *fn++); rn++) {
613 >                if (*rn == '/')
614                          dp = rn;
615                  else if (*rn == '.')
616                          tp = rn;
617 +        }
618          if (tp != NULL && tp > dp)
619                  *tp = '\0';
620   }
621  
622  
623 < badvalue(vc)                    /* report bad variable value and exit */
624 < int     vc;
623 > static void
624 > badvalue(                       /* report bad variable value and exit */
625 >        int     vc
626 > )
627   {
628          sprintf(errmsg, "bad value for variable '%s'", vnam(vc));
629          error(USER, errmsg);
630   }
631  
632  
633 + void
634   eputs(s)                        /* put error message to stderr */
635   register char  *s;
636   {
# Line 585 | Line 650 | register char  *s;
650   }
651  
652  
653 < wputs(s)                        /* put warning string to stderr */
589 < char    *s;
590 < {
591 <        if (!nowarn)
592 <                eputs(s);
593 < }
594 <
595 <
653 > void
654   quit(ec)                        /* exit program gracefully */
655   int     ec;
656   {
657          int     status = 0;
658  
659 <        if (hdlist[0] != NULL) {        /* flush holodeck */
660 <                if (ncprocs > 0) {
661 <                        done_packets(flush_queue());
662 <                        status = end_rtrace();  /* close rtrace */
663 <                        hdflush(NULL);
664 <                        if (vdef(REPORT)) {
665 <                                long    fsiz, fuse;
666 <                                report(0);
609 <                                fsiz = hdfilen(hdlist[0]->fd);
610 <                                fuse = hdfiluse(hdlist[0]->fd, 1);
611 <                                fprintf(stderr,
659 >        if (hdlist[0] != NULL) {        /* close holodeck */
660 >                if (nprocs > 0)
661 >                        status = done_rtrace();         /* calls hdsync() */
662 >                if ((ncprocs > 0) & (force >= 0) && vdef(REPORT)) {
663 >                        off_t   fsiz, fuse;
664 >                        fsiz = hdfilen(hdlist[0]->fd);
665 >                        fuse = hdfiluse(hdlist[0]->fd);
666 >                        fprintf(stderr,
667                          "%s: %.1f Mbyte holodeck file, %.1f%% fragmentation\n",
668 <                                                hdkfile, fsiz/(1024.*1024.),
669 <                                                100.*(fsiz-fuse)/fsiz);
670 <                        }
616 <                } else
617 <                        hdflush(NULL);
668 >                                        hdkfile, fsiz/(1024.*1024.),
669 >                                        100.*(fsiz-fuse)/fsiz);
670 >                }
671          }
672          if (orig_mode >= 0)             /* reset holodeck access mode */
673                  fchmod(hdlist[0]->fd, orig_mode);

Diff Legend

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