--- ray/src/ot/oconv.c 1989/06/21 08:56:35 1.5 +++ ray/src/ot/oconv.c 1990/08/29 13:01:16 1.10 @@ -39,7 +39,9 @@ int nfiles = 0; /* number of object files */ double mincusize; /* minimum cube size from resolu */ +int (*addobjnotify[])() = {NULL}; /* new object notifier functions */ + main(argc, argv) /* convert object files to an octree */ int argc; char **argv; @@ -177,15 +179,15 @@ char *s; eputs(s) /* put string to stderr */ register char *s; { - static int inline = 0; + static int inln = 0; - if (!inline++) { + if (!inln++) { fputs(progname, stderr); fputs(": ", stderr); } fputs(s, stderr); if (*s && s[strlen(s)-1] == '\n') - inline = 0; + inln = 0; } @@ -193,15 +195,18 @@ addobject(cu, obj) /* add an object to a cube */ register CUBE *cu; OBJECT obj; { +#define nexti(n) ((ndx += cnt*cnt++)%(n)) + static unsigned long ndx; + static unsigned int cnt; CUBE cukid; OCTREE ot; OBJECT oset[MAXSET+1]; - int in; + int in, k; register int i, j; in = (*ofun[objptr(obj)->otype].funp)(objptr(obj), cu); - if (!in) + if (in == O_MISS) return; /* no intersection */ if (istree(cu->cutree)) { @@ -228,7 +233,7 @@ OBJECT obj; objset(oset, cu->cutree); cukid.cusize = cu->cusize * 0.5; - if (in == 2 || oset[0] < objlim || cukid.cusize < mincusize) { + if (in==O_IN || oset[0] < objlim || cukid.cusize < mincusize) { /* add to set */ if (oset[0] >= MAXSET) { sprintf(errmsg, @@ -250,12 +255,24 @@ OBJECT obj; if ((1<otype)) + addobject(&cukid, oset[j]); + /* then this object */ addobject(&cukid, obj); + /* volumes last */ + k = nexti(oset[0]); /* random start */ + for (j = k+1; j <= oset[0]; j++) + if (isvolume(objptr(oset[j])->otype)) + addobject(&cukid, oset[j]); + for (j = 1; j <= k; j++) + if (isvolume(objptr(oset[j])->otype)) + addobject(&cukid, oset[j]); octkid(ot, i) = cukid.cutree; } cu->cutree = ot; } } +#undef nexti }