| 361 |
|
|
| 362 |
|
extern void |
| 363 |
|
raycontrib( /* compute (cumulative) ray contribution */ |
| 364 |
< |
double rc[3], |
| 364 |
> |
RREAL rc[3], |
| 365 |
|
const RAY *r, |
| 366 |
|
int flags |
| 367 |
|
) |
| 520 |
|
if (r->rmax > FTINY) { /* except aft plane if one */ |
| 521 |
|
r->ro = &Aftplane; |
| 522 |
|
r->rot = r->rmax; |
| 523 |
< |
for (i = 0; i < 3; i++) |
| 524 |
< |
r->rop[i] = r->rorg[i] + r->rot*r->rdir[i]; |
| 523 |
> |
VSUM(r->rop, r->rorg, r->rdir, r->rot); |
| 524 |
|
} |
| 525 |
|
/* find global cube entrance point */ |
| 526 |
|
t = 0.0; |
| 543 |
|
if (t >= r->rot) /* clipped already */ |
| 544 |
|
return(0); |
| 545 |
|
/* advance position */ |
| 546 |
< |
for (i = 0; i < 3; i++) |
| 548 |
< |
curpos[i] += r->rdir[i]*t; |
| 546 |
> |
VSUM(curpos, curpos, r->rdir, t); |
| 547 |
|
|
| 548 |
|
if (!incube(scene, curpos)) /* non-intersecting ray */ |
| 549 |
|
return(0); |
| 634 |
|
ax = 2; |
| 635 |
|
} |
| 636 |
|
} |
| 637 |
< |
pos[0] += r->rdir[0]*t; |
| 640 |
< |
pos[1] += r->rdir[1]*t; |
| 641 |
< |
pos[2] += r->rdir[2]*t; |
| 637 |
> |
VSUM(pos, pos, r->rdir, t); |
| 638 |
|
return(ax); |
| 639 |
|
} |
| 640 |
|
|
