63 |
|
static void setoutput(char *vs); |
64 |
|
extern void tranotify(OBJECT obj); |
65 |
|
static void bogusray(void); |
66 |
< |
static void rad(FVECT org, FVECT dir, double dmax); |
67 |
< |
static void irrad(FVECT org, FVECT dir); |
66 |
> |
static void raycast(RAY *r); |
67 |
> |
static void rayirrad(RAY *r); |
68 |
> |
static void rtcompute(FVECT org, FVECT dir, double dmax); |
69 |
|
static int printvals(RAY *r); |
70 |
|
static int getvec(FVECT vec, int fmt, FILE *fp); |
71 |
|
static void tabin(RAY *r); |
109 |
|
|
110 |
|
extern void |
111 |
|
rtrace( /* trace rays from file */ |
112 |
< |
char *fname |
112 |
> |
char *fname, |
113 |
> |
int nproc |
114 |
|
) |
115 |
|
{ |
116 |
|
unsigned long vcount = (hresolu > 1) ? (unsigned long)hresolu*vresolu |
130 |
|
SET_FILE_BINARY(fp); |
131 |
|
/* set up output */ |
132 |
|
setoutput(outvals); |
133 |
+ |
if (imm_irrad) |
134 |
+ |
castonly = 0; |
135 |
+ |
else if (castonly) |
136 |
+ |
nproc = 1; /* don't bother multiprocessing */ |
137 |
|
switch (outform) { |
138 |
|
case 'a': putreal = puta; break; |
139 |
|
case 'f': putreal = putf; break; |
145 |
|
default: |
146 |
|
error(CONSISTENCY, "botched output format"); |
147 |
|
} |
148 |
< |
if (ray_pnprocs > 1) |
148 |
> |
if (nproc > 1) { /* start multiprocessing */ |
149 |
> |
ray_popen(nproc); |
150 |
|
ray_fifo_out = printvals; |
151 |
+ |
} |
152 |
|
if (hresolu > 0) { |
153 |
|
if (vresolu > 0) |
154 |
|
fprtresolu(hresolu, vresolu, stdout); |
160 |
|
|
161 |
|
d = normalize(direc); |
162 |
|
if (d == 0.0) { /* zero ==> flush */ |
163 |
< |
if (ray_pnprocs > 1 && ray_fifo_flush() < 0) |
163 |
> |
if (nproc > 1 && ray_fifo_flush() < 0) |
164 |
|
error(USER, "lost children"); |
165 |
|
bogusray(); |
166 |
|
if (--nextflush <= 0 || !vcount) { |
168 |
|
nextflush = hresolu; |
169 |
|
} |
170 |
|
} else { /* compute and print */ |
171 |
< |
if (imm_irrad) |
164 |
< |
irrad(orig, direc); |
165 |
< |
else |
166 |
< |
rad(orig, direc, lim_dist ? d : 0.0); |
171 |
> |
rtcompute(orig, direc, lim_dist ? d : 0.0); |
172 |
|
/* flush if time */ |
173 |
|
if (!--nextflush) { |
174 |
< |
if (ray_pnprocs > 1 && ray_fifo_flush() < 0) |
174 |
> |
if (nproc > 1 && ray_fifo_flush() < 0) |
175 |
|
error(USER, "lost children"); |
176 |
|
fflush(stdout); |
177 |
|
nextflush = hresolu; |
182 |
|
if (vcount && !--vcount) /* check for end */ |
183 |
|
break; |
184 |
|
} |
185 |
< |
if (ray_pnprocs > 1 && ray_fifo_flush() < 0) |
186 |
< |
error(USER, "unable to complete processing"); |
185 |
> |
if (nproc > 1) { /* clean up children */ |
186 |
> |
if (ray_fifo_flush() < 0) |
187 |
> |
error(USER, "unable to complete processing"); |
188 |
> |
ray_pclose(0); |
189 |
> |
} |
190 |
|
if (fflush(stdout) < 0) |
191 |
|
error(SYSTEM, "write error"); |
192 |
|
if (vcount) |
301 |
|
|
302 |
|
|
303 |
|
static void |
304 |
< |
rad( /* compute and print ray value(s) */ |
305 |
< |
FVECT org, |
298 |
< |
FVECT dir, |
299 |
< |
double dmax |
304 |
> |
raycast( /* compute first ray intersection only */ |
305 |
> |
RAY *r |
306 |
|
) |
307 |
|
{ |
308 |
< |
VCOPY(thisray.rorg, org); |
309 |
< |
VCOPY(thisray.rdir, dir); |
310 |
< |
thisray.rmax = dmax; |
311 |
< |
if (!castonly && ray_pnprocs > 1) { |
312 |
< |
if (ray_fifo_in(&thisray) < 0) |
313 |
< |
error(USER, "lost children"); |
308 |
< |
return; |
308 |
> |
if (!localhit(r, &thescene)) { |
309 |
> |
if (r->ro == &Aftplane) { /* clipped */ |
310 |
> |
r->ro = NULL; |
311 |
> |
r->rot = FHUGE; |
312 |
> |
} else |
313 |
> |
sourcehit(r); |
314 |
|
} |
310 |
– |
rayorigin(&thisray, PRIMARY, NULL, NULL); |
311 |
– |
if (castonly) { |
312 |
– |
if (!localhit(&thisray, &thescene)) { |
313 |
– |
if (thisray.ro == &Aftplane) { /* clipped */ |
314 |
– |
thisray.ro = NULL; |
315 |
– |
thisray.rot = FHUGE; |
316 |
– |
} else |
317 |
– |
sourcehit(&thisray); |
318 |
– |
} |
319 |
– |
} else |
320 |
– |
ray_trace(&thisray); |
321 |
– |
printvals(&thisray); |
315 |
|
} |
316 |
|
|
317 |
|
|
318 |
|
static void |
319 |
< |
irrad( /* compute immediate irradiance value */ |
319 |
> |
rayirrad( /* compute irradiance rather than radiance */ |
320 |
> |
RAY *r |
321 |
> |
) |
322 |
> |
{ |
323 |
> |
r->rot = 1e-5; /* pretend we hit surface */ |
324 |
> |
VSUM(r->rop, r->rorg, r->rdir, r->rot); |
325 |
> |
r->ron[0] = -r->rdir[0]; |
326 |
> |
r->ron[1] = -r->rdir[1]; |
327 |
> |
r->ron[2] = -r->rdir[2]; |
328 |
> |
r->rod = 1.0; |
329 |
> |
/* compute result */ |
330 |
> |
(*ofun[Lamb.otype].funp)(&Lamb, r); |
331 |
> |
} |
332 |
> |
|
333 |
> |
|
334 |
> |
static void |
335 |
> |
rtcompute( /* compute and print ray value(s) */ |
336 |
|
FVECT org, |
337 |
< |
FVECT dir |
337 |
> |
FVECT dir, |
338 |
> |
double dmax |
339 |
|
) |
340 |
|
{ |
341 |
< |
VSUM(thisray.rorg, org, dir, 1.1e-4); |
342 |
< |
thisray.rdir[0] = -dir[0]; |
343 |
< |
thisray.rdir[1] = -dir[1]; |
344 |
< |
thisray.rdir[2] = -dir[2]; |
345 |
< |
thisray.rmax = 0.0; |
341 |
> |
if (imm_irrad) { /* set up ray */ |
342 |
> |
VSUM(thisray.rorg, org, dir, 1.1e-4); |
343 |
> |
thisray.rdir[0] = -dir[0]; |
344 |
> |
thisray.rdir[1] = -dir[1]; |
345 |
> |
thisray.rdir[2] = -dir[2]; |
346 |
> |
thisray.rmax = 0.0; |
347 |
> |
} else { |
348 |
> |
VCOPY(thisray.rorg, org); |
349 |
> |
VCOPY(thisray.rdir, dir); |
350 |
> |
thisray.rmax = dmax; |
351 |
> |
} |
352 |
|
rayorigin(&thisray, PRIMARY, NULL, NULL); |
353 |
< |
/* pretend we hit surface */ |
354 |
< |
thisray.rot = 1e-5; |
355 |
< |
thisray.rod = 1.0; |
356 |
< |
VCOPY(thisray.ron, dir); |
357 |
< |
VSUM(thisray.rop, org, dir, 1e-4); |
358 |
< |
/* compute and print */ |
359 |
< |
(*ofun[Lamb.otype].funp)(&Lamb, &thisray); |
353 |
> |
/* special case evaluators */ |
354 |
> |
if (castonly) |
355 |
> |
thisray.revf = raycast; |
356 |
> |
else if (imm_irrad) |
357 |
> |
thisray.revf = rayirrad; |
358 |
> |
|
359 |
> |
if (ray_pnprocs > 1) { /* multiprocessing FIFO? */ |
360 |
> |
if (ray_fifo_in(&thisray) < 0) |
361 |
> |
error(USER, "lost children"); |
362 |
> |
return; |
363 |
> |
} |
364 |
> |
samplendx++; /* else do it ourselves */ |
365 |
> |
rayvalue(&thisray); |
366 |
|
printvals(&thisray); |
367 |
|
} |
368 |
|
|