| 1 | greg | 1.1 | #ifndef lint | 
| 2 | greg | 2.13 | static const char       RCSid[] = "$Id$"; | 
| 3 | greg | 1.1 | #endif | 
| 4 |  |  | /* | 
| 5 |  |  | *  readoct.c - routines to read octree information. | 
| 6 | greg | 2.13 | */ | 
| 7 |  |  |  | 
| 8 |  |  | /* ==================================================================== | 
| 9 |  |  | * The Radiance Software License, Version 1.0 | 
| 10 |  |  | * | 
| 11 |  |  | * Copyright (c) 1990 - 2002 The Regents of the University of California, | 
| 12 |  |  | * through Lawrence Berkeley National Laboratory.   All rights reserved. | 
| 13 |  |  | * | 
| 14 |  |  | * Redistribution and use in source and binary forms, with or without | 
| 15 |  |  | * modification, are permitted provided that the following conditions | 
| 16 |  |  | * are met: | 
| 17 |  |  | * | 
| 18 |  |  | * 1. Redistributions of source code must retain the above copyright | 
| 19 |  |  | *         notice, this list of conditions and the following disclaimer. | 
| 20 | greg | 1.1 | * | 
| 21 | greg | 2.13 | * 2. Redistributions in binary form must reproduce the above copyright | 
| 22 |  |  | *       notice, this list of conditions and the following disclaimer in | 
| 23 |  |  | *       the documentation and/or other materials provided with the | 
| 24 |  |  | *       distribution. | 
| 25 |  |  | * | 
| 26 |  |  | * 3. The end-user documentation included with the redistribution, | 
| 27 |  |  | *           if any, must include the following acknowledgment: | 
| 28 |  |  | *             "This product includes Radiance software | 
| 29 |  |  | *                 (http://radsite.lbl.gov/) | 
| 30 |  |  | *                 developed by the Lawrence Berkeley National Laboratory | 
| 31 |  |  | *               (http://www.lbl.gov/)." | 
| 32 |  |  | *       Alternately, this acknowledgment may appear in the software itself, | 
| 33 |  |  | *       if and wherever such third-party acknowledgments normally appear. | 
| 34 |  |  | * | 
| 35 |  |  | * 4. The names "Radiance," "Lawrence Berkeley National Laboratory" | 
| 36 |  |  | *       and "The Regents of the University of California" must | 
| 37 |  |  | *       not be used to endorse or promote products derived from this | 
| 38 |  |  | *       software without prior written permission. For written | 
| 39 |  |  | *       permission, please contact [email protected]. | 
| 40 |  |  | * | 
| 41 |  |  | * 5. Products derived from this software may not be called "Radiance", | 
| 42 |  |  | *       nor may "Radiance" appear in their name, without prior written | 
| 43 |  |  | *       permission of Lawrence Berkeley National Laboratory. | 
| 44 |  |  | * | 
| 45 |  |  | * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESSED OR IMPLIED | 
| 46 |  |  | * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES | 
| 47 |  |  | * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE | 
| 48 |  |  | * DISCLAIMED.   IN NO EVENT SHALL Lawrence Berkeley National Laboratory OR | 
| 49 |  |  | * ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, | 
| 50 |  |  | * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT | 
| 51 |  |  | * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF | 
| 52 |  |  | * USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND | 
| 53 |  |  | * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, | 
| 54 |  |  | * OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT | 
| 55 |  |  | * OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF | 
| 56 |  |  | * SUCH DAMAGE. | 
| 57 |  |  | * ==================================================================== | 
| 58 |  |  | * | 
| 59 |  |  | * This software consists of voluntary contributions made by many | 
| 60 |  |  | * individuals on behalf of Lawrence Berkeley National Laboratory.   For more | 
| 61 |  |  | * information on Lawrence Berkeley National Laboratory, please see | 
| 62 |  |  | * <http://www.lbl.gov/>. | 
| 63 | greg | 1.1 | */ | 
| 64 |  |  |  | 
| 65 |  |  | #include  "standard.h" | 
| 66 |  |  |  | 
| 67 |  |  | #include  "octree.h" | 
| 68 |  |  |  | 
| 69 |  |  | #include  "object.h" | 
| 70 |  |  |  | 
| 71 |  |  | #include  "otypes.h" | 
| 72 |  |  |  | 
| 73 | greg | 2.4 | static double  ogetflt(); | 
| 74 |  |  | static long  ogetint(); | 
| 75 |  |  | static char  *ogetstr(); | 
| 76 | greg | 2.12 | static int  nonsurfinset(); | 
| 77 | greg | 2.10 | static int  getobj(), octerror(), skiptree(); | 
| 78 | greg | 1.7 | static OCTREE  getfullnode(), gettree(); | 
| 79 | greg | 1.1 |  | 
| 80 | greg | 2.13 | static char  *infn;                     /* input file specification */ | 
| 81 | greg | 1.1 | static FILE  *infp;                     /* input file stream */ | 
| 82 | greg | 1.14 | static int  objsize;                    /* size of stored OBJECT's */ | 
| 83 | greg | 1.1 | static OBJECT  objorig;                 /* zeroeth object */ | 
| 84 | greg | 2.12 | static OBJECT  fnobjects;               /* number of objects in this file */ | 
| 85 | greg | 1.1 | static short  otypmap[NUMOTYPE+8];      /* object type map */ | 
| 86 |  |  |  | 
| 87 |  |  |  | 
| 88 |  |  | int | 
| 89 | greg | 2.13 | readoct(inpspec, load, scene, ofn)      /* read in octree file or stream */ | 
| 90 |  |  | char  *inpspec; | 
| 91 | greg | 1.1 | int  load; | 
| 92 |  |  | CUBE  *scene; | 
| 93 |  |  | char  *ofn[]; | 
| 94 |  |  | { | 
| 95 | greg | 1.11 | char  sbuf[512]; | 
| 96 | greg | 1.1 | int  nf; | 
| 97 |  |  | register int  i; | 
| 98 | greg | 1.15 | long  m; | 
| 99 | greg | 1.1 |  | 
| 100 | greg | 2.13 | if (inpspec == NULL) { | 
| 101 | greg | 1.1 | infn = "standard input"; | 
| 102 |  |  | infp = stdin; | 
| 103 | greg | 2.13 | } else if (inpspec[0] == '!') { | 
| 104 |  |  | infn = inpspec; | 
| 105 |  |  | if ((infp = popen(inpspec+1, "r")) == NULL) { | 
| 106 |  |  | sprintf(errmsg, "cannot execute \"%s\"", inpspec); | 
| 107 |  |  | error(SYSTEM, errmsg); | 
| 108 |  |  | } | 
| 109 | greg | 1.1 | } else { | 
| 110 | greg | 2.13 | infn = inpspec; | 
| 111 |  |  | if ((infp = fopen(inpspec, "r")) == NULL) { | 
| 112 | greg | 1.1 | sprintf(errmsg, "cannot open octree file \"%s\"", | 
| 113 | greg | 2.13 | inpspec); | 
| 114 | greg | 1.1 | error(SYSTEM, errmsg); | 
| 115 |  |  | } | 
| 116 |  |  | } | 
| 117 | greg | 2.7 | #ifdef MSDOS | 
| 118 |  |  | setmode(fileno(infp), O_BINARY); | 
| 119 |  |  | #endif | 
| 120 | greg | 1.1 | /* get header */ | 
| 121 | greg | 2.11 | if (checkheader(infp, OCTFMT, load&IO_INFO ? stdout : (FILE *)NULL) < 0) | 
| 122 | greg | 1.12 | octerror(USER, "not an octree"); | 
| 123 | greg | 1.1 | /* check format */ | 
| 124 | greg | 2.4 | if ((objsize = ogetint(2)-OCTMAGIC) <= 0 || | 
| 125 | greg | 1.14 | objsize > MAXOBJSIZ || objsize > sizeof(long)) | 
| 126 | greg | 1.12 | octerror(USER, "incompatible octree format"); | 
| 127 | greg | 1.1 | /* get boundaries */ | 
| 128 |  |  | if (load & IO_BOUNDS) { | 
| 129 |  |  | for (i = 0; i < 3; i++) | 
| 130 | greg | 2.4 | scene->cuorg[i] = atof(ogetstr(sbuf)); | 
| 131 |  |  | scene->cusize = atof(ogetstr(sbuf)); | 
| 132 | greg | 1.1 | } else { | 
| 133 |  |  | for (i = 0; i < 4; i++) | 
| 134 | greg | 2.4 | ogetstr(sbuf); | 
| 135 | greg | 1.1 | } | 
| 136 |  |  | objorig = nobjects;             /* set object offset */ | 
| 137 |  |  | nf = 0;                         /* get object files */ | 
| 138 | greg | 2.4 | while (*ogetstr(sbuf)) { | 
| 139 | greg | 1.1 | if (load & IO_SCENE) | 
| 140 |  |  | readobj(sbuf); | 
| 141 |  |  | if (load & IO_FILES) | 
| 142 |  |  | ofn[nf] = savqstr(sbuf); | 
| 143 |  |  | nf++; | 
| 144 |  |  | } | 
| 145 |  |  | if (load & IO_FILES) | 
| 146 |  |  | ofn[nf] = NULL; | 
| 147 | greg | 1.8 | /* get number of objects */ | 
| 148 | greg | 2.4 | fnobjects = m = ogetint(objsize); | 
| 149 | greg | 1.15 | if (fnobjects != m) | 
| 150 |  |  | octerror(USER, "too many objects"); | 
| 151 | greg | 1.1 |  | 
| 152 | greg | 2.13 | if (load & IO_TREE)             /* get the octree */ | 
| 153 | greg | 1.1 | scene->cutree = gettree(); | 
| 154 | greg | 2.5 | else if (load & IO_SCENE && nf == 0) | 
| 155 |  |  | skiptree(); | 
| 156 |  |  |  | 
| 157 |  |  | if (load & IO_SCENE)            /* get the scene */ | 
| 158 |  |  | if (nf == 0) { | 
| 159 |  |  | for (i = 0; *ogetstr(sbuf); i++) | 
| 160 |  |  | if ((otypmap[i] = otype(sbuf)) < 0) { | 
| 161 |  |  | sprintf(errmsg, "unknown type \"%s\"", sbuf); | 
| 162 |  |  | octerror(WARNING, errmsg); | 
| 163 |  |  | } | 
| 164 |  |  | while (getobj() != OVOID) | 
| 165 |  |  | ; | 
| 166 |  |  | } else {                    /* consistency checks */ | 
| 167 |  |  | /* check object count */ | 
| 168 |  |  | if (nobjects != objorig+fnobjects) | 
| 169 |  |  | octerror(USER, "bad object count; octree stale?"); | 
| 170 |  |  | /* check for non-surfaces */ | 
| 171 | greg | 2.12 | if (dosets(nonsurfinset)) | 
| 172 | greg | 2.5 | octerror(USER, "modifier in tree; octree stale?"); | 
| 173 |  |  | } | 
| 174 | greg | 2.13 | /* close the input */ | 
| 175 |  |  | if (infn[0] == '!') | 
| 176 |  |  | pclose(infp); | 
| 177 |  |  | else | 
| 178 |  |  | fclose(infp); | 
| 179 | greg | 1.1 | return(nf); | 
| 180 |  |  | } | 
| 181 |  |  |  | 
| 182 |  |  |  | 
| 183 |  |  | static char * | 
| 184 | greg | 2.4 | ogetstr(s)                      /* get null-terminated string */ | 
| 185 | greg | 1.1 | char  *s; | 
| 186 |  |  | { | 
| 187 | greg | 2.4 | extern char  *getstr(); | 
| 188 | greg | 1.1 |  | 
| 189 | greg | 2.4 | if (getstr(s, infp) == NULL) | 
| 190 |  |  | octerror(USER, "truncated octree"); | 
| 191 |  |  | return(s); | 
| 192 | greg | 1.1 | } | 
| 193 |  |  |  | 
| 194 |  |  |  | 
| 195 |  |  | static OCTREE | 
| 196 |  |  | getfullnode()                   /* get a set, return fullnode */ | 
| 197 |  |  | { | 
| 198 | greg | 2.7 | OBJECT  set[MAXSET+1]; | 
| 199 | greg | 1.1 | register int  i; | 
| 200 | greg | 1.14 | register long  m; | 
| 201 | greg | 1.1 |  | 
| 202 | greg | 2.4 | if ((set[0] = ogetint(objsize)) > MAXSET) | 
| 203 | greg | 1.1 | octerror(USER, "bad set in getfullnode"); | 
| 204 | greg | 1.14 | for (i = 1; i <= set[0]; i++) { | 
| 205 | greg | 2.4 | m = ogetint(objsize) + objorig; | 
| 206 | greg | 1.15 | if ((set[i] = m) != m) | 
| 207 | greg | 1.14 | octerror(USER, "too many objects"); | 
| 208 |  |  | } | 
| 209 | greg | 1.1 | return(fullnode(set)); | 
| 210 |  |  | } | 
| 211 |  |  |  | 
| 212 |  |  |  | 
| 213 |  |  | static long | 
| 214 | greg | 2.4 | ogetint(siz)                    /* get a siz-byte integer */ | 
| 215 |  |  | int  siz; | 
| 216 | greg | 1.1 | { | 
| 217 | greg | 2.4 | extern long  getint(); | 
| 218 | greg | 1.3 | register long  r; | 
| 219 | greg | 1.1 |  | 
| 220 | greg | 2.4 | r = getint(siz, infp); | 
| 221 |  |  | if (feof(infp)) | 
| 222 |  |  | octerror(USER, "truncated octree"); | 
| 223 | greg | 1.1 | return(r); | 
| 224 |  |  | } | 
| 225 |  |  |  | 
| 226 |  |  |  | 
| 227 |  |  | static double | 
| 228 | greg | 2.4 | ogetflt()                       /* get a floating point number */ | 
| 229 | greg | 1.1 | { | 
| 230 | greg | 2.4 | extern double  getflt(); | 
| 231 | greg | 2.7 | double  r; | 
| 232 | greg | 1.1 |  | 
| 233 | greg | 2.4 | r = getflt(infp); | 
| 234 |  |  | if (feof(infp)) | 
| 235 |  |  | octerror(USER, "truncated octree"); | 
| 236 |  |  | return(r); | 
| 237 | greg | 1.1 | } | 
| 238 |  |  |  | 
| 239 |  |  |  | 
| 240 |  |  | static OCTREE | 
| 241 |  |  | gettree()                       /* get a pre-ordered octree */ | 
| 242 |  |  | { | 
| 243 | greg | 2.7 | register OCTREE  ot; | 
| 244 | greg | 1.1 | register int  i; | 
| 245 |  |  |  | 
| 246 |  |  | switch (getc(infp)) { | 
| 247 |  |  | case OT_EMPTY: | 
| 248 |  |  | return(EMPTY); | 
| 249 |  |  | case OT_FULL: | 
| 250 |  |  | return(getfullnode()); | 
| 251 |  |  | case OT_TREE: | 
| 252 |  |  | if ((ot = octalloc()) == EMPTY) | 
| 253 |  |  | octerror(SYSTEM, "out of tree space in gettree"); | 
| 254 |  |  | for (i = 0; i < 8; i++) | 
| 255 |  |  | octkid(ot, i) = gettree(); | 
| 256 |  |  | return(ot); | 
| 257 | greg | 2.5 | case EOF: | 
| 258 |  |  | octerror(USER, "truncated octree"); | 
| 259 |  |  | default: | 
| 260 |  |  | octerror(USER, "damaged octree"); | 
| 261 |  |  | } | 
| 262 | greg | 2.12 | } | 
| 263 |  |  |  | 
| 264 |  |  |  | 
| 265 |  |  | static int | 
| 266 |  |  | nonsurfinset(os)                        /* check set for modifier */ | 
| 267 |  |  | register OBJECT  *os; | 
| 268 |  |  | { | 
| 269 |  |  | register OBJECT  s; | 
| 270 |  |  | register int  i; | 
| 271 |  |  |  | 
| 272 |  |  | for (i = *os; i-- > 0; ) | 
| 273 |  |  | if ((s = *++os) >= objorig && s < objorig+fnobjects && | 
| 274 |  |  | ismodifier(objptr(s)->otype)) | 
| 275 |  |  | return(1); | 
| 276 |  |  | return(0); | 
| 277 | greg | 2.5 | } | 
| 278 |  |  |  | 
| 279 |  |  |  | 
| 280 |  |  | static | 
| 281 |  |  | skiptree()                              /* skip octree on input */ | 
| 282 |  |  | { | 
| 283 |  |  | register int  i; | 
| 284 |  |  |  | 
| 285 |  |  | switch (getc(infp)) { | 
| 286 |  |  | case OT_EMPTY: | 
| 287 |  |  | return; | 
| 288 |  |  | case OT_FULL: | 
| 289 | greg | 2.6 | for (i = ogetint(objsize)*objsize; i-- > 0; ) | 
| 290 |  |  | if (getc(infp) == EOF) | 
| 291 |  |  | octerror(USER, "truncated octree"); | 
| 292 | greg | 2.5 | return; | 
| 293 |  |  | case OT_TREE: | 
| 294 |  |  | for (i = 0; i < 8; i++) | 
| 295 |  |  | skiptree(); | 
| 296 |  |  | return; | 
| 297 | greg | 1.1 | case EOF: | 
| 298 |  |  | octerror(USER, "truncated octree"); | 
| 299 |  |  | default: | 
| 300 |  |  | octerror(USER, "damaged octree"); | 
| 301 |  |  | } | 
| 302 |  |  | } | 
| 303 |  |  |  | 
| 304 |  |  |  | 
| 305 |  |  | static | 
| 306 |  |  | getobj()                                /* get next object */ | 
| 307 |  |  | { | 
| 308 |  |  | char  sbuf[MAXSTR]; | 
| 309 |  |  | int  obj; | 
| 310 |  |  | register int  i; | 
| 311 | greg | 1.14 | register long  m; | 
| 312 | greg | 2.7 | register OBJREC  *objp; | 
| 313 | greg | 1.1 |  | 
| 314 | greg | 2.4 | i = ogetint(1); | 
| 315 | greg | 1.3 | if (i == -1) | 
| 316 | greg | 1.1 | return(OVOID);          /* terminator */ | 
| 317 |  |  | if ((obj = newobject()) == OVOID) | 
| 318 |  |  | error(SYSTEM, "out of object space"); | 
| 319 |  |  | objp = objptr(obj); | 
| 320 |  |  | if ((objp->otype = otypmap[i]) < 0) | 
| 321 |  |  | octerror(USER, "reference to unknown type"); | 
| 322 | greg | 2.4 | if ((m = ogetint(objsize)) != OVOID) { | 
| 323 | greg | 1.14 | m += objorig; | 
| 324 |  |  | if ((OBJECT)m != m) | 
| 325 |  |  | octerror(USER, "too many objects"); | 
| 326 |  |  | } | 
| 327 |  |  | objp->omod = m; | 
| 328 | greg | 2.4 | objp->oname = savqstr(ogetstr(sbuf)); | 
| 329 |  |  | if (objp->oargs.nsargs = ogetint(2)) { | 
| 330 | greg | 2.13 | objp->oargs.sarg = (char **)malloc | 
| 331 | greg | 1.1 | (objp->oargs.nsargs*sizeof(char *)); | 
| 332 |  |  | if (objp->oargs.sarg == NULL) | 
| 333 |  |  | goto memerr; | 
| 334 |  |  | for (i = 0; i < objp->oargs.nsargs; i++) | 
| 335 | greg | 2.4 | objp->oargs.sarg[i] = savestr(ogetstr(sbuf)); | 
| 336 | greg | 1.1 | } else | 
| 337 |  |  | objp->oargs.sarg = NULL; | 
| 338 | greg | 2.7 | #ifdef  IARGS | 
| 339 | greg | 2.4 | if (objp->oargs.niargs = ogetint(2)) { | 
| 340 | greg | 2.13 | objp->oargs.iarg = (long *)malloc | 
| 341 | greg | 1.1 | (objp->oargs.niargs*sizeof(long)); | 
| 342 |  |  | if (objp->oargs.iarg == NULL) | 
| 343 |  |  | goto memerr; | 
| 344 |  |  | for (i = 0; i < objp->oargs.niargs; i++) | 
| 345 | greg | 2.4 | objp->oargs.iarg[i] = ogetint(4); | 
| 346 | greg | 1.1 | } else | 
| 347 |  |  | objp->oargs.iarg = NULL; | 
| 348 |  |  | #endif | 
| 349 | greg | 2.4 | if (objp->oargs.nfargs = ogetint(2)) { | 
| 350 | greg | 2.13 | objp->oargs.farg = (FLOAT *)malloc | 
| 351 | greg | 1.15 | (objp->oargs.nfargs*sizeof(FLOAT)); | 
| 352 | greg | 1.1 | if (objp->oargs.farg == NULL) | 
| 353 |  |  | goto memerr; | 
| 354 |  |  | for (i = 0; i < objp->oargs.nfargs; i++) | 
| 355 | greg | 2.4 | objp->oargs.farg[i] = ogetflt(); | 
| 356 | greg | 1.1 | } else | 
| 357 |  |  | objp->oargs.farg = NULL; | 
| 358 |  |  | /* initialize */ | 
| 359 |  |  | objp->os = NULL; | 
| 360 |  |  | /* insert */ | 
| 361 |  |  | insertobject(obj); | 
| 362 |  |  | return(obj); | 
| 363 |  |  | memerr: | 
| 364 |  |  | error(SYSTEM, "out of memory in getobj"); | 
| 365 |  |  | } | 
| 366 |  |  |  | 
| 367 |  |  |  | 
| 368 |  |  | static | 
| 369 |  |  | octerror(etyp, msg)                     /* octree error */ | 
| 370 |  |  | int  etyp; | 
| 371 |  |  | char  *msg; | 
| 372 |  |  | { | 
| 373 |  |  | char  msgbuf[128]; | 
| 374 |  |  |  | 
| 375 |  |  | sprintf(msgbuf, "(%s): %s", infn, msg); | 
| 376 |  |  | error(etyp, msgbuf); | 
| 377 |  |  | } |