src/common/readoct.c

#ifndef lint
static const char       RCSid[] = "$Id: readoct.c,v 2.32 2023/02/07 20:28:16 greg Exp $";
#endif
/*
 *  readoct.c - routines to read octree information.
 */

#include "copyright.h"

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

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

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