src/rt/rmain.c

/* Copyright (c) 1991 Regents of the University of California */

#ifndef lint
static char SCCSid[] = "$SunId$ LBL";
#endif

/*
 *  rmain.c - main for ray tracing programs
 *
 *     3/24/87
 */

                                        /* for flaky pre-processors */
#ifndef  RPICT
#define  RPICT          0
#endif
#ifndef  RTRACE
#define  RTRACE         0
#endif
#ifndef  RVIEW
#define  RVIEW          0
#endif

#include  "ray.h"

#include  "resolu.h"

#include  "octree.h"

#include  <signal.h>

#include  "view.h"

#ifndef  DEFPATH
#define  DEFPATH        ":/usr/local/lib/ray"
#endif
#ifndef  ULIBVAR
#define  ULIBVAR        "RAYPATH"
#endif

char  *progname;                        /* argv[0] */

char  *octname;                         /* octree name */

char  *libpath;                         /* library directory list */

char  *sigerr[NSIG];                    /* signal error messages */

extern char  VersionID[];               /* version ID string */

extern int  stderr_v();                 /* standard error output */
int  (*errvec)() = stderr_v;            /* error output vector */
int  (*wrnvec)() = stderr_v;            /* warning output vector */
int  (*cmdvec)() = NULL;                /* command error vector */

int  (*trace)() = NULL;                 /* trace call */
int  do_irrad = 0;                      /* compute irradiance? */

extern long  time();
long  tstart;                           /* start time */

extern int  ambnotify();                /* new object notify functions */
int  (*addobjnotify[])() = {ambnotify, NULL};

CUBE  thescene;                         /* our scene */

extern int  imm_irrad;                  /* calculate immediate irradiance? */

extern int  ralrm;                      /* seconds between reports */

extern int  greyscale;                  /* map colors to brightness? */
extern char  *devname;                  /* output device name */

extern int  inform;                     /* input format */
extern int  outform;                    /* output format */
extern char  *outvals;                  /* output values */

extern VIEW  ourview;                   /* viewing parameters */

extern int  hresolu;                    /* horizontal resolution */
extern int  vresolu;                    /* vertical resolution */
extern double  pixaspect;               /* pixel aspect ratio */

extern int  psample;                    /* pixel sample size */
extern double  maxdiff;                 /* max. sample difference */
extern double  dstrpix;                 /* square pixel distribution */

extern double  dstrsrc;                 /* square source distribution */
extern double  shadthresh;              /* shadow threshold */
extern double  shadcert;                /* shadow testing certainty */
extern int  directrelay;                /* number of source relays */
extern int  vspretest;                  /* virtual source pretest density */
extern int  directinvis;                /* light sources invisible to eye? */
extern double  srcsizerat;              /* maximum source size/dist. ratio */

extern double  specthresh;              /* specular sampling threshold */
extern double  specjitter;              /* specular sampling jitter */

extern int  maxdepth;                   /* maximum recursion depth */
extern double  minweight;               /* minimum ray weight */

extern COLOR  ambval;                   /* ambient value */
extern double  ambacc;                  /* ambient accuracy */
extern int  ambres;                     /* ambient resolution */
extern int  ambdiv;                     /* ambient divisions */
extern int  ambssamp;                   /* ambient super-samples */
extern int  ambounce;                   /* ambient bounces */
extern char  *amblist[];                /* ambient include/exclude list */
extern int  ambincl;                    /* include == 1, exclude == 0 */


main(argc, argv)
int  argc;
char  *argv[];
{
#define  check(ol,al)           if (argv[i][ol] || \
                                badarg(argc-i-1,argv+i+1,al)) \
                                goto badopt
#define  bool(olen,var)         switch (argv[i][olen]) { \
                                case '\0': var = !var; break; \
                                case 'y': case 'Y': case 't': case 'T': \
                                case '+': case '1': var = 1; break; \
                                case 'n': case 'N': case 'f': case 'F': \
                                case '-': case '0': var = 0; break; \
                                default: goto badopt; }
        extern char  *getenv();
        char  *err;
        char  *recover = NULL;
        char  *zfile = NULL;
        char  *errfile = NULL;
        char  *ambfile = NULL;
        char  **amblp = amblist;
        int  loadflags = ~IO_FILES;
        int  rval, gotvfile = 0;
        RESOLU  rs;
        int  i;
                                        /* record start time */
        tstart = time((long *)0);
                                        /* global program name */
        progname = argv[0];
                                        /* get library path */
        if ((libpath = getenv(ULIBVAR)) == NULL)
                libpath = DEFPATH;
                                        /* initialize object types */
        initotypes();
                                        /* initialize urand */
        initurand(2048);
                                        /* option city */
        for (i = 1; i < argc && argv[i][0] == '-'; i++) {
                if (!strcmp(argv[i], "-version")) {
                        puts(VersionID);
                        quit(0);
                }
                if (!strcmp(argv[i], "-defaults") ||
                                !strcmp(argv[i], "-help")) {
                        printdefaults();
                        quit(0);
                }
#if  RVIEW
                if (!strcmp(argv[i], "-devices")) {
                        printdevices();
                        quit(0);
                }
#endif
#if  RPICT|RVIEW
                rval = getviewopt(&ourview, argc-i, argv+i);
                if (rval >= 0) {
                        i += rval;
                        continue;
                }
#endif
                switch (argv[i][1]) {
#if  RPICT|RVIEW
                case 'v':                               /* view file */
                        if (argv[i][2] != 'f')
                                goto badopt;
                        check(3,"s");
                        rval = viewfile(argv[++i], &ourview, NULL);
                        if (rval < 0) {
                                sprintf(errmsg,
                                "cannot open view file \"%s\"",
                                                argv[i]);
                                error(SYSTEM, errmsg);
                        } else if (rval == 0) {
                                sprintf(errmsg,
                                        "bad view file \"%s\"",
                                                argv[i]);
                                error(USER, errmsg);
                        } else
                                gotvfile += rval;
                        break;
#endif
                case 'd':                               /* direct */
                        switch (argv[i][2]) {
                        case 't':                               /* threshold */
                                check(3,"f");
                                shadthresh = atof(argv[++i]);
                                break;
                        case 'c':                               /* certainty */
                                check(3,"f");
                                shadcert = atof(argv[++i]);
                                break;
                        case 'j':                               /* jitter */
                                check(3,"f");
                                dstrsrc = atof(argv[++i]);
                                break;
                        case 'r':                               /* relays */
                                check(3,"i");
                                directrelay = atoi(argv[++i]);
                                break;
                        case 'p':                               /* pretest */
                                check(3,"i");
                                vspretest = atoi(argv[++i]);
                                break;
                        case 'i':                               /* invis. */
                                bool(3,directinvis);
                                break;
                        case 's':                               /* size */
                                check(3,"f");
                                srcsizerat = atof(argv[++i]);
                                break;
                        default:
                                goto badopt;
                        }
                        break;
                case 's':                               /* specular */
                        switch (argv[i][2]) {
                        case 't':                               /* threshold */
                                check(3,"f");
                                specthresh = atof(argv[++i]);
                                break;
                        case 'j':                               /* jitter */
                                check(3,"f");
                                specjitter = atof(argv[++i]);
                                break;
                        default:
                                goto badopt;
                        }
                        break;
#if  RPICT|RVIEW
                case 'p':                               /* pixel */
                        switch (argv[i][2]) {
                        case 's':                               /* sample */
                                check(3,"i");
                                psample = atoi(argv[++i]);
                                break;
                        case 't':                               /* threshold */
                                check(3,"f");
                                maxdiff = atof(argv[++i]);
                                break;
#if  RPICT
                        case 'j':                               /* jitter */
                                check(3,"f");
                                dstrpix = atof(argv[++i]);
                                break;
                        case 'a':                               /* aspect */
                                check(3,"f");
                                pixaspect = atof(argv[++i]);
                                break;
#endif
                        default:
                                goto badopt;
                        }
                        break;
#endif
#if  RPICT|RTRACE
                case 'x':                               /* x resolution */
                        check(2,"i");
                        hresolu = atoi(argv[++i]);
                        break;
                case 'y':                               /* y resolution */
                        check(2,"i");
                        vresolu = atoi(argv[++i]);
                        break;
#endif
                case 'w':                               /* warnings */
                        rval = wrnvec != NULL;
                        bool(2,rval);
                        if (rval) wrnvec = stderr_v;
                        else wrnvec = NULL;
                        break;
                case 'e':                               /* error file */
                        check(2,"s");
                        errfile = argv[++i];
                        break;
                case 'l':                               /* limit */
                        switch (argv[i][2]) {
                        case 'r':                               /* recursion */
                                check(3,"i");
                                maxdepth = atoi(argv[++i]);
                                break;
                        case 'w':                               /* weight */
                                check(3,"f");
                                minweight = atof(argv[++i]);
                                break;
                        default:
                                goto badopt;
                        }
                        break;
                case 'i':                               /* irradiance */
                        bool(2,do_irrad);
                        break;
#if  RPICT
                case 'z':                               /* z file */
                        check(2,"s");
                        zfile = argv[++i];
                        break;
                case 'r':                               /* recover file */
                        check(2,"s");
                        recover = argv[++i];
                        rval = viewfile(recover, &ourview, &rs);
                        if (rval <= 0 || rs.or != PIXSTANDARD) {
                                sprintf(errmsg,
                        "cannot recover view parameters from \"%s\"", recover);
                                error(WARNING, errmsg);
                        } else {
                                gotvfile += rval;
                                pixaspect = 0.0;
                                hresolu = rs.xr;
                                vresolu = rs.yr;
                        }
                        break;
                case 't':                               /* timer */
                        check(2,"i");
                        ralrm = atoi(argv[++i]);
                        break;
#endif
                case 'a':                               /* ambient */
                        switch (argv[i][2]) {
                        case 'v':                               /* value */
                                check(3,"fff");
                                setcolor(ambval, atof(argv[i+1]),
                                                atof(argv[i+2]),
                                                atof(argv[i+3]));
                                i += 3;
                                break;
                        case 'a':                               /* accuracy */
                                check(3,"f");
                                ambacc = atof(argv[++i]);
                                break;
                        case 'r':                               /* resolution */
                                check(3,"i");
                                ambres = atoi(argv[++i]);
                                break;
                        case 'd':                               /* divisions */
                                check(3,"i");
                                ambdiv = atoi(argv[++i]);
                                break;
                        case 's':                               /* super-samp */
                                check(3,"i");
                                ambssamp = atoi(argv[++i]);
                                break;
                        case 'b':                               /* bounces */
                                check(3,"i");
                                ambounce = atoi(argv[++i]);
                                break;
                        case 'i':                               /* include */
                        case 'I':
                                check(3,"s");
                                if (ambincl != 1) {
                                        ambincl = 1;
                                        amblp = amblist;
                                }
                                if (argv[i][2] == 'I') {        /* file */
                                        rval = wordfile(amblp, argv[++i]);
                                        if (rval < 0) {
                                                sprintf(errmsg,
                                "cannot open ambient include file \"%s\"",
                                                                argv[i]);
                                                error(SYSTEM, errmsg);
                                        }
                                        amblp += rval;
                                } else
                                        *amblp++ = argv[++i];
                                break;
                        case 'e':                               /* exclude */
                        case 'E':
                                check(3,"s");
                                if (ambincl != 0) {
                                        ambincl = 0;
                                        amblp = amblist;
                                }
                                if (argv[i][2] == 'E') {        /* file */
                                        rval = wordfile(amblp, argv[++i]);
                                        if (rval < 0) {
                                                sprintf(errmsg,
                                "cannot open ambient exclude file \"%s\"",
                                                                argv[i]);
                                                error(SYSTEM, errmsg);
                                        }
                                        amblp += rval;
                                } else
                                        *amblp++ = argv[++i];
                                break;
                        case 'f':                               /* file */
                                check(3,"s");
                                ambfile= argv[++i];
                                break;
                        default:
                                goto badopt;
                        }
                        break;
#if  RTRACE
                case 'I':                               /* immed. irradiance */
                        bool(2,imm_irrad);
                        break;
                case 'f':                               /* format i/o */
                        switch (argv[i][2]) {
                        case 'a':                               /* ascii */
                        case 'f':                               /* float */
                        case 'd':                               /* double */
                                inform = argv[i][2];
                                break;
                        default:
                                goto badopt;
                        }
                        switch (argv[i][3]) {
                        case '\0':
                                outform = inform;
                                break;
                        case 'a':                               /* ascii */
                        case 'f':                               /* float */
                        case 'd':                               /* double */
                                check(4,"");
                                outform = argv[i][3];
                                break;
                        default:
                                goto badopt;
                        }
                        break;
                case 'o':                               /* output */
                        outvals = argv[i]+2;
                        break;
                case 'h':                               /* header output */
                        rval = loadflags & IO_INFO;
                        bool(2,rval);
                        loadflags = rval ? loadflags | IO_INFO :
                                        loadflags & ~IO_INFO;
                        break;
#endif
#if  RVIEW
                case 'b':                               /* black and white */
                        bool(2,greyscale);
                        break;
                case 'o':                               /* output device */
                        check(2,"s");
                        devname = argv[++i];
                        break;
#endif
                default:
                        goto badopt;
                }
        }
        *amblp = NULL;
#if  RPICT|RVIEW
        err = setview(&ourview);        /* set viewing parameters */
        if (err != NULL)
                error(USER, err);
#endif
#if  RPICT
        normaspect(viewaspect(&ourview), &pixaspect, &hresolu, &vresolu);
#endif
                                        /* set up signal handling */
        sigdie(SIGINT, "Interrupt");
        sigdie(SIGHUP, "Hangup");
        sigdie(SIGTERM, "Terminate");
        sigdie(SIGPIPE, "Broken pipe");
        sigdie(SIGALRM, "Alarm clock");
#ifdef  SIGXCPU
        sigdie(SIGXCPU, "CPU limit exceeded");
        sigdie(SIGXFSZ, "File size exceeded");
#endif
                                        /* open error file */
        if (errfile != NULL) {
                if (freopen(errfile, "a", stderr) == NULL)
                        quit(2);
                fprintf(stderr, "**************\n*** PID %5d: ",
                                getpid());
                printargs(argc, argv, stderr);
                fputs("\n", stderr);
                fflush(stderr);
        }
#ifdef  NICE
        nice(NICE);                     /* lower priority */
#endif
                                        /* get octree */
#if  RVIEW
        loadflags &= ~IO_INFO;
#endif
        if (i == argc)
                octname = NULL;
        else if (i == argc-1)
                octname = argv[i];
        else
                goto badopt;
#if  RVIEW|RTRACE
        if (octname == NULL)
                error(USER, "missing octree argument");
#endif
        readoct(octname, loadflags, &thescene, NULL);

        if (loadflags & IO_INFO) {      /* print header */
                printargs(i, argv, stdout);
                printf("SOFTWARE= %s\n", VersionID);
#if  RPICT
                if (gotvfile) {
                        printf(VIEWSTR);
                        fprintview(&ourview, stdout);
                        printf("\n");
                }
                if (pixaspect < .99 || pixaspect > 1.01)
                        fputaspect(pixaspect, stdout);
                fputformat(COLRFMT, stdout);
#endif
#if  RTRACE
                fputformat(     outform=='a' ? "ascii" :
                                outform=='f' ? "float" :
                                "double",       stdout  );
#endif
                printf("\n");
        }

        marksources();                  /* find and mark sources */

        setambient(ambfile);            /* initialize ambient calculation */

#if  RPICT
        render(zfile, recover);         /* render the scene */
#endif
#if  RTRACE
        rtrace(NULL);                   /* trace rays from stdin */
#endif
#if  RVIEW
        rview();                        /* go */
#endif
        quit(0);

badopt:
        sprintf(errmsg, "command line error at '%s'", argv[i]);
        error(USER, errmsg);

#undef  check
#undef  bool
}


eputs(s)                                /* error output */
char  *s;
{
        if (errvec != NULL)
                (*errvec)(s);
}


wputs(s)                                /* warning output */
char  *s;
{
        int  lasterrno = errno;         /* save errno */

        if (wrnvec != NULL)
                (*wrnvec)(s);

        errno = lasterrno;
}


cputs(s)                                /* command error output */
char  *s;
{
        if (cmdvec != NULL)
                (*cmdvec)(s);
}


stderr_v(s)                             /* put string to stderr */
register char  *s;
{
        static int  midline = 0;

        if (!*s)
                return;
        if (!midline++) {
                fputs(progname, stderr);
                fputs(": ", stderr);
        }
        fputs(s, stderr);
        if (s[strlen(s)-1] == '\n') {
                fflush(stderr);
                midline = 0;
        }
}


onsig(signo)                            /* fatal signal */
int  signo;
{
        static int  gotsig = 0;

        if (gotsig++)                   /* two signals and we're gone! */
                _exit(signo);

        alarm(15);                      /* allow 15 seconds to clean up */
        signal(SIGALRM, SIG_DFL);       /* make certain we do die */
        eputs("signal - ");
        eputs(sigerr[signo]);
        eputs("\n");
        quit(3);
}


sigdie(signo, msg)                      /* set fatal signal */
int  signo;
char  *msg;
{
        if (signal(signo, onsig) == SIG_IGN)
                signal(signo, SIG_IGN);
        sigerr[signo] = msg;
}


printdefaults()                 /* print default values to stdout */
{
        register char  *cp;

#if  RTRACE
        if (imm_irrad)
                printf("-I+\t\t\t\t# immediate irradiance on\n");
        else
#endif
        printf(do_irrad ? "-i+\t\t\t\t# irradiance calculation on\n" :
                        "-i-\t\t\t\t# irradiance calculation off\n");
#if  RVIEW
        printf(greyscale ? "-b+\t\t\t\t# greyscale on\n" :
                        "-b-\t\t\t\t# greyscale off\n");
#endif
#if  RPICT|RVIEW
        printf("-vt%c\t\t\t\t# view type %s\n", ourview.type,
                        ourview.type==VT_PER ? "perspective" :
                        ourview.type==VT_PAR ? "parallel" :
                        ourview.type==VT_HEM ? "hemispherical" :
                        ourview.type==VT_ANG ? "angular" :
                        "unknown");
        printf("-vp %f %f %f\t# view point\n",
                        ourview.vp[0], ourview.vp[1], ourview.vp[2]);
        printf("-vd %f %f %f\t# view direction\n", 
                        ourview.vdir[0], ourview.vdir[1], ourview.vdir[2]);
        printf("-vu %f %f %f\t# view up\n", 
                        ourview.vup[0], ourview.vup[1], ourview.vup[2]);
        printf("-vh %f\t\t\t# view horizontal size\n", ourview.horiz);
        printf("-vv %f\t\t\t# view vertical size\n", ourview.vert);
        printf("-vs %f\t\t\t# view shift\n", ourview.hoff);
        printf("-vl %f\t\t\t# view lift\n", ourview.voff);
#endif
#if  RPICT|RTRACE
        printf("-x  %-9d\t\t\t# x resolution\n", hresolu);
        printf("-y  %-9d\t\t\t# y resolution\n", vresolu);
#endif
#if  RPICT
        printf("-pa %f\t\t\t# pixel aspect ratio\n", pixaspect);
        printf("-pj %f\t\t\t# pixel jitter\n", dstrpix);
#endif
#if  RPICT|RVIEW
        printf("-ps %-9d\t\t\t# pixel sample\n", psample);
        printf("-pt %f\t\t\t# pixel threshold\n", maxdiff);
#endif
        printf("-dt %f\t\t\t# direct threshold\n", shadthresh);
        printf("-dc %f\t\t\t# direct certainty\n", shadcert);
        printf("-dj %f\t\t\t# direct jitter\n", dstrsrc);
        printf("-ds %f\t\t\t# direct sampling\n", srcsizerat);
        printf("-dr %-9d\t\t\t# direct relays\n", directrelay);
        printf("-dp %-9d\t\t\t# direct pretest density\n", vspretest);
        printf(directinvis ? "-di+\t\t\t\t# direct invisibility on\n" :
                        "-di-\t\t\t\t# direct invisibility off\n");
        printf("-sj %f\t\t\t# specular jitter\n", specjitter);
        printf("-st %f\t\t\t# specular threshold\n", specthresh);
        printf("-av %f %f %f\t# ambient value\n", colval(ambval,RED),
                        colval(ambval,GRN), colval(ambval, BLU));
        printf("-ab %-9d\t\t\t# ambient bounces\n", ambounce);
        printf("-aa %f\t\t\t# ambient accuracy\n", ambacc);
        printf("-ar %-9d\t\t\t# ambient resolution\n", ambres);
        printf("-ad %-9d\t\t\t# ambient divisions\n", ambdiv);
        printf("-as %-9d\t\t\t# ambient super-samples\n", ambssamp);
        printf("-lr %-9d\t\t\t# limit reflection\n", maxdepth);
        printf("-lw %f\t\t\t# limit weight\n", minweight);
#if  RPICT
        printf("-t  %-9d\t\t\t# time between reports\n", ralrm);
#endif
#if  RVIEW
        printf("-o %s\t\t\t\t# output device\n", devname);
#endif
#if  RTRACE
        printf("-f%c%c\t\t\t\t# format input/output = %s/%s\n",
                        inform, outform, 
                        inform=='a' ? "ascii" :
                        inform=='f' ? "float" : "double",
                        outform=='a' ? "ascii" :
                        outform=='f' ? "float" : "double");
        printf("-o%s\t\t\t\t# output", outvals);
        for (cp = outvals; *cp; cp++)
                switch (*cp) {
                case 't': printf(" trace"); break;
                case 'o': printf(" origin"); break;
                case 'd': printf(" direction"); break;
                case 'v': printf(" value"); break;
                case 'l': printf(" length"); break;
                case 'p': printf(" point"); break;
                case 'n': printf(" normal"); break;
                case 's': printf(" surface"); break;
                case 'w': printf(" weight"); break;
                case 'm': printf(" modifier"); break;
                }
        printf("\n");
#endif
        printf(wrnvec != NULL ? "-w+\t\t\t\t# warning messages on\n" :
                        "-w-\t\t\t\t# warning messages off\n");
}
Revision:	2.6
Committed:	Thu May 21 10:52:17 1992 UTC (33 years ago) by greg
Content type:	text/plain
Branch:	MAIN
Changes since 2.5:	+2 -1 lines
Log Message:	fixed complaint by IBM/RT
#	User	Rev	Content
1	greg	1.34	/* Copyright (c) 1991 Regents of the University of California */
2	greg	1.1
3			#ifndef lint
4			static char SCCSid[] = "$SunId$ LBL";
5			#endif
6
7			/*
8			* rmain.c - main for ray tracing programs
9			*
10			* 3/24/87
11			*/
12
13			/* for flaky pre-processors */
14			#ifndef RPICT
15			#define RPICT 0
16			#endif
17			#ifndef RTRACE
18			#define RTRACE 0
19			#endif
20			#ifndef RVIEW
21			#define RVIEW 0
22			#endif
23
24			#include "ray.h"
25
26	greg	1.50	#include "resolu.h"
27
28	greg	1.1	#include "octree.h"
29
30			#include <signal.h>
31
32			#include "view.h"
33
34			#ifndef DEFPATH
35			#define DEFPATH ":/usr/local/lib/ray"
36			#endif
37	greg	1.17	#ifndef ULIBVAR
38	greg	1.1	#define ULIBVAR "RAYPATH"
39	greg	1.17	#endif
40	greg	1.1
41			char progname; / argv[0] */
42
43	greg	1.19	char octname; / octree name */
44
45	greg	1.1	char libpath; / library directory list */
46
47			char sigerr[NSIG]; / signal error messages */
48
49	greg	1.46	extern char VersionID[]; /* version ID string */
50
51	greg	1.1	extern int stderr_v(); /* standard error output */
52			int (errvec)() = stderr_v; / error output vector */
53			int (wrnvec)() = stderr_v; / warning output vector */
54			int (cmdvec)() = NULL; / command error vector */
55
56			int (trace)() = NULL; / trace call */
57	greg	1.36	int do_irrad = 0; /* compute irradiance? */
58	greg	1.1
59	greg	1.34	extern long time();
60			long tstart; /* start time */
61
62	greg	1.24	extern int ambnotify(); /* new object notify functions */
63			int (*addobjnotify[])() = {ambnotify, NULL};
64
65	greg	1.1	CUBE thescene; /* our scene */
66
67	greg	1.36	extern int imm_irrad; /* calculate immediate irradiance? */
68
69	greg	1.1	extern int ralrm; /* seconds between reports */
70
71			extern int greyscale; /* map colors to brightness? */
72			extern char devname; / output device name */
73
74			extern int inform; /* input format */
75			extern int outform; /* output format */
76			extern char outvals; / output values */
77
78			extern VIEW ourview; /* viewing parameters */
79
80			extern int hresolu; /* horizontal resolution */
81			extern int vresolu; /* vertical resolution */
82	greg	1.22	extern double pixaspect; /* pixel aspect ratio */
83	greg	1.1
84			extern int psample; /* pixel sample size */
85			extern double maxdiff; /* max. sample difference */
86	greg	1.11	extern double dstrpix; /* square pixel distribution */
87	greg	1.1
88			extern double dstrsrc; /* square source distribution */
89	greg	1.11	extern double shadthresh; /* shadow threshold */
90	greg	1.15	extern double shadcert; /* shadow testing certainty */
91	greg	1.39	extern int directrelay; /* number of source relays */
92	greg	1.40	extern int vspretest; /* virtual source pretest density */
93	greg	1.42	extern int directinvis; /* light sources invisible to eye? */
94	greg	1.47	extern double srcsizerat; /* maximum source size/dist. ratio */
95	greg	1.1
96	greg	2.4	extern double specthresh; /* specular sampling threshold */
97			extern double specjitter; /* specular sampling jitter */
98
99	greg	1.1	extern int maxdepth; /* maximum recursion depth */
100			extern double minweight; /* minimum ray weight */
101
102			extern COLOR ambval; /* ambient value */
103			extern double ambacc; /* ambient accuracy */
104			extern int ambres; /* ambient resolution */
105			extern int ambdiv; /* ambient divisions */
106			extern int ambssamp; /* ambient super-samples */
107			extern int ambounce; /* ambient bounces */
108			extern char amblist[]; / ambient include/exclude list */
109			extern int ambincl; /* include == 1, exclude == 0 */
110
111
112			main(argc, argv)
113			int argc;
114			char *argv[];
115			{
116	greg	1.49	#define check(ol,al) if (argv[i][ol] \|\| \
117			badarg(argc-i-1,argv+i+1,al)) \
118			goto badopt
119	greg	1.43	#define bool(olen,var) switch (argv[i][olen]) { \
120			case '\0': var = !var; break; \
121			case 'y': case 'Y': case 't': case 'T': \
122			case '+': case '1': var = 1; break; \
123			case 'n': case 'N': case 'f': case 'F': \
124			case '-': case '0': var = 0; break; \
125			default: goto badopt; }
126	greg	1.45	extern char *getenv();
127	greg	1.1	char *err;
128			char *recover = NULL;
129	greg	1.20	char *zfile = NULL;
130	greg	1.1	char *errfile = NULL;
131			char *ambfile = NULL;
132			char **amblp = amblist;
133			int loadflags = ~IO_FILES;
134	greg	1.4	int rval, gotvfile = 0;
135	greg	1.50	RESOLU rs;
136	greg	1.1	int i;
137	greg	1.34	/* record start time */
138			tstart = time((long *)0);
139	greg	1.1	/* global program name */
140			progname = argv[0];
141			/* get library path */
142			if ((libpath = getenv(ULIBVAR)) == NULL)
143			libpath = DEFPATH;
144	greg	1.31	/* initialize object types */
145			initotypes();
146	greg	1.37	/* initialize urand */
147	greg	1.41	initurand(2048);
148	greg	1.10	/* option city */
149	greg	1.1	for (i = 1; i < argc && argv[i][0] == '-'; i++) {
150	greg	1.46	if (!strcmp(argv[i], "-version")) {
151			puts(VersionID);
152			quit(0);
153			}
154	greg	1.18	if (!strcmp(argv[i], "-defaults") \|\|
155			!strcmp(argv[i], "-help")) {
156	greg	1.1	printdefaults();
157			quit(0);
158			}
159			#if RVIEW
160			if (!strcmp(argv[i], "-devices")) {
161			printdevices();
162			quit(0);
163			}
164			#endif
165	greg	1.21	#if RPICT\|RVIEW
166			rval = getviewopt(&ourview, argc-i, argv+i);
167			if (rval >= 0) {
168			i += rval;
169			continue;
170			}
171			#endif
172	greg	1.1	switch (argv[i][1]) {
173			#if RPICT\|RVIEW
174	greg	1.21	case 'v': /* view file */
175			if (argv[i][2] != 'f')
176	greg	1.10	goto badopt;
177	greg	1.49	check(3,"s");
178	greg	1.50	rval = viewfile(argv[++i], &ourview, NULL);
179	greg	1.21	if (rval < 0) {
180			sprintf(errmsg,
181			"cannot open view file \"%s\"",
182			argv[i]);
183			error(SYSTEM, errmsg);
184			} else if (rval == 0) {
185			sprintf(errmsg,
186			"bad view file \"%s\"",
187			argv[i]);
188			error(USER, errmsg);
189			} else
190			gotvfile += rval;
191	greg	1.1	break;
192			#endif
193	greg	1.11	case 'd': /* direct */
194	greg	1.1	switch (argv[i][2]) {
195	greg	1.13	case 't': /* threshold */
196	greg	1.49	check(3,"f");
197	greg	1.11	shadthresh = atof(argv[++i]);
198	greg	1.1	break;
199	greg	1.15	case 'c': /* certainty */
200	greg	1.49	check(3,"f");
201	greg	1.15	shadcert = atof(argv[++i]);
202			break;
203	greg	1.11	case 'j': /* jitter */
204	greg	1.49	check(3,"f");
205	greg	1.1	dstrsrc = atof(argv[++i]);
206			break;
207	greg	1.47	case 'r': /* relays */
208	greg	1.49	check(3,"i");
209	greg	1.39	directrelay = atoi(argv[++i]);
210			break;
211	greg	1.47	case 'p': /* pretest */
212	greg	1.49	check(3,"i");
213	greg	1.40	vspretest = atoi(argv[++i]);
214	greg	1.42	break;
215	greg	1.47	case 'i': /* invis. */
216	greg	1.43	bool(3,directinvis);
217	greg	1.40	break;
218	greg	1.47	case 's': /* size */
219	greg	1.49	check(3,"f");
220	greg	1.47	srcsizerat = atof(argv[++i]);
221			break;
222	greg	1.1	default:
223	greg	1.10	goto badopt;
224	greg	1.1	}
225			break;
226	greg	2.4	case 's': /* specular */
227			switch (argv[i][2]) {
228			case 't': /* threshold */
229			check(3,"f");
230			specthresh = atof(argv[++i]);
231			break;
232			case 'j': /* jitter */
233			check(3,"f");
234			specjitter = atof(argv[++i]);
235			break;
236			default:
237			goto badopt;
238			}
239			break;
240	greg	1.1	#if RPICT\|RVIEW
241	greg	2.4	case 'p': /* pixel */
242	greg	1.1	switch (argv[i][2]) {
243	greg	2.4	case 's': /* sample */
244	greg	1.49	check(3,"i");
245	greg	1.1	psample = atoi(argv[++i]);
246			break;
247	greg	1.13	case 't': /* threshold */
248	greg	1.49	check(3,"f");
249	greg	1.1	maxdiff = atof(argv[++i]);
250			break;
251	greg	1.11	#if RPICT
252			case 'j': /* jitter */
253	greg	1.49	check(3,"f");
254	greg	1.11	dstrpix = atof(argv[++i]);
255			break;
256	greg	2.4	case 'a': /* aspect */
257			check(3,"f");
258			pixaspect = atof(argv[++i]);
259			break;
260	greg	1.11	#endif
261	greg	1.1	default:
262	greg	1.10	goto badopt;
263	greg	1.1	}
264			break;
265			#endif
266	greg	1.22	#if RPICT\|RTRACE
267	greg	1.1	case 'x': /* x resolution */
268	greg	1.49	check(2,"i");
269	greg	1.1	hresolu = atoi(argv[++i]);
270			break;
271			case 'y': /* y resolution */
272	greg	1.49	check(2,"i");
273	greg	1.1	vresolu = atoi(argv[++i]);
274			break;
275	greg	1.22	#endif
276	greg	1.1	case 'w': /* warnings */
277	greg	1.43	rval = wrnvec != NULL;
278			bool(2,rval);
279	greg	2.6	if (rval) wrnvec = stderr_v;
280			else wrnvec = NULL;
281	greg	1.1	break;
282			case 'e': /* error file */
283	greg	1.49	check(2,"s");
284	greg	1.1	errfile = argv[++i];
285			break;
286			case 'l': /* limit */
287			switch (argv[i][2]) {
288			case 'r': /* recursion */
289	greg	1.49	check(3,"i");
290	greg	1.1	maxdepth = atoi(argv[++i]);
291			break;
292			case 'w': /* weight */
293	greg	1.49	check(3,"f");
294	greg	1.1	minweight = atof(argv[++i]);
295			break;
296			default:
297	greg	1.10	goto badopt;
298	greg	1.1	}
299			break;
300	greg	1.36	case 'i': /* irradiance */
301	greg	1.43	bool(2,do_irrad);
302	greg	1.36	break;
303	greg	1.1	#if RPICT
304	greg	1.20	case 'z': /* z file */
305	greg	1.49	check(2,"s");
306	greg	1.20	zfile = argv[++i];
307			break;
308	greg	1.1	case 'r': /* recover file */
309	greg	1.49	check(2,"s");
310	greg	1.1	recover = argv[++i];
311	greg	1.50	rval = viewfile(recover, &ourview, &rs);
312			if (rval <= 0 \|\| rs.or != PIXSTANDARD) {
313	greg	1.5	sprintf(errmsg,
314			"cannot recover view parameters from \"%s\"", recover);
315			error(WARNING, errmsg);
316	greg	1.25	} else {
317	greg	1.4	gotvfile += rval;
318	greg	1.25	pixaspect = 0.0;
319	greg	1.50	hresolu = rs.xr;
320			vresolu = rs.yr;
321	greg	1.25	}
322	greg	1.1	break;
323			case 't': /* timer */
324	greg	1.49	check(2,"i");
325	greg	1.1	ralrm = atoi(argv[++i]);
326			break;
327			#endif
328			case 'a': /* ambient */
329			switch (argv[i][2]) {
330			case 'v': /* value */
331	greg	1.49	check(3,"fff");
332	greg	1.1	setcolor(ambval, atof(argv[i+1]),
333			atof(argv[i+2]),
334			atof(argv[i+3]));
335			i += 3;
336			break;
337			case 'a': /* accuracy */
338	greg	1.49	check(3,"f");
339	greg	1.1	ambacc = atof(argv[++i]);
340			break;
341			case 'r': /* resolution */
342	greg	1.49	check(3,"i");
343	greg	1.1	ambres = atoi(argv[++i]);
344			break;
345			case 'd': /* divisions */
346	greg	1.49	check(3,"i");
347	greg	1.1	ambdiv = atoi(argv[++i]);
348			break;
349			case 's': /* super-samp */
350	greg	1.49	check(3,"i");
351	greg	1.1	ambssamp = atoi(argv[++i]);
352			break;
353			case 'b': /* bounces */
354	greg	1.49	check(3,"i");
355	greg	1.1	ambounce = atoi(argv[++i]);
356			break;
357			case 'i': /* include */
358	greg	2.5	case 'I':
359	greg	1.49	check(3,"s");
360	greg	1.1	if (ambincl != 1) {
361			ambincl = 1;
362			amblp = amblist;
363			}
364	greg	2.5	if (argv[i][2] == 'I') { /* file */
365			rval = wordfile(amblp, argv[++i]);
366			if (rval < 0) {
367			sprintf(errmsg,
368			"cannot open ambient include file \"%s\"",
369			argv[i]);
370			error(SYSTEM, errmsg);
371			}
372			amblp += rval;
373			} else
374			*amblp++ = argv[++i];
375	greg	1.1	break;
376			case 'e': /* exclude */
377	greg	2.5	case 'E':
378	greg	1.49	check(3,"s");
379	greg	1.1	if (ambincl != 0) {
380			ambincl = 0;
381			amblp = amblist;
382			}
383	greg	2.5	if (argv[i][2] == 'E') { /* file */
384			rval = wordfile(amblp, argv[++i]);
385			if (rval < 0) {
386			sprintf(errmsg,
387			"cannot open ambient exclude file \"%s\"",
388			argv[i]);
389			error(SYSTEM, errmsg);
390			}
391			amblp += rval;
392			} else
393			*amblp++ = argv[++i];
394	greg	1.1	break;
395			case 'f': /* file */
396	greg	1.49	check(3,"s");
397	greg	1.1	ambfile= argv[++i];
398			break;
399			default:
400	greg	1.10	goto badopt;
401	greg	1.1	}
402			break;
403			#if RTRACE
404	greg	1.36	case 'I': /* immed. irradiance */
405	greg	1.43	bool(2,imm_irrad);
406	greg	1.36	break;
407	greg	1.1	case 'f': /* format i/o */
408			switch (argv[i][2]) {
409			case 'a': /* ascii */
410			case 'f': /* float */
411			case 'd': /* double */
412			inform = argv[i][2];
413			break;
414			default:
415	greg	1.10	goto badopt;
416	greg	1.1	}
417			switch (argv[i][3]) {
418			case '\0':
419			outform = inform;
420			break;
421			case 'a': /* ascii */
422			case 'f': /* float */
423			case 'd': /* double */
424	greg	1.49	check(4,"");
425	greg	1.1	outform = argv[i][3];
426			break;
427			default:
428	greg	1.10	goto badopt;
429	greg	1.1	}
430			break;
431			case 'o': /* output */
432			outvals = argv[i]+2;
433			break;
434	greg	1.44	case 'h': /* header output */
435			rval = loadflags & IO_INFO;
436			bool(2,rval);
437			loadflags = rval ? loadflags \| IO_INFO :
438			loadflags & ~IO_INFO;
439	greg	1.1	break;
440			#endif
441			#if RVIEW
442			case 'b': /* black and white */
443	greg	1.49	bool(2,greyscale);
444	greg	1.1	break;
445			case 'o': /* output device */
446	greg	1.49	check(2,"s");
447	greg	1.1	devname = argv[++i];
448			break;
449			#endif
450			default:
451	greg	1.19	goto badopt;
452	greg	1.1	}
453			}
454	greg	1.24	*amblp = NULL;
455	greg	1.1	#if RPICT\|RVIEW
456			err = setview(&ourview); /* set viewing parameters */
457			if (err != NULL)
458			error(USER, err);
459			#endif
460	greg	1.22	#if RPICT
461	greg	1.23	normaspect(viewaspect(&ourview), &pixaspect, &hresolu, &vresolu);
462	greg	1.22	#endif
463	greg	1.1	/* set up signal handling */
464			sigdie(SIGINT, "Interrupt");
465			sigdie(SIGHUP, "Hangup");
466			sigdie(SIGTERM, "Terminate");
467			sigdie(SIGPIPE, "Broken pipe");
468	greg	1.45	sigdie(SIGALRM, "Alarm clock");
469	greg	1.1	#ifdef SIGXCPU
470			sigdie(SIGXCPU, "CPU limit exceeded");
471			sigdie(SIGXFSZ, "File size exceeded");
472			#endif
473			/* open error file */
474			if (errfile != NULL) {
475			if (freopen(errfile, "a", stderr) == NULL)
476	greg	1.14	quit(2);
477	greg	1.1	fprintf(stderr, "************\n* PID %5d: ",
478			getpid());
479			printargs(argc, argv, stderr);
480			fputs("\n", stderr);
481			fflush(stderr);
482			}
483			#ifdef NICE
484			nice(NICE); /* lower priority */
485			#endif
486	greg	1.19	/* get octree */
487	greg	1.1	#if RVIEW
488			loadflags &= ~IO_INFO;
489			#endif
490			if (i == argc)
491	greg	1.19	octname = NULL;
492			else if (i == argc-1)
493			octname = argv[i];
494	greg	1.1	else
495	greg	1.19	goto badopt;
496			#if RVIEW\|RTRACE
497			if (octname == NULL)
498			error(USER, "missing octree argument");
499	greg	1.1	#endif
500	greg	1.19	readoct(octname, loadflags, &thescene, NULL);
501	greg	1.1
502			if (loadflags & IO_INFO) { /* print header */
503			printargs(i, argv, stdout);
504	greg	1.46	printf("SOFTWARE= %s\n", VersionID);
505	greg	1.1	#if RPICT
506			if (gotvfile) {
507			printf(VIEWSTR);
508			fprintview(&ourview, stdout);
509			printf("\n");
510			}
511	greg	1.23	if (pixaspect < .99 \|\| pixaspect > 1.01)
512	greg	1.22	fputaspect(pixaspect, stdout);
513	greg	1.35	fputformat(COLRFMT, stdout);
514			#endif
515			#if RTRACE
516			fputformat( outform=='a' ? "ascii" :
517			outform=='f' ? "float" :
518			"double", stdout );
519	greg	1.1	#endif
520			printf("\n");
521			}
522
523			marksources(); /* find and mark sources */
524
525	greg	1.24	setambient(ambfile); /* initialize ambient calculation */
526	greg	1.1
527			#if RPICT
528	greg	1.20	render(zfile, recover); /* render the scene */
529	greg	1.1	#endif
530			#if RTRACE
531			rtrace(NULL); /* trace rays from stdin */
532			#endif
533			#if RVIEW
534			rview(); /* go */
535			#endif
536			quit(0);
537	greg	1.19
538			badopt:
539			sprintf(errmsg, "command line error at '%s'", argv[i]);
540			error(USER, errmsg);
541
542	greg	1.10	#undef check
543	greg	1.43	#undef bool
544	greg	1.1	}
545
546
547			eputs(s) /* error output */
548			char *s;
549			{
550			if (errvec != NULL)
551			(*errvec)(s);
552			}
553
554
555			wputs(s) /* warning output */
556			char *s;
557			{
558	greg	2.3	int lasterrno = errno; /* save errno */
559
560	greg	1.1	if (wrnvec != NULL)
561			(*wrnvec)(s);
562	greg	2.3
563			errno = lasterrno;
564	greg	1.1	}
565
566
567			cputs(s) /* command error output */
568			char *s;
569			{
570			if (cmdvec != NULL)
571			(*cmdvec)(s);
572			}
573
574
575			stderr_v(s) /* put string to stderr */
576			register char *s;
577			{
578	greg	1.26	static int midline = 0;
579	greg	1.1
580	greg	1.27	if (!*s)
581			return;
582	greg	1.26	if (!midline++) {
583	greg	1.1	fputs(progname, stderr);
584			fputs(": ", stderr);
585			}
586			fputs(s, stderr);
587	greg	1.27	if (s[strlen(s)-1] == '\n') {
588	greg	1.1	fflush(stderr);
589	greg	1.26	midline = 0;
590	greg	1.1	}
591			}
592
593
594			onsig(signo) /* fatal signal */
595			int signo;
596			{
597	greg	1.8	static int gotsig = 0;
598
599			if (gotsig++) /* two signals and we're gone! */
600	greg	1.9	_exit(signo);
601	greg	1.8
602	greg	1.48	alarm(15); /* allow 15 seconds to clean up */
603			signal(SIGALRM, SIG_DFL); /* make certain we do die */
604	greg	1.1	eputs("signal - ");
605			eputs(sigerr[signo]);
606			eputs("\n");
607	greg	1.14	quit(3);
608	greg	1.1	}
609
610
611			sigdie(signo, msg) /* set fatal signal */
612			int signo;
613			char *msg;
614			{
615			if (signal(signo, onsig) == SIG_IGN)
616			signal(signo, SIG_IGN);
617			sigerr[signo] = msg;
618			}
619
620
621			printdefaults() /* print default values to stdout */
622			{
623			register char *cp;
624
625	greg	1.43	#if RTRACE
626			if (imm_irrad)
627			printf("-I+\t\t\t\t# immediate irradiance on\n");
628			else
629			#endif
630			printf(do_irrad ? "-i+\t\t\t\t# irradiance calculation on\n" :
631			"-i-\t\t\t\t# irradiance calculation off\n");
632	greg	1.49	#if RVIEW
633			printf(greyscale ? "-b+\t\t\t\t# greyscale on\n" :
634			"-b-\t\t\t\t# greyscale off\n");
635			#endif
636	greg	1.1	#if RPICT\|RVIEW
637	greg	1.7	printf("-vt%c\t\t\t\t# view type %s\n", ourview.type,
638	greg	1.1	ourview.type==VT_PER ? "perspective" :
639			ourview.type==VT_PAR ? "parallel" :
640	greg	1.28	ourview.type==VT_HEM ? "hemispherical" :
641			ourview.type==VT_ANG ? "angular" :
642	greg	1.1	"unknown");
643			printf("-vp %f %f %f\t# view point\n",
644			ourview.vp[0], ourview.vp[1], ourview.vp[2]);
645			printf("-vd %f %f %f\t# view direction\n",
646			ourview.vdir[0], ourview.vdir[1], ourview.vdir[2]);
647			printf("-vu %f %f %f\t# view up\n",
648			ourview.vup[0], ourview.vup[1], ourview.vup[2]);
649			printf("-vh %f\t\t\t# view horizontal size\n", ourview.horiz);
650			printf("-vv %f\t\t\t# view vertical size\n", ourview.vert);
651	greg	1.22	printf("-vs %f\t\t\t# view shift\n", ourview.hoff);
652			printf("-vl %f\t\t\t# view lift\n", ourview.voff);
653	greg	1.1	#endif
654	greg	1.22	#if RPICT\|RTRACE
655	greg	1.1	printf("-x %-9d\t\t\t# x resolution\n", hresolu);
656			printf("-y %-9d\t\t\t# y resolution\n", vresolu);
657	greg	1.22	#endif
658			#if RPICT
659	greg	2.4	printf("-pa %f\t\t\t# pixel aspect ratio\n", pixaspect);
660			printf("-pj %f\t\t\t# pixel jitter\n", dstrpix);
661	greg	1.22	#endif
662	greg	1.1	#if RPICT\|RVIEW
663	greg	2.4	printf("-ps %-9d\t\t\t# pixel sample\n", psample);
664			printf("-pt %f\t\t\t# pixel threshold\n", maxdiff);
665	greg	1.22	#endif
666	greg	1.29	printf("-dt %f\t\t\t# direct threshold\n", shadthresh);
667			printf("-dc %f\t\t\t# direct certainty\n", shadcert);
668			printf("-dj %f\t\t\t# direct jitter\n", dstrsrc);
669	greg	1.47	printf("-ds %f\t\t\t# direct sampling\n", srcsizerat);
670	greg	1.39	printf("-dr %-9d\t\t\t# direct relays\n", directrelay);
671	greg	1.40	printf("-dp %-9d\t\t\t# direct pretest density\n", vspretest);
672	greg	1.43	printf(directinvis ? "-di+\t\t\t\t# direct invisibility on\n" :
673			"-di-\t\t\t\t# direct invisibility off\n");
674	greg	2.4	printf("-sj %f\t\t\t# specular jitter\n", specjitter);
675			printf("-st %f\t\t\t# specular threshold\n", specthresh);
676	greg	1.1	printf("-av %f %f %f\t# ambient value\n", colval(ambval,RED),
677			colval(ambval,GRN), colval(ambval, BLU));
678			printf("-ab %-9d\t\t\t# ambient bounces\n", ambounce);
679			printf("-aa %f\t\t\t# ambient accuracy\n", ambacc);
680			printf("-ar %-9d\t\t\t# ambient resolution\n", ambres);
681			printf("-ad %-9d\t\t\t# ambient divisions\n", ambdiv);
682			printf("-as %-9d\t\t\t# ambient super-samples\n", ambssamp);
683			printf("-lr %-9d\t\t\t# limit reflection\n", maxdepth);
684			printf("-lw %f\t\t\t# limit weight\n", minweight);
685			#if RPICT
686			printf("-t %-9d\t\t\t# time between reports\n", ralrm);
687			#endif
688			#if RVIEW
689			printf("-o %s\t\t\t\t# output device\n", devname);
690			#endif
691			#if RTRACE
692			printf("-f%c%c\t\t\t\t# format input/output = %s/%s\n",
693			inform, outform,
694			inform=='a' ? "ascii" :
695			inform=='f' ? "float" : "double",
696			outform=='a' ? "ascii" :
697			outform=='f' ? "float" : "double");
698			printf("-o%s\t\t\t\t# output", outvals);
699			for (cp = outvals; *cp; cp++)
700			switch (*cp) {
701			case 't': printf(" trace"); break;
702			case 'o': printf(" origin"); break;
703			case 'd': printf(" direction"); break;
704			case 'v': printf(" value"); break;
705			case 'l': printf(" length"); break;
706			case 'p': printf(" point"); break;
707			case 'n': printf(" normal"); break;
708			case 's': printf(" surface"); break;
709			case 'w': printf(" weight"); break;
710			case 'm': printf(" modifier"); break;
711			}
712			printf("\n");
713			#endif
714	greg	1.43	printf(wrnvec != NULL ? "-w+\t\t\t\t# warning messages on\n" :
715			"-w-\t\t\t\t# warning messages off\n");
716	greg	1.1	}