| 13 |
|
|
| 14 |
|
#include "ray.h" |
| 15 |
|
#include "rpaint.h" |
| 16 |
+ |
#include "otypes.h" |
| 17 |
+ |
#include "otspecial.h" |
| 18 |
+ |
#include "source.h" |
| 19 |
|
#include "random.h" |
| 20 |
|
|
| 21 |
|
#ifndef WFLUSH |
| 22 |
< |
#define WFLUSH 64 /* flush after this many rays */ |
| 22 |
> |
#define WFLUSH 64 /* flush after this many primary rays */ |
| 23 |
|
#endif |
| 24 |
+ |
#ifndef WFLUSH1 |
| 25 |
+ |
#define WFLUSH1 512 /* or this many total rays */ |
| 26 |
+ |
#endif |
| 27 |
|
|
| 28 |
+ |
|
| 29 |
|
#ifdef SMLFLT |
| 30 |
|
#define sscanvec(s,v) (sscanf(s,"%f %f %f",v,v+1,v+2)==3) |
| 31 |
|
#else |
| 32 |
|
#define sscanvec(s,v) (sscanf(s,"%lf %lf %lf",v,v+1,v+2)==3) |
| 33 |
|
#endif |
| 34 |
|
|
| 35 |
< |
static int niflush; /* flushes since newimage() */ |
| 35 |
> |
extern int ray_pnprocs; |
| 36 |
|
|
| 37 |
+ |
static RNUMBER niflush; /* flushes since newimage() */ |
| 38 |
+ |
|
| 39 |
|
int |
| 40 |
|
getrect( /* get a box */ |
| 41 |
|
char *s, |
| 117 |
|
} |
| 118 |
|
if (!direc || ourview.type == VT_PAR) { |
| 119 |
|
rayorigin(&thisray, PRIMARY, NULL, NULL); |
| 120 |
< |
if (!localhit(&thisray, &thescene)) { |
| 120 |
> |
while (localhit(&thisray, &thescene)) { |
| 121 |
> |
OBJREC *m = findmaterial(thisray.ro); |
| 122 |
> |
if (m != NULL && !istransp(m->otype) && |
| 123 |
> |
!isBSDFproxy(m) && |
| 124 |
> |
(thisray.clipset == NULL || |
| 125 |
> |
!inset(thisray.clipset, |
| 126 |
> |
thisray.ro->omod))) |
| 127 |
> |
break; /* found something */ |
| 128 |
> |
VCOPY(thisray.rorg, thisray.rop); |
| 129 |
> |
rayclear(&thisray); /* skip invisible */ |
| 130 |
> |
} |
| 131 |
> |
if (thisray.ro == NULL) { |
| 132 |
|
error(COMMAND, "not a local object"); |
| 133 |
|
return(-1); |
| 134 |
|
} |
| 153 |
|
} |
| 154 |
|
|
| 155 |
|
|
| 156 |
< |
float * /* keep consistent with COLOR typedef */ |
| 156 |
> |
COLORV * |
| 157 |
|
greyof( /* convert color to greyscale */ |
| 158 |
|
COLOR col |
| 159 |
|
) |
| 166 |
|
return(gcol); |
| 167 |
|
} |
| 168 |
|
|
| 169 |
+ |
static void |
| 170 |
+ |
recolor( /* recolor the given node */ |
| 171 |
+ |
PNODE *p |
| 172 |
+ |
) |
| 173 |
+ |
{ |
| 174 |
+ |
while (p->kid != NULL) { /* need to propogate down */ |
| 175 |
+ |
int mx = (p->xmin + p->xmax) >> 1; |
| 176 |
+ |
int my = (p->ymin + p->ymax) >> 1; |
| 177 |
+ |
int ki; |
| 178 |
+ |
if (p->x >= mx) |
| 179 |
+ |
ki = (p->y >= my) ? UR : DR; |
| 180 |
+ |
else |
| 181 |
+ |
ki = (p->y >= my) ? UL : DL; |
| 182 |
+ |
pcopy(p, p->kid+ki); |
| 183 |
+ |
p = p->kid + ki; |
| 184 |
+ |
} |
| 185 |
|
|
| 186 |
+ |
(*dev->paintr)(greyscale?greyof(p->v):p->v, |
| 187 |
+ |
p->xmin, p->ymin, p->xmax, p->ymax); |
| 188 |
+ |
} |
| 189 |
+ |
|
| 190 |
|
int |
| 191 |
|
paint( /* compute and paint a rectangle */ |
| 192 |
|
PNODE *p |
| 193 |
|
) |
| 194 |
|
{ |
| 155 |
– |
extern int ray_pnprocs; |
| 156 |
– |
static unsigned long lastflush = 0; |
| 195 |
|
static RAY thisray; |
| 196 |
|
double h, v; |
| 197 |
|
|
| 198 |
< |
if (p->xmax - p->xmin <= 0 || p->ymax - p->ymin <= 0) { /* empty */ |
| 198 |
> |
if ((p->xmax <= p->xmin) | (p->ymax <= p->ymin)) { /* empty */ |
| 199 |
|
p->x = p->xmin; |
| 200 |
|
p->y = p->ymin; |
| 201 |
|
setcolor(p->v, 0.0, 0.0, 0.0); |
| 208 |
|
if ((thisray.rmax = viewray(thisray.rorg, thisray.rdir, &ourview, |
| 209 |
|
h/hresolu, v/vresolu)) < -FTINY) { |
| 210 |
|
setcolor(thisray.rcol, 0.0, 0.0, 0.0); |
| 211 |
< |
} else { |
| 212 |
< |
int rval; |
| 211 |
> |
} else if (!ray_pnprocs) { /* immediate mode */ |
| 212 |
> |
ray_trace(&thisray); |
| 213 |
> |
} else { /* queuing mode */ |
| 214 |
> |
int rval; |
| 215 |
|
rayorigin(&thisray, PRIMARY, NULL, NULL); |
| 216 |
< |
thisray.rno = (unsigned long)p; |
| 216 |
> |
thisray.rno = (RNUMBER)p; |
| 217 |
|
rval = ray_pqueue(&thisray); |
| 218 |
|
if (!rval) |
| 219 |
|
return(0); |
| 220 |
|
if (rval < 0) |
| 221 |
|
return(-1); |
| 222 |
+ |
/* get node for returned ray */ |
| 223 |
|
p = (PNODE *)thisray.rno; |
| 224 |
|
} |
| 225 |
|
|
| 226 |
|
copycolor(p->v, thisray.rcol); |
| 227 |
|
scalecolor(p->v, exposure); |
| 228 |
|
|
| 229 |
< |
(*dev->paintr)(greyscale?greyof(p->v):p->v, |
| 189 |
< |
p->xmin, p->ymin, p->xmax, p->ymax); |
| 229 |
> |
recolor(p); /* paint it */ |
| 230 |
|
|
| 231 |
< |
if (dev->flush != NULL && raynum - lastflush >= ray_pnprocs * |
| 232 |
< |
(ambounce > 0 && niflush < WFLUSH ? niflush : WFLUSH)) { |
| 233 |
< |
lastflush = raynum; |
| 234 |
< |
(*dev->flush)(); |
| 235 |
< |
niflush++; |
| 231 |
> |
if (dev->flush != NULL) { /* shall we check for input? */ |
| 232 |
> |
static RNUMBER lastflush = 0; |
| 233 |
> |
RNUMBER counter = raynum; |
| 234 |
> |
int flushintvl; |
| 235 |
> |
if (!ray_pnprocs) { |
| 236 |
> |
counter = nrays; |
| 237 |
> |
flushintvl = WFLUSH1; |
| 238 |
> |
} else if (ambounce == 0) |
| 239 |
> |
flushintvl = ray_pnprocs*WFLUSH; |
| 240 |
> |
else if (niflush < WFLUSH) |
| 241 |
> |
flushintvl = ray_pnprocs*niflush/(ambounce*(ambounce>0)+1); |
| 242 |
> |
else |
| 243 |
> |
flushintvl = ray_pnprocs*WFLUSH/(ambounce*(ambounce>0)+1); |
| 244 |
> |
if (lastflush > counter) |
| 245 |
> |
lastflush = 0; /* counter wrapped */ |
| 246 |
> |
|
| 247 |
> |
if (counter - lastflush >= flushintvl) { |
| 248 |
> |
lastflush = counter; |
| 249 |
> |
(*dev->flush)(); |
| 250 |
> |
niflush++; |
| 251 |
> |
} |
| 252 |
|
} |
| 253 |
|
return(1); |
| 254 |
|
} |
| 255 |
|
|
| 256 |
|
|
| 257 |
|
int |
| 258 |
< |
waitrays(void) /* finish up pending rays */ |
| 258 |
> |
waitrays(void) /* finish up pending rays */ |
| 259 |
|
{ |
| 260 |
|
int nwaited = 0; |
| 261 |
|
int rval; |
| 262 |
|
RAY raydone; |
| 263 |
< |
|
| 263 |
> |
|
| 264 |
> |
if (!ray_pnprocs) /* immediate mode? */ |
| 265 |
> |
return(0); |
| 266 |
|
while ((rval = ray_presult(&raydone, 0)) > 0) { |
| 267 |
|
PNODE *p = (PNODE *)raydone.rno; |
| 268 |
|
copycolor(p->v, raydone.rcol); |
| 269 |
|
scalecolor(p->v, exposure); |
| 270 |
< |
(*dev->paintr)(greyscale?greyof(p->v):p->v, |
| 213 |
< |
p->xmin, p->ymin, p->xmax, p->ymax); |
| 270 |
> |
recolor(p); |
| 271 |
|
nwaited++; |
| 272 |
|
} |
| 273 |
|
if (rval < 0) |
| 281 |
|
char *s |
| 282 |
|
) |
| 283 |
|
{ |
| 284 |
< |
int newnp; |
| 285 |
< |
/* change in nproc? */ |
| 286 |
< |
if (s != NULL && sscanf(s, "%d", &newnp) == 1 && newnp > 0) { |
| 287 |
< |
if (!newparam) { |
| 231 |
< |
if (newnp < nproc) |
| 232 |
< |
ray_pclose(nproc - newnp); |
| 233 |
< |
else |
| 234 |
< |
ray_popen(newnp - nproc); |
| 235 |
< |
} |
| 236 |
< |
nproc = newnp; |
| 237 |
< |
} |
| 284 |
> |
int newnp = 0; |
| 285 |
> |
/* # rendering procs arg? */ |
| 286 |
> |
if (s != NULL) |
| 287 |
> |
sscanf(s, "%d", &newnp); |
| 288 |
|
/* free old image */ |
| 289 |
|
freepkids(&ptrunk); |
| 290 |
|
/* compute resolution */ |
| 299 |
|
(*dev->clear)(hresolu, vresolu); |
| 300 |
|
|
| 301 |
|
if (newparam) { /* (re)start rendering procs */ |
| 302 |
< |
ray_pclose(0); |
| 303 |
< |
ray_popen(nproc); |
| 302 |
> |
if (ray_pnprocs) |
| 303 |
> |
ray_pclose(0); /* should already be closed */ |
| 304 |
> |
if (newnp > 0) |
| 305 |
> |
nproc = newnp; |
| 306 |
> |
if (nproc > 1) |
| 307 |
> |
ray_popen(nproc); |
| 308 |
|
newparam = 0; |
| 309 |
+ |
} else if ((newnp > 0) & (newnp != nproc)) { |
| 310 |
+ |
if (newnp == 1) /* change # rendering procs */ |
| 311 |
+ |
ray_pclose(0); |
| 312 |
+ |
else if (newnp < ray_pnprocs) |
| 313 |
+ |
ray_pclose(ray_pnprocs - newnp); |
| 314 |
+ |
else |
| 315 |
+ |
ray_popen(newnp - ray_pnprocs); |
| 316 |
+ |
nproc = newnp; |
| 317 |
|
} |
| 318 |
|
niflush = 0; /* get first value */ |
| 319 |
|
paint(&ptrunk); |
| 414 |
|
|
| 415 |
|
|
| 416 |
|
void |
| 417 |
+ |
compavg( /* recompute averages */ |
| 418 |
+ |
PNODE *p |
| 419 |
+ |
) |
| 420 |
+ |
{ |
| 421 |
+ |
int i, navg; |
| 422 |
+ |
|
| 423 |
+ |
if (p->kid == NULL) |
| 424 |
+ |
return; |
| 425 |
+ |
|
| 426 |
+ |
setcolor(p->v, .0, .0, .0); |
| 427 |
+ |
navg = 0; |
| 428 |
+ |
for (i = 0; i < 4; i++) { |
| 429 |
+ |
if (p->kid[i].xmin >= p->kid[i].xmax) continue; |
| 430 |
+ |
if (p->kid[i].ymin >= p->kid[i].ymax) continue; |
| 431 |
+ |
compavg(p->kid+i); |
| 432 |
+ |
addcolor(p->v, p->kid[i].v); |
| 433 |
+ |
navg++; |
| 434 |
+ |
} |
| 435 |
+ |
if (navg > 1) |
| 436 |
+ |
scalecolor(p->v, 1./navg); |
| 437 |
+ |
} |
| 438 |
+ |
|
| 439 |
+ |
|
| 440 |
+ |
void |
| 441 |
|
scalepict( /* scale picture values */ |
| 442 |
|
PNODE *p, |
| 443 |
|
double sf |
| 531 |
|
) |
| 532 |
|
{ |
| 533 |
|
double d; |
| 448 |
– |
FVECT v1; |
| 534 |
|
VIEW nv = ourview; |
| 535 |
|
int i; |
| 536 |
|
|
| 537 |
|
spinvector(nv.vdir, ourview.vdir, ourview.vup, angle*(PI/180.)); |
| 538 |
< |
if (elev != 0.0) { |
| 539 |
< |
fcross(v1, ourview.vup, nv.vdir); |
| 540 |
< |
normalize(v1); |
| 456 |
< |
spinvector(nv.vdir, nv.vdir, v1, elev*(PI/180.)); |
| 457 |
< |
} |
| 538 |
> |
if (elev != 0.0) |
| 539 |
> |
geodesic(nv.vdir, nv.vdir, nv.vup, elev*(-PI/180.), GEOD_RAD); |
| 540 |
> |
|
| 541 |
|
if (nv.type == VT_PAR) { |
| 542 |
|
nv.horiz /= mag; |
| 543 |
|
nv.vert /= mag; |
