extern int
o_instance( /* compute ray intersection with octree */
OBJREC *o,
register RAY *r
)
{
RAY rcont;
double d;
register INSTANCE *ins;
register int i;
/* get the octree */
ins = getinstance(o, IO_ALL);
/* copy and transform ray */
rcont = *r;
multp3(rcont.rorg, r->rorg, ins->x.b.xfm);
multv3(rcont.rdir, r->rdir, ins->x.b.xfm);
for (i = 0; i < 3; i++)
rcont.rdir[i] /= ins->x.b.sca;
rcont.rmax *= ins->x.b.sca;
/* clear and trace it */
rayclear(&rcont);
if (!localhit(&rcont, &ins->obj->scube))
return(0); /* missed */
if (rcont.rot * ins->x.f.sca >= r->rot)
return(0); /* not close enough */
if (o->omod != OVOID) { /* if we have modifier, use it */
r->ro = o;
r->rox = NULL;
} else { /* else use theirs */
r->ro = rcont.ro;
if (rcont.rox != NULL) {
newrayxf(r); /* allocate transformation */
/* NOTE: r->rox may equal rcont.rox! */
multmat4(r->rox->f.xfm, rcont.rox->f.xfm, ins->x.f.xfm);
r->rox->f.sca = rcont.rox->f.sca * ins->x.f.sca;
multmat4(r->rox->b.xfm, ins->x.b.xfm, rcont.rox->b.xfm);
r->rox->b.sca = ins->x.b.sca * rcont.rox->b.sca;
} else
r->rox = &ins->x;
}
/* transform it back */
r->rot = rcont.rot * ins->x.f.sca;
multp3(r->rop, rcont.rop, ins->x.f.xfm);
multv3(r->ron, rcont.ron, ins->x.f.xfm);
multv3(r->pert, rcont.pert, ins->x.f.xfm);
d = 1./ins->x.f.sca;
for (i = 0; i < 3; i++) {
r->ron[i] *= d;
r->pert[i] *= d;
}
r->rod = rcont.rod;
r->uv[0] = rcont.uv[0];
r->uv[1] = rcont.uv[1];
/* return hit */
return(1);
}
DLLEXPORT int DLLCALL tjDestroy(tjhandle handle)
{
getinstance(handle);
if(setjmp(this->jerr.setjmp_buffer)) return -1;
if(this->init&COMPRESS) jpeg_destroy_compress(cinfo);
if(this->init&DECOMPRESS) jpeg_destroy_decompress(dinfo);
free(this);
return 0;
}
int
o_instance( /* determine if instance intersects */
OBJREC *o,
CUBE *cu
)
{
INSTANCE *ins;
/* get octree bounds */
ins = getinstance(o, IO_BOUNDS);
/* call o_cube to do the work */
return(o_cube(&ins->obj->scube, &ins->x, cu));
}
void
add2bbox( /* expand bounding box to fit object */
register OBJREC *o,
FVECT bbmin,
FVECT bbmax
)
{
CONE *co;
FACE *fo;
INSTANCE *io;
MESHINST *mi;
FVECT v;
register int i, j;
switch (o->otype) {
case OBJ_SPHERE:
case OBJ_BUBBLE:
if (o->oargs.nfargs != 4)
objerror(o, USER, "bad arguments");
for (i = 0; i < 3; i++) {
VCOPY(v, o->oargs.farg);
v[i] -= o->oargs.farg[3];
point2bbox(v, bbmin, bbmax);
v[i] += 2.0 * o->oargs.farg[3];
point2bbox(v, bbmin, bbmax);
}
break;
case OBJ_FACE:
fo = getface(o);
j = fo->nv;
while (j--)
point2bbox(VERTEX(fo,j), bbmin, bbmax);
break;
case OBJ_CONE:
case OBJ_CUP:
case OBJ_CYLINDER:
case OBJ_TUBE:
case OBJ_RING:
co = getcone(o, 0);
if (o->otype != OBJ_RING)
circle2bbox(CO_P0(co), co->ad, CO_R0(co), bbmin, bbmax);
circle2bbox(CO_P1(co), co->ad, CO_R1(co), bbmin, bbmax);
break;
case OBJ_INSTANCE:
io = getinstance(o, IO_BOUNDS);
for (j = 0; j < 8; j++) {
for (i = 0; i < 3; i++) {
v[i] = io->obj->scube.cuorg[i];
if (j & 1<<i)
v[i] += io->obj->scube.cusize;
}
multp3(v, v, io->x.f.xfm);
point2bbox(v, bbmin, bbmax);
}
break;
case OBJ_MESH:
mi = getmeshinst(o, IO_BOUNDS);
for (j = 0; j < 8; j++) {
for (i = 0; i < 3; i++) {
v[i] = mi->msh->mcube.cuorg[i];
if (j & 1<<i)
v[i] += mi->msh->mcube.cusize;
}
multp3(v, v, mi->x.f.xfm);
point2bbox(v, bbmin, bbmax);
}
break;
}
}
extern int
load_os( /* load associated data for object */
register OBJREC *op
)
{
DATARRAY *dp;
switch (op->otype) {
case OBJ_FACE: /* polygon */
getface(op);
return(1);
case OBJ_CONE: /* cone */
case OBJ_RING: /* disk */
case OBJ_CYLINDER: /* cylinder */
case OBJ_CUP: /* inverted cone */
case OBJ_TUBE: /* inverted cylinder */
getcone(op, 1);
return(1);
case OBJ_INSTANCE: /* octree instance */
getinstance(op, IO_ALL);
return(1);
case OBJ_MESH: /* mesh instance */
getmeshinst(op, IO_ALL);
return(1);
case PAT_CPICT: /* color picture */
if (op->oargs.nsargs < 4)
goto sargerr;
getpict(op->oargs.sarg[3]);
getfunc(op, 4, 0x3<<5, 0);
return(1);
case PAT_CDATA: /* color data */
dp = getdata(op->oargs.sarg[3]);
getdata(op->oargs.sarg[4]);
getdata(op->oargs.sarg[5]);
getfunc(op, 6, ((1<<dp->nd)-1)<<7, 0);
return(1);
case PAT_BDATA: /* brightness data */
if (op->oargs.nsargs < 2)
goto sargerr;
dp = getdata(op->oargs.sarg[1]);
getfunc(op, 2, ((1<<dp->nd)-1)<<3, 0);
return(1);
case PAT_BFUNC: /* brightness function */
getfunc(op, 1, 0x1, 0);
return(1);
case PAT_CFUNC: /* color function */
getfunc(op, 3, 0x7, 0);
return(1);
case TEX_DATA: /* texture data */
if (op->oargs.nsargs < 6)
goto sargerr;
dp = getdata(op->oargs.sarg[3]);
getdata(op->oargs.sarg[4]);
getdata(op->oargs.sarg[5]);
getfunc(op, 6, ((1<<dp->nd)-1)<<7, 1);
return(1);
case TEX_FUNC: /* texture function */
getfunc(op, 3, 0x7, 1);
return(1);
case MIX_DATA: /* mixture data */
dp = getdata(op->oargs.sarg[3]);
getfunc(op, 4, ((1<<dp->nd)-1)<<5, 0);
return(1);
case MIX_PICT: /* mixture picture */
getpict(op->oargs.sarg[3]);
getfunc(op, 4, 0x3<<5, 0);
return(1);
case MIX_FUNC: /* mixture function */
getfunc(op, 3, 0x4, 0);
return(1);
case MAT_PLASTIC2: /* anisotropic plastic */
case MAT_METAL2: /* anisotropic metal */
getfunc(op, 3, 0x7, 1);
return(1);
case MAT_BRTDF: /* BRDTfunc material */
getfunc(op, 9, 0x3f, 0);
return(1);
case MAT_BSDF: /* BSDF material */
if (op->oargs.nsargs < 6)
goto sargerr;
getfunc(op, 5, 0x1d, 1);
loadBSDF(op->oargs.sarg[1]);
return(1);
case MAT_PDATA: /* plastic BRDF data */
case MAT_MDATA: /* metal BRDF data */
case MAT_TDATA: /* trans BRDF data */
if (op->oargs.nsargs < 2)
goto sargerr;
getdata(op->oargs.sarg[1]);
getfunc(op, 2, 0, 0);
return(1);
case MAT_PFUNC: /* plastic BRDF func */
case MAT_MFUNC: /* metal BRDF func */
case MAT_TFUNC: /* trans BRDF func */
getfunc(op, 1, 0, 0);
return(1);
case MAT_DIRECT1: /* prism1 material */
getfunc(op, 4, 0xf, 1);
return(1);
case MAT_DIRECT2: /* prism2 material */
getfunc(op, 8, 0xff, 1);
return(1);
}
/* nothing to load for the remaining types */
return(0);
sargerr:
objerror(op, USER, "too few string arguments");
return 0; /* pro forma return */
}
