static void
addface( /* add a face to a cube */
register CUBE *cu,
OBJECT obj
)
{
if (o_face(objptr(obj), cu) == O_MISS)
return;
if (istree(cu->cutree)) {
CUBE cukid; /* do children */
int i, j;
cukid.cusize = cu->cusize * 0.5;
for (i = 0; i < 8; i++) {
cukid.cutree = octkid(cu->cutree, i);
for (j = 0; j < 3; j++) {
cukid.cuorg[j] = cu->cuorg[j];
if ((1<<j) & i)
cukid.cuorg[j] += cukid.cusize;
}
addface(&cukid, obj);
octkid(cu->cutree, i) = cukid.cutree;
}
return;
}
if (isempty(cu->cutree)) {
OBJECT oset[2]; /* singular set */
oset[0] = 1; oset[1] = obj;
cu->cutree = fullnode(oset);
return;
}
/* add to full node */
add2full(cu, obj);
}
static void
puttree( /* write octree to stdout in pre-order form */
register OCTREE ot
)
{
register int i;
if (istree(ot)) {
putc(OT_TREE, stdout); /* indicate tree */
for (i = 0; i < 8; i++) /* write tree */
puttree(octkid(ot, i));
} else if (isfull(ot)) {
putc(OT_FULL, stdout); /* indicate fullnode */
putfullnode(ot); /* write fullnode */
} else
putc(OT_EMPTY, stdout); /* indicate empty */
}
static int
nonsurfintree(OCTREE ot) /* check tree for modifiers */
{
OBJECT set[MAXSET+1];
register int i;
if (isempty(ot))
return(0);
if (istree(ot)) {
for (i = 0; i < 8; i++)
if (nonsurfintree(octkid(ot, i)))
return(1);
return(0);
}
objset(set, ot);
for (i = set[0]; i > 0; i-- )
if (ismodifier(objptr(set[i])->otype))
return(1);
return(0);
}
static void
puttree( /* write octree to fp in pre-order form */
register OCTREE ot,
FILE *fp
)
{
if (istree(ot)) {
register int i;
putc(OT_TREE, fp); /* indicate tree */
for (i = 0; i < 8; i++) /* write tree */
puttree(octkid(ot, i), fp);
return;
}
if (isfull(ot)) {
putc(OT_FULL, fp); /* indicate fullnode */
putfullnode(ot, fp); /* write fullnode */
return;
}
putc(OT_EMPTY, fp); /* indicate empty */
}
static void
addobject( /* add an object to a cube */
register CUBE *cu,
OBJECT obj
)
{
int inc;
inc = (*ofun[objptr(obj)->otype].funp)(objptr(obj), cu);
if (inc == O_MISS)
return; /* no intersection */
if (istree(cu->cutree)) {
CUBE cukid; /* do children */
int i, j;
cukid.cusize = cu->cusize * 0.5;
for (i = 0; i < 8; i++) {
cukid.cutree = octkid(cu->cutree, i);
for (j = 0; j < 3; j++) {
cukid.cuorg[j] = cu->cuorg[j];
if ((1<<j) & i)
cukid.cuorg[j] += cukid.cusize;
}
addobject(&cukid, obj);
octkid(cu->cutree, i) = cukid.cutree;
}
return;
}
if (isempty(cu->cutree)) {
OBJECT oset[2]; /* singular set */
oset[0] = 1; oset[1] = obj;
cu->cutree = fullnode(oset);
return;
}
/* add to full node */
add2full(cu, obj, inc);
}
static int
raymove( /* check for hit as we move */
FVECT pos, /* current position, modified herein */
OBJECT *cxs, /* checked objects, modified by checkhit */
int dirf, /* direction indicators to speed tests */
RAY *r,
CUBE *cu
)
{
int ax;
double dt, t;
if (istree(cu->cutree)) { /* recurse on subcubes */
CUBE cukid;
int br, sgn;
cukid.cusize = cu->cusize * 0.5; /* find subcube */
VCOPY(cukid.cuorg, cu->cuorg);
br = 0;
if (pos[0] >= cukid.cuorg[0]+cukid.cusize) {
cukid.cuorg[0] += cukid.cusize;
br |= 1;
}
if (pos[1] >= cukid.cuorg[1]+cukid.cusize) {
cukid.cuorg[1] += cukid.cusize;
br |= 2;
}
if (pos[2] >= cukid.cuorg[2]+cukid.cusize) {
cukid.cuorg[2] += cukid.cusize;
br |= 4;
}
for ( ; ; ) {
cukid.cutree = octkid(cu->cutree, br);
if ((ax = raymove(pos,cxs,dirf,r,&cukid)) == RAYHIT)
return(RAYHIT);
sgn = 1 << ax;
if (sgn & dirf) /* positive axis? */
if (sgn & br)
return(ax); /* overflow */
else {
cukid.cuorg[ax] += cukid.cusize;
br |= sgn;
}
else
if (sgn & br) {
cukid.cuorg[ax] -= cukid.cusize;
br &= ~sgn;
} else
return(ax); /* underflow */
}
/*NOTREACHED*/
}
if (isfull(cu->cutree)) {
if (checkhit(r, cu, cxs))
return(RAYHIT);
} else if (r->ro == &Aftplane && incube(cu, r->rop))
return(RAYHIT);
/* advance to next cube */
if (dirf&0x11) {
dt = dirf&1 ? cu->cuorg[0] + cu->cusize : cu->cuorg[0];
t = (dt - pos[0])/r->rdir[0];
ax = 0;
} else
t = FHUGE;
if (dirf&0x22) {
dt = dirf&2 ? cu->cuorg[1] + cu->cusize : cu->cuorg[1];
dt = (dt - pos[1])/r->rdir[1];
if (dt < t) {
t = dt;
ax = 1;
}
}
if (dirf&0x44) {
dt = dirf&4 ? cu->cuorg[2] + cu->cusize : cu->cuorg[2];
dt = (dt - pos[2])/r->rdir[2];
if (dt < t) {
t = dt;
ax = 2;
}
}
VSUM(pos, pos, r->rdir, t);
return(ax);
}
