9 |
|
*/ |
10 |
|
|
11 |
|
#include "rholo.h" |
12 |
< |
#include "paths.h" |
12 |
> |
#include "random.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 |
|
|
28 |
– |
double expval = 1.; /* global exposure value */ |
29 |
– |
|
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; |
61 |
– |
int force = 0; |
63 |
|
/* mark start time */ |
64 |
|
starttime = time(NULL); |
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)) |
105 |
< |
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); |
206 |
|
if (!vdef(TIME) || vflt(TIME) <= FTINY) |
207 |
|
endtime = 0; |
208 |
|
else |
209 |
< |
endtime = starttime + vflt(TIME)*3600.; |
209 |
> |
endtime = starttime + vflt(TIME)*3600. + .5; |
210 |
|
/* set up memory cache */ |
211 |
|
if (outdev == NULL) |
212 |
|
hdcachesize = 0; /* manual flushing */ |
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); |
233 |
|
goto memerr; |
234 |
|
freepacks[--i].nr = 0; |
235 |
|
freepacks[i].next = NULL; |
236 |
< |
if (!vbool(OBSTRUCTIONS)) { |
236 |
> |
if (!vdef(OBSTRUCTIONS) || !vbool(OBSTRUCTIONS)) { |
237 |
|
freepacks[i].offset = (float *)bmalloc( |
238 |
|
RPACKSIZ*sizeof(float)*(i+1) ); |
239 |
|
if (freepacks[i].offset == NULL) |
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; |
298 |
< |
/* are we out of stuff to do? */ |
299 |
< |
if (idle && outdev == NULL) |
300 |
< |
return(0); |
301 |
< |
/* else process packets */ |
312 |
> |
/* process packets */ |
313 |
|
done_packets(do_packets(pl)); |
314 |
|
return(1); /* and continue */ |
315 |
|
} |
316 |
|
|
317 |
|
|
307 |
– |
report(t) /* report progress so far */ |
308 |
– |
time_t t; |
309 |
– |
{ |
310 |
– |
if (t == 0) |
311 |
– |
t = time(NULL); |
312 |
– |
fprintf(stderr, "%s: %ld packets (%ld rays) done after %.2f hours\n", |
313 |
– |
progname, npacksdone, nraysdone, (t-starttime)/3600.); |
314 |
– |
fflush(stderr); |
315 |
– |
if (vdef(REPORT)) |
316 |
– |
reporttime = t + (time_t)(vflt(REPORT)*60.+.5); |
317 |
– |
} |
318 |
– |
|
319 |
– |
|
318 |
|
setdefaults(gp) /* set default values */ |
319 |
|
register HDGRID *gp; |
320 |
|
{ |
331 |
|
sprintf(errmsg, "ignoring all but first %s", vnam(SECTION)); |
332 |
|
error(WARNING, errmsg); |
333 |
|
} |
334 |
+ |
if (!vdef(OCTREE)) { |
335 |
+ |
if ((vval(OCTREE) = bmalloc(strlen(froot)+5)) == NULL) |
336 |
+ |
error(SYSTEM, "out of memory"); |
337 |
+ |
sprintf(vval(OCTREE), "%s.oct", froot); |
338 |
+ |
vdef(OCTREE)++; |
339 |
+ |
} |
340 |
+ |
if (!vdef(VDIST)) { |
341 |
+ |
vval(VDIST) = "F"; |
342 |
+ |
vdef(VDIST)++; |
343 |
+ |
} |
344 |
+ |
if (!vdef(OCCUPANCY)) { |
345 |
+ |
vval(OCCUPANCY) = "U"; |
346 |
+ |
vdef(OCCUPANCY)++; |
347 |
+ |
} |
348 |
+ |
/* append rendering options */ |
349 |
+ |
if (vdef(RENDER)) |
350 |
+ |
rtargc += wordstring(rtargv+rtargc, vval(RENDER)); |
351 |
+ |
|
352 |
+ |
if (gp == NULL) /* already initialized? */ |
353 |
+ |
return; |
354 |
+ |
/* set grid parameters */ |
355 |
+ |
gp->grid[0] = gp->grid[1] = gp->grid[2] = 0; |
356 |
|
if (sscanf(vval(SECTION), |
357 |
< |
"%lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf", |
357 |
> |
"%lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %hd %hd %hd", |
358 |
|
&gp->orig[0], &gp->orig[1], &gp->orig[2], |
359 |
|
&gp->xv[0][0], &gp->xv[0][1], &gp->xv[0][2], |
360 |
|
&gp->xv[1][0], &gp->xv[1][1], &gp->xv[1][2], |
361 |
< |
&gp->xv[2][0], &gp->xv[2][1], &gp->xv[2][2]) != 12) |
361 |
> |
&gp->xv[2][0], &gp->xv[2][1], &gp->xv[2][2], |
362 |
> |
&gp->grid[0], &gp->grid[1], &gp->grid[2]) < 12) |
363 |
|
badvalue(SECTION); |
364 |
|
maxlen = 0.; |
365 |
|
for (i = 0; i < 3; i++) |
373 |
|
if ((d = vflt(GRID)) <= FTINY) |
374 |
|
badvalue(GRID); |
375 |
|
for (i = 0; i < 3; i++) |
376 |
< |
gp->grid[i] = len[i]/d + (1.-FTINY); |
377 |
< |
if (!vdef(EXPOSURE)) { |
357 |
< |
sprintf(errmsg, "%s must be defined", vnam(EXPOSURE)); |
358 |
< |
error(USER, errmsg); |
359 |
< |
} |
360 |
< |
expval = vval(EXPOSURE)[0] == '-' || vval(EXPOSURE)[0] == '+' ? |
361 |
< |
pow(2., vflt(EXPOSURE)) : vflt(EXPOSURE); |
362 |
< |
if (!vdef(OCTREE)) { |
363 |
< |
if ((vval(OCTREE) = bmalloc(strlen(froot)+5)) == NULL) |
364 |
< |
error(SYSTEM, "out of memory"); |
365 |
< |
sprintf(vval(OCTREE), "%s.oct", froot); |
366 |
< |
vdef(OCTREE)++; |
367 |
< |
} |
368 |
< |
if (!vdef(OBSTRUCTIONS)) { |
369 |
< |
vval(OBSTRUCTIONS) = "T"; |
370 |
< |
vdef(OBSTRUCTIONS)++; |
371 |
< |
} |
372 |
< |
if (!vdef(OCCUPANCY)) { |
373 |
< |
vval(OCCUPANCY) = "U"; |
374 |
< |
vdef(OCCUPANCY)++; |
375 |
< |
} |
376 |
< |
/* append rendering options */ |
377 |
< |
if (vdef(RENDER)) |
378 |
< |
rtargc += wordstring(rtargv+rtargc, vval(RENDER)); |
376 |
> |
if (gp->grid[i] <= 0) |
377 |
> |
gp->grid[i] = len[i]/d + (1.-FTINY); |
378 |
|
} |
379 |
|
|
380 |
|
|
482 |
|
p->next = freepacks; |
483 |
|
freepacks = p; |
484 |
|
} |
485 |
< |
if (n2flush > 512*RPACKSIZ*ncprocs) { |
485 |
> |
if (n2flush > 512*RPACKSIZ*nprocs) { |
486 |
|
hdflush(NULL); /* flush holodeck buffers */ |
487 |
|
n2flush = 0; |
488 |
|
} |
489 |
|
} |
490 |
|
|
491 |
|
|
493 |
– |
getradfile(rfargs) /* run rad and get needed variables */ |
494 |
– |
char *rfargs; |
495 |
– |
{ |
496 |
– |
static short mvar[] = {OCTREE,EXPOSURE,-1}; |
497 |
– |
static char tf1[] = TEMPLATE; |
498 |
– |
char tf2[64]; |
499 |
– |
char combuf[256]; |
500 |
– |
char *pippt; |
501 |
– |
register int i; |
502 |
– |
register char *cp; |
503 |
– |
/* create rad command */ |
504 |
– |
mktemp(tf1); |
505 |
– |
sprintf(tf2, "%s.rif", tf1); |
506 |
– |
sprintf(combuf, |
507 |
– |
"rad -v 0 -s -e -w %s OPTFILE=%s | egrep '^[ \t]*(NOMATCH", |
508 |
– |
rfargs, tf1); |
509 |
– |
cp = combuf; |
510 |
– |
while (*cp){ |
511 |
– |
if (*cp == '|') pippt = cp; |
512 |
– |
cp++; |
513 |
– |
} /* match unset variables */ |
514 |
– |
for (i = 0; mvar[i] >= 0; i++) |
515 |
– |
if (!vdef(mvar[i])) { |
516 |
– |
*cp++ = '|'; |
517 |
– |
strcpy(cp, vnam(mvar[i])); |
518 |
– |
while (*cp) cp++; |
519 |
– |
pippt = NULL; |
520 |
– |
} |
521 |
– |
if (pippt != NULL) |
522 |
– |
strcpy(pippt, "> /dev/null"); /* nothing to match */ |
523 |
– |
else |
524 |
– |
sprintf(cp, ")[ \t]*=' > %s", tf2); |
525 |
– |
if (system(combuf)) { |
526 |
– |
error(SYSTEM, "cannot execute rad command"); |
527 |
– |
unlink(tf2); /* clean up */ |
528 |
– |
unlink(tf1); |
529 |
– |
quit(1); |
530 |
– |
} |
531 |
– |
if (pippt == NULL) { |
532 |
– |
loadvars(tf2); /* load variables */ |
533 |
– |
unlink(tf2); |
534 |
– |
} |
535 |
– |
rtargc += wordfile(rtargv+rtargc, tf1); /* get rtrace options */ |
536 |
– |
unlink(tf1); /* clean up */ |
537 |
– |
} |
538 |
– |
|
539 |
– |
|
492 |
|
rootname(rn, fn) /* remove tail from end of fn */ |
493 |
|
register char *rn, *fn; |
494 |
|
{ |
495 |
|
char *tp, *dp; |
496 |
|
|
497 |
|
for (tp = NULL, dp = rn; *rn = *fn++; rn++) |
498 |
< |
if (ISDIRSEP(*rn)) |
498 |
> |
if (*rn == '/') |
499 |
|
dp = rn; |
500 |
|
else if (*rn == '.') |
501 |
|
tp = rn; |
544 |
|
{ |
545 |
|
int status = 0; |
546 |
|
|
547 |
< |
if (hdlist[0] != NULL) { /* flush holodeck */ |
548 |
< |
if (ncprocs > 0) { |
549 |
< |
done_packets(flush_queue()); |
550 |
< |
status = end_rtrace(); /* close rtrace */ |
551 |
< |
hdflush(NULL); |
552 |
< |
if (vdef(REPORT)) { |
553 |
< |
long fsiz, fuse; |
554 |
< |
report(0); |
603 |
< |
fsiz = hdfilen(hdlist[0]->fd); |
604 |
< |
fuse = hdfiluse(hdlist[0]->fd, 1); |
605 |
< |
fprintf(stderr, |
547 |
> |
if (hdlist[0] != NULL) { /* close holodeck */ |
548 |
> |
if (nprocs > 0) |
549 |
> |
status = done_rtrace(); /* calls hdsync() */ |
550 |
> |
if (ncprocs > 0 && vdef(REPORT)) { |
551 |
> |
long fsiz, fuse; |
552 |
> |
fsiz = hdfilen(hdlist[0]->fd); |
553 |
> |
fuse = hdfiluse(hdlist[0]->fd, 1); |
554 |
> |
fprintf(stderr, |
555 |
|
"%s: %.1f Mbyte holodeck file, %.1f%% fragmentation\n", |
556 |
< |
hdkfile, fsiz/(1024.*1024.), |
557 |
< |
100.*(fsiz-fuse)/fsiz); |
558 |
< |
} |
610 |
< |
} else |
611 |
< |
hdflush(NULL); |
556 |
> |
hdkfile, fsiz/(1024.*1024.), |
557 |
> |
100.*(fsiz-fuse)/fsiz); |
558 |
> |
} |
559 |
|
} |
560 |
|
if (orig_mode >= 0) /* reset holodeck access mode */ |
561 |
|
fchmod(hdlist[0]->fd, orig_mode); |