src/common/readoct.c

/* Copyright (c) 1992 Regents of the University of California */

#ifndef lint
static char SCCSid[] = "$SunId$ LBL";
#endif

/*
 *  readoct.c - routines to read octree information.
 *
 *     7/30/85
 */

#include  "standard.h"

#include  "octree.h"

#include  "object.h"

#include  "otypes.h"

static double  ogetflt();
static long  ogetint();
static char  *ogetstr();
static int  nonsurfinset();
static int  getobj(), octerror(), skiptree();
static OCTREE  getfullnode(), gettree();

static char  *infn;                     /* input file name */
static FILE  *infp;                     /* input file stream */
static int  objsize;                    /* size of stored OBJECT's */
static OBJECT  objorig;                 /* zeroeth object */
static OBJECT  fnobjects;               /* number of objects in this file */
static short  otypmap[NUMOTYPE+8];      /* object type map */


int
readoct(fname, load, scene, ofn)        /* read in octree from file */
char  *fname;
int  load;
CUBE  *scene;
char  *ofn[];
{
        char  sbuf[512];
        int  nf;
        register int  i;
        long  m;
        
        if (fname == NULL) {
                infn = "standard input";
                infp = stdin;
        } else {
                infn = fname;
                if ((infp = fopen(fname, "r")) == NULL) {
                        sprintf(errmsg, "cannot open octree file \"%s\"",
                                        fname);
                        error(SYSTEM, errmsg);
                }
        }
#ifdef MSDOS
        setmode(fileno(infp), O_BINARY);
#endif
                                        /* get header */
        if (checkheader(infp, OCTFMT, load&IO_INFO ? stdout : (FILE *)NULL) < 0)
                octerror(USER, "not an octree");
                                        /* check format */
        if ((objsize = ogetint(2)-OCTMAGIC) <= 0 ||
                        objsize > MAXOBJSIZ || objsize > sizeof(long))
                octerror(USER, "incompatible octree format");
                                        /* get boundaries */
        if (load & IO_BOUNDS) {
                for (i = 0; i < 3; i++)
                        scene->cuorg[i] = atof(ogetstr(sbuf));
                scene->cusize = atof(ogetstr(sbuf));
        } else {
                for (i = 0; i < 4; i++)
                        ogetstr(sbuf);
        }
        objorig = nobjects;             /* set object offset */
        nf = 0;                         /* get object files */
        while (*ogetstr(sbuf)) {
                if (load & IO_SCENE)
                        readobj(sbuf);
                if (load & IO_FILES)
                        ofn[nf] = savqstr(sbuf);
                nf++;
        }
        if (load & IO_FILES)
                ofn[nf] = NULL;
                                        /* get number of objects */
        fnobjects = m = ogetint(objsize);
        if (fnobjects != m)
                octerror(USER, "too many objects");

        if (load & IO_TREE)                     /* get the octree */
                scene->cutree = gettree();
        else if (load & IO_SCENE && nf == 0)
                skiptree();
                
        if (load & IO_SCENE)            /* get the scene */
            if (nf == 0) {
                for (i = 0; *ogetstr(sbuf); i++)
                        if ((otypmap[i] = otype(sbuf)) < 0) {
                                sprintf(errmsg, "unknown type \"%s\"", sbuf);
                                octerror(WARNING, errmsg);
                        }
                while (getobj() != OVOID)
                        ;
            } else {                    /* consistency checks */
                                /* check object count */
                if (nobjects != objorig+fnobjects)
                        octerror(USER, "bad object count; octree stale?");
                                /* check for non-surfaces */
                if (dosets(nonsurfinset))
                        octerror(USER, "modifier in tree; octree stale?");
            }
        fclose(infp);
        return(nf);
}


static char *
ogetstr(s)                      /* get null-terminated string */
char  *s;
{
        extern char  *getstr();

        if (getstr(s, infp) == NULL)
                octerror(USER, "truncated octree");
        return(s);
}


static OCTREE
getfullnode()                   /* get a set, return fullnode */
{
        OBJECT  set[MAXSET+1];
        register int  i;
        register long  m;

        if ((set[0] = ogetint(objsize)) > MAXSET)
                octerror(USER, "bad set in getfullnode");
        for (i = 1; i <= set[0]; i++) {
                m = ogetint(objsize) + objorig;
                if ((set[i] = m) != m)
                        octerror(USER, "too many objects");
        }
        return(fullnode(set));
}       


static long
ogetint(siz)                    /* get a siz-byte integer */
int  siz;
{
        extern long  getint();
        register long  r;

        r = getint(siz, infp);
        if (feof(infp))
                octerror(USER, "truncated octree");
        return(r);
}


static double
ogetflt()                       /* get a floating point number */
{
        extern double  getflt();
        double  r;

        r = getflt(infp);
        if (feof(infp))
                octerror(USER, "truncated octree");
        return(r);
}
        

static OCTREE
gettree()                       /* get a pre-ordered octree */
{
        register OCTREE  ot;
        register int  i;
        
        switch (getc(infp)) {
        case OT_EMPTY:
                return(EMPTY);
        case OT_FULL:
                return(getfullnode());
        case OT_TREE:
                if ((ot = octalloc()) == EMPTY)
                        octerror(SYSTEM, "out of tree space in gettree");
                for (i = 0; i < 8; i++)
                        octkid(ot, i) = gettree();
                return(ot);
        case EOF:
                octerror(USER, "truncated octree");
        default:
                octerror(USER, "damaged octree");
        }
}


static int
nonsurfinset(os)                        /* check set for modifier */
register OBJECT  *os;
{
        register OBJECT  s;
        register int  i;

        for (i = *os; i-- > 0; )
                if ((s = *++os) >= objorig && s < objorig+fnobjects &&
                                ismodifier(objptr(s)->otype))
                        return(1);
        return(0);
}


static
skiptree()                              /* skip octree on input */
{
        register int  i;
        
        switch (getc(infp)) {
        case OT_EMPTY:
                return;
        case OT_FULL:
                for (i = ogetint(objsize)*objsize; i-- > 0; )
                        if (getc(infp) == EOF)
                                octerror(USER, "truncated octree");
                return;
        case OT_TREE:
                for (i = 0; i < 8; i++)
                        skiptree();
                return;
        case EOF:
                octerror(USER, "truncated octree");
        default:
                octerror(USER, "damaged octree");
        }
}


static
getobj()                                /* get next object */
{
        char  sbuf[MAXSTR];
        int  obj;
        register int  i;
        register long  m;
        register OBJREC  *objp;

        i = ogetint(1);
        if (i == -1)
                return(OVOID);          /* terminator */
        if ((obj = newobject()) == OVOID)
                error(SYSTEM, "out of object space");
        objp = objptr(obj);
        if ((objp->otype = otypmap[i]) < 0)
                octerror(USER, "reference to unknown type");
        if ((m = ogetint(objsize)) != OVOID) {
                m += objorig;
                if ((OBJECT)m != m)
                        octerror(USER, "too many objects");
        }
        objp->omod = m;
        objp->oname = savqstr(ogetstr(sbuf));
        if (objp->oargs.nsargs = ogetint(2)) {
                objp->oargs.sarg = (char **)bmalloc
                                (objp->oargs.nsargs*sizeof(char *));
                if (objp->oargs.sarg == NULL)
                        goto memerr;
                for (i = 0; i < objp->oargs.nsargs; i++)
                        objp->oargs.sarg[i] = savestr(ogetstr(sbuf));
        } else
                objp->oargs.sarg = NULL;
#ifdef  IARGS
        if (objp->oargs.niargs = ogetint(2)) {
                objp->oargs.iarg = (long *)bmalloc
                                (objp->oargs.niargs*sizeof(long));
                if (objp->oargs.iarg == NULL)
                        goto memerr;
                for (i = 0; i < objp->oargs.niargs; i++)
                        objp->oargs.iarg[i] = ogetint(4);
        } else
                objp->oargs.iarg = NULL;
#endif
        if (objp->oargs.nfargs = ogetint(2)) {
                objp->oargs.farg = (FLOAT *)bmalloc
                                (objp->oargs.nfargs*sizeof(FLOAT));
                if (objp->oargs.farg == NULL)
                        goto memerr;
                for (i = 0; i < objp->oargs.nfargs; i++)
                        objp->oargs.farg[i] = ogetflt();
        } else
                objp->oargs.farg = NULL;
                                                /* initialize */
        objp->os = NULL;
                                                /* insert */
        insertobject(obj);
        return(obj);
memerr:
        error(SYSTEM, "out of memory in getobj");
}


static
octerror(etyp, msg)                     /* octree error */
int  etyp;
char  *msg;
{
        char  msgbuf[128];

        sprintf(msgbuf, "(%s): %s", infn, msg);
        error(etyp, msgbuf);
}
Revision:	2.12
Committed:	Fri Mar 7 15:45:26 1997 UTC (27 years, 1 month ago) by greg
Content type:	text/plain
Branch:	MAIN
Changes since 2.11:	+18 -2 lines
Log Message:	added new set routines
#	User	Rev	Content
1	greg	2.5	/* Copyright (c) 1992 Regents of the University of California */
2	greg	1.1
3			#ifndef lint
4			static char SCCSid[] = "$SunId$ LBL";
5			#endif
6
7			/*
8			* readoct.c - routines to read octree information.
9			*
10			* 7/30/85
11			*/
12
13			#include "standard.h"
14
15			#include "octree.h"
16
17			#include "object.h"
18
19			#include "otypes.h"
20
21	greg	2.4	static double ogetflt();
22			static long ogetint();
23			static char *ogetstr();
24	greg	2.12	static int nonsurfinset();
25	greg	2.10	static int getobj(), octerror(), skiptree();
26	greg	1.7	static OCTREE getfullnode(), gettree();
27	greg	1.1
28			static char infn; / input file name */
29			static FILE infp; / input file stream */
30	greg	1.14	static int objsize; /* size of stored OBJECT's */
31	greg	1.1	static OBJECT objorig; /* zeroeth object */
32	greg	2.12	static OBJECT fnobjects; /* number of objects in this file */
33	greg	1.1	static short otypmap[NUMOTYPE+8]; /* object type map */
34
35
36			int
37			readoct(fname, load, scene, ofn) /* read in octree from file */
38			char *fname;
39			int load;
40			CUBE *scene;
41			char *ofn[];
42			{
43	greg	1.11	char sbuf[512];
44	greg	1.1	int nf;
45			register int i;
46	greg	1.15	long m;
47	greg	1.1
48			if (fname == NULL) {
49			infn = "standard input";
50			infp = stdin;
51			} else {
52			infn = fname;
53			if ((infp = fopen(fname, "r")) == NULL) {
54			sprintf(errmsg, "cannot open octree file \"%s\"",
55			fname);
56			error(SYSTEM, errmsg);
57			}
58			}
59	greg	2.7	#ifdef MSDOS
60			setmode(fileno(infp), O_BINARY);
61			#endif
62	greg	1.1	/* get header */
63	greg	2.11	if (checkheader(infp, OCTFMT, load&IO_INFO ? stdout : (FILE *)NULL) < 0)
64	greg	1.12	octerror(USER, "not an octree");
65	greg	1.1	/* check format */
66	greg	2.4	if ((objsize = ogetint(2)-OCTMAGIC) <= 0 \|\|
67	greg	1.14	objsize > MAXOBJSIZ \|\| objsize > sizeof(long))
68	greg	1.12	octerror(USER, "incompatible octree format");
69	greg	1.1	/* get boundaries */
70			if (load & IO_BOUNDS) {
71			for (i = 0; i < 3; i++)
72	greg	2.4	scene->cuorg[i] = atof(ogetstr(sbuf));
73			scene->cusize = atof(ogetstr(sbuf));
74	greg	1.1	} else {
75			for (i = 0; i < 4; i++)
76	greg	2.4	ogetstr(sbuf);
77	greg	1.1	}
78			objorig = nobjects; /* set object offset */
79			nf = 0; /* get object files */
80	greg	2.4	while (*ogetstr(sbuf)) {
81	greg	1.1	if (load & IO_SCENE)
82			readobj(sbuf);
83			if (load & IO_FILES)
84			ofn[nf] = savqstr(sbuf);
85			nf++;
86			}
87			if (load & IO_FILES)
88			ofn[nf] = NULL;
89	greg	1.8	/* get number of objects */
90	greg	2.4	fnobjects = m = ogetint(objsize);
91	greg	1.15	if (fnobjects != m)
92			octerror(USER, "too many objects");
93	greg	1.1
94	greg	2.5	if (load & IO_TREE) /* get the octree */
95	greg	1.1	scene->cutree = gettree();
96	greg	2.5	else if (load & IO_SCENE && nf == 0)
97			skiptree();
98
99			if (load & IO_SCENE) /* get the scene */
100			if (nf == 0) {
101			for (i = 0; *ogetstr(sbuf); i++)
102			if ((otypmap[i] = otype(sbuf)) < 0) {
103			sprintf(errmsg, "unknown type \"%s\"", sbuf);
104			octerror(WARNING, errmsg);
105			}
106			while (getobj() != OVOID)
107			;
108			} else { /* consistency checks */
109			/* check object count */
110			if (nobjects != objorig+fnobjects)
111			octerror(USER, "bad object count; octree stale?");
112			/* check for non-surfaces */
113	greg	2.12	if (dosets(nonsurfinset))
114	greg	2.5	octerror(USER, "modifier in tree; octree stale?");
115			}
116	greg	1.1	fclose(infp);
117			return(nf);
118			}
119
120
121			static char *
122	greg	2.4	ogetstr(s) /* get null-terminated string */
123	greg	1.1	char *s;
124			{
125	greg	2.4	extern char *getstr();
126	greg	1.1
127	greg	2.4	if (getstr(s, infp) == NULL)
128			octerror(USER, "truncated octree");
129			return(s);
130	greg	1.1	}
131
132
133			static OCTREE
134			getfullnode() /* get a set, return fullnode */
135			{
136	greg	2.7	OBJECT set[MAXSET+1];
137	greg	1.1	register int i;
138	greg	1.14	register long m;
139	greg	1.1
140	greg	2.4	if ((set[0] = ogetint(objsize)) > MAXSET)
141	greg	1.1	octerror(USER, "bad set in getfullnode");
142	greg	1.14	for (i = 1; i <= set[0]; i++) {
143	greg	2.4	m = ogetint(objsize) + objorig;
144	greg	1.15	if ((set[i] = m) != m)
145	greg	1.14	octerror(USER, "too many objects");
146			}
147	greg	1.1	return(fullnode(set));
148			}
149
150
151			static long
152	greg	2.4	ogetint(siz) /* get a siz-byte integer */
153			int siz;
154	greg	1.1	{
155	greg	2.4	extern long getint();
156	greg	1.3	register long r;
157	greg	1.1
158	greg	2.4	r = getint(siz, infp);
159			if (feof(infp))
160			octerror(USER, "truncated octree");
161	greg	1.1	return(r);
162			}
163
164
165			static double
166	greg	2.4	ogetflt() /* get a floating point number */
167	greg	1.1	{
168	greg	2.4	extern double getflt();
169	greg	2.7	double r;
170	greg	1.1
171	greg	2.4	r = getflt(infp);
172			if (feof(infp))
173			octerror(USER, "truncated octree");
174			return(r);
175	greg	1.1	}
176
177
178			static OCTREE
179			gettree() /* get a pre-ordered octree */
180			{
181	greg	2.7	register OCTREE ot;
182	greg	1.1	register int i;
183
184			switch (getc(infp)) {
185			case OT_EMPTY:
186			return(EMPTY);
187			case OT_FULL:
188			return(getfullnode());
189			case OT_TREE:
190			if ((ot = octalloc()) == EMPTY)
191			octerror(SYSTEM, "out of tree space in gettree");
192			for (i = 0; i < 8; i++)
193			octkid(ot, i) = gettree();
194			return(ot);
195	greg	2.5	case EOF:
196			octerror(USER, "truncated octree");
197			default:
198			octerror(USER, "damaged octree");
199			}
200	greg	2.12	}
201
202
203			static int
204			nonsurfinset(os) /* check set for modifier */
205			register OBJECT *os;
206			{
207			register OBJECT s;
208			register int i;
209
210			for (i = *os; i-- > 0; )
211			if ((s = *++os) >= objorig && s < objorig+fnobjects &&
212			ismodifier(objptr(s)->otype))
213			return(1);
214			return(0);
215	greg	2.5	}
216
217
218			static
219			skiptree() /* skip octree on input */
220			{
221			register int i;
222
223			switch (getc(infp)) {
224			case OT_EMPTY:
225			return;
226			case OT_FULL:
227	greg	2.6	for (i = ogetint(objsize)*objsize; i-- > 0; )
228			if (getc(infp) == EOF)
229			octerror(USER, "truncated octree");
230	greg	2.5	return;
231			case OT_TREE:
232			for (i = 0; i < 8; i++)
233			skiptree();
234			return;
235	greg	1.1	case EOF:
236			octerror(USER, "truncated octree");
237			default:
238			octerror(USER, "damaged octree");
239			}
240			}
241
242
243			static
244			getobj() /* get next object */
245			{
246			char sbuf[MAXSTR];
247			int obj;
248			register int i;
249	greg	1.14	register long m;
250	greg	2.7	register OBJREC *objp;
251	greg	1.1
252	greg	2.4	i = ogetint(1);
253	greg	1.3	if (i == -1)
254	greg	1.1	return(OVOID); /* terminator */
255			if ((obj = newobject()) == OVOID)
256			error(SYSTEM, "out of object space");
257			objp = objptr(obj);
258			if ((objp->otype = otypmap[i]) < 0)
259			octerror(USER, "reference to unknown type");
260	greg	2.4	if ((m = ogetint(objsize)) != OVOID) {
261	greg	1.14	m += objorig;
262			if ((OBJECT)m != m)
263			octerror(USER, "too many objects");
264			}
265			objp->omod = m;
266	greg	2.4	objp->oname = savqstr(ogetstr(sbuf));
267			if (objp->oargs.nsargs = ogetint(2)) {
268	greg	1.1	objp->oargs.sarg = (char **)bmalloc
269			(objp->oargs.nsargssizeof(char ));
270			if (objp->oargs.sarg == NULL)
271			goto memerr;
272			for (i = 0; i < objp->oargs.nsargs; i++)
273	greg	2.4	objp->oargs.sarg[i] = savestr(ogetstr(sbuf));
274	greg	1.1	} else
275			objp->oargs.sarg = NULL;
276	greg	2.7	#ifdef IARGS
277	greg	2.4	if (objp->oargs.niargs = ogetint(2)) {
278	greg	1.1	objp->oargs.iarg = (long *)bmalloc
279			(objp->oargs.niargs*sizeof(long));
280			if (objp->oargs.iarg == NULL)
281			goto memerr;
282			for (i = 0; i < objp->oargs.niargs; i++)
283	greg	2.4	objp->oargs.iarg[i] = ogetint(4);
284	greg	1.1	} else
285			objp->oargs.iarg = NULL;
286			#endif
287	greg	2.4	if (objp->oargs.nfargs = ogetint(2)) {
288	greg	1.15	objp->oargs.farg = (FLOAT *)bmalloc
289			(objp->oargs.nfargs*sizeof(FLOAT));
290	greg	1.1	if (objp->oargs.farg == NULL)
291			goto memerr;
292			for (i = 0; i < objp->oargs.nfargs; i++)
293	greg	2.4	objp->oargs.farg[i] = ogetflt();
294	greg	1.1	} else
295			objp->oargs.farg = NULL;
296			/* initialize */
297			objp->os = NULL;
298			/* insert */
299			insertobject(obj);
300			return(obj);
301			memerr:
302			error(SYSTEM, "out of memory in getobj");
303			}
304
305
306			static
307			octerror(etyp, msg) /* octree error */
308			int etyp;
309			char *msg;
310			{
311			char msgbuf[128];
312
313			sprintf(msgbuf, "(%s): %s", infn, msg);
314			error(etyp, msgbuf);
315			}