9 |
|
*/ |
10 |
|
|
11 |
|
#include "rholo.h" |
12 |
+ |
#include "random.h" |
13 |
|
#include "paths.h" |
14 |
|
#include <signal.h> |
15 |
|
#include <sys/types.h> |
26 |
|
|
27 |
|
int nowarn = 0; /* turn warnings off? */ |
28 |
|
|
28 |
– |
double expval = 1.; /* global exposure value */ |
29 |
– |
|
29 |
|
int ncprocs = 0; /* desired number of compute processes */ |
30 |
|
|
31 |
|
char *outdev = NULL; /* output device name */ |
32 |
|
|
33 |
+ |
int readinp = 0; /* read commands from stdin */ |
34 |
+ |
|
35 |
|
time_t starttime; /* time we got started */ |
36 |
|
time_t endtime; /* time we should end by */ |
37 |
|
time_t reporttime; /* time for next report */ |
62 |
|
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; |
118 |
|
creatholo(&hdg); |
119 |
|
} else { /* else load holodeck */ |
120 |
|
loadholo(); |
121 |
< |
if (vdef(RIF)) /* load RIF if any */ |
121 |
> |
/* check settings */ |
122 |
> |
checkvalues(); |
123 |
> |
/* load RIF if any */ |
124 |
> |
if (vdef(RIF)) |
125 |
|
getradfile(vval(RIF)); |
126 |
+ |
/* set defaults */ |
127 |
+ |
setdefaults(NULL); |
128 |
|
} |
129 |
|
/* initialize */ |
130 |
|
initrholo(); |
131 |
< |
/* run */ |
131 |
> |
/* main loop */ |
132 |
|
while (rholo()) |
133 |
|
; |
134 |
|
/* done */ |
192 |
|
{ |
193 |
|
extern int global_packet(); |
194 |
|
register int i; |
195 |
+ |
/* close holodeck on exec() */ |
196 |
+ |
fcntl(hdlist[0]->fd, F_SETFD, FD_CLOEXEC); |
197 |
|
|
198 |
|
if (outdev != NULL) /* open output device */ |
199 |
|
disp_open(outdev); |
345 |
|
sprintf(errmsg, "ignoring all but first %s", vnam(SECTION)); |
346 |
|
error(WARNING, errmsg); |
347 |
|
} |
336 |
– |
if (sscanf(vval(SECTION), |
337 |
– |
"%lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf", |
338 |
– |
&gp->orig[0], &gp->orig[1], &gp->orig[2], |
339 |
– |
&gp->xv[0][0], &gp->xv[0][1], &gp->xv[0][2], |
340 |
– |
&gp->xv[1][0], &gp->xv[1][1], &gp->xv[1][2], |
341 |
– |
&gp->xv[2][0], &gp->xv[2][1], &gp->xv[2][2]) != 12) |
342 |
– |
badvalue(SECTION); |
343 |
– |
maxlen = 0.; |
344 |
– |
for (i = 0; i < 3; i++) |
345 |
– |
if ((len[i] = VLEN(gp->xv[i])) > maxlen) |
346 |
– |
maxlen = len[i]; |
347 |
– |
if (!vdef(GRID)) { |
348 |
– |
sprintf(buf, "%.4f", maxlen/8.); |
349 |
– |
vval(GRID) = savqstr(buf); |
350 |
– |
vdef(GRID)++; |
351 |
– |
} |
352 |
– |
if ((d = vflt(GRID)) <= FTINY) |
353 |
– |
badvalue(GRID); |
354 |
– |
for (i = 0; i < 3; i++) |
355 |
– |
gp->grid[i] = len[i]/d + (1.-FTINY); |
356 |
– |
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); |
348 |
|
if (!vdef(OCTREE)) { |
349 |
|
if ((vval(OCTREE) = bmalloc(strlen(froot)+5)) == NULL) |
350 |
|
error(SYSTEM, "out of memory"); |
355 |
|
vval(OBSTRUCTIONS) = "T"; |
356 |
|
vdef(OBSTRUCTIONS)++; |
357 |
|
} |
358 |
+ |
if (!vdef(VDIST)) { |
359 |
+ |
vval(VDIST) = "F"; |
360 |
+ |
vdef(VDIST)++; |
361 |
+ |
} |
362 |
|
if (!vdef(OCCUPANCY)) { |
363 |
|
vval(OCCUPANCY) = "U"; |
364 |
|
vdef(OCCUPANCY)++; |
366 |
|
/* append rendering options */ |
367 |
|
if (vdef(RENDER)) |
368 |
|
rtargc += wordstring(rtargv+rtargc, vval(RENDER)); |
369 |
+ |
|
370 |
+ |
if (gp == NULL) /* already initialized? */ |
371 |
+ |
return; |
372 |
+ |
/* set grid parameters */ |
373 |
+ |
gp->grid[0] = gp->grid[1] = gp->grid[2] = 0; |
374 |
+ |
if (sscanf(vval(SECTION), |
375 |
+ |
"%lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %hd %hd %hd", |
376 |
+ |
&gp->orig[0], &gp->orig[1], &gp->orig[2], |
377 |
+ |
&gp->xv[0][0], &gp->xv[0][1], &gp->xv[0][2], |
378 |
+ |
&gp->xv[1][0], &gp->xv[1][1], &gp->xv[1][2], |
379 |
+ |
&gp->xv[2][0], &gp->xv[2][1], &gp->xv[2][2], |
380 |
+ |
&gp->grid[0], &gp->grid[1], &gp->grid[2]) < 12) |
381 |
+ |
badvalue(SECTION); |
382 |
+ |
maxlen = 0.; |
383 |
+ |
for (i = 0; i < 3; i++) |
384 |
+ |
if ((len[i] = VLEN(gp->xv[i])) > maxlen) |
385 |
+ |
maxlen = len[i]; |
386 |
+ |
if (!vdef(GRID)) { |
387 |
+ |
sprintf(buf, "%.4f", maxlen/8.); |
388 |
+ |
vval(GRID) = savqstr(buf); |
389 |
+ |
vdef(GRID)++; |
390 |
+ |
} |
391 |
+ |
if ((d = vflt(GRID)) <= FTINY) |
392 |
+ |
badvalue(GRID); |
393 |
+ |
for (i = 0; i < 3; i++) |
394 |
+ |
if (gp->grid[i] <= 0) |
395 |
+ |
gp->grid[i] = len[i]/d + (1.-FTINY); |
396 |
|
} |
397 |
|
|
398 |
|
|
480 |
|
done_packets(pl) /* handle finished packets */ |
481 |
|
PACKET *pl; |
482 |
|
{ |
483 |
< |
static int nunflushed = 0; |
483 |
> |
static int n2flush = 0; |
484 |
|
register PACKET *p; |
485 |
|
|
486 |
|
while (pl != NULL) { |
491 |
|
p->nr*sizeof(RAYVAL)); |
492 |
|
if (outdev != NULL) /* display it */ |
493 |
|
disp_packet((PACKHEAD *)p); |
494 |
< |
else |
495 |
< |
nunflushed += p->nr; |
494 |
> |
if (hdcachesize <= 0) /* manual flushing */ |
495 |
> |
n2flush += p->nr; |
496 |
|
nraysdone += p->nr; |
497 |
|
npacksdone++; |
498 |
|
} |
500 |
|
p->next = freepacks; |
501 |
|
freepacks = p; |
502 |
|
} |
503 |
< |
if (nunflushed >= 256*RPACKSIZ) { |
503 |
> |
if (n2flush > 512*RPACKSIZ*ncprocs) { |
504 |
|
hdflush(NULL); /* flush holodeck buffers */ |
505 |
< |
nunflushed = 0; |
505 |
> |
n2flush = 0; |
506 |
|
} |
507 |
|
} |
508 |
|
|
510 |
|
getradfile(rfargs) /* run rad and get needed variables */ |
511 |
|
char *rfargs; |
512 |
|
{ |
513 |
< |
static short mvar[] = {OCTREE,EXPOSURE,-1}; |
513 |
> |
static short mvar[] = {OCTREE,-1}; |
514 |
|
static char tf1[] = TEMPLATE; |
515 |
|
char tf2[64]; |
516 |
|
char combuf[256]; |