cv/mgflib/mgf2inv.c

#ifndef lint
static const char       RCSid[] = "$Id: mgf2inv.c,v 1.9 2003/07/27 22:12:02 schorsch Exp $";
#endif
/*
 * Convert MGF to Inventor file.
 *
 *      December 1995   Greg Ward
 */

#include <stdio.h>

#include <stdlib.h>

#include <math.h>

#include <ctype.h>

#include <string.h>

#include "parser.h"

#include "lookup.h"

#define O_INV1          1       /* Inventor 1.0 output */
#define O_INV2          2       /* Inventor 2.0 output */
#define O_VRML1         3       /* VRML 1.0 output */

#define MAXID           48      /* maximum identifier length */

#define VERTFMT         "%+16.9e %+16.9e %+16.9e\n%+6.3f %+6.3f %+6.3f"
#define VZVECT          "+0.000 +0.000 +0.000"
#define VFSEPPOS        50      /* position of newline in above */
#define VFLEN           72      /* total vertex string length */
#define MAXVERT         10240   /* maximum cached vertices */

#define setvkey(k,v)    sprintf(k,VERTFMT,(v)->p[0],(v)->p[1],(v)->p[2],\
                                        (v)->n[0],(v)->n[1],(v)->n[2]);

char    vlist[MAXVERT][VFLEN];  /* our vertex cache */
int     nverts;                 /* current cache size */

LUTAB   vert_tab = LU_SINIT(NULL,NULL);

struct face {
        struct face     *next;          /* next face in list */
        short           nv;             /* number of vertices */
        short           vl[3];          /* vertex index list (variable) */
}       *flist, *flast;         /* our face cache */

#define newface(n)      (struct face *)malloc(sizeof(struct face) + \
                                ((n) > 3 ? (n)-3 : 0)*sizeof(short))
#define freeface(f)     free((MEM_PTR)f)

#define TABSTOP         8       /* assumed number of characters per tab */
#define SHIFTW          2       /* nesting shift width */
#define MAXIND          15      /* maximum indent level */

char    tabs[MAXIND*SHIFTW+1];  /* current tab-in string */

#define curmatname      (c_cmname == NULL ? "mat" : to_id(c_cmname))

int     outtype = O_INV2;       /* output format */

extern int      i_comment(), i_object(), i_xf(),
                i_cyl(), i_face(), i_sph();

extern char     *to_id();


main(argc, argv)
int     argc;
char    *argv[];
{
        int     i;
                                /* initialize dispatch table */
        mg_ehand[MG_E_COMMENT] = i_comment;     /* we pass comments */
        mg_ehand[MG_E_COLOR] = c_hcolor;        /* they get color */
        mg_ehand[MG_E_CMIX] = c_hcolor;         /* they mix colors */
        mg_ehand[MG_E_CSPEC] = c_hcolor;        /* they get spectra */
        mg_ehand[MG_E_CXY] = c_hcolor;          /* they get chromaticities */
        mg_ehand[MG_E_CCT] = c_hcolor;          /* they get color temp's */
        mg_ehand[MG_E_CYL] = i_cyl;             /* we do cylinders */
        mg_ehand[MG_E_ED] = c_hmaterial;        /* they get emission */
        mg_ehand[MG_E_FACE] = i_face;           /* we do faces */
        mg_ehand[MG_E_MATERIAL] = c_hmaterial;  /* they get materials */
        mg_ehand[MG_E_NORMAL] = c_hvertex;      /* they get normals */
        mg_ehand[MG_E_OBJECT] = i_object;       /* we track object names */
        mg_ehand[MG_E_POINT] = c_hvertex;       /* they get points */
        mg_ehand[MG_E_RD] = c_hmaterial;        /* they get diffuse refl. */
        mg_ehand[MG_E_RS] = c_hmaterial;        /* they get specular refl. */
        mg_ehand[MG_E_SIDES] = c_hmaterial;     /* they get # sides */
        mg_ehand[MG_E_SPH] = i_sph;             /* we do spheres */
        mg_ehand[MG_E_TD] = c_hmaterial;        /* they get diffuse trans. */
        mg_ehand[MG_E_TS] = c_hmaterial;        /* they get specular trans. */
        mg_ehand[MG_E_VERTEX] = c_hvertex;      /* they get vertices */
        mg_ehand[MG_E_XF] = i_xf;               /* we track transforms */
        mg_init();              /* initialize the parser */
                                /* get options and print format line */
        for (i = 1; i < argc && argv[i][0] == '-'; i++)
                if (!strcmp(argv[i], "-vrml"))
                        outtype = O_VRML1;
                else if (!strcmp(argv[i], "-1"))
                        outtype = O_INV1;
                else if (!strcmp(argv[i], "-2"))
                        outtype = O_INV2;
                else
                        goto userr;
        switch (outtype) {
        case O_INV1:
                printf("#Inventor V1.0 ascii\n");
                break;
        case O_INV2:
                printf("#Inventor V2.0 ascii\n");
                break;
        case O_VRML1:
                printf("#VRML V1.0 ascii\n");
                break;
        }
        printf("## Translated from MGF Version %d.%d\n", MG_VMAJOR, MG_VMINOR);
        printf("Separator {\n");                /* begin root node */
                                                /* general properties */
        printf("MaterialBinding { value OVERALL }\n");
        printf("NormalBinding { value PER_VERTEX_INDEXED }\n");
        if (outtype != O_INV1) {
                printf("ShapeHints {\n");
                printf("\tvertexOrdering CLOCKWISE\n");
                printf("\tfaceType UNKNOWN_FACE_TYPE\n");
                printf("}\n");
        }
        if (i == argc) {        /* load standard input */
                if (mg_load(NULL) != MG_OK)
                        exit(1);
                if (mg_nunknown)
                        printf("## %s: %u unknown entities\n",
                                        argv[0], mg_nunknown);
        }
                                /* load MGF files */
        for ( ; i < argc; i++) {
                printf("## %s %s ##############################\n",
                                argv[0], argv[i]);
                mg_nunknown = 0;
                if (mg_load(argv[i]) != MG_OK)
                        exit(1);
                if (mg_nunknown)
                        printf("## %s %s: %u unknown entities\n",
                                        argv[0], argv[i], mg_nunknown);
        }
        flush_cache();          /* flush face cache, just in case */
        printf("}\n");          /* close root node */
        exit(0);
userr:
        fprintf(stderr, "%s: [-1|-2|-vrml] [file] ..\n", argv[0]);
        exit(1);
}


indent(deeper)                          /* indent in or out */
int     deeper;
{
        static int      level;          /* current nesting level */
        register int    i;
        register char   *cp;

        if (deeper) level++;            /* in or out? */
        else if (level > 0) level--;
                                        /* compute actual shift */
        if ((i = level) > MAXIND) i = MAXIND;
        cp = tabs;
        for (i *= SHIFTW; i >= TABSTOP; i -= TABSTOP)
                *cp++ = '\t';
        while (i--)
                *cp++ = ' ';
        *cp = '\0';
}


int
i_comment(ac, av)                       /* transfer comment as is */
int     ac;
char    **av;
{
        fputs(tabs, stdout);
        putchar('#');                   /* Inventor comment character */
        while (--ac > 0) {
                putchar(' ');
                fputs(*++av, stdout);
        }
        putchar('\n');
        return(MG_OK);
}


int
i_object(ac, av)                        /* group object name */
int     ac;
char    **av;
{
        static int      objnest;

        flush_cache();                  /* flush cached objects */
        if (ac == 2) {                          /* start group */
                printf("%sDEF %s Group {\n", tabs, to_id(av[1]));
                indent(1);
                objnest++;
                return(MG_OK);
        }
        if (ac == 1) {                          /* end group */
                if (--objnest < 0)
                        return(MG_ECNTXT);
                indent(0);
                fputs(tabs, stdout);
                fputs("}\n", stdout);
                return(MG_OK);
        }
        return(MG_EARGC);
}


int
i_xf(ac, av)                            /* transform object(s) */
int     ac;
char    **av;
{
        static long     xfid;
        register XF_SPEC        *spec;

        flush_cache();                  /* flush cached objects */
        if (ac == 1) {                  /* end of transform */
                if ((spec = xf_context) == NULL)
                        return(MG_ECNTXT);
                indent(0);                      /* close original segment */
                printf("%s}\n", tabs);
                indent(0);
                printf("%s}\n", tabs);
                if (spec->xarr != NULL) {       /* check for iteration */
                        register struct xf_array        *ap = spec->xarr;
                        register int    n;

                        ap->aarg[ap->ndim-1].i = 1;     /* iterate array */
                        for ( ; ; ) {
                                n = ap->ndim-1;
                                while (ap->aarg[n].i < ap->aarg[n].n) {
                                        sprintf(ap->aarg[n].arg, "%d",
                                                        ap->aarg[n].i);
                                        printf("%sSeparator {\n", tabs);
                                        indent(1);
                                        (void)put_xform(spec);
                                        printf("%sUSE _xf%ld\n", tabs,
                                                        spec->xid);
                                        indent(0);
                                        printf("%s}\n", tabs);
                                        ++ap->aarg[n].i;
                                }
                                ap->aarg[n].i = 0;
                                (void)strcpy(ap->aarg[n].arg, "0");
                                while (n-- && ++ap->aarg[n].i >= ap->aarg[n].n) {
                                        ap->aarg[n].i = 0;
                                        (void)strcpy(ap->aarg[n].arg, "0");
                                }
                                if (n < 0)
                                        break;
                                sprintf(ap->aarg[n].arg, "%d", ap->aarg[n].i);
                        }
                }
                                                /* pop transform */
                xf_context = spec->prev;
                free_xf(spec);
                return(MG_OK);
        }
                                        /* else allocate new transform */
        if ((spec = new_xf(ac-1, av+1)) == NULL)
                return(MG_EMEM);
        spec->xid = ++xfid;             /* assign unique ID */
        spec->prev = xf_context;        /* push onto stack */
        xf_context = spec;
                                        /* translate xf specification */
        printf("%sSeparator {\n", tabs);
        indent(1);
        if (put_xform(spec) < 0)
                return(MG_ETYPE);
        printf("%sDEF _xf%ld Group {\n", tabs, spec->xid);      /* begin */
        indent(1);
        return(MG_OK);
}


int
put_xform(spec)                 /* translate and print transform */
register XF_SPEC        *spec;
{
        register char   **av;
        register int    n;

        n = xf_ac(spec) - xf_ac(spec->prev);
        if (xf(&spec->xf, n, av=xf_av(spec)) != n)
                return(-1);
        printf("%sMatrixTransform {\n", tabs);
        indent(1);
        printf("%s# xf", tabs);         /* put out original as comment */
        while (n--) {
                putchar(' ');
                fputs(*av++, stdout);
        }
        putchar('\n');                  /* put out computed matrix */
        printf("%smatrix %13.9g %13.9g %13.9g %13.9g\n", tabs,
                        spec->xf.xfm[0][0], spec->xf.xfm[0][1],
                        spec->xf.xfm[0][2], spec->xf.xfm[0][3]);
        for (n = 1; n < 4; n++)
                printf("%s       %13.9g %13.9g %13.9g %13.9g\n", tabs,
                                spec->xf.xfm[n][0], spec->xf.xfm[n][1],
                                spec->xf.xfm[n][2], spec->xf.xfm[n][3]);
        indent(0);
        printf("%s}\n", tabs);
        return(0);
}


int
put_material()                  /* put out current material */
{
        char    *mname = curmatname;
        float   rgbval[3];

        if (!c_cmaterial->clock) {      /* current, just use it */
                printf("%sUSE %s\n", tabs, mname);
                return(0);
        }
                                /* else update definition */
        printf("%sDEF %s Group {\n", tabs, mname);
        indent(1);
        printf("%sMaterial {\n", tabs);
        indent(1);
        mgf2rgb(&c_cmaterial->rd_c, c_cmaterial->rd, rgbval);
        printf("%sambientColor %.4f %.4f %.4f\n", tabs,
                        rgbval[0], rgbval[1], rgbval[2]);
        printf("%sdiffuseColor %.4f %.4f %.4f\n", tabs,
                        rgbval[0], rgbval[1], rgbval[2]);
        if (c_cmaterial->rs > FTINY) {
                mgf2rgb(&c_cmaterial->rs_c, c_cmaterial->rs, rgbval);
                printf("%sspecularColor %.4f %.4f %.4f\n", tabs,
                                rgbval[0], rgbval[1], rgbval[2]);
                printf("%sshininess %.3f\n", tabs, 1.-sqrt(c_cmaterial->rs_a));
        }
        if (c_cmaterial->ed > FTINY) {
                mgf2rgb(&c_cmaterial->ed_c, 1.0, rgbval);
                printf("%semissiveColor %.4f %.4f %.4f\n", tabs,
                                rgbval[0], rgbval[1], rgbval[2]);
        }
        if (c_cmaterial->ts > FTINY)
                printf("%stransparency %.4f\n", tabs,
                                c_cmaterial->ts + c_cmaterial->td);
        indent(0);
        printf("%s}\n", tabs);
        if (outtype != O_INV1)
                printf("%sShapeHints { shapeType %s faceType UNKNOWN_FACE_TYPE }\n",
                        tabs,
                        c_cmaterial->sided ? "SOLID" : "UNKNOWN_SHAPE_TYPE");
        indent(0);
        printf("%s}\n", tabs);
        c_cmaterial->clock = 0;
        return(0);
}


int
i_face(ac, av)                  /* translate an N-sided face */
int     ac;
char    **av;
{
        static char     lastmat[MAXID];
        struct face     *newf;
        register C_VERTEX       *vp;
        register LUENT  *lp;
        register int    i;

        if (ac < 4)
                return(MG_EARGC);
        if ( strcmp(lastmat, curmatname) || c_cmaterial->clock ||
                        nverts == 0 || nverts+ac-1 >= MAXVERT) {
                flush_cache();                  /* new cache */
                lu_init(&vert_tab, MAXVERT);
                printf("%sSeparator {\n", tabs);
                indent(1);
                if (put_material() < 0)         /* put out material */
                        return(MG_EBADMAT);
                (void)strcpy(lastmat, curmatname);
        }
                                /* allocate new face */
        if ((newf = newface(ac-1)) == NULL)
                return(MG_EMEM);
        newf->nv = ac-1;
                                /* get vertex references */
        for (i = 0; i < newf->nv; i++) {
                if ((vp = c_getvert(av[i+1])) == NULL)
                        return(MG_EUNDEF);
                setvkey(vlist[nverts], vp);
                lp = lu_find(&vert_tab, vlist[nverts]);
                if (lp == NULL)
                        return(MG_EMEM);
                if (lp->key == NULL)
                        lp->key = (char *)vlist[nverts++];
                newf->vl[i] = ((char (*)[VFLEN])lp->key - vlist);
        }
                                /* add to face list */
        newf->next = NULL;
        if (flist == NULL)
                flist = newf;
        else
                flast->next = newf;
        flast = newf;
        return(MG_OK);          /* we'll actually put it out later */
}


int
i_sph(ac, av)                   /* translate sphere description */
int     ac;
char    **av;
{
        register C_VERTEX       *cent;

        if (ac != 3)
                return(MG_EARGC);
        flush_cache();          /* flush vertex cache */
        printf("%sSeparator {\n", tabs);
        indent(1);
                                /* put out current material */
        if (put_material() < 0)
                return(MG_EBADMAT);
                                /* get center */
        if ((cent = c_getvert(av[1])) == NULL)
                return(MG_EUNDEF);
                                /* get radius */
        if (!isflt(av[2]))
                return(MG_ETYPE);
        printf("%sTranslation { translation %13.9g %13.9g %13.9g }\n", tabs,
                        cent->p[0], cent->p[1], cent->p[2]);
        printf("%sSphere { radius %s }\n", tabs, av[2]);
        indent(0);
        printf("%s}\n", tabs);
        return(MG_OK);
}


int
i_cyl(ac, av)                   /* translate a cylinder description */
int     ac;
char    **av;
{
        register C_VERTEX       *v1, *v2;
        FVECT   va;
        double  length, angle;

        if (ac != 4)
                return(MG_EARGC);
        flush_cache();          /* flush vertex cache */
        printf("%sSeparator {\n", tabs);
        indent(1);
                                /* put out current material */
        if (put_material() < 0)
                return(MG_EBADMAT);
                                /* get endpoints */
        if (((v1 = c_getvert(av[1])) == NULL) | ((v2 = c_getvert(av[3])) == NULL))
                return(MG_EUNDEF);
                                /* get radius */
        if (!isflt(av[2]))
                return(MG_ETYPE);
                                /* compute transform */
        va[0] = v2->p[0] - v1->p[0];
        va[1] = v2->p[1] - v1->p[1];
        va[2] = v2->p[2] - v1->p[2];
        length = sqrt(DOT(va,va));
        if (va[1] >= length)
                angle = 0.;
        else if (va[1] <= -length)
                angle = PI;
        else
                angle = acos(va[1]/length);
        printf("%sTranslation { translation %13.9g %13.9g %13.9g }\n", tabs,
                        .5*(v1->p[0]+v2->p[0]), .5*(v1->p[1]+v2->p[1]),
                        .5*(v1->p[2]+v2->p[2]));
        printf("%sRotation { rotation %.9g %.9g %.9g %.9g }\n", tabs,
                        va[2], 0., -va[0], angle);
                                /* open-ended */
        printf("%sCylinder { parts SIDES height %13.9g radius %s }\n", tabs,
                        length, av[2]);
        indent(0);
        printf("%s}\n", tabs);
        return(MG_OK);
}


char *
to_id(name)                     /* make sure a name is a valid Inventor ID */
register char   *name;
{
        static char     id[MAXID];
        register char   *cp;

        for (cp = id; *name && cp < MAXID-1+id; name++)
                if (isalnum(*name) || *name == '_')
                        *cp++ = *name;
                else
                        *cp++ = '_';
        *cp = '\0';
        return(id);
}


flush_cache()                   /* put out cached faces */
{
        int     donorms = 0;
        register struct face    *f;
        register int    i;

        if (nverts == 0)
                return;
                                        /* put out coordinates */
        printf("%sCoordinate3 {\n", tabs);
        indent(1);
        vlist[0][VFSEPPOS] = '\0';
        printf("%spoint [ %s", tabs, vlist[0]);
        for (i = 1; i < nverts; i++) {
                vlist[i][VFSEPPOS] = '\0';
                printf(",\n%s        %s", tabs, vlist[i]);
                if (strcmp(VFSEPPOS+1+vlist[i], VZVECT))
                        donorms++;
        }
        indent(0);
        printf(" ]\n%s}\n", tabs);
        if (donorms) {                  /* put out normals */
                printf("%sNormal {\n", tabs);
                indent(1);
                printf("%svector [ %s", tabs, VFSEPPOS+1+vlist[0]);
                for (i = 1; i < nverts; i++)
                        printf(",\n%s         %s", tabs, VFSEPPOS+1+vlist[i]);
                indent(0);
                printf(" ]\n%s}\n", tabs);
        }
                                        /* put out faces */
        printf("%sIndexedFaceSet {\n", tabs);
        indent(1);
        f = flist;                      /* coordinate indices */
        printf("%scoordIndex [ %d", tabs, f->vl[0]);
        for (i = 1; i < f->nv; i++)
                printf(", %d", f->vl[i]);
        for (f = f->next; f != NULL; f = f->next) {
                printf(", -1,\n%s             %d", tabs, f->vl[0]);
                for (i = 1; i < f->nv; i++)
                        printf(", %d", f->vl[i]);
        }
        printf(" ]\n");
        if (donorms) {
                f = flist;                      /* normal indices */
                printf("%snormalIndex [ %d", tabs, f->vl[0]);
                for (i = 1; i < f->nv; i++)
                        printf(", %d", f->vl[i]);
                for (f = f->next; f != NULL; f = f->next) {
                        printf(", -1,\n%s              %d", tabs, f->vl[0]);
                        for (i = 1; i < f->nv; i++)
                                printf(", %d", f->vl[i]);
                }
                printf(" ]\n");
        }
        indent(0);                      /* close IndexedFaceSet */
        printf("%s}\n", tabs);
        indent(0);                      /* close face group */
        printf("%s}\n", tabs);
        while ((f = flist) != NULL) {   /* free face list */
                flist = f->next;
                freeface(f);
        }
        lu_done(&vert_tab);             /* clear lookup table */
        nverts = 0;
}
Revision:	1.10
Committed:	Tue Aug 5 20:40:16 2003 UTC (20 years, 9 months ago) by greg
Content type:	text/plain
Branch:	MAIN
Changes since 1.9:	+3 -2 lines
Log Message:	Fixed minor bug in mgf2inv where faceType was not specified as unknown
#	User	Rev	Content
1	greg	1.1	#ifndef lint
2	greg	1.10	static const char RCSid[] = "$Id: mgf2inv.c,v 1.9 2003/07/27 22:12:02 schorsch Exp $";
3	greg	1.1	#endif
4			/*
5			* Convert MGF to Inventor file.
6			*
7			* December 1995 Greg Ward
8			*/
9
10			#include <stdio.h>
11	greg	1.8
12			#include <stdlib.h>
13	greg	1.1
14			#include <math.h>
15
16	greg	1.2	#include <ctype.h>
17
18	greg	1.6	#include <string.h>
19
20	greg	1.1	#include "parser.h"
21
22	greg	1.3	#include "lookup.h"
23
24	greg	1.2	#define O_INV1 1 /* Inventor 1.0 output */
25			#define O_INV2 2 /* Inventor 2.0 output */
26	gregl	1.7	#define O_VRML1 3 /* VRML 1.0 output */
27	greg	1.1
28	greg	1.3	#define MAXID 48 /* maximum identifier length */
29
30			#define VERTFMT "%+16.9e %+16.9e %+16.9e\n%+6.3f %+6.3f %+6.3f"
31			#define VZVECT "+0.000 +0.000 +0.000"
32			#define VFSEPPOS 50 /* position of newline in above */
33			#define VFLEN 72 /* total vertex string length */
34			#define MAXVERT 10240 /* maximum cached vertices */
35
36			#define setvkey(k,v) sprintf(k,VERTFMT,(v)->p[0],(v)->p[1],(v)->p[2],\
37			(v)->n[0],(v)->n[1],(v)->n[2]);
38
39			char vlist[MAXVERT][VFLEN]; /* our vertex cache */
40			int nverts; /* current cache size */
41
42			LUTAB vert_tab = LU_SINIT(NULL,NULL);
43
44			struct face {
45			struct face next; / next face in list */
46			short nv; /* number of vertices */
47			short vl[3]; /* vertex index list (variable) */
48			} flist, flast; /* our face cache */
49
50			#define newface(n) (struct face *)malloc(sizeof(struct face) + \
51			((n) > 3 ? (n)-3 : 0)*sizeof(short))
52			#define freeface(f) free((MEM_PTR)f)
53
54	greg	1.1	#define TABSTOP 8 /* assumed number of characters per tab */
55			#define SHIFTW 2 /* nesting shift width */
56			#define MAXIND 15 /* maximum indent level */
57
58	greg	1.2	char tabs[MAXINDSHIFTW+1]; / current tab-in string */
59	greg	1.1
60	greg	1.3	#define curmatname (c_cmname == NULL ? "mat" : to_id(c_cmname))
61
62	greg	1.2	int outtype = O_INV2; /* output format */
63	greg	1.1
64	greg	1.2	extern int i_comment(), i_object(), i_xf(),
65			i_cyl(), i_face(), i_sph();
66	greg	1.1
67	greg	1.2	extern char *to_id();
68	greg	1.1
69
70			main(argc, argv)
71			int argc;
72			char *argv[];
73			{
74			int i;
75			/* initialize dispatch table */
76			mg_ehand[MG_E_COMMENT] = i_comment; /* we pass comments */
77			mg_ehand[MG_E_COLOR] = c_hcolor; /* they get color */
78			mg_ehand[MG_E_CMIX] = c_hcolor; /* they mix colors */
79			mg_ehand[MG_E_CSPEC] = c_hcolor; /* they get spectra */
80			mg_ehand[MG_E_CXY] = c_hcolor; /* they get chromaticities */
81			mg_ehand[MG_E_CCT] = c_hcolor; /* they get color temp's */
82			mg_ehand[MG_E_CYL] = i_cyl; /* we do cylinders */
83			mg_ehand[MG_E_ED] = c_hmaterial; /* they get emission */
84			mg_ehand[MG_E_FACE] = i_face; /* we do faces */
85			mg_ehand[MG_E_MATERIAL] = c_hmaterial; /* they get materials */
86			mg_ehand[MG_E_NORMAL] = c_hvertex; /* they get normals */
87			mg_ehand[MG_E_OBJECT] = i_object; /* we track object names */
88			mg_ehand[MG_E_POINT] = c_hvertex; /* they get points */
89			mg_ehand[MG_E_RD] = c_hmaterial; /* they get diffuse refl. */
90			mg_ehand[MG_E_RS] = c_hmaterial; /* they get specular refl. */
91			mg_ehand[MG_E_SIDES] = c_hmaterial; /* they get # sides */
92			mg_ehand[MG_E_SPH] = i_sph; /* we do spheres */
93			mg_ehand[MG_E_TD] = c_hmaterial; /* they get diffuse trans. */
94			mg_ehand[MG_E_TS] = c_hmaterial; /* they get specular trans. */
95			mg_ehand[MG_E_VERTEX] = c_hvertex; /* they get vertices */
96			mg_ehand[MG_E_XF] = i_xf; /* we track transforms */
97			mg_init(); /* initialize the parser */
98	greg	1.2	/* get options and print format line */
99			for (i = 1; i < argc && argv[i][0] == '-'; i++)
100			if (!strcmp(argv[i], "-vrml"))
101	gregl	1.7	outtype = O_VRML1;
102	greg	1.2	else if (!strcmp(argv[i], "-1"))
103			outtype = O_INV1;
104			else if (!strcmp(argv[i], "-2"))
105			outtype = O_INV2;
106			else
107			goto userr;
108			switch (outtype) {
109			case O_INV1:
110			printf("#Inventor V1.0 ascii\n");
111			break;
112			case O_INV2:
113			printf("#Inventor V2.0 ascii\n");
114			break;
115	gregl	1.7	case O_VRML1:
116			printf("#VRML V1.0 ascii\n");
117	greg	1.2	break;
118			}
119	greg	1.1	printf("## Translated from MGF Version %d.%d\n", MG_VMAJOR, MG_VMINOR);
120			printf("Separator {\n"); /* begin root node */
121			/* general properties */
122			printf("MaterialBinding { value OVERALL }\n");
123	greg	1.2	printf("NormalBinding { value PER_VERTEX_INDEXED }\n");
124			if (outtype != O_INV1) {
125			printf("ShapeHints {\n");
126			printf("\tvertexOrdering CLOCKWISE\n");
127			printf("\tfaceType UNKNOWN_FACE_TYPE\n");
128			printf("}\n");
129			}
130	greg	1.1	if (i == argc) { /* load standard input */
131			if (mg_load(NULL) != MG_OK)
132			exit(1);
133			if (mg_nunknown)
134			printf("## %s: %u unknown entities\n",
135			argv[0], mg_nunknown);
136			}
137			/* load MGF files */
138			for ( ; i < argc; i++) {
139			printf("## %s %s ##############################\n",
140			argv[0], argv[i]);
141			mg_nunknown = 0;
142			if (mg_load(argv[i]) != MG_OK)
143			exit(1);
144			if (mg_nunknown)
145			printf("## %s %s: %u unknown entities\n",
146			argv[0], argv[i], mg_nunknown);
147			}
148	greg	1.3	flush_cache(); /* flush face cache, just in case */
149			printf("}\n"); /* close root node */
150	greg	1.1	exit(0);
151	greg	1.2	userr:
152			fprintf(stderr, "%s: [-1\|-2\|-vrml] [file] ..\n", argv[0]);
153			exit(1);
154	greg	1.1	}
155
156
157			indent(deeper) /* indent in or out */
158			int deeper;
159			{
160			static int level; /* current nesting level */
161			register int i;
162			register char *cp;
163
164			if (deeper) level++; /* in or out? */
165			else if (level > 0) level--;
166			/* compute actual shift */
167			if ((i = level) > MAXIND) i = MAXIND;
168			cp = tabs;
169			for (i *= SHIFTW; i >= TABSTOP; i -= TABSTOP)
170			*cp++ = '\t';
171			while (i--)
172			*cp++ = ' ';
173			*cp = '\0';
174			}
175
176
177			int
178			i_comment(ac, av) /* transfer comment as is */
179			int ac;
180			char **av;
181			{
182			fputs(tabs, stdout);
183			putchar('#'); /* Inventor comment character */
184			while (--ac > 0) {
185			putchar(' ');
186			fputs(*++av, stdout);
187			}
188			putchar('\n');
189			return(MG_OK);
190			}
191
192
193			int
194			i_object(ac, av) /* group object name */
195			int ac;
196			char **av;
197			{
198			static int objnest;
199
200	greg	1.3	flush_cache(); /* flush cached objects */
201	greg	1.1	if (ac == 2) { /* start group */
202	greg	1.2	printf("%sDEF %s Group {\n", tabs, to_id(av[1]));
203	greg	1.1	indent(1);
204			objnest++;
205			return(MG_OK);
206			}
207			if (ac == 1) { /* end group */
208			if (--objnest < 0)
209			return(MG_ECNTXT);
210			indent(0);
211			fputs(tabs, stdout);
212			fputs("}\n", stdout);
213			return(MG_OK);
214			}
215			return(MG_EARGC);
216			}
217
218
219			int
220			i_xf(ac, av) /* transform object(s) */
221			int ac;
222			char **av;
223			{
224			static long xfid;
225			register XF_SPEC *spec;
226
227	greg	1.3	flush_cache(); /* flush cached objects */
228	greg	1.1	if (ac == 1) { /* end of transform */
229			if ((spec = xf_context) == NULL)
230			return(MG_ECNTXT);
231			indent(0); /* close original segment */
232			printf("%s}\n", tabs);
233			indent(0);
234			printf("%s}\n", tabs);
235			if (spec->xarr != NULL) { /* check for iteration */
236			register struct xf_array *ap = spec->xarr;
237			register int n;
238
239			ap->aarg[ap->ndim-1].i = 1; /* iterate array */
240			for ( ; ; ) {
241			n = ap->ndim-1;
242			while (ap->aarg[n].i < ap->aarg[n].n) {
243			sprintf(ap->aarg[n].arg, "%d",
244			ap->aarg[n].i);
245			printf("%sSeparator {\n", tabs);
246			indent(1);
247			(void)put_xform(spec);
248			printf("%sUSE _xf%ld\n", tabs,
249			spec->xid);
250			indent(0);
251			printf("%s}\n", tabs);
252			++ap->aarg[n].i;
253			}
254			ap->aarg[n].i = 0;
255			(void)strcpy(ap->aarg[n].arg, "0");
256			while (n-- && ++ap->aarg[n].i >= ap->aarg[n].n) {
257			ap->aarg[n].i = 0;
258			(void)strcpy(ap->aarg[n].arg, "0");
259			}
260			if (n < 0)
261			break;
262			sprintf(ap->aarg[n].arg, "%d", ap->aarg[n].i);
263			}
264			}
265			/* pop transform */
266			xf_context = spec->prev;
267			free_xf(spec);
268			return(MG_OK);
269			}
270			/* else allocate new transform */
271			if ((spec = new_xf(ac-1, av+1)) == NULL)
272			return(MG_EMEM);
273			spec->xid = ++xfid; /* assign unique ID */
274			spec->prev = xf_context; /* push onto stack */
275			xf_context = spec;
276			/* translate xf specification */
277			printf("%sSeparator {\n", tabs);
278			indent(1);
279			if (put_xform(spec) < 0)
280			return(MG_ETYPE);
281			printf("%sDEF _xf%ld Group {\n", tabs, spec->xid); /* begin */
282			indent(1);
283			return(MG_OK);
284			}
285
286
287			int
288			put_xform(spec) /* translate and print transform */
289			register XF_SPEC *spec;
290			{
291			register char **av;
292			register int n;
293
294			n = xf_ac(spec) - xf_ac(spec->prev);
295			if (xf(&spec->xf, n, av=xf_av(spec)) != n)
296			return(-1);
297			printf("%sMatrixTransform {\n", tabs);
298			indent(1);
299			printf("%s# xf", tabs); /* put out original as comment */
300			while (n--) {
301			putchar(' ');
302			fputs(*av++, stdout);
303			}
304			putchar('\n'); /* put out computed matrix */
305			printf("%smatrix %13.9g %13.9g %13.9g %13.9g\n", tabs,
306			spec->xf.xfm[0][0], spec->xf.xfm[0][1],
307			spec->xf.xfm[0][2], spec->xf.xfm[0][3]);
308			for (n = 1; n < 4; n++)
309			printf("%s %13.9g %13.9g %13.9g %13.9g\n", tabs,
310			spec->xf.xfm[n][0], spec->xf.xfm[n][1],
311			spec->xf.xfm[n][2], spec->xf.xfm[n][3]);
312			indent(0);
313			printf("%s}\n", tabs);
314			return(0);
315			}
316
317
318			int
319			put_material() /* put out current material */
320			{
321	greg	1.3	char *mname = curmatname;
322	greg	1.1	float rgbval[3];
323
324			if (!c_cmaterial->clock) { /* current, just use it */
325	greg	1.3	printf("%sUSE %s\n", tabs, mname);
326	greg	1.1	return(0);
327			}
328			/* else update definition */
329	greg	1.3	printf("%sDEF %s Group {\n", tabs, mname);
330	greg	1.1	indent(1);
331			printf("%sMaterial {\n", tabs);
332			indent(1);
333			mgf2rgb(&c_cmaterial->rd_c, c_cmaterial->rd, rgbval);
334	greg	1.2	printf("%sambientColor %.4f %.4f %.4f\n", tabs,
335	greg	1.1	rgbval[0], rgbval[1], rgbval[2]);
336	greg	1.2	printf("%sdiffuseColor %.4f %.4f %.4f\n", tabs,
337	greg	1.1	rgbval[0], rgbval[1], rgbval[2]);
338			if (c_cmaterial->rs > FTINY) {
339			mgf2rgb(&c_cmaterial->rs_c, c_cmaterial->rs, rgbval);
340	greg	1.2	printf("%sspecularColor %.4f %.4f %.4f\n", tabs,
341	greg	1.1	rgbval[0], rgbval[1], rgbval[2]);
342	greg	1.4	printf("%sshininess %.3f\n", tabs, 1.-sqrt(c_cmaterial->rs_a));
343	greg	1.1	}
344			if (c_cmaterial->ed > FTINY) {
345			mgf2rgb(&c_cmaterial->ed_c, 1.0, rgbval);
346			printf("%semissiveColor %.4f %.4f %.4f\n", tabs,
347			rgbval[0], rgbval[1], rgbval[2]);
348			}
349			if (c_cmaterial->ts > FTINY)
350			printf("%stransparency %.4f\n", tabs,
351			c_cmaterial->ts + c_cmaterial->td);
352			indent(0);
353			printf("%s}\n", tabs);
354	greg	1.2	if (outtype != O_INV1)
355	greg	1.10	printf("%sShapeHints { shapeType %s faceType UNKNOWN_FACE_TYPE }\n",
356			tabs,
357	greg	1.1	c_cmaterial->sided ? "SOLID" : "UNKNOWN_SHAPE_TYPE");
358			indent(0);
359			printf("%s}\n", tabs);
360			c_cmaterial->clock = 0;
361			return(0);
362			}
363
364
365			int
366			i_face(ac, av) /* translate an N-sided face */
367			int ac;
368			char **av;
369			{
370	greg	1.3	static char lastmat[MAXID];
371			struct face *newf;
372			register C_VERTEX *vp;
373			register LUENT *lp;
374	greg	1.1	register int i;
375
376			if (ac < 4)
377			return(MG_EARGC);
378	greg	1.3	if ( strcmp(lastmat, curmatname) \|\| c_cmaterial->clock \|\|
379			nverts == 0 \|\| nverts+ac-1 >= MAXVERT) {
380			flush_cache(); /* new cache */
381			lu_init(&vert_tab, MAXVERT);
382			printf("%sSeparator {\n", tabs);
383			indent(1);
384			if (put_material() < 0) /* put out material */
385			return(MG_EBADMAT);
386			(void)strcpy(lastmat, curmatname);
387			}
388			/* allocate new face */
389			if ((newf = newface(ac-1)) == NULL)
390			return(MG_EMEM);
391			newf->nv = ac-1;
392	greg	1.1	/* get vertex references */
393	greg	1.3	for (i = 0; i < newf->nv; i++) {
394			if ((vp = c_getvert(av[i+1])) == NULL)
395	greg	1.1	return(MG_EUNDEF);
396	greg	1.3	setvkey(vlist[nverts], vp);
397			lp = lu_find(&vert_tab, vlist[nverts]);
398			if (lp == NULL)
399			return(MG_EMEM);
400			if (lp->key == NULL)
401			lp->key = (char *)vlist[nverts++];
402			newf->vl[i] = ((char (*)[VFLEN])lp->key - vlist);
403	greg	1.1	}
404	greg	1.3	/* add to face list */
405			newf->next = NULL;
406			if (flist == NULL)
407			flist = newf;
408			else
409			flast->next = newf;
410			flast = newf;
411			return(MG_OK); /* we'll actually put it out later */
412	greg	1.1	}
413
414
415			int
416			i_sph(ac, av) /* translate sphere description */
417			int ac;
418			char **av;
419			{
420			register C_VERTEX *cent;
421
422			if (ac != 3)
423			return(MG_EARGC);
424	greg	1.3	flush_cache(); /* flush vertex cache */
425	greg	1.1	printf("%sSeparator {\n", tabs);
426			indent(1);
427			/* put out current material */
428			if (put_material() < 0)
429			return(MG_EBADMAT);
430			/* get center */
431			if ((cent = c_getvert(av[1])) == NULL)
432			return(MG_EUNDEF);
433			/* get radius */
434			if (!isflt(av[2]))
435			return(MG_ETYPE);
436			printf("%sTranslation { translation %13.9g %13.9g %13.9g }\n", tabs,
437			cent->p[0], cent->p[1], cent->p[2]);
438			printf("%sSphere { radius %s }\n", tabs, av[2]);
439			indent(0);
440			printf("%s}\n", tabs);
441			return(MG_OK);
442			}
443
444
445			int
446			i_cyl(ac, av) /* translate a cylinder description */
447			int ac;
448			char **av;
449			{
450			register C_VERTEX v1, v2;
451			FVECT va;
452			double length, angle;
453
454			if (ac != 4)
455			return(MG_EARGC);
456	greg	1.3	flush_cache(); /* flush vertex cache */
457	greg	1.1	printf("%sSeparator {\n", tabs);
458			indent(1);
459			/* put out current material */
460			if (put_material() < 0)
461			return(MG_EBADMAT);
462			/* get endpoints */
463	schorsch	1.9	if (((v1 = c_getvert(av[1])) == NULL) \| ((v2 = c_getvert(av[3])) == NULL))
464	greg	1.1	return(MG_EUNDEF);
465			/* get radius */
466			if (!isflt(av[2]))
467			return(MG_ETYPE);
468			/* compute transform */
469			va[0] = v2->p[0] - v1->p[0];
470			va[1] = v2->p[1] - v1->p[1];
471			va[2] = v2->p[2] - v1->p[2];
472			length = sqrt(DOT(va,va));
473	greg	1.5	if (va[1] >= length)
474			angle = 0.;
475			else if (va[1] <= -length)
476			angle = PI;
477			else
478			angle = acos(va[1]/length);
479	greg	1.1	printf("%sTranslation { translation %13.9g %13.9g %13.9g }\n", tabs,
480			.5(v1->p[0]+v2->p[0]), .5(v1->p[1]+v2->p[1]),
481			.5*(v1->p[2]+v2->p[2]));
482	greg	1.2	printf("%sRotation { rotation %.9g %.9g %.9g %.9g }\n", tabs,
483			va[2], 0., -va[0], angle);
484	greg	1.1	/* open-ended */
485			printf("%sCylinder { parts SIDES height %13.9g radius %s }\n", tabs,
486			length, av[2]);
487			indent(0);
488			printf("%s}\n", tabs);
489			return(MG_OK);
490	greg	1.2	}
491
492
493			char *
494			to_id(name) /* make sure a name is a valid Inventor ID */
495			register char *name;
496			{
497	greg	1.3	static char id[MAXID];
498	greg	1.2	register char *cp;
499
500	greg	1.3	for (cp = id; *name && cp < MAXID-1+id; name++)
501	greg	1.2	if (isalnum(name) \|\| name == '_')
502			cp++ = name;
503			else
504			*cp++ = '_';
505			*cp = '\0';
506			return(id);
507	greg	1.3	}
508
509
510			flush_cache() /* put out cached faces */
511			{
512			int donorms = 0;
513			register struct face *f;
514			register int i;
515
516			if (nverts == 0)
517			return;
518			/* put out coordinates */
519			printf("%sCoordinate3 {\n", tabs);
520			indent(1);
521			vlist[0][VFSEPPOS] = '\0';
522			printf("%spoint [ %s", tabs, vlist[0]);
523			for (i = 1; i < nverts; i++) {
524			vlist[i][VFSEPPOS] = '\0';
525			printf(",\n%s %s", tabs, vlist[i]);
526			if (strcmp(VFSEPPOS+1+vlist[i], VZVECT))
527			donorms++;
528			}
529			indent(0);
530			printf(" ]\n%s}\n", tabs);
531			if (donorms) { /* put out normals */
532			printf("%sNormal {\n", tabs);
533			indent(1);
534			printf("%svector [ %s", tabs, VFSEPPOS+1+vlist[0]);
535			for (i = 1; i < nverts; i++)
536			printf(",\n%s %s", tabs, VFSEPPOS+1+vlist[i]);
537			indent(0);
538			printf(" ]\n%s}\n", tabs);
539			}
540			/* put out faces */
541			printf("%sIndexedFaceSet {\n", tabs);
542			indent(1);
543			f = flist; /* coordinate indices */
544			printf("%scoordIndex [ %d", tabs, f->vl[0]);
545			for (i = 1; i < f->nv; i++)
546			printf(", %d", f->vl[i]);
547			for (f = f->next; f != NULL; f = f->next) {
548			printf(", -1,\n%s %d", tabs, f->vl[0]);
549			for (i = 1; i < f->nv; i++)
550			printf(", %d", f->vl[i]);
551			}
552			printf(" ]\n");
553			if (donorms) {
554			f = flist; /* normal indices */
555			printf("%snormalIndex [ %d", tabs, f->vl[0]);
556			for (i = 1; i < f->nv; i++)
557			printf(", %d", f->vl[i]);
558			for (f = f->next; f != NULL; f = f->next) {
559			printf(", -1,\n%s %d", tabs, f->vl[0]);
560			for (i = 1; i < f->nv; i++)
561			printf(", %d", f->vl[i]);
562			}
563			printf(" ]\n");
564			}
565			indent(0); /* close IndexedFaceSet */
566			printf("%s}\n", tabs);
567			indent(0); /* close face group */
568			printf("%s}\n", tabs);
569			while ((f = flist) != NULL) { /* free face list */
570			flist = f->next;
571			freeface(f);
572			}
573			lu_done(&vert_tab); /* clear lookup table */
574			nverts = 0;
575	greg	1.1	}