src/common/readoct.c

#ifndef lint
static const char       RCSid[] = "$Id: readoct.c,v 2.27 2004/12/15 22:02:04 greg Exp $";
#endif
/*
 *  readoct.c - routines to read octree information.
 */

#include "copyright.h"

#include  <stdio.h>
#include  <time.h>

#include  "platform.h"
#include  "rtprocess.h"
#include  "standard.h"
#include  "octree.h"
#include  "object.h"
#include  "otypes.h"
#include  "resolu.h"

#ifdef getc_unlocked            /* avoid horrendous overhead of flockfile */
#undef getc
#define getc    getc_unlocked
#endif

static double  ogetflt(void);
static long  ogetint(int);
static char  *ogetstr(char *);
static int  nonsurfintree(OCTREE ot);
static void  octerror(int  etyp, char  *msg);
static void  skiptree(void);
static OCTREE  getfullnode(void), gettree(void);

static char  *infn;                     /* input file specification */
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 */


int
readoct(                                /* read in octree file or stream */
        char  *inpspec,
        int  load,
        CUBE  *scene,
        char  *ofn[]
)
{
        char  sbuf[512];
        int  nf;
        int  i;
        long  m;
        
        if (inpspec == NULL) {
                infn = "standard input";
                infp = stdin;
        } else if (inpspec[0] == '!') {
                infn = inpspec;
                if ((infp = popen(inpspec+1, "r")) == NULL) {
                        sprintf(errmsg, "cannot execute \"%s\"", inpspec);
                        error(SYSTEM, errmsg);
                }
        } else {
                infn = inpspec;
                if ((infp = fopen(inpspec, "r")) == NULL) {
                        sprintf(errmsg, "cannot open octree file \"%s\"",
                                        inpspec);
                        error(SYSTEM, errmsg);
                }
        }
        SET_FILE_BINARY(infp);
                                        /* 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) {
                                        /* load binary scene data */
                readscene(infp, objsize);

            } else {                    /* consistency checks */
                                /* check object count */
                if (nobjects != objorig+fnobjects)
                        octerror(USER, "bad object count; octree stale?");
                                /* check for non-surfaces */
                if (nonsurfintree(scene->cutree))
                        octerror(USER, "modifier in tree; octree stale?");
            }
        }
                                /* close the input */
        if (infn[0] == '!')
                pclose(infp);
        else
                fclose(infp);
        return(nf);
}


static char *
ogetstr(char *s)                /* get null-terminated string */
{
        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(int siz)                /* get a siz-byte integer */
{
        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");
        }
        return EMPTY; /* pro forma return */
}


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
skiptree(void)                          /* 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 void
octerror(int etyp, char *msg)           /* octree error */
{
        char  msgbuf[128];

        sprintf(msgbuf, "(%s): %s", infn, msg);
        error(etyp, msgbuf);
}
Revision:	2.28
Committed:	Thu Dec 16 19:31:36 2004 UTC (20 years, 10 months ago) by greg
Content type:	text/plain
Branch:	MAIN
CVS Tags:	rad5R0, rad4R2, rad3R7P2, rad3R7P1, rad4R1, rad4R0, rad3R8, rad3R9, rad4R2P1, rad4R2P2
Changes since 2.27:	+2 -2 lines
Log Message:	Fixed problem with last bug-fix (serious)
#	User	Rev	Content
1	greg	1.1	#ifndef lint
2	greg	2.28	static const char RCSid[] = "$Id: readoct.c,v 2.27 2004/12/15 22:02:04 greg Exp $";
3	greg	1.1	#endif
4			/*
5			* readoct.c - routines to read octree information.
6	greg	2.13	*/
7
8	greg	2.14	#include "copyright.h"
9	greg	1.1
10	schorsch	2.17	#include <stdio.h>
11	greg	2.19	#include <time.h>
12	schorsch	2.17
13	greg	2.23	#include "platform.h"
14	schorsch	2.24	#include "rtprocess.h"
15	greg	1.1	#include "standard.h"
16			#include "octree.h"
17			#include "object.h"
18			#include "otypes.h"
19	schorsch	2.17	#include "resolu.h"
20	greg	1.1
21	greg	2.25	#ifdef getc_unlocked /* avoid horrendous overhead of flockfile */
22	greg	2.26	#undef getc
23	greg	2.25	#define getc getc_unlocked
24			#endif
25
26	schorsch	2.20	static double ogetflt(void);
27			static long ogetint(int);
28			static char ogetstr(char );
29	greg	2.27	static int nonsurfintree(OCTREE ot);
30	schorsch	2.17	static void octerror(int etyp, char *msg);
31			static void skiptree(void);
32	schorsch	2.20	static OCTREE getfullnode(void), gettree(void);
33	greg	1.1
34	greg	2.13	static char infn; / input file specification */
35	greg	1.1	static FILE infp; / input file stream */
36	greg	1.14	static int objsize; /* size of stored OBJECT's */
37	greg	1.1	static OBJECT objorig; /* zeroeth object */
38	greg	2.12	static OBJECT fnobjects; /* number of objects in this file */
39	greg	1.1
40
41			int
42	greg	2.27	readoct( /* read in octree file or stream */
43			char *inpspec,
44			int load,
45			CUBE *scene,
46			char *ofn[]
47			)
48	greg	1.1	{
49	greg	1.11	char sbuf[512];
50	greg	1.1	int nf;
51	greg	2.15	int i;
52	greg	1.15	long m;
53	greg	1.1
54	greg	2.13	if (inpspec == NULL) {
55	greg	1.1	infn = "standard input";
56			infp = stdin;
57	greg	2.13	} else if (inpspec[0] == '!') {
58			infn = inpspec;
59			if ((infp = popen(inpspec+1, "r")) == NULL) {
60			sprintf(errmsg, "cannot execute \"%s\"", inpspec);
61			error(SYSTEM, errmsg);
62			}
63	greg	1.1	} else {
64	greg	2.13	infn = inpspec;
65			if ((infp = fopen(inpspec, "r")) == NULL) {
66	greg	1.1	sprintf(errmsg, "cannot open octree file \"%s\"",
67	greg	2.13	inpspec);
68	greg	1.1	error(SYSTEM, errmsg);
69			}
70			}
71	schorsch	2.16	SET_FILE_BINARY(infp);
72	greg	1.1	/* get header */
73	greg	2.11	if (checkheader(infp, OCTFMT, load&IO_INFO ? stdout : (FILE *)NULL) < 0)
74	greg	1.12	octerror(USER, "not an octree");
75	greg	1.1	/* check format */
76	greg	2.4	if ((objsize = ogetint(2)-OCTMAGIC) <= 0 \|\|
77	greg	1.14	objsize > MAXOBJSIZ \|\| objsize > sizeof(long))
78	greg	1.12	octerror(USER, "incompatible octree format");
79	greg	1.1	/* get boundaries */
80			if (load & IO_BOUNDS) {
81			for (i = 0; i < 3; i++)
82	greg	2.4	scene->cuorg[i] = atof(ogetstr(sbuf));
83			scene->cusize = atof(ogetstr(sbuf));
84	greg	1.1	} else {
85			for (i = 0; i < 4; i++)
86	greg	2.4	ogetstr(sbuf);
87	greg	1.1	}
88			objorig = nobjects; /* set object offset */
89			nf = 0; /* get object files */
90	greg	2.4	while (*ogetstr(sbuf)) {
91	greg	1.1	if (load & IO_SCENE)
92			readobj(sbuf);
93			if (load & IO_FILES)
94			ofn[nf] = savqstr(sbuf);
95			nf++;
96			}
97			if (load & IO_FILES)
98			ofn[nf] = NULL;
99	greg	1.8	/* get number of objects */
100	greg	2.4	fnobjects = m = ogetint(objsize);
101	greg	1.15	if (fnobjects != m)
102			octerror(USER, "too many objects");
103	greg	1.1
104	greg	2.13	if (load & IO_TREE) /* get the octree */
105	greg	1.1	scene->cutree = gettree();
106	greg	2.5	else if (load & IO_SCENE && nf == 0)
107			skiptree();
108
109	schorsch	2.21	if (load & IO_SCENE) { /* get the scene */
110	greg	2.5	if (nf == 0) {
111	greg	2.15	/* load binary scene data */
112			readscene(infp, objsize);
113
114	greg	2.5	} else { /* consistency checks */
115			/* check object count */
116			if (nobjects != objorig+fnobjects)
117			octerror(USER, "bad object count; octree stale?");
118			/* check for non-surfaces */
119	greg	2.27	if (nonsurfintree(scene->cutree))
120	greg	2.5	octerror(USER, "modifier in tree; octree stale?");
121			}
122	schorsch	2.21	}
123	greg	2.13	/* close the input */
124			if (infn[0] == '!')
125			pclose(infp);
126			else
127			fclose(infp);
128	greg	1.1	return(nf);
129			}
130
131
132			static char *
133	greg	2.27	ogetstr(char s) / get null-terminated string */
134	greg	1.1	{
135	greg	2.4	extern char *getstr();
136	greg	1.1
137	greg	2.4	if (getstr(s, infp) == NULL)
138			octerror(USER, "truncated octree");
139			return(s);
140	greg	1.1	}
141
142
143			static OCTREE
144			getfullnode() /* get a set, return fullnode */
145			{
146	greg	2.7	OBJECT set[MAXSET+1];
147	greg	1.1	register int i;
148	greg	1.14	register long m;
149	greg	1.1
150	greg	2.4	if ((set[0] = ogetint(objsize)) > MAXSET)
151	greg	1.1	octerror(USER, "bad set in getfullnode");
152	greg	1.14	for (i = 1; i <= set[0]; i++) {
153	greg	2.4	m = ogetint(objsize) + objorig;
154	greg	1.15	if ((set[i] = m) != m)
155	greg	1.14	octerror(USER, "too many objects");
156			}
157	greg	1.1	return(fullnode(set));
158			}
159
160
161			static long
162	greg	2.27	ogetint(int siz) /* get a siz-byte integer */
163	greg	1.1	{
164	greg	2.4	extern long getint();
165	greg	1.3	register long r;
166	greg	1.1
167	greg	2.4	r = getint(siz, infp);
168			if (feof(infp))
169			octerror(USER, "truncated octree");
170	greg	1.1	return(r);
171			}
172
173
174			static double
175	greg	2.4	ogetflt() /* get a floating point number */
176	greg	1.1	{
177	greg	2.4	extern double getflt();
178	greg	2.7	double r;
179	greg	1.1
180	greg	2.4	r = getflt(infp);
181			if (feof(infp))
182			octerror(USER, "truncated octree");
183			return(r);
184	greg	1.1	}
185
186
187			static OCTREE
188			gettree() /* get a pre-ordered octree */
189			{
190	greg	2.7	register OCTREE ot;
191	greg	1.1	register int i;
192
193			switch (getc(infp)) {
194			case OT_EMPTY:
195			return(EMPTY);
196			case OT_FULL:
197			return(getfullnode());
198			case OT_TREE:
199			if ((ot = octalloc()) == EMPTY)
200			octerror(SYSTEM, "out of tree space in gettree");
201			for (i = 0; i < 8; i++)
202			octkid(ot, i) = gettree();
203			return(ot);
204	greg	2.5	case EOF:
205			octerror(USER, "truncated octree");
206			default:
207			octerror(USER, "damaged octree");
208			}
209	greg	2.22	return EMPTY; /* pro forma return */
210	greg	2.12	}
211
212
213			static int
214	greg	2.27	nonsurfintree(OCTREE ot) /* check tree for modifiers */
215	greg	2.12	{
216	greg	2.27	OBJECT set[MAXSET+1];
217	greg	2.12	register int i;
218
219	greg	2.27	if (isempty(ot))
220			return(0);
221			if (istree(ot)) {
222			for (i = 0; i < 8; i++)
223			if (nonsurfintree(octkid(ot, i)))
224			return(1);
225			return(0);
226			}
227			objset(set, ot);
228	greg	2.28	for (i = set[0]; i > 0; i-- )
229	greg	2.27	if (ismodifier(objptr(set[i])->otype))
230	greg	2.12	return(1);
231			return(0);
232	greg	2.5	}
233
234
235	schorsch	2.17	static void
236			skiptree(void) /* skip octree on input */
237	greg	2.5	{
238			register int i;
239
240			switch (getc(infp)) {
241			case OT_EMPTY:
242			return;
243			case OT_FULL:
244	greg	2.6	for (i = ogetint(objsize)*objsize; i-- > 0; )
245			if (getc(infp) == EOF)
246			octerror(USER, "truncated octree");
247	greg	2.5	return;
248			case OT_TREE:
249			for (i = 0; i < 8; i++)
250			skiptree();
251			return;
252	greg	1.1	case EOF:
253			octerror(USER, "truncated octree");
254			default:
255			octerror(USER, "damaged octree");
256			}
257			}
258
259
260	schorsch	2.17	static void
261	greg	2.27	octerror(int etyp, char msg) / octree error */
262	greg	1.1	{
263			char msgbuf[128];
264
265			sprintf(msgbuf, "(%s): %s", infn, msg);
266			error(etyp, msgbuf);
267			}