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.2 by gregl, Fri Oct 31 11:44:09 1997 UTC vs.
Revision 3.74 by greg, Fri Dec 5 00:22:33 2008 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 <stdio.h>
9 + #include <time.h>
10 + #include <signal.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 "paths.h"
13 < #include <sys/types.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 19 | 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  
24 int     nowarn = 0;             /* turn warnings off? */
25
26 double  expval = 1.;            /* global exposure value */
27
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 + off_t   maxdisk;                /* maximum file space (bytes) */
56 +
57   int     rtargc = 1;             /* rtrace command */
58   char    *rtargv[128] = {"rtrace", NULL};
59  
60 + int     orig_mode = -1;         /* original file mode (-1 if unchanged) */
61 +
62   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 < extern time_t   time();
68 > char  *sigerr[NSIG];            /* signal error messages */
69  
70 + extern int      nowarn;         /* turn warnings off? */
71  
72 < main(argc, argv)
73 < int     argc;
74 < char    *argv[];
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 >
85 > int
86 > main(
87 >        int     argc,
88 >        char    *argv[]
89 > )
90   {
51        HDGRID  hdg;
91          int     i;
53        int     force = 0;
92  
93          progname = argv[0];                     /* get arguments */
94          for (i = 1; i < argc && argv[i][0] == '-'; i++)
# Line 59 | 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 74 | Line 118 | char   *argv[];
118                  default:
119                          goto userr;
120                  }
77                                                /* do we have a job? */
78        if (outdev == NULL && ncprocs <= 0)
79                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();
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 < initrholo()                     /* get our holodeck running */
174 > static void
175 > onsig(                          /* fatal signal */
176 >        int  signo
177 > )
178   {
179 +        static int  gotsig = 0;
180 +
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(300);                     /* allow 5 minutes to clean up */
188 +        eputs("signal - ");
189 +        eputs(sigerr[signo]);
190 +        eputs("\n");
191 +        quit(3);
192 + }
193 +
194 +
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);
203 +        sigerr[signo] = msg;
204 + }
205 +
206 +
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 & 0777))
219 +                return(-1);             /* already set */
220 +                                        /* else change it */
221 +        if (fchmod(fd, mod) < 0) {
222 +                error(WARNING, "cannot change holodeck file access mode");
223 +                return(-1);
224 +        }
225 +        return(stbuf.st_mode);          /* return original mode */
226 + }
227 +
228 +
229 + static void
230 + initrholo(void)                 /* get our holodeck running */
231 + {
232          extern int      global_packet();
233          register int    i;
234 <                                                /* check output device */
235 <        if (outdev != NULL)
236 <                open_display(outdev);
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);
239          else if (ncprocs > 0)                   /* else use global ray feed */
240                  init_global();
241 <                                                /* record the time */
242 <        starttime = time(NULL);
243 <        if (!vdef(TIME) || vflt(TIME) <= FTINY)
134 <                endtime = 0;
241 >                                                /* record disk space limit */
242 >        if (!vdef(DISKSPACE))
243 >                maxdisk = ((off_t)1<<(sizeof(off_t)*8-2)) - 1024;
244          else
245 <                endtime = starttime + vflt(TIME)*3600.;
245 >                maxdisk = 1024.*1024.*vflt(DISKSPACE);
246                                                  /* set up memory cache */
247 <        hdcachesize = 1024.*1024.*vflt(CACHE);
247 >        if (outdev == NULL)
248 >                hdcachesize = 0;                /* manual flushing */
249 >        else if (vdef(CACHE))
250 >                hdcachesize = 1024.*1024.*vflt(CACHE);
251                                                  /* open report file */
252          if (vdef(REPORT)) {
253                  register char   *s = sskip2(vval(REPORT), 1);
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 157 | 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 171 | Line 290 | initrholo()                    /* get our holodeck running */
290                          freepacks[i].next = &freepacks[i+1];
291                  }
292          }
293 +                                        /* set up signal handling */
294 +        sigdie(SIGINT, "Interrupt");
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) & (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;
186        long    l;
329                                          /* check display */
330 <        if (outdev != NULL && !disp_check(idle))
331 <                return(0);
332 <                                        /* display only? */
333 <        if (ncprocs <= 0)
334 <                return(1);
335 <                                        /* check file size */
336 <        if ((l = 1024.*1024.*vflt(DISKSPACE)) > 0 &&
337 <                        hdfiluse(hdlist[0]->fd, 0) + hdmemuse(0) >= l)
338 <                return(0);
339 <                                        /* check time */
340 <        if (endtime > 0 || vdef(REPORT))
341 <                t = time(NULL);
342 <        if (endtime > 0 && t >= endtime)
343 <                return(0);
344 <        if (vdef(REPORT) && t >= reporttime)
330 >        if (nprocs <= 0)
331 >                idle = 1;
332 >        if (outdev != NULL) {
333 >                if (!disp_check(idle))
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 > #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(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;
216 <                                        /* are we out of stuff to do? */
217 <        if (idle && outdev == NULL)
218 <                return(0);
219 <                                        /* 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   {
228        if (t == 0)
229                t = time(NULL);
230        fprintf(stderr, "%s: %ld packets (%ld rays) done after %.2f hours\n",
231                        progname, npacksdone, nraysdone, (t-starttime)/3600.);
232        fflush(stderr);
233        reporttime = t + (time_t)(vflt(REPORT)*60.);
234 }
235
236
237 setdefaults(gp)                 /* set default values */
238 register HDGRID *gp;
239 {
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          }
249        if (vdef(SECTION) > 1) {
250                sprintf(errmsg, "ignoring all but first %s", vnam(SECTION));
251                error(WARNING, errmsg);
252        }
253        if (sscanf(vval(SECTION),
254                        "%lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf",
255                        &gp->orig[0], &gp->orig[1], &gp->orig[2],
256                        &gp->xv[0][0], &gp->xv[0][1], &gp->xv[0][2],
257                        &gp->xv[1][0], &gp->xv[1][1], &gp->xv[1][2],
258                        &gp->xv[2][0], &gp->xv[2][1], &gp->xv[2][2]) != 12)
259                badvalue(SECTION);
260        maxlen = 0.;
261        for (i = 0; i < 3; i++)
262                if ((len[i] = VLEN(gp->xv[i])) > maxlen)
263                        maxlen = len[i];
264        if (!vdef(GRID)) {
265                sprintf(buf, "%.4f", maxlen/8.);
266                vval(GRID) = savqstr(buf);
267                vdef(GRID)++;
268        }
269        if ((d = vflt(GRID)) <= FTINY)
270                badvalue(GRID);
271        for (i = 0; i < 3; i++)
272                gp->grid[i] = len[i]/d + (1.-FTINY);
273        if (!vdef(EXPOSURE)) {
274                sprintf(errmsg, "%s must be defined", vnam(EXPOSURE));
275                error(USER, errmsg);
276        }
277        expval = vval(EXPOSURE)[0] == '-' || vval(EXPOSURE)[0] == '+' ?
278                        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          }
412 <        if (!vdef(DISKSPACE)) {
413 <                sprintf(errmsg,
414 <                        "no %s setting, assuming 100 Mbytes available",
288 <                                vnam(DISKSPACE));
289 <                error(WARNING, errmsg);
290 <                vval(DISKSPACE) = "100";
291 <                vdef(DISKSPACE)++;
412 >        if (!vdef(VDIST)) {
413 >                vval(VDIST) = "F";
414 >                vdef(VDIST)++;
415          }
293        if (!vdef(CACHE)) {
294                sprintf(errmsg,
295                        "no %s setting, assuming 10 Mbytes available",
296                                vnam(CACHE));
297                error(WARNING, errmsg);
298                vval(CACHE) = "10";
299                vdef(CACHE)++;
300        }
301        if (!vdef(OBSTRUCTIONS)) {
302                vval(OBSTRUCTIONS) = "T";
303                vdef(OBSTRUCTIONS)++;
304        }
305        if (!vdef(OCCUPANCY)) {
306                vval(OCCUPANCY) = "U";
307                vdef(OCCUPANCY)++;
308        }
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 <        long    endloc = 0;
455 >        extern char     VersionID[];
456 >        int32   lastloc, nextloc;
457 >        int     n;
458 >        int     fd;
459          FILE    *fp;
460                                          /* open & truncate file */
461          if ((fp = fopen(hdkfile, "w+")) == NULL) {
# Line 324 | Line 464 | HDGRID *gp;
464          }
465                                          /* write information header */
466          newheader("RADIANCE", fp);
467 +        fprintf(fp, "SOFTWARE= %s\n", VersionID);
468          printvars(fp);
469          fputformat(HOLOFMT, fp);
470          fputc('\n', fp);
471 <        putw(HOLOMAGIC, fp);            /* put magic number & terminus */
472 <        fwrite(&endloc, sizeof(long), 1, fp);
473 <        fflush(fp);                     /* flush buffer */
474 <        hdinit(fileno(fp), gp);         /* allocate and initialize index */
475 <                        /* we're dropping fp here.... */
471 >        putw(HOLOMAGIC, fp);            /* put magic number */
472 >        fd = dup(fileno(fp));
473 >        fclose(fp);                     /* flush and close stdio stream */
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, (off_t)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, (off_t)(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 348 | 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 356 | 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   {
522          FILE    *fp;
523 <        long    endloc;
524 <                                        /* open input file */
525 <        if ((fp = fopen(hdkfile, "r+")) == NULL) {
526 <                sprintf(errmsg, "cannot open \"%s\" for appending", hdkfile);
527 <                error(SYSTEM, errmsg);
523 >        int     fd;
524 >        int     n;
525 >        int32   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 376 | Line 550 | loadholo()                     /* start loading a holodeck from fname */
550                                  hdkfile);
551                  error(USER, errmsg);
552          }
553 <        fread(&endloc, sizeof(long), 1, fp);
554 <        if (endloc != 0)
555 <                error(WARNING, "ignoring multiple sections in holodeck file");
556 <        fseek(fp, 0L, 1);                       /* align system file pointer */
557 <        hdinit(fileno(fp), NULL);               /* allocate and load index */
558 <                        /* we're dropping fp here.... */
553 >        nextloc = ftell(fp);                    /* get stdio position */
554 >        fd = dup(fileno(fp));
555 >        fclose(fp);                             /* done with stdio */
556 >        for (n = 0; nextloc > 0L; n++) {        /* initialize each section */
557 >                lseek(fd, (off_t)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;
576  
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(p);
584 >                                disp_packet((PACKHEAD *)p);
585 >                        if (hdcachesize <= 0)
586 >                                n2flush++;
587 >                        nraysdone += p->nr;
588 >                        npacksdone++;
589 >                        p->nr = 0;
590                  }
591 <                nraysdone += p->nr;
403 <                npacksdone++;
404 <                p->nr = 0;                      /* push onto free list */
405 <                p->next = freepacks;
591 >                p->next = freepacks;            /* push onto free list */
592                  freepacks = p;
593          }
594 < }
595 <
596 <
597 < getradfile(rfargs)              /* run rad and get needed variables */
598 < char    *rfargs;
599 < {
414 <        static short    mvar[] = {VIEW,OCTREE,EXPOSURE,-1};
415 <        static char     tf1[] = TEMPLATE;
416 <        char    tf2[64];
417 <        char    combuf[256];
418 <        register int    i;
419 <        register char   *cp;
420 <                                        /* create rad command */
421 <        mktemp(tf1);
422 <        sprintf(tf2, "%s.rif", tf1);
423 <        sprintf(combuf,
424 <                "rad -v 0 -s -e -w %s OPTFILE=%s | egrep '^[ \t]*(NOMATCH",
425 <                        rfargs, tf1);
426 <        cp = combuf;
427 <        while (*cp) cp++;               /* match unset variables */
428 <        for (i = 0; mvar[i] >= 0; i++)
429 <                if (!vdef(mvar[i])) {
430 <                        *cp++ = '|';
431 <                        strcpy(cp, vnam(mvar[i]));
432 <                        while (*cp) cp++;
433 <                }
434 <        sprintf(cp, ")[ \t]*=' > %s", tf2);
435 <        if (system(combuf)) {
436 <                error(SYSTEM, "cannot execute rad command");
437 <                unlink(tf2);                    /* clean up */
438 <                unlink(tf1);
439 <                quit(1);
594 >        if (n2flush >= RTFLUSH) {
595 >                if (outdev != NULL)
596 >                        hdsync(NULL, 1);
597 >                else
598 >                        hdflush(NULL);
599 >                n2flush = 0;
600          }
441        loadvars(tf2);                  /* load variables */
442        rtargc += wordfile(rtargv+rtargc, tf1); /* get rtrace options */
443        unlink(tf2);                    /* clean up */
444        unlink(tf1);
601   }
602  
603  
604 < rootname(rn, fn)                /* remove tail from end of fn */
605 < register char   *rn, *fn;
604 > static void
605 > rootname(               /* remove tail from end of fn */
606 >        register char   *rn,
607 >        register char   *fn
608 > )
609   {
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 487 | Line 650 | register char  *s;
650   }
651  
652  
653 < wputs(s)                        /* put warning string to stderr */
491 < char    *s;
492 < {
493 <        if (!nowarn)
494 <                eputs(s);
495 < }
496 <
497 <
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);
667 <                                fsiz = lseek(hdlist[0]->fd, 0L, 2);
668 <                                fuse = hdfiluse(hdlist[0]->fd, 1);
669 <                                fprintf(stderr,
670 <                        "%s: %.1f Mbyte holodeck file, %.2f%% fragmentation\n",
515 <                                                hdkfile, fsiz/(1024.*1024.),
516 <                                                100.*(fsiz-fuse)/fsiz);
517 <                        }
518 <                } else
519 <                        hdflush(NULL);
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 >                }
671          }
672 +        if (orig_mode >= 0)             /* reset holodeck access mode */
673 +                fchmod(hdlist[0]->fd, orig_mode);
674 +        if (outdev != NULL)             /* close display */
675 +                disp_close();
676          exit(ec ? ec : status);         /* exit */
677   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines