src/rt/rmain.c

/* Copyright (c) 1992 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"

#include  "paths.h"

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 */
OBJECT  nsceneobjs;                     /* number of objects in 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 char  *formstr();                /* string from format */
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  *outfile = NULL;
        char  *zfile = NULL;
        char  *errfile = NULL;
        char  *ambfile = NULL;
        char  **amblp = amblist;
        int  loadflags = ~IO_FILES;
        int  seqstart = 0;
        int  rval;
        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);
                        }
                        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 'S':                               /* slave index */
                        check(2,"i");
                        seqstart = atoi(argv[++i]);
                        break;
                case 'o':                               /* output file */
                        check(2,"s");
                        outfile = argv[++i];
                        break;
                case 'z':                               /* z file */
                        check(2,"s");
                        zfile = argv[++i];
                        break;
                case 'r':                               /* recover file */
                        check(2,"s");
                        recover = argv[++i];
                        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 */
                        case 'c':                               /* color */
                                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
                                        /* 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 (
#if  RPICT
                seqstart > 0 &&
#endif
                        octname == NULL)
                error(USER, "missing octree argument");
#if  RPICT
        if (outfile != NULL)
                openheader();
#endif
#ifdef  MSDOS
#if  RTRACE
        if (outform != 'a')
#endif
        setmode(fileno(stdout), O_BINARY);
        if (octname == NULL)
                setmode(fileno(stdin), O_BINARY);
#endif
        readoct(octname, loadflags, &thescene, NULL);
        nsceneobjs = nobjects;

        if (loadflags & IO_INFO) {      /* print header */
                printargs(i, argv, stdout);
                printf("SOFTWARE= %s\n", VersionID);
#if  RTRACE
                fputformat(formstr(outform), stdout);
                putchar('\n');
#endif
        }

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

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

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