src/gen/xform.c

#ifndef lint
static const char RCSid[] = "$Id: xform.c,v 2.50 2016/03/04 00:21:21 greg Exp $";
#endif
/*
 *  xform.c - program to transform object files.
 *              Transformations must preserve aspect ratio.
 *
 *     10/19/85
 *     11/6/86          Finally added error checking!
 */

#include  <ctype.h>

#include  "platform.h"
#include  "paths.h"
#include  "rtio.h"
#include  "rtmath.h"
#include  "object.h"
#include  "otypes.h"

int  xac;                               /* global xform argument count */
char  **xav;                            /* global xform argument pointer */
int  xfa;                               /* start of xf arguments */

XF  tot;                                /* total transformation */
int  reverse;                           /* boolean true if scene mirrored */

int  invert;                            /* boolean true to invert surfaces */

int  expand;                            /* boolean true to expand commands */

char  *newmod;                          /* new modifier for surfaces */

char  *idprefix;                        /* prefix for object identifiers */

FUN  ofun[NUMOTYPE] = INIT_OTYPE;       /* default types and actions */

short  tinvers[NUMOTYPE];               /* inverse types for surfaces */

int  nrept = 1;                         /* number of array repetitions */

int stdinused = 0;                      /* stdin has been used by -f option? */

char  mainfn[PATH_MAX];                 /* main file name */
FILE  *mainfp = NULL;                   /* main file pointer */

#define  progname  (xav[0])

static int doargf(int ac, char **av, int fi);
static int doarray(int ac, char **av, int ai);
static void xform(char *name, FILE *fin);
static void xfcomm(char *fname, FILE *fin);
static void xfobject(char *fname, FILE *fin);
static int addxform(FILE *fin);
static int alias(FILE *fin);
void initotypes(void); /* XXX conflict with otypes.h */
static void openmain(char *iname);


int
main(           /* get transform options and transform file */
        int  argc,
        char  *argv[]
)
{
        int  mal_prefix = 0;
        int  a;
                                        /* check for argument list file */
        for (a = 1; a < argc; a++)
                if (!strcmp(argv[a], "-f"))
                        return(doargf(argc, argv, a));
                                        /* check for regular array */
        for (a = 1; a < argc; a++)
                if (!strcmp(argv[a], "-a"))
                        return(doarray(argc, argv, a));

        initotypes();                   /* initialize */
        invert = 0;
        expand = 1;
        newmod = NULL;
        idprefix = NULL;

        for (a = 1; a < argc; a++) {
                if (argv[a][0] == '-')
                        switch (argv[a][1]) {
                        case 'm':
                                if (argv[a][2] | (a+1 >= argc))
                                        break;
                                a++;
                                if (newmod == NULL)
                                        newmod = argv[a];
                                continue;
                        case 'n':
                                if (argv[a][2] | (a+1 >= argc))
                                        break;
                                a++;
                                if (idprefix == NULL)
                                        idprefix = argv[a];
                                else {
                                        char    *newp;
                                        newp = (char *)malloc(strlen(idprefix)+
                                                        strlen(argv[a])+2);
                                        if (newp == NULL)
                                                exit(2);
                                        sprintf(newp, "%s.%s",
                                                        idprefix, argv[a]);
                                        if (mal_prefix++)
                                                free(idprefix);
                                        idprefix = newp;
                                }
                                continue;
                        case 'c':
                                if (argv[a][2])
                                        break;
                                expand = 0;
                                continue;
                        case 'e':
                                if (argv[a][2])
                                        break;
                                expand = 1;
                                continue;
                        case 'I':
                                if (argv[a][2])
                                        break;
                                invert = !invert;
                                continue;
                        }
                break;
        }

        xav = argv;
        xfa = a;

        a += xf(&tot, argc-a, argv+a);

        if ( (reverse = (tot.sca < 0.0)) )
                tot.sca = -tot.sca;
        if (invert)
                reverse = !reverse;

        if (a < argc && argv[a][0] == '-') {
                fprintf(stderr, "%s: command line error at '%s'\n",
                                argv[0], argv[a]);
                exit(1);
        }

        xac = a;
                                        /* simple header */
        putchar('#');
        for (a = 0; a < xac; a++)
                printf(" %s", xav[a]);
        putchar('\n');
                                        /* transform input */
        if (xac == argc) {
                if (stdinused) {
                        fprintf(stderr, "%s: cannot use stdin more than once\n",
                                        argv[0]);
                        exit(1);
                }
                openmain(NULL);
                xform(mainfn, mainfp);
        } else
                for (a = xac; a < argc; a++) {
                        openmain(argv[a]);
                        xform(mainfn, mainfp);
                }

        if (mal_prefix)
                free(idprefix);
        return(0);
}


int
doargf(                 /* take argument list from file */
        int ac,
        char  **av,
        int fi
)
{
        int  inquote;
        char  *newav[256], **avp;
        char  argbuf[2048];
        char  *newid, newidbuf[128];
        char  *oldid;
        char    *cp;
        FILE    *argfp;
        int  n, i, k, newac, err;
        
        if (fi >= ac-1 || (av[fi+1][0] == '-' && av[fi+1][1] != '\0')) {
                fprintf(stderr, "%s: missing file for -f option\n", av[0]);
                exit(1);
        }
        if (av[fi+1][0] == '-' && av[fi+1][1] == '\0') {
                if (stdinused++) {
                        fprintf(stderr,
                                "%s: cannot use stdin more than once\n",
                                        av[0]);
                        exit(1);
                }
                argfp = stdin;
                n = 100;                /* we just don't know! */
        } else {
                if ((argfp = fopen(av[fi+1], "r")) == NULL) {
                        fprintf(stderr,
                                "%s: cannot open argument file \"%s\"\n",
                                        av[0], av[fi+1]);
                        exit(1);
                }
                n = 0;                  /* count number of lines in file */
                while (fgetline(argbuf,sizeof(argbuf),argfp) != NULL)
                        n += (argbuf[0] != '\0') & (argbuf[0] != '#');
                if (!n) {
                        fprintf(stderr, "%s: empty argument file \"%s\"\n",
                                        av[0], av[fi+1]);
                        exit(1);
                }
                rewind(argfp);
        }
        nrept *= n;
        err = 0; k = 0;                 /* read each arg list and call main */
        while (fgetline(argbuf,sizeof(argbuf),argfp) != NULL) {
                if (!argbuf[0] | (argbuf[0] == '#')) {
                        printf("%s\n", argbuf);
                        continue;
                }
                avp = newav+2;
                avp[0] = av[0];
                for (i = 1; i < fi; i++)
                        avp[i] = av[i];
                newac = i;
                cp = argbuf;            /* parse new words */
                if (*cp == '!') cp++;
                if (!strncmp(cp, "xform ", 6)) cp += 6;
                inquote = 0;
                for ( ; ; ) {
                skipspaces:
                        while (isspace(*cp))    /* nullify spaces */
                                *cp++ = '\0';
                        if ((*cp == '"') | (*cp == '\''))
                                inquote = *cp++;
                        if (!*cp)               /* all done? */
                                break;
                        if (cp[0] == '\\' && cp[1])
                                if (*++cp == '\n')
                                        goto skipspaces;
                        avp[newac++] = cp;      /* add argument to list */
                        if (inquote) {
                                while (*++cp)
                                        if (*cp == inquote) {
                                                *cp++ = '\0';
                                                break;
                                        }
                        } else {
                                while (*++cp && !isspace(*cp))
                                        ;
                        }
                }
                for (i = fi+2; i < ac; i++)
                        avp[newac++] = av[i];
                avp[newac] = NULL;
                newid = newidbuf;
                oldid = NULL;
                for (i = 2; i < newac; i++)
                        if (!strcmp(avp[i-1], "-n")) {
                                oldid = avp[i];
                                if (strlen(oldid)+32 > sizeof(newidbuf)) {
                                        newid = (char *)malloc(strlen(oldid)+32);
                                        if (newid == NULL)
                                                exit(2);
                                }
                                avp[i] = newid;
                                break;
                        }
                if (oldid == NULL) {
                        newav[0] = av[0];
                        newav[1] = "-n";
                        newav[2] = newid;
                        avp = newav;
                        newac += 2;
                }
                if (oldid == NULL)
                        sprintf(newid, "i%d", k);
                else
                        sprintf(newid, "%s.%d", oldid, k);
                err |= main(newac, avp);
                if (newid != newidbuf)
                        free(newid);
                k++;
        }
        fclose(argfp);
        return(err);
}


int
doarray(                        /* make array */
        int  ac,
        char  **av,
        int  ai
)
{
        char  *newav[256], **avp;
        char  newid[128], repts[32];
        char  *oldid = NULL;
        int  n, i, err;
        
        if (ai >= ac-1 || (n = atoi(av[ai+1])) <= 0) {
                fprintf(stderr, "%s: missing count for -a option\n", av[0]);
                exit(1);
        }
        nrept *= n;
        avp = newav+2;
        avp[0] = av[0];
        for (i = 1; i < ac; i++)
                if (!strcmp(av[i-1], "-n")) {
                        oldid = av[i];
                        avp[i] = newid;
                } else
                        avp[i] = av[i];
        avp[ai] = "-i";
        avp[ai+1] = repts;
        avp[i] = NULL;
        if (oldid == NULL) {
                newav[0] = av[0];
                newav[1] = "-n";
                newav[2] = newid;
                avp = newav;
                ac += 2;
        }
        err = 0;
        for (i = 0; i < n; i++) {
                if (oldid == NULL)
                        sprintf(newid, "a%d", i);
                else
                        sprintf(newid, "%s.%d", oldid, i);
                sprintf(repts, "%d", i);
                err |= main(ac, avp);
        }
        return(err);
}


void
xform(                  /* transform stream by tot.xfm */
        char  *name,
        FILE  *fin
)
{
        int  nobjs = 0;
        int  c;

        while ((c = getc(fin)) != EOF) {
                if (isspace(c))                         /* blank */
                        continue;
                if (c == '#') {                         /* comment */
                        putchar(c);
                        do {
                                if ((c = getc(fin)) == EOF)
                                        return;
                                putchar(c);
                        } while (c != '\n');
                } else if (c == '!') {                  /* command */
                        ungetc(c, fin);
                        xfcomm(name, fin);
                        nobjs++;
                } else {                                /* object */
                        ungetc(c, fin);
                        xfobject(name, fin);
                        nobjs++;
                }
        }
        if (nobjs == 0)
                fprintf(stderr, "%s: (%s): warning - empty file\n",
                                progname, name);
}


void
xfcomm(                 /* transform a command */
        char  *fname,
        FILE  *fin
)
{
        FILE  *pin;
        char  buf[2048];
        int  i;

        fgetline(buf, sizeof(buf), fin);
        if (expand) {
                if ((pin = popen(buf+1, "r")) == NULL) {
                        fprintf(stderr, "%s: (%s): cannot execute \"%s\"\n",
                                        progname, fname, buf);
                        exit(1);
                }
                xform(buf, pin);
                pclose(pin);
        } else {
                printf("\n%s", buf);
                if (xac > 1) {
                        printf(" | %s -e", xav[0]);     /* expand next time */
                        for (i = 1; i < xac; i++)
                                if (i >= xfa || strcmp(xav[i], "-c"))
                                        printf(" %s", xav[i]);
                }
                putchar('\n');
        }
}


void
xfobject(                               /* transform an object */
        char  *fname,
        FILE  *fin
)
{
        char  typ[16], nam[MAXSTR];
        int  fn;
                                                /* modifier and type */
        strcpy(typ, "EOF");
        fgetword(nam, sizeof(nam), fin);
        fgetword(typ, sizeof(typ), fin);
        if ((fn = otype(typ)) < 0) {
                fprintf(stderr, "%s: (%s): unknown object type \"%s\"\n",
                                progname, fname, typ);
                exit(1);
        }
        putchar('\n');
        if (ismodifier(fn)) {
                fputword(nam, stdout);
                printf(" %s ", typ);
        } else {
                fputword(newmod != NULL ? newmod : nam, stdout);
                printf(" %s ", invert ? ofun[tinvers[fn]].funame : typ);
        }
                                                /* object name */
        fgetword(nam, sizeof(nam), fin);
        if (idprefix == NULL || ismodifier(fn))
                fputword(nam, stdout);
        else {
                char    nnam[MAXSTR];
                sprintf(nnam, "%s.%s", idprefix, nam);
                fputword(nnam, stdout);
        }
        putchar('\n');
                                                /* transform arguments */
        if ((*ofun[fn].funp)(fin) < 0) {
                fprintf(stderr, "%s: (%s): bad %s \"%s\"\n",
                                progname, fname, ofun[fn].funame, nam);
                exit(1);
        }
}


int
o_default(                      /* pass on arguments unchanged */
        FILE  *fin
)
{
        int  i;
        FUNARGS  fa;

        if (readfargs(&fa, fin) != 1)
                return(-1);
                                        /* string arguments */
        printf("%d", fa.nsargs);
        for (i = 0; i < fa.nsargs; i++) {
                fputc(' ', stdout);
                fputword(fa.sarg[i], stdout);
        }
        printf("\n");
#ifdef  IARGS
                                        /* integer arguments */
        printf("%d", fa.niargs);
        for (i = 0; i < fa.niargs; i++)
                printf(" %d", fa.iarg[i]);
        printf("\n");
#else
        printf("0\n");
#endif
                                        /* float arguments */
        printf("%d", fa.nfargs);
        for (i = 0; i < fa.nfargs; i++)
                printf(" %18.12g", fa.farg[i]);
        printf("\n");
        freefargs(&fa);
        return(0);
}


int
addxform(                       /* add xf arguments to strings */
        FILE  *fin
)
{
        int  i;
        int  resetarr = 0;
        FUNARGS  fa;

        if (readfargs(&fa, fin) != 1)
                return(-1);
                                        /* string arguments */
        if (xac > xfa && strcmp(xav[xfa], "-i"))
                resetarr = 2;
        printf("%d", fa.nsargs + resetarr + xac-xfa);
        for (i = 0; i < fa.nsargs; i++) {
                fputc(' ', stdout);
                fputword(fa.sarg[i], stdout);
        }
        if (resetarr)
                printf(" -i 1");
        for (i = xfa; i < xac; i++)     /* add xf arguments */
                printf(" %s", xav[i]);
        printf("\n");
#ifdef  IARGS
                                        /* integer arguments */
        printf("%d", fa.niargs);
        for (i = 0; i < fa.niargs; i++)
                printf(" %d", fa.iarg[i]);
        printf("\n");
#else
        printf("0\n");
#endif
                                        /* float arguments */
        printf("%d", fa.nfargs);
        for (i = 0; i < fa.nfargs; i++)
                printf(" %18.12g", fa.farg[i]);
        printf("\n");
        freefargs(&fa);
        return(0);
}


int
alias(                  /* transfer alias */
        FILE  *fin
)
{
        char  aliasnm[MAXSTR];

        if (fgetword(aliasnm, MAXSTR, fin) == NULL)
                return(-1);
        printf("\t%s\n", aliasnm);
        return(0);
}


int
m_glow(                 /* transform arguments for proximity light */
        FILE  *fin
)
{
        FUNARGS  fa;

        if (readfargs(&fa, fin) != 1)
                return(-1);
        if ((fa.nsargs != 0) | (fa.nfargs != 4))
                return(-1);
        printf("0\n0\n4");
        printf(" %18.12g %18.12g %18.12g",
                        fa.farg[0], fa.farg[1], fa.farg[2]);
        printf(" %18.12g\n", fa.farg[3] * tot.sca);
        freefargs(&fa);
        return(0);
}


int
m_spot(                 /* transform arguments for spotlight */
        FILE  *fin
)
{
        FVECT  v;
        FUNARGS  fa;

        if (readfargs(&fa, fin) != 1)
                return(-1);
        if ((fa.nsargs != 0) | ( fa.nfargs != 7))
                return(-1);
        printf("0\n0\n7");
        printf(" %18.12g %18.12g %18.12g %18.12g\n",
                        fa.farg[0], fa.farg[1], fa.farg[2], fa.farg[3]);
        multv3(v, fa.farg+4, tot.xfm);
        printf("\t%18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
        freefargs(&fa);
        return(0);
}


int
m_mist(         /* transform arguments for mist */
        FILE  *fin
)
{
        FUNARGS  fa;
        int     i;

        if (readfargs(&fa, fin) != 1)
                return(-1);
        if (fa.nfargs > 7)
                return(-1);
        printf("%d", fa.nsargs);
        if (idprefix == NULL)
                for (i = 0; i < fa.nsargs; i++)
                        printf(" %s", fa.sarg[i]);
        else
                for (i = 0; i < fa.nsargs; i++) {
                        char    sname[256], *sp;
                        char    *cp1, *cp2 = sname;
                                                        /* add idprefix */
                        for (sp = fa.sarg[i]; *sp; sp = cp1) {
                                for (cp1 = idprefix; *cp1; )
                                        *cp2++ = *cp1++;
                                *cp2++ = '.';
                                for (cp1 = sp; *cp1 &&
                                                (*cp2++ = *cp1++) != '>'; )
                                        ;
                        }
                        *cp2 = '\0';
                        printf(" %s", sname);
                }
        printf("\n0\n%d", fa.nfargs);
        if (fa.nfargs > 2)
                printf(" %12.6g %12.6g %12.6g", fa.farg[0]/tot.sca,
                                fa.farg[1]/tot.sca, fa.farg[2]/tot.sca);
        for (i = 3; i < fa.nfargs; i++)
                printf(" %12.6g", fa.farg[i]);
        printf("\n");
        freefargs(&fa);
        return(0);
}


int
m_dielectric(           /* transform arguments for dielectric */
        FILE  *fin
)
{
        FUNARGS  fa;

        if (readfargs(&fa, fin) != 1)
                return(-1);
        if (fa.nsargs != 0  || fa.nfargs != 5)
                return(-1);
        printf("0\n0\n5");
        printf(" %12.6g %12.6g %12.6g",
                pow(fa.farg[0], 1.0/tot.sca),
                pow(fa.farg[1], 1.0/tot.sca),
                pow(fa.farg[2], 1.0/tot.sca));
        printf(" %12.6g %12.6g\n", fa.farg[3], fa.farg[4]);
        freefargs(&fa);
        return(0);
}


int
m_interface(            /* transform arguments for interface */
        FILE  *fin
)
{
        FUNARGS  fa;

        if (readfargs(&fa, fin) != 1)
                return(-1);
        if (fa.nsargs != 0  || fa.nfargs != 8)
                return(-1);
        printf("0\n0\n8\n");
        printf("%12.6g %12.6g %12.6g",
                pow(fa.farg[0], 1.0/tot.sca),
                pow(fa.farg[1], 1.0/tot.sca),
                pow(fa.farg[2], 1.0/tot.sca));
        printf(" %12.6g\n", fa.farg[3]);
        printf("%12.6g %12.6g %12.6g",
                pow(fa.farg[4], 1.0/tot.sca),
                pow(fa.farg[5], 1.0/tot.sca),
                pow(fa.farg[6], 1.0/tot.sca));
        printf(" %12.6g\n", fa.farg[7]);
        freefargs(&fa);
        return(0);
}


int
text(                   /* transform text arguments */
        FILE  *fin
)
{
        int  i;
        FVECT  v;
        FUNARGS  fa;

        if (readfargs(&fa, fin) != 1)
                return(-1);
        if (fa.nfargs < 9)
                return(-1);
                                        /* string arguments */
        printf("%d", fa.nsargs);
        for (i = 0; i < fa.nsargs; i++) {
                fputc(' ', stdout);
                fputword(fa.sarg[i], stdout);
        }
        printf("\n0\n%d\n", fa.nfargs);
                                        /* anchor point */
        multp3(v, fa.farg, tot.xfm);
        printf(" %18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
                                        /* right vector */
        multv3(v, fa.farg+3, tot.xfm);
        printf(" %18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
                                        /* down vector */
        multv3(v, fa.farg+6, tot.xfm);
        printf(" %18.12g %18.12g %18.12g", v[0], v[1], v[2]);
                                        /* remaining arguments */
        for (i = 9; i < fa.nfargs; i++) {
                if (i%3 == 0)
                        putchar('\n');
                printf(" %18.12g", fa.farg[i]);
        }
        putchar('\n');
        freefargs(&fa);
        return(0);
}


int
o_source(                       /* transform source arguments */
        FILE  *fin
)
{
        FVECT  dv;
        FUNARGS  fa;

        if (readfargs(&fa, fin) != 1)
                return(-1);
        if (fa.nsargs != 0  || fa.nfargs != 4)
                return(-1);
                                        /* transform direction vector */
        multv3(dv, fa.farg, tot.xfm);
                                        /* output */
        printf("0\n0\n4");
        printf(" %18.12g %18.12g %18.12g %18.12g\n",
                        dv[0], dv[1], dv[2], fa.farg[3]);
        freefargs(&fa);
        return(0);
}


int
o_sphere(                       /* transform sphere arguments */
        FILE  *fin
)
{
        FVECT  cent;
        double  rad;
        FUNARGS  fa;

        if (readfargs(&fa, fin) != 1)
                return(-1);
        if (fa.nsargs != 0  || fa.nfargs != 4)
                return(-1);
        
        multp3(cent, fa.farg, tot.xfm); /* transform center */
        
        rad = fa.farg[3] * tot.sca;             /* scale radius */
        
        printf("0\n0\n4");
        printf(" %18.12g %18.12g %18.12g %18.12g\n",
                        cent[0], cent[1], cent[2], rad);
        freefargs(&fa);
        return(0);
}


int
o_face(                 /* transform face arguments */
        FILE  *fin
)
{
        FVECT  p;
        int  i;
        FUNARGS  fa;

        if (readfargs(&fa, fin) != 1)
                return(-1);
        if (fa.nsargs != 0  || fa.nfargs % 3)
                return(-1);
        
        printf("0\n0\n%d\n", fa.nfargs);
        
        for (i = 0; i < fa.nfargs; i += 3) {
                if (reverse)
                        multp3(p, fa.farg+(fa.nfargs-i-3), tot.xfm);
                else
                        multp3(p, fa.farg+i, tot.xfm);
                printf(" %18.12g %18.12g %18.12g\n", p[0], p[1], p[2]);
        }
        freefargs(&fa);
        return(0);
}


int
o_cone(                 /* transform cone and cup arguments */
        FILE  *fin
)
{
        FVECT  p0, p1;
        double  r0, r1;
        FUNARGS  fa;

        if (readfargs(&fa, fin) != 1)
                return(-1);
        if (fa.nsargs != 0  || fa.nfargs != 8)
                return(-1);

        printf("0\n0\n8\n");

        multp3(p0, fa.farg, tot.xfm);
        multp3(p1, fa.farg+3, tot.xfm);
        r0 = fa.farg[6] * tot.sca;
        r1 = fa.farg[7] * tot.sca;
        printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
        printf(" %18.12g %18.12g %18.12g\n", p1[0], p1[1], p1[2]);
        printf(" %18.12g %18.12g\n", r0, r1);

        freefargs(&fa);
        return(0);
}


int
o_cylinder(                     /* transform cylinder and tube arguments */
        FILE  *fin
)
{
        FVECT  p0, p1;
        double  rad;
        FUNARGS  fa;

        if (readfargs(&fa, fin) != 1)
                return(-1);
        if (fa.nsargs != 0  || fa.nfargs != 7)
                return(-1);

        printf("0\n0\n7\n");

        multp3(p0, fa.farg, tot.xfm);
        multp3(p1, fa.farg+3, tot.xfm);
        rad = fa.farg[6] * tot.sca;
        printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
        printf(" %18.12g %18.12g %18.12g\n", p1[0], p1[1], p1[2]);
        printf(" %18.12g\n", rad);
        freefargs(&fa);
        return(0);
}


int
o_ring(                 /* transform ring arguments */
        FILE  *fin
)
{
        FVECT  p0, pd;
        double  r0, r1;
        FUNARGS  fa;

        if (readfargs(&fa, fin) != 1)
                return(-1);
        if (fa.nsargs != 0  || fa.nfargs != 8)
                return(-1);

        printf("0\n0\n8\n");

        multp3(p0, fa.farg, tot.xfm);
        multv3(pd, fa.farg+3, tot.xfm);
        if (invert) {
                pd[0] = -pd[0];
                pd[1] = -pd[1];
                pd[2] = -pd[2];
        }
        r0 = fa.farg[6] * tot.sca;
        r1 = fa.farg[7] * tot.sca;
        printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
        printf(" %18.12g %18.12g %18.12g\n", pd[0], pd[1], pd[2]);
        printf(" %18.12g %18.12g\n", r0, r1);
        freefargs(&fa);
        return(0);
}


void
initotypes(void)                        /* initialize ofun[] array */
{
        int  i;

        if (ofun[OBJ_SOURCE].funp == o_source)
                return;                 /* done already */
                                        /* functions get new transform */
        for (i = 0; i < NUMOTYPE; i++)
                if (hasfunc(i))
                        ofun[i].funp = addxform;
                                        /* special cases */
        ofun[OBJ_SOURCE].funp = o_source;
        ofun[OBJ_SPHERE].funp =
        ofun[OBJ_BUBBLE].funp = o_sphere;
        ofun[OBJ_FACE].funp = o_face;
        ofun[OBJ_CONE].funp =
        ofun[OBJ_CUP].funp = o_cone;
        ofun[OBJ_CYLINDER].funp =
        ofun[OBJ_TUBE].funp = o_cylinder;
        ofun[OBJ_RING].funp = o_ring;
        ofun[OBJ_INSTANCE].funp =
        ofun[OBJ_MESH].funp = addxform;
        ofun[MAT_GLOW].funp = m_glow;
        ofun[MAT_SPOT].funp = m_spot;
        ofun[MAT_DIELECTRIC].funp = m_dielectric;
        ofun[MAT_INTERFACE].funp = m_interface;
        ofun[MAT_MIST].funp = m_mist;
        ofun[PAT_CTEXT].funp =
        ofun[PAT_BTEXT].funp =
        ofun[MIX_TEXT].funp = text;
        ofun[MOD_ALIAS].funp = alias;
                                        /* surface inverses */
        tinvers[OBJ_FACE] = OBJ_FACE;
        tinvers[OBJ_SOURCE] = OBJ_SOURCE;
        tinvers[OBJ_CONE] = OBJ_CUP;
        tinvers[OBJ_CUP] = OBJ_CONE;
        tinvers[OBJ_SPHERE] = OBJ_BUBBLE;
        tinvers[OBJ_BUBBLE] = OBJ_SPHERE;
        tinvers[OBJ_RING] = OBJ_RING;
        tinvers[OBJ_CYLINDER] = OBJ_TUBE;
        tinvers[OBJ_TUBE] = OBJ_CYLINDER;
        tinvers[OBJ_INSTANCE] = OBJ_INSTANCE;   /* oh, well */
        tinvers[OBJ_MESH] = OBJ_MESH;           /* ditto */
}


#ifdef  OLDXFORM
void
openmain(
        char  *fname
)
{
        if (fname == NULL) {
                strcpy(mainfn, "standard input");
                mainfp = stdin;
                return;
        }
        if (mainfp != NULL) {
                if (!strcmp(fname, mainfn)) {
                        rewind(mainfp);
                        return;
                }
                fclose(mainfp);
        }
        if ((mainfp = fopen(fname, "r")) == NULL) {
                fprintf(stderr, "%s: cannot open file \"%s\"\n",
                                progname, fname);
                exit(1);
        }
        strcpy(mainfn, fname);
}
#else
void
openmain(               /* open input, changing directory for file */
        char  *iname
)
{
        static char  origdir[PATH_MAX];
        static char  curfn[PATH_MAX];
        static int  diffdir;
        char  *fpath;

        if (iname == NULL) {                    /* standard input */
                if (mainfp == NULL) {
                        int  c;
                        strcpy(mainfn, "standard input");
                        if (nrept <= 1) {
                                mainfp = stdin;
                                return;                 /* just read once */
                        }
                                                        /* else copy */
                        if ((mainfp = tmpfile()) == NULL) {
                                fprintf(stderr,
                                        "%s: cannot create temporary file\n",
                                                progname);
                                exit(1);
                        }
                        while ((c = getc(stdin)) != EOF)
                                putc(c, mainfp);
                }
                rewind(mainfp);                 /* rewind copy */
                return;
        }
        if (mainfp == NULL) {                   /* first call, initialize */
                getcwd(origdir, sizeof(origdir));
        } else if (!strcmp(iname, curfn)) {     /* just need to rewind? */
                rewind(mainfp);
                return;
        } else {                                /* else close old stream */
                fclose(mainfp);
                mainfp = NULL;
                if (diffdir) {                  /* return to our directory */
                        chdir(origdir);
                        diffdir = 0;
                }
        }
        strcpy(curfn, iname);                   /* remember input name */
                                                /* get full path for file */
        if ((fpath = getpath(iname, getrlibpath(), R_OK)) == NULL) {
                fprintf(stderr, "%s: cannot find file \"%s\"\n",
                                progname, iname);
                exit(1);
        }
        if (fpath[0] == '.' && ISDIRSEP(fpath[1]))      /* remove leading ./ */
                fpath += 2;
                                                /* record path name */
        strcpy(mainfn, fpath);
        if (expand) {                           /* change to local directory */
                char  *cp = fpath + strlen(fpath);      /* get dir. */
                while (cp > fpath) {
                        cp--;
                        if (ISDIRSEP(*cp)) {
                                if (cp == fpath)
                                        cp++;   /* root special case */
                                break;
                        }
                }
                *cp = '\0';
                if (fpath[0]) {                 /* change to new directory? */
                        if (chdir(fpath) < 0) {
                                fprintf(stderr,
                                "%s: cannot change directory to \"%s\"\n",
                                                progname, fpath);
                                exit(1);
                        }
                        diffdir++;
                }
                                                /* get final path component */
                for (fpath = iname+strlen(iname);
                                fpath > iname && !ISDIRSEP(fpath[-1]); fpath--)
                        ;
        }
                                                /* finally, open the file */
        if ((mainfp = fopen(fpath, "r")) == NULL) {
                fprintf(stderr, "%s: cannot open file \"%s\"\n",
                                progname, mainfn);
                exit(1);
        }
}
#endif
Revision:	2.51
Committed:	Sun Dec 2 16:30:30 2018 UTC (5 years, 5 months ago) by greg
Content type:	text/plain
Branch:	MAIN
CVS Tags:	rad5R4, rad5R3, HEAD
Changes since 2.50:	+8 -6 lines
Log Message:	Echo empty lines and comments in output when using -f option