src/cv/mgf2meta.c

#ifndef lint
static const char       RCSid[] = "$Id: mgf2meta.c,v 2.7 2003/06/30 14:59:11 schorsch Exp $";
#endif
/*
 * Convert MGF (Materials and Geometry Format) to Metafile 2-d graphics
 */

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include "random.h"
#include "mgflib/parser.h"

#define MSIZE   ((1<<14)-1)
#define MX(v)   (int)(MSIZE*(v)[(proj_axis+1)%3])
#define MY(v)   (int)(MSIZE*(v)[(proj_axis+2)%3])

int     r_face(int ac, char **av);
int     proj_axis;
double  limit[3][2];
int     layer;
long    rthresh = 1;

extern int      mg_nqcdivs;


main(argc, argv)                /* convert files to stdout */
int     argc;
char    *argv[];
{
        int     i;
                                /* initialize dispatch table */
        mg_ehand[MG_E_FACE] = r_face;
        mg_ehand[MG_E_POINT] = c_hvertex;
        mg_ehand[MG_E_VERTEX] = c_hvertex;
        mg_ehand[MG_E_XF] = xf_handler;
        mg_nqcdivs = 3;         /* reduce object subdivision */
        mg_init();              /* initialize the parser */
                                        /* get arguments */
        if (argc > 9 && !strcmp(argv[1], "-t")) {
                rthresh = atof(argv[2])*MSIZE + 0.5;
                rthresh *= rthresh;
                argv += 2;
                argc -= 2;
        }
        if (argc < 8 || (proj_axis = argv[1][0]-'x') < 0 || proj_axis > 2)
                goto userr;
        limit[0][0] = atof(argv[2]); limit[0][1] = atof(argv[3]);
        limit[1][0] = atof(argv[4]); limit[1][1] = atof(argv[5]);
        limit[2][0] = atof(argv[6]); limit[2][1] = atof(argv[7]);
        
        if (argc == 8) {                /* convert stdin */
                if (mg_load(NULL) != MG_OK)
                        exit(1);
        } else                          /* convert each file */
                for (i = 8; i < argc; i++) {
                        if (mg_load(argv[i]) != MG_OK)
                                exit(1);
                        newlayer();
                }
        mendpage();                     /* print page */
        mdone();                        /* close output */
        exit(0);
userr:
        fputs("Usage: mgf2meta [-t thresh] {x|y|z} xmin xmax ymin ymax zmin zmax [file.mgf] ..\n",
                        stderr);
        exit(1);
}


int
r_face(ac, av)                  /* convert a face */
int     ac;
char    **av;
{
        static FVECT    bbmin = {0,0,0}, bbmax = {1,1,1};
        register int    i, j;
        register C_VERTEX       *cv;
        FVECT   v1, v2, vo;

        if (ac < 4)
                return(MG_EARGC);
                                                /* connect to last point */
        if ((cv = c_getvert(av[ac-1])) == NULL)
                return(MG_EUNDEF);
        xf_xfmpoint(vo, cv->p);
        for (j = 0; j < 3; j++)
                vo[j] = (vo[j] - limit[j][0])/(limit[j][1]-limit[j][0]);
        for (i = 1; i < ac; i++) {              /* go around face */
                if ((cv = c_getvert(av[i])) == NULL)
                        return(MG_EUNDEF);
                xf_xfmpoint(v2, cv->p);
                for (j = 0; j < 3; j++)
                        v2[j] = (v2[j] - limit[j][0])/(limit[j][1]-limit[j][0]);
                VCOPY(v1, vo);
                VCOPY(vo, v2);
                if (clip(v1, v2, bbmin, bbmax))
                        doline(MX(v1), MY(v1), MX(v2), MY(v2));
        }
        return(MG_OK);
}


#define  HTBLSIZ        16381           /* prime hash table size */

short   hshtab[HTBLSIZ][4];             /* done line segments */

#define  hash(mx1,my1,mx2,my2)  ((long)(mx1)<<15 ^ (long)(my1)<<10 ^ \
                                        (long)(mx2)<<5 ^ (long)(my2))


newlayer()                              /* start a new layer */
{
        (void)memset((char *)hshtab, '\0', sizeof(hshtab));
        if (++layer >= 16) {
                mendpage();
                layer = 0;
        }
}


int
doline(v1x, v1y, v2x, v2y)              /* draw line conditionally */
int     v1x, v1y, v2x, v2y;
{
        register int    h;

        if (v1x > v2x || (v1x == v2x && v1y > v2y)) {   /* sort endpoints */
                h=v1x; v1x=v2x; v2x=h;
                h=v1y; v1y=v2y; v2y=h;
        }
        h = hash(v1x, v1y, v2x, v2y) % HTBLSIZ;
        if (hshtab[h][0] == v1x && hshtab[h][1] == v1y &&
                        hshtab[h][2] == v2x && hshtab[h][3] == v2y)
                return(0);
        hshtab[h][0] = v1x; hshtab[h][1] = v1y;
        hshtab[h][2] = v2x; hshtab[h][3] = v2y;
        if ((long)(v2x-v1x)*(v2x-v1x) + (long)(v2y-v1y)*(v2y-v1y)
                        <= random() % rthresh)
                return(0);
        mline(v1x, v1y, layer/4, 0, layer%4);
        mdraw(v2x, v2y);
        return(1);
}
Revision:	2.8
Committed:	Tue Oct 21 19:19:28 2003 UTC (21 years ago) by schorsch
Content type:	text/plain
Branch:	MAIN
Changes since 2.7:	+2 -2 lines
Log Message:	Various platform compatibility fixes.
#	User	Rev	Content
1	greg	2.1	#ifndef lint
2	schorsch	2.8	static const char RCSid[] = "$Id: mgf2meta.c,v 2.7 2003/06/30 14:59:11 schorsch Exp $";
3	greg	2.1	#endif
4			/*
5	greg	2.2	* Convert MGF (Materials and Geometry Format) to Metafile 2-d graphics
6	greg	2.1	*/
7
8			#include <stdio.h>
9	greg	2.6	#include <stdlib.h>
10	schorsch	2.7	#include <string.h>
11	greg	2.1	#include <math.h>
12	greg	2.4	#include "random.h"
13	greg	2.1	#include "mgflib/parser.h"
14
15	greg	2.5	#define MSIZE ((1<<14)-1)
16			#define MX(v) (int)(MSIZE*(v)[(proj_axis+1)%3])
17			#define MY(v) (int)(MSIZE*(v)[(proj_axis+2)%3])
18
19	schorsch	2.8	int r_face(int ac, char **av);
20	greg	2.1	int proj_axis;
21			double limit[3][2];
22	greg	2.3	int layer;
23	greg	2.5	long rthresh = 1;
24	greg	2.1
25			extern int mg_nqcdivs;
26
27
28			main(argc, argv) /* convert files to stdout */
29			int argc;
30			char *argv[];
31			{
32			int i;
33			/* initialize dispatch table */
34			mg_ehand[MG_E_FACE] = r_face;
35			mg_ehand[MG_E_POINT] = c_hvertex;
36			mg_ehand[MG_E_VERTEX] = c_hvertex;
37			mg_ehand[MG_E_XF] = xf_handler;
38			mg_nqcdivs = 3; /* reduce object subdivision */
39			mg_init(); /* initialize the parser */
40			/* get arguments */
41	greg	2.5	if (argc > 9 && !strcmp(argv[1], "-t")) {
42			rthresh = atof(argv[2])*MSIZE + 0.5;
43			rthresh *= rthresh;
44			argv += 2;
45			argc -= 2;
46			}
47	greg	2.1	if (argc < 8 \|\| (proj_axis = argv[1][0]-'x') < 0 \|\| proj_axis > 2)
48			goto userr;
49			limit[0][0] = atof(argv[2]); limit[0][1] = atof(argv[3]);
50			limit[1][0] = atof(argv[4]); limit[1][1] = atof(argv[5]);
51			limit[2][0] = atof(argv[6]); limit[2][1] = atof(argv[7]);
52
53			if (argc == 8) { /* convert stdin */
54			if (mg_load(NULL) != MG_OK)
55			exit(1);
56			} else /* convert each file */
57	greg	2.2	for (i = 8; i < argc; i++) {
58	greg	2.1	if (mg_load(argv[i]) != MG_OK)
59			exit(1);
60	greg	2.4	newlayer();
61	greg	2.2	}
62	greg	2.3	mendpage(); /* print page */
63	greg	2.2	mdone(); /* close output */
64	greg	2.1	exit(0);
65			userr:
66	greg	2.5	fputs("Usage: mgf2meta [-t thresh] {x\|y\|z} xmin xmax ymin ymax zmin zmax [file.mgf] ..\n",
67			stderr);
68	greg	2.1	exit(1);
69			}
70
71
72			int
73			r_face(ac, av) /* convert a face */
74			int ac;
75			char **av;
76			{
77			static FVECT bbmin = {0,0,0}, bbmax = {1,1,1};
78			register int i, j;
79			register C_VERTEX *cv;
80			FVECT v1, v2, vo;
81
82			if (ac < 4)
83			return(MG_EARGC);
84			/* connect to last point */
85			if ((cv = c_getvert(av[ac-1])) == NULL)
86			return(MG_EUNDEF);
87			xf_xfmpoint(vo, cv->p);
88			for (j = 0; j < 3; j++)
89			vo[j] = (vo[j] - limit[j][0])/(limit[j][1]-limit[j][0]);
90			for (i = 1; i < ac; i++) { /* go around face */
91			if ((cv = c_getvert(av[i])) == NULL)
92			return(MG_EUNDEF);
93			xf_xfmpoint(v2, cv->p);
94			for (j = 0; j < 3; j++)
95			v2[j] = (v2[j] - limit[j][0])/(limit[j][1]-limit[j][0]);
96			VCOPY(v1, vo);
97			VCOPY(vo, v2);
98	greg	2.4	if (clip(v1, v2, bbmin, bbmax))
99			doline(MX(v1), MY(v1), MX(v2), MY(v2));
100	greg	2.1	}
101			return(MG_OK);
102	greg	2.4	}
103
104
105			#define HTBLSIZ 16381 /* prime hash table size */
106
107			short hshtab[HTBLSIZ][4]; /* done line segments */
108
109			#define hash(mx1,my1,mx2,my2) ((long)(mx1)<<15 ^ (long)(my1)<<10 ^ \
110			(long)(mx2)<<5 ^ (long)(my2))
111
112
113			newlayer() /* start a new layer */
114			{
115	schorsch	2.7	(void)memset((char *)hshtab, '\0', sizeof(hshtab));
116	greg	2.4	if (++layer >= 16) {
117			mendpage();
118			layer = 0;
119			}
120			}
121
122
123			int
124			doline(v1x, v1y, v2x, v2y) /* draw line conditionally */
125			int v1x, v1y, v2x, v2y;
126			{
127			register int h;
128
129			if (v1x > v2x \|\| (v1x == v2x && v1y > v2y)) { /* sort endpoints */
130			h=v1x; v1x=v2x; v2x=h;
131			h=v1y; v1y=v2y; v2y=h;
132			}
133			h = hash(v1x, v1y, v2x, v2y) % HTBLSIZ;
134			if (hshtab[h][0] == v1x && hshtab[h][1] == v1y &&
135			hshtab[h][2] == v2x && hshtab[h][3] == v2y)
136			return(0);
137			hshtab[h][0] = v1x; hshtab[h][1] = v1y;
138			hshtab[h][2] = v2x; hshtab[h][3] = v2y;
139			if ((long)(v2x-v1x)(v2x-v1x) + (long)(v2y-v1y)(v2y-v1y)
140	greg	2.5	<= random() % rthresh)
141	greg	2.4	return(0);
142			mline(v1x, v1y, layer/4, 0, layer%4);
143			mdraw(v2x, v2y);
144			return(1);
145	greg	2.1	}