src/common/readoct.c

#ifndef lint
static const char       RCSid[] = "$Id: readoct.c,v 2.24 2003/10/27 10:19:31 schorsch 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 */
#define getc    getc_unlocked
#endif

static double  ogetflt(void);
static long  ogetint(int);
static char  *ogetstr(char *);
static int  nonsurfinset(OBJECT *);
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(inpspec, load, scene, ofn)      /* 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 (dosets(nonsurfinset))
                        octerror(USER, "modifier in tree; octree stale?");
            }
        }
                                /* close the input */
        if (infn[0] == '!')
                pclose(infp);
        else
                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");
        }
        return EMPTY; /* pro forma return */
}


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

        sprintf(msgbuf, "(%s): %s", infn, msg);
        error(etyp, msgbuf);
}
Revision:	2.25
Committed:	Thu Apr 29 14:36:49 2004 UTC (20 years ago) by greg
Content type:	text/plain
Branch:	MAIN
Changes since 2.24:	+5 -1 lines
Log Message:	Added macros to avoid flockfile(3) overhead
#	User	Rev	Content
1	greg	1.1	#ifndef lint
2	greg	2.25	static const char RCSid[] = "$Id: readoct.c,v 2.24 2003/10/27 10:19:31 schorsch 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			#define getc getc_unlocked
23			#endif
24
25	schorsch	2.20	static double ogetflt(void);
26			static long ogetint(int);
27			static char ogetstr(char );
28			static int nonsurfinset(OBJECT *);
29	schorsch	2.17	static void octerror(int etyp, char *msg);
30			static void skiptree(void);
31	schorsch	2.20	static OCTREE getfullnode(void), gettree(void);
32	greg	1.1
33	greg	2.13	static char infn; / input file specification */
34	greg	1.1	static FILE infp; / input file stream */
35	greg	1.14	static int objsize; /* size of stored OBJECT's */
36	greg	1.1	static OBJECT objorig; /* zeroeth object */
37	greg	2.12	static OBJECT fnobjects; /* number of objects in this file */
38	greg	1.1
39
40			int
41	greg	2.13	readoct(inpspec, load, scene, ofn) /* read in octree file or stream */
42			char *inpspec;
43	greg	1.1	int load;
44			CUBE *scene;
45			char *ofn[];
46			{
47	greg	1.11	char sbuf[512];
48	greg	1.1	int nf;
49	greg	2.15	int i;
50	greg	1.15	long m;
51	greg	1.1
52	greg	2.13	if (inpspec == NULL) {
53	greg	1.1	infn = "standard input";
54			infp = stdin;
55	greg	2.13	} else if (inpspec[0] == '!') {
56			infn = inpspec;
57			if ((infp = popen(inpspec+1, "r")) == NULL) {
58			sprintf(errmsg, "cannot execute \"%s\"", inpspec);
59			error(SYSTEM, errmsg);
60			}
61	greg	1.1	} else {
62	greg	2.13	infn = inpspec;
63			if ((infp = fopen(inpspec, "r")) == NULL) {
64	greg	1.1	sprintf(errmsg, "cannot open octree file \"%s\"",
65	greg	2.13	inpspec);
66	greg	1.1	error(SYSTEM, errmsg);
67			}
68			}
69	schorsch	2.16	SET_FILE_BINARY(infp);
70	greg	1.1	/* get header */
71	greg	2.11	if (checkheader(infp, OCTFMT, load&IO_INFO ? stdout : (FILE *)NULL) < 0)
72	greg	1.12	octerror(USER, "not an octree");
73	greg	1.1	/* check format */
74	greg	2.4	if ((objsize = ogetint(2)-OCTMAGIC) <= 0 \|\|
75	greg	1.14	objsize > MAXOBJSIZ \|\| objsize > sizeof(long))
76	greg	1.12	octerror(USER, "incompatible octree format");
77	greg	1.1	/* get boundaries */
78			if (load & IO_BOUNDS) {
79			for (i = 0; i < 3; i++)
80	greg	2.4	scene->cuorg[i] = atof(ogetstr(sbuf));
81			scene->cusize = atof(ogetstr(sbuf));
82	greg	1.1	} else {
83			for (i = 0; i < 4; i++)
84	greg	2.4	ogetstr(sbuf);
85	greg	1.1	}
86			objorig = nobjects; /* set object offset */
87			nf = 0; /* get object files */
88	greg	2.4	while (*ogetstr(sbuf)) {
89	greg	1.1	if (load & IO_SCENE)
90			readobj(sbuf);
91			if (load & IO_FILES)
92			ofn[nf] = savqstr(sbuf);
93			nf++;
94			}
95			if (load & IO_FILES)
96			ofn[nf] = NULL;
97	greg	1.8	/* get number of objects */
98	greg	2.4	fnobjects = m = ogetint(objsize);
99	greg	1.15	if (fnobjects != m)
100			octerror(USER, "too many objects");
101	greg	1.1
102	greg	2.13	if (load & IO_TREE) /* get the octree */
103	greg	1.1	scene->cutree = gettree();
104	greg	2.5	else if (load & IO_SCENE && nf == 0)
105			skiptree();
106
107	schorsch	2.21	if (load & IO_SCENE) { /* get the scene */
108	greg	2.5	if (nf == 0) {
109	greg	2.15	/* load binary scene data */
110			readscene(infp, objsize);
111
112	greg	2.5	} else { /* consistency checks */
113			/* check object count */
114			if (nobjects != objorig+fnobjects)
115			octerror(USER, "bad object count; octree stale?");
116			/* check for non-surfaces */
117	greg	2.12	if (dosets(nonsurfinset))
118	greg	2.5	octerror(USER, "modifier in tree; octree stale?");
119			}
120	schorsch	2.21	}
121	greg	2.13	/* close the input */
122			if (infn[0] == '!')
123			pclose(infp);
124			else
125			fclose(infp);
126	greg	1.1	return(nf);
127			}
128
129
130			static char *
131	greg	2.4	ogetstr(s) /* get null-terminated string */
132	greg	1.1	char *s;
133			{
134	greg	2.4	extern char *getstr();
135	greg	1.1
136	greg	2.4	if (getstr(s, infp) == NULL)
137			octerror(USER, "truncated octree");
138			return(s);
139	greg	1.1	}
140
141
142			static OCTREE
143			getfullnode() /* get a set, return fullnode */
144			{
145	greg	2.7	OBJECT set[MAXSET+1];
146	greg	1.1	register int i;
147	greg	1.14	register long m;
148	greg	1.1
149	greg	2.4	if ((set[0] = ogetint(objsize)) > MAXSET)
150	greg	1.1	octerror(USER, "bad set in getfullnode");
151	greg	1.14	for (i = 1; i <= set[0]; i++) {
152	greg	2.4	m = ogetint(objsize) + objorig;
153	greg	1.15	if ((set[i] = m) != m)
154	greg	1.14	octerror(USER, "too many objects");
155			}
156	greg	1.1	return(fullnode(set));
157			}
158
159
160			static long
161	greg	2.4	ogetint(siz) /* get a siz-byte integer */
162			int siz;
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			nonsurfinset(os) /* check set for modifier */
215			register OBJECT *os;
216			{
217			register OBJECT s;
218			register int i;
219
220			for (i = *os; i-- > 0; )
221			if ((s = *++os) >= objorig && s < objorig+fnobjects &&
222			ismodifier(objptr(s)->otype))
223			return(1);
224			return(0);
225	greg	2.5	}
226
227
228	schorsch	2.17	static void
229			skiptree(void) /* skip octree on input */
230	greg	2.5	{
231			register int i;
232
233			switch (getc(infp)) {
234			case OT_EMPTY:
235			return;
236			case OT_FULL:
237	greg	2.6	for (i = ogetint(objsize)*objsize; i-- > 0; )
238			if (getc(infp) == EOF)
239			octerror(USER, "truncated octree");
240	greg	2.5	return;
241			case OT_TREE:
242			for (i = 0; i < 8; i++)
243			skiptree();
244			return;
245	greg	1.1	case EOF:
246			octerror(USER, "truncated octree");
247			default:
248			octerror(USER, "damaged octree");
249			}
250			}
251
252
253	schorsch	2.17	static void
254	greg	1.1	octerror(etyp, msg) /* octree error */
255			int etyp;
256			char *msg;
257			{
258			char msgbuf[128];
259
260			sprintf(msgbuf, "(%s): %s", infn, msg);
261			error(etyp, msgbuf);
262			}