31 |
|
char *octname; /* octree name */ |
32 |
|
char *sigerr[NSIG]; /* signal error messages */ |
33 |
|
char *errfile = NULL; /* error output file */ |
34 |
< |
|
34 |
> |
#ifdef SSKIPOPT |
35 |
> |
int *ssndx_map = NULL; /* loaded source index map */ |
36 |
> |
float *scorr_map = NULL; /* loaded source correction map */ |
37 |
> |
#endif |
38 |
|
extern time_t time(); |
39 |
|
extern time_t tstart; /* start time */ |
40 |
|
|
98 |
|
char *recover = NULL; |
99 |
|
char *outfile = NULL; |
100 |
|
char *zfile = NULL; |
101 |
+ |
#ifdef SSKIPOPT |
102 |
+ |
char *sskipspec = NULL; |
103 |
+ |
char *scorrfile = NULL; |
104 |
+ |
#endif |
105 |
|
int loadflags = ~IO_FILES; |
106 |
|
int seqstart = 0; |
107 |
|
int persist = 0; |
248 |
|
check(2,"s"); |
249 |
|
recover = argv[++i]; |
250 |
|
break; |
251 |
< |
case 't': /* timer */ |
252 |
< |
check(2,"i"); |
253 |
< |
ralrm = atoi(argv[++i]); |
251 |
> |
#ifdef SSKIPOPT |
252 |
> |
case 'd': |
253 |
> |
switch (argv[i][2]) { |
254 |
> |
case 'S': /* direct skip */ |
255 |
> |
check(3,"s"); |
256 |
> |
sskipspec = argv[++i]; |
257 |
> |
break; |
258 |
> |
case 'C': /* direct correction */ |
259 |
> |
check(3,"s"); |
260 |
> |
scorrfile = argv[++i]; |
261 |
> |
break; |
262 |
> |
default: |
263 |
> |
goto badopt; |
264 |
> |
} |
265 |
|
break; |
266 |
+ |
#endif |
267 |
|
#ifdef PERSIST |
268 |
|
case 'P': /* persist file */ |
269 |
|
if (argv[i][2] == 'P') { |
276 |
|
persistfile(argv[++i]); |
277 |
|
break; |
278 |
|
#endif |
279 |
+ |
case 't': /* timer */ |
280 |
+ |
check(2,"i"); |
281 |
+ |
ralrm = atoi(argv[++i]); |
282 |
+ |
break; |
283 |
|
case 'w': /* warnings */ |
284 |
|
rval = erract[WARNING].pf != NULL; |
285 |
|
check_bool(2,rval); |
377 |
|
ray_init_pmap(); /* PMAP: set up & load photon maps */ |
378 |
|
|
379 |
|
marksources(); /* find and mark sources */ |
380 |
< |
|
380 |
> |
#ifdef SSKIPOPT |
381 |
> |
if (sskipspec != NULL) { /* load source skip maps? */ |
382 |
> |
if ((seqstart > 0) | (persist > 0)) |
383 |
> |
error(USER, |
384 |
> |
"source skip map incompatible with -S/-P/-PP"); |
385 |
> |
if (srcskip_open(sskipspec, scorrfile) <= 0) |
386 |
> |
quit(1); |
387 |
> |
if ((ssndx_map = srcskip_ndxmap()) == NULL) |
388 |
> |
quit(1); |
389 |
> |
if (scorrfile != NULL && |
390 |
> |
(scorr_map = srcskip_corrmap()) == NULL) |
391 |
> |
quit(1); |
392 |
> |
srcskip_close(); /* leaves maps loaded */ |
393 |
> |
} |
394 |
> |
#endif |
395 |
|
setambient(); /* initialize ambient calculation */ |
396 |
|
|
397 |
|
fflush(stdout); /* in case we're duplicating header */ |
462 |
|
#undef check_bool |
463 |
|
} |
464 |
|
|
465 |
+ |
#ifdef SSKIPOPT |
466 |
+ |
void /* source skip parameters for primary */ |
467 |
+ |
sskip_ray(RAY *r, double h, double v) |
468 |
+ |
{ |
469 |
+ |
size_t ssi; |
470 |
+ |
|
471 |
+ |
if (ssndx_map == NULL) |
472 |
+ |
return; /* nothing to do */ |
473 |
+ |
|
474 |
+ |
/* rows counted from the top! */ |
475 |
+ |
ssi = (size_t)((1.-v)*sskip_dim[1])*sskip_dim[0] + |
476 |
+ |
(size_t)(h*sskip_dim[0]); |
477 |
+ |
|
478 |
+ |
if (ssi >= sskip_dim[0]*sskip_dim[1]) |
479 |
+ |
error(CONSISTENCY, "bad index in sskip_ray()"); |
480 |
+ |
|
481 |
+ |
r->rsrc = ssndx_map[ssi]; |
482 |
+ |
if (scorr_map != NULL) |
483 |
+ |
r->scorr = scorr_map[ssi]; |
484 |
+ |
} |
485 |
+ |
#endif |
486 |
|
|
487 |
|
void |
488 |
|
wputs( /* warning output function */ |