10 |
|
|
11 |
|
#include "rholo.h" |
12 |
|
#include "random.h" |
13 |
– |
#include "paths.h" |
13 |
|
#include <signal.h> |
14 |
|
#include <sys/types.h> |
15 |
|
#include <sys/stat.h> |
21 |
|
|
22 |
|
char *progname; /* our program name */ |
23 |
|
char *hdkfile; /* holodeck file name */ |
24 |
< |
char froot[MAXPATH]; /* root file name */ |
24 |
> |
char froot[256]; /* root file name */ |
25 |
|
|
26 |
|
int nowarn = 0; /* turn warnings off? */ |
27 |
|
|
29 |
– |
double expval = 1.; /* global exposure value */ |
30 |
– |
|
28 |
|
int ncprocs = 0; /* desired number of compute processes */ |
29 |
|
|
30 |
|
char *outdev = NULL; /* output device name */ |
31 |
|
|
32 |
+ |
int readinp = 0; /* read commands from stdin */ |
33 |
+ |
|
34 |
+ |
int force = 0; /* allow overwrite of holodeck */ |
35 |
+ |
|
36 |
|
time_t starttime; /* time we got started */ |
37 |
|
time_t endtime; /* time we should end by */ |
38 |
|
time_t reporttime; /* time for next report */ |
60 |
|
{ |
61 |
|
HDGRID hdg; |
62 |
|
int i; |
62 |
– |
int force = 0; |
63 |
|
/* mark start time */ |
64 |
|
starttime = time(NULL); |
65 |
< |
initurand(10240); /* initialize urand */ |
65 |
> |
initurand(16384); /* initialize urand */ |
66 |
|
progname = argv[0]; /* get arguments */ |
67 |
|
for (i = 1; i < argc && argv[i][0] == '-'; i++) |
68 |
|
switch (argv[i][1]) { |
72 |
|
case 'f': /* force overwrite */ |
73 |
|
force++; |
74 |
|
break; |
75 |
+ |
case 'i': /* read input from stdin */ |
76 |
+ |
readinp++; |
77 |
+ |
break; |
78 |
|
case 'n': /* compute processes */ |
79 |
|
if (i >= argc-2) |
80 |
|
goto userr; |
106 |
|
/* check settings */ |
107 |
|
checkvalues(); |
108 |
|
/* load RIF if any */ |
109 |
< |
if (vdef(RIF)) |
107 |
< |
getradfile(vval(RIF)); |
109 |
> |
getradfile(); |
110 |
|
/* set defaults */ |
111 |
|
setdefaults(&hdg); |
112 |
|
/* holodeck exists? */ |
117 |
|
creatholo(&hdg); |
118 |
|
} else { /* else load holodeck */ |
119 |
|
loadholo(); |
120 |
< |
if (vdef(RIF)) /* load RIF if any */ |
121 |
< |
getradfile(vval(RIF)); |
120 |
> |
/* check settings */ |
121 |
> |
checkvalues(); |
122 |
> |
/* load RIF if any */ |
123 |
> |
getradfile(); |
124 |
> |
/* set defaults */ |
125 |
> |
setdefaults(NULL); |
126 |
|
} |
127 |
|
/* initialize */ |
128 |
|
initrholo(); |
129 |
< |
/* run */ |
129 |
> |
/* main loop */ |
130 |
|
while (rholo()) |
131 |
|
; |
132 |
|
/* done */ |
190 |
|
{ |
191 |
|
extern int global_packet(); |
192 |
|
register int i; |
193 |
+ |
/* close holodeck on exec() */ |
194 |
+ |
fcntl(hdlist[0]->fd, F_SETFD, FD_CLOEXEC); |
195 |
|
|
196 |
|
if (outdev != NULL) /* open output device */ |
197 |
|
disp_open(outdev); |
222 |
|
if (ncprocs > 0) { |
223 |
|
i = start_rtrace(); |
224 |
|
if (i < 1) |
225 |
< |
error(USER, "cannot start rtrace process"); |
225 |
> |
error(USER, "cannot start rtrace process(es)"); |
226 |
|
if (vdef(REPORT)) { /* make first report */ |
227 |
|
printargs(rtargc, rtargv, stderr); |
228 |
|
report(0); |
277 |
|
if (!disp_check(idle)) |
278 |
|
return(0); |
279 |
|
/* display only? */ |
280 |
< |
if (ncprocs <= 0) |
281 |
< |
return(1); |
280 |
> |
if (nprocs <= 0) |
281 |
> |
return(outdev != NULL); |
282 |
|
/* check file size */ |
283 |
|
if (maxdisk > 0 && hdfilen(hdlist[0]->fd) >= maxdisk) { |
284 |
|
error(WARNING, "file limit exceeded"); |
285 |
< |
return(0); |
285 |
> |
done_rtrace(); |
286 |
> |
idle = 1; |
287 |
> |
return(1); /* comes back */ |
288 |
|
} |
289 |
|
/* check time */ |
290 |
|
if (endtime > 0 || reporttime > 0) |
291 |
|
t = time(NULL); |
292 |
|
if (endtime > 0 && t >= endtime) { |
293 |
|
error(WARNING, "time limit exceeded"); |
294 |
< |
return(0); |
294 |
> |
done_rtrace(); |
295 |
> |
idle = 1; |
296 |
> |
return(1); /* comes back */ |
297 |
|
} |
298 |
|
if (reporttime > 0 && t >= reporttime) |
299 |
|
report(t); |
302 |
|
p = freepacks; freepacks = p->next; p->next = NULL; |
303 |
|
if (!next_packet(p)) { |
304 |
|
p->next = freepacks; freepacks = p; |
305 |
+ |
idle = 1; |
306 |
|
break; |
307 |
|
} |
308 |
|
if (pl == NULL) pl = p; |
309 |
|
else plend->next = p; |
310 |
|
plend = p; |
311 |
|
} |
312 |
< |
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 */ |
312 |
> |
/* process packets */ |
313 |
|
done_packets(do_packets(pl)); |
314 |
|
return(1); /* and continue */ |
315 |
|
} |
318 |
|
report(t) /* report progress so far */ |
319 |
|
time_t t; |
320 |
|
{ |
321 |
< |
if (t == 0) |
321 |
> |
static time_t seconds2go = 1000000; |
322 |
> |
|
323 |
> |
if (t == 0L) |
324 |
|
t = time(NULL); |
325 |
< |
fprintf(stderr, "%s: %ld packets (%ld rays) done after %.2f hours\n", |
326 |
< |
progname, npacksdone, nraysdone, (t-starttime)/3600.); |
327 |
< |
fflush(stderr); |
328 |
< |
if (vdef(REPORT)) |
329 |
< |
reporttime = t + (time_t)(vflt(REPORT)*60.+.5); |
325 |
> |
sprintf(errmsg, "%ld packets (%ld rays) done after %.2f hours\n", |
326 |
> |
npacksdone, nraysdone, (t-starttime)/3600.); |
327 |
> |
eputs(errmsg); |
328 |
> |
if (seconds2go == 1000000) |
329 |
> |
seconds2go = vdef(REPORT) ? (long)(vflt(REPORT)*60. + .5) : 0L; |
330 |
> |
if (seconds2go) |
331 |
> |
reporttime = t + seconds2go; |
332 |
|
} |
333 |
|
|
334 |
|
|
348 |
|
sprintf(errmsg, "ignoring all but first %s", vnam(SECTION)); |
349 |
|
error(WARNING, errmsg); |
350 |
|
} |
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); |
351 |
|
if (!vdef(OCTREE)) { |
352 |
|
if ((vval(OCTREE) = bmalloc(strlen(froot)+5)) == NULL) |
353 |
|
error(SYSTEM, "out of memory"); |
369 |
|
/* append rendering options */ |
370 |
|
if (vdef(RENDER)) |
371 |
|
rtargc += wordstring(rtargv+rtargc, vval(RENDER)); |
372 |
+ |
|
373 |
+ |
if (gp == NULL) /* already initialized? */ |
374 |
+ |
return; |
375 |
+ |
/* set grid parameters */ |
376 |
+ |
gp->grid[0] = gp->grid[1] = gp->grid[2] = 0; |
377 |
+ |
if (sscanf(vval(SECTION), |
378 |
+ |
"%lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %hd %hd %hd", |
379 |
+ |
&gp->orig[0], &gp->orig[1], &gp->orig[2], |
380 |
+ |
&gp->xv[0][0], &gp->xv[0][1], &gp->xv[0][2], |
381 |
+ |
&gp->xv[1][0], &gp->xv[1][1], &gp->xv[1][2], |
382 |
+ |
&gp->xv[2][0], &gp->xv[2][1], &gp->xv[2][2], |
383 |
+ |
&gp->grid[0], &gp->grid[1], &gp->grid[2]) < 12) |
384 |
+ |
badvalue(SECTION); |
385 |
+ |
maxlen = 0.; |
386 |
+ |
for (i = 0; i < 3; i++) |
387 |
+ |
if ((len[i] = VLEN(gp->xv[i])) > maxlen) |
388 |
+ |
maxlen = len[i]; |
389 |
+ |
if (!vdef(GRID)) { |
390 |
+ |
sprintf(buf, "%.4f", maxlen/8.); |
391 |
+ |
vval(GRID) = savqstr(buf); |
392 |
+ |
vdef(GRID)++; |
393 |
+ |
} |
394 |
+ |
if ((d = vflt(GRID)) <= FTINY) |
395 |
+ |
badvalue(GRID); |
396 |
+ |
for (i = 0; i < 3; i++) |
397 |
+ |
if (gp->grid[i] <= 0) |
398 |
+ |
gp->grid[i] = len[i]/d + (1.-FTINY); |
399 |
|
} |
400 |
|
|
401 |
|
|
503 |
|
p->next = freepacks; |
504 |
|
freepacks = p; |
505 |
|
} |
506 |
< |
if (n2flush > 512*RPACKSIZ*ncprocs) { |
506 |
> |
if (n2flush > 512*RPACKSIZ*nprocs) { |
507 |
|
hdflush(NULL); /* flush holodeck buffers */ |
508 |
|
n2flush = 0; |
509 |
|
} |
510 |
|
} |
511 |
|
|
512 |
|
|
499 |
– |
getradfile(rfargs) /* run rad and get needed variables */ |
500 |
– |
char *rfargs; |
501 |
– |
{ |
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 |
– |
|
513 |
|
rootname(rn, fn) /* remove tail from end of fn */ |
514 |
|
register char *rn, *fn; |
515 |
|
{ |
516 |
|
char *tp, *dp; |
517 |
|
|
518 |
|
for (tp = NULL, dp = rn; *rn = *fn++; rn++) |
519 |
< |
if (ISDIRSEP(*rn)) |
519 |
> |
if (*rn == '/') |
520 |
|
dp = rn; |
521 |
|
else if (*rn == '.') |
522 |
|
tp = rn; |
566 |
|
int status = 0; |
567 |
|
|
568 |
|
if (hdlist[0] != NULL) { /* flush holodeck */ |
569 |
< |
if (ncprocs > 0) { |
570 |
< |
done_packets(flush_queue()); |
571 |
< |
status = end_rtrace(); /* close rtrace */ |
572 |
< |
hdflush(NULL); |
573 |
< |
if (vdef(REPORT)) { |
574 |
< |
long fsiz, fuse; |
575 |
< |
report(0); |
576 |
< |
fsiz = hdfilen(hdlist[0]->fd); |
610 |
< |
fuse = hdfiluse(hdlist[0]->fd, 1); |
611 |
< |
fprintf(stderr, |
569 |
> |
if (nprocs > 0) |
570 |
> |
status = done_rtrace(); |
571 |
> |
hdflush(NULL); |
572 |
> |
if (ncprocs > 0 && vdef(REPORT)) { |
573 |
> |
long fsiz, fuse; |
574 |
> |
fsiz = hdfilen(hdlist[0]->fd); |
575 |
> |
fuse = hdfiluse(hdlist[0]->fd, 1); |
576 |
> |
fprintf(stderr, |
577 |
|
"%s: %.1f Mbyte holodeck file, %.1f%% fragmentation\n", |
578 |
< |
hdkfile, fsiz/(1024.*1024.), |
579 |
< |
100.*(fsiz-fuse)/fsiz); |
580 |
< |
} |
616 |
< |
} else |
617 |
< |
hdflush(NULL); |
578 |
> |
hdkfile, fsiz/(1024.*1024.), |
579 |
> |
100.*(fsiz-fuse)/fsiz); |
580 |
> |
} |
581 |
|
} |
582 |
|
if (orig_mode >= 0) /* reset holodeck access mode */ |
583 |
|
fchmod(hdlist[0]->fd, orig_mode); |