src/gen/mkillum.c

#ifndef lint
static const char RCSid[] = "$Id: mkillum.c,v 2.23 2004/10/28 00:50:47 greg Exp $";
#endif
/*
 * Make illum sources for optimizing rendering process
 */

#include  <signal.h>
#include  <ctype.h>

#include  "platform.h"
#include  "mkillum.h"
#include  "random.h"

                                /* default parameters */
#define  SAMPDENS       48              /* points per projected steradian */
#define  NSAMPS         32              /* samples per point */
#define  DFLMAT         "illum_mat"     /* material name */
#define  DFLDAT         "illum"         /* data file name */
                                /* selection options */
#define  S_NONE         0               /* select none */
#define  S_ELEM         1               /* select specified element */
#define  S_COMPL        2               /* select all but element */
#define  S_ALL          3               /* select all */

                                /* rtrace command and defaults */
char  *rtargv[64] = { "rtrace", "-dj", ".25", "-dr", "3", "-dv-",
                "-ab", "2", "-ad", "1024", "-as", "512", "-aa", ".1", };
int  rtargc = 14;
                                /* overriding rtrace options */
char  *myrtopts[] = { "-I-", "-i-", "-ld-", "-ov", "-h-",
                        "-fff", "-y", "0", NULL };

struct rtproc   rt0;            /* head of rtrace process list */

struct illum_args  thisillum = {        /* our illum and default values */
                0,
                DFLMAT,
                DFLDAT,
                0,
                VOIDID,
                SAMPDENS,
                NSAMPS,
                0.,
        };

char    matcheck[MAXSTR];       /* current material to include or exclude */
int     matselect = S_ALL;      /* selection criterion */

FUN     ofun[NUMOTYPE] = INIT_OTYPE;    /* object types */

char    persistfn[] = "pfXXXXXX";       /* persist file name */

int     gargc;                  /* global argc */
char    **gargv;                /* global argv */
#define  progname       gargv[0]

int     doneheader = 0;         /* printed header yet? */
#define  checkhead()    if (!doneheader++) printhead(gargc,gargv)

int     warnings = 1;           /* print warnings? */

int done_rprocs(struct rtproc *rtp);
void init(int np);
void filter(register FILE       *infp, char     *name);
void xoptions(char      *s, char        *nm);
void printopts(void);
void printhead(register int  ac, register char  **av);
void xobject(FILE  *fp, char  *nm);


int
main(           /* compute illum distributions using rtrace */
        int     argc,
        char    *argv[]
)
{
        int     nprocs = 1;
        char    *rtpath;
        FILE    *fp;
        register int    i;
                                /* set global arguments */
        gargv = argv;
                                /* check for -n option */
        if (!strcmp(argv[1], "-n")) {
                nprocs = atoi(argv[2]);
                if (nprocs <= 0)
                        error(USER, "illegal number of processes");
                i = 3;
        } else
                i = 1;
                                /* set up rtrace command */
        for ( ; i < argc; i++) {
                if (argv[i][0] == '<' && argv[i][1] == '\0')
                        break;
                rtargv[rtargc++] = argv[i];
                if (argv[i][0] == '-' && argv[i][1] == 'w')
                        switch (argv[i][2]) {
                        case '\0':
                                warnings = !warnings;
                                break;
                        case '+':
                        case 'T': case 't':
                        case 'Y': case 'y':
                        case '1':
                                warnings = 1;
                                break;
                        case '-':
                        case 'F': case 'f':
                        case 'N': case 'n':
                        case '0':
                                warnings = 0;
                                break;
                        }
        }
        gargc = i;
        if (!strcmp(rtargv[--rtargc], "-defaults"))
                nprocs = 0;
        if (nprocs > 1) {       /* add persist file if parallel invocation */
                rtargv[rtargc++] = "-PP";
                rtargv[rtargc++] = mktemp(persistfn);
        }
                                /* add "mandatory" rtrace options */
        for (i = 0; myrtopts[i] != NULL; i++)
                rtargv[rtargc++] = myrtopts[i];
                                /* finally, put back final argument */
        rtargv[rtargc++] = argv[gargc-1];
        rtargv[rtargc] = NULL;
        if (!nprocs) {          /* just asking for defaults? */
                printopts(); fflush(stdout);
                rtpath = getpath(rtargv[0], getenv("PATH"), X_OK);
                if (rtpath == NULL) {
                        eputs(rtargv[0]);
                        eputs(": command not found\n");
                        exit(1);
                }
                execv(rtpath, rtargv);
                perror(rtpath); /* execv() should not return */
                exit(1);
        }
        if (gargc < 2 || argv[gargc-1][0] == '-')
                error(USER, "missing octree argument");
                                /* else initialize and run our calculation */
        init(nprocs);
        if (gargc+1 < argc)
                for (i = gargc+1; i < argc; i++) {
                        if ((fp = fopen(argv[i], "r")) == NULL) {
                                sprintf(errmsg,
                                "cannot open scene file \"%s\"", argv[i]);
                                error(SYSTEM, errmsg);
                        }
                        filter(fp, argv[i]);
                        fclose(fp);
                }
        else
                filter(stdin, "standard input");
        quit(0);
}

static void
killpersist(void)                       /* kill persistent rtrace process */
{
        FILE    *fp = fopen(persistfn, "r");
        int     pid;

        if (fp == NULL)
                return;
        if (fscanf(fp, "%*s %d", &pid) != 1 || kill(pid, SIGALRM) < 0)
                unlink(persistfn);
        fclose(fp);
}

int
done_rprocs(struct rtproc *rtp)
{
        int     st0, st1 = 0;

        if (rtp->next != NULL) {        /* close last opened first! */
                st1 = done_rprocs(rtp->next);
                free((void *)rtp->next);
                rtp->next = NULL;
        }
        st0 = close_process(&rtp->pd);
        if (st0 < 0)
                error(WARNING, "unknown return status from rtrace process");
        else if (st0 > 0)
                return(st0);
        return(st1);
}

void
quit(int status)                        /* exit with status */
{
        int     rtstat;

        if (rt0.next != NULL)           /* terminate persistent rtrace */
                killpersist();
                                        /* clean up rtrace process(es) */
        rtstat = done_rprocs(&rt0);
        if (status == 0)
                status = rtstat;
        exit(status);
}

void
init(int np)                            /* start rtrace and set up buffers */
{
        struct rtproc   *rtp;
        int     i;
        int     maxbytes;
                                        /* set up object functions */
        ofun[OBJ_FACE].funp = o_face;
        ofun[OBJ_SPHERE].funp = o_sphere;
        ofun[OBJ_RING].funp = o_ring;
                                        /* set up signal handling */
#ifdef SIGPIPE /* not present on Windows */
        signal(SIGPIPE, quit);
#endif
        rtp = &rt0;                     /* start rtrace process(es) */
        for (i = 0; i++ < np; ) {
                errno = 0;
                maxbytes = open_process(&rtp->pd, rtargv);
                if (maxbytes == 0) {
                        eputs(rtargv[0]);
                        eputs(": command not found\n");
                        exit(1);
                }
                if (maxbytes < 0)
                        error(SYSTEM, "cannot start rtrace process");
                rtp->bsiz = maxbytes/(6*sizeof(float));
                rtp->buf = (float *)malloc(6*sizeof(float)*rtp->bsiz--);
                rtp->dest = (float **)calloc(rtp->bsiz, sizeof(float *));
                if (rtp->buf == NULL || rtp->dest == NULL)
                        error(SYSTEM, "out of memory in init");
                rtp->nrays = 0;
                if (i == np)            /* last process? */
                        break;
                if (np > 1)
                        sleep(2);       /* wait for persist file */
                rtp->next = (struct rtproc *)malloc(sizeof(struct rtproc));
                if (rtp->next == NULL)
                        error(SYSTEM, "out of memory in init");
                rtp = rtp->next;
        }
        rtp->next = NULL;
                                        /* set up urand */
        initurand(16384);
}


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

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


void
wputs(s)                        /* print warning if enabled */
char  *s;
{
        if (warnings)
                eputs(s);
}


void
filter(         /* process stream */
        register FILE   *infp,
        char    *name
)
{
        char    buf[512];
        FILE    *pfp;
        register int    c;

        while ((c = getc(infp)) != EOF) {
                if (isspace(c))
                        continue;
                if (c == '#') {                         /* comment/options */
                        buf[0] = c;
                        fgets(buf+1, sizeof(buf)-1, infp);
                        xoptions(buf, name);
                } else if (c == '!') {                  /* command */
                        buf[0] = c;
                        fgetline(buf+1, sizeof(buf)-1, infp);
                        if ((pfp = popen(buf+1, "r")) == NULL) {
                                sprintf(errmsg, "cannot execute \"%s\"", buf);
                                error(SYSTEM, errmsg);
                        }
                        filter(pfp, buf);
                        pclose(pfp);
                } else {                                /* object */
                        ungetc(c, infp);
                        xobject(infp, name);
                }
        }
}


void
xoptions(                       /* process options in string s */
        char    *s,
        char    *nm
)
{
        extern FILE     *freopen();
        char    buf[64];
        int     nerrs = 0;
        register char   *cp;

        if (strncmp(s, "#@mkillum", 9) || !isspace(s[9])) {
                fputs(s, stdout);               /* not for us */
                return;
        }
        cp = s+10;
        while (*cp) {
                switch (*cp) {
                case ' ':
                case '\t':
                case '\n':
                case '\r':
                case '\f':
                        cp++;
                        continue;
                case 'm':                       /* material name */
                        if (*++cp != '=')
                                break;
                        if (!*++cp || isspace(*cp))
                                break;
                        atos(thisillum.matname, MAXSTR, cp);
                        cp = sskip(cp);
                        if (!(thisillum.flags & IL_DATCLB)) {
                                strcpy(thisillum.datafile, thisillum.matname);
                                thisillum.dfnum = 0;
                        }
                        continue;
                case 'f':                       /* data file name */
                        if (*++cp != '=')
                                break;
                        if (!*++cp || isspace(*cp)) {
                                strcpy(thisillum.datafile,thisillum.matname);
                                thisillum.dfnum = 0;
                                thisillum.flags &= ~IL_DATCLB;
                                continue;
                        }
                        atos(thisillum.datafile, MAXSTR, cp);
                        cp = sskip(cp);
                        thisillum.dfnum = 0;
                        thisillum.flags |= IL_DATCLB;
                        continue;
                case 'i':                       /* include material */
                case 'e':                       /* exclude material */
                        if (cp[1] != '=')
                                break;
                        matselect = (*cp == 'i') ? S_ELEM : S_COMPL;
                        cp += 2;
                        atos(matcheck, MAXSTR, cp);
                        cp = sskip(cp);
                        continue;
                case 'a':                       /* use everything */
                        cp = sskip(cp);
                        matselect = S_ALL;
                        continue;
                case 'n':                       /* use nothing (passive) */
                        cp = sskip(cp);
                        matselect = S_NONE;
                        continue;
                case 'c':                       /* color calculation */
                        if (*++cp != '=')
                                break;
                        switch (*++cp) {
                        case 'a':                       /* average */
                                thisillum.flags = (thisillum.flags|IL_COLAVG)
                                                        & ~IL_COLDST;
                                break;
                        case 'd':                       /* distribution */
                                thisillum.flags |= (IL_COLDST|IL_COLAVG);
                                break;
                        case 'n':                       /* none */
                                thisillum.flags &= ~(IL_COLAVG|IL_COLDST);
                                break;
                        default:
                                goto opterr;
                        }
                        cp = sskip(cp);
                        continue;
                case 'd':                       /* point sample density */
                        if (*++cp != '=')
                                break;
                        if (!isintd(++cp, " \t\n\r"))
                                break;
                        thisillum.sampdens = atoi(cp);
                        cp = sskip(cp);
                        continue;
                case 's':                       /* point super-samples */
                        if (*++cp != '=')
                                break;
                        if (!isintd(++cp, " \t\n\r"))
                                break;
                        thisillum.nsamps = atoi(cp);
                        cp = sskip(cp);
                        continue;
                case 'l':                       /* light sources */
                        cp++;
                        if (*cp == '+')
                                thisillum.flags |= IL_LIGHT;
                        else if (*cp == '-')
                                thisillum.flags &= ~IL_LIGHT;
                        else
                                break;
                        cp++;
                        continue;
                case 'b':                       /* brightness */
                        if (*++cp != '=')
                                break;
                        if (!isfltd(++cp, " \t\n\r"))
                                break;
                        thisillum.minbrt = atof(cp);
                        if (thisillum.minbrt < 0.)
                                thisillum.minbrt = 0.;
                        cp = sskip(cp);
                        continue;
                case 'o':                       /* output file */
                        if (*++cp != '=')
                                break;
                        if (!*++cp || isspace(*cp))
                                break;
                        atos(buf, sizeof(buf), cp);
                        cp = sskip(cp);
                        if (freopen(buf, "w", stdout) == NULL) {
                                sprintf(errmsg,
                                "cannot open output file \"%s\"", buf);
                                error(SYSTEM, errmsg);
                        }
                        doneheader = 0;
                        continue;
                case '!':                       /* processed file! */
                        sprintf(errmsg, "(%s): already processed!", nm);
                        error(WARNING, errmsg);
                        matselect = S_NONE;
                        return;
                }
        opterr:                                 /* skip faulty option */
                while (*cp && !isspace(*cp))
                        cp++;
                nerrs++;
        }
                                                /* print header? */
        checkhead();
                                                /* issue warnings? */
        if (nerrs) {
                sprintf(errmsg, "(%s): %d error(s) in option line:",
                                nm, nerrs);
                error(WARNING, errmsg);
                wputs(s);
                printf("# %s: the following option line has %d error(s):\n",
                                progname, nerrs);
        }
                                                /* print pure comment */
        printf("# %s", s+2);
}

void
printopts(void)                 /* print out option default values */
{
        printf("m=%-15s\t\t# material name\n", thisillum.matname);
        printf("f=%-15s\t\t# data file name\n", thisillum.datafile);
        if (thisillum.flags & IL_COLAVG)
                if (thisillum.flags & IL_COLDST)
                        printf("c=d\t\t\t\t# color distribution\n");
                else
                        printf("c=a\t\t\t\t# color average\n");
        else
                printf("c=n\t\t\t\t# color none\n");
        if (thisillum.flags & IL_LIGHT)
                printf("l+\t\t\t\t# light type on\n");
        else
                printf("l-\t\t\t\t# light type off\n");
        printf("d=%d\t\t\t\t# density of points\n", thisillum.sampdens);
        printf("s=%d\t\t\t\t# samples per point\n", thisillum.nsamps);
        printf("b=%f\t\t\t# minimum average brightness\n", thisillum.minbrt);
}


void
printhead(                      /* print out header */
        register int  ac,
        register char  **av
)
{
        putchar('#');
        while (ac-- > 0) {
                putchar(' ');
                fputs(*av++, stdout);
        }
        fputs("\n#@mkillum !\n", stdout);
}


void
xobject(                                /* translate an object from fp */
        FILE  *fp,
        char  *nm
)
{
        OBJREC  thisobj;
        char  str[MAXSTR];
        int  doit;
                                        /* read the object */
        if (fgetword(thisillum.altmat, MAXSTR, fp) == NULL)
                goto readerr;
        if (fgetword(str, MAXSTR, fp) == NULL)
                goto readerr;
                                        /* is it an alias? */
        if (!strcmp(str, ALIASKEY)) {
                if (fgetword(str, MAXSTR, fp) == NULL)
                        goto readerr;
                printf("\n%s %s %s", thisillum.altmat, ALIASKEY, str);
                if (fgetword(str, MAXSTR, fp) == NULL)
                        goto readerr;
                printf("\t%s\n", str);
                return;
        }
        thisobj.omod = OVOID;           /* unused field */
        if ((thisobj.otype = otype(str)) < 0) {
                sprintf(errmsg, "(%s): unknown type \"%s\"", nm, str);
                error(USER, errmsg);
        }
        if (fgetword(str, MAXSTR, fp) == NULL)
                goto readerr;
        thisobj.oname = str;
        if (readfargs(&thisobj.oargs, fp) != 1)
                goto readerr;
        thisobj.os = NULL;
                                        /* check for translation */
        switch (matselect) {
        case S_NONE:
                doit = 0;
                break;
        case S_ALL:
                doit = 1;
                break;
        case S_ELEM:
                doit = !strcmp(thisillum.altmat, matcheck);
                break;
        case S_COMPL:
                doit = strcmp(thisillum.altmat, matcheck);
                break;
        }
        doit = doit && issurface(thisobj.otype);
                                                /* print header? */
        checkhead();
                                                /* process object */
        if (doit)
                (*ofun[thisobj.otype].funp)(&thisobj, &thisillum, &rt0, nm);
        else
                printobj(thisillum.altmat, &thisobj);
                                                /* free arguments */
        freefargs(&thisobj.oargs);
        return;
readerr:
        sprintf(errmsg, "(%s): error reading scene", nm);
        error(USER, errmsg);
}
Revision:	2.24
Committed:	Wed May 25 04:44:25 2005 UTC (18 years, 11 months ago) by greg
Content type:	text/plain
Branch:	MAIN
Changes since 2.23:	+3 -4 lines
Log Message:	Created rtcontrib program for computing ray contributions and coefficients
#	User	Rev	Content
1	greg	1.1	#ifndef lint
2	greg	2.24	static const char RCSid[] = "$Id: mkillum.c,v 2.23 2004/10/28 00:50:47 greg Exp $";
3	greg	1.1	#endif
4			/*
5			* Make illum sources for optimizing rendering process
6			*/
7
8	greg	1.3	#include <signal.h>
9	schorsch	2.13	#include <ctype.h>
10	greg	1.3
11	greg	2.18	#include "platform.h"
12	greg	2.15	#include "mkillum.h"
13	schorsch	2.19	#include "random.h"
14	greg	1.1
15			/* default parameters */
16	greg	1.8	#define SAMPDENS 48 /* points per projected steradian */
17	greg	1.1	#define NSAMPS 32 /* samples per point */
18			#define DFLMAT "illum_mat" /* material name */
19	greg	1.10	#define DFLDAT "illum" /* data file name */
20	greg	1.1	/* selection options */
21			#define S_NONE 0 /* select none */
22			#define S_ELEM 1 /* select specified element */
23			#define S_COMPL 2 /* select all but element */
24			#define S_ALL 3 /* select all */
25
26			/* rtrace command and defaults */
27	greg	2.5	char *rtargv[64] = { "rtrace", "-dj", ".25", "-dr", "3", "-dv-",
28	greg	2.17	"-ab", "2", "-ad", "1024", "-as", "512", "-aa", ".1", };
29	greg	1.4	int rtargc = 14;
30	greg	1.1	/* overriding rtrace options */
31	gregl	2.10	char *myrtopts[] = { "-I-", "-i-", "-ld-", "-ov", "-h-",
32			"-fff", "-y", "0", NULL };
33	greg	1.1
34	greg	2.20	struct rtproc rt0; /* head of rtrace process list */
35	greg	1.1
36	greg	1.2	struct illum_args thisillum = { /* our illum and default values */
37			0,
38			DFLMAT,
39	greg	1.10	DFLDAT,
40	greg	1.2	0,
41			VOIDID,
42	greg	1.3	SAMPDENS,
43	greg	1.2	NSAMPS,
44	greg	1.12	0.,
45	greg	1.2	};
46	greg	1.1
47			char matcheck[MAXSTR]; /* current material to include or exclude */
48			int matselect = S_ALL; /* selection criterion */
49
50	greg	1.2	FUN ofun[NUMOTYPE] = INIT_OTYPE; /* object types */
51	greg	1.1
52	greg	2.20	char persistfn[] = "pfXXXXXX"; /* persist file name */
53
54	greg	1.1	int gargc; /* global argc */
55			char *gargv; / global argv */
56			#define progname gargv[0]
57
58			int doneheader = 0; /* printed header yet? */
59	greg	1.2	#define checkhead() if (!doneheader++) printhead(gargc,gargv)
60	greg	1.1
61	greg	1.2	int warnings = 1; /* print warnings? */
62
63	greg	2.23	int done_rprocs(struct rtproc *rtp);
64	greg	2.20	void init(int np);
65	schorsch	2.19	void filter(register FILE infp, char name);
66			void xoptions(char s, char nm);
67			void printopts(void);
68			void printhead(register int ac, register char **av);
69			void xobject(FILE fp, char nm);
70	greg	1.1
71	schorsch	2.19
72			int
73			main( /* compute illum distributions using rtrace */
74			int argc,
75			char *argv[]
76			)
77	greg	1.1	{
78	greg	2.20	int nprocs = 1;
79	greg	1.1	char *rtpath;
80	greg	1.3	FILE *fp;
81	greg	1.1	register int i;
82			/* set global arguments */
83	greg	1.3	gargv = argv;
84	greg	2.20	/* check for -n option */
85			if (!strcmp(argv[1], "-n")) {
86			nprocs = atoi(argv[2]);
87			if (nprocs <= 0)
88			error(USER, "illegal number of processes");
89			i = 3;
90			} else
91			i = 1;
92	greg	1.1	/* set up rtrace command */
93	greg	2.20	for ( ; i < argc; i++) {
94	greg	2.6	if (argv[i][0] == '<' && argv[i][1] == '\0')
95	greg	1.3	break;
96	greg	1.1	rtargv[rtargc++] = argv[i];
97	greg	1.2	if (argv[i][0] == '-' && argv[i][1] == 'w')
98	greg	2.6	switch (argv[i][2]) {
99			case '\0':
100			warnings = !warnings;
101			break;
102			case '+':
103			case 'T': case 't':
104			case 'Y': case 'y':
105			case '1':
106			warnings = 1;
107			break;
108			case '-':
109			case 'F': case 'f':
110			case 'N': case 'n':
111			case '0':
112			warnings = 0;
113			break;
114			}
115	greg	1.2	}
116	greg	2.3	gargc = i;
117	greg	2.20	if (!strcmp(rtargv[--rtargc], "-defaults"))
118			nprocs = 0;
119			if (nprocs > 1) { /* add persist file if parallel invocation */
120			rtargv[rtargc++] = "-PP";
121			rtargv[rtargc++] = mktemp(persistfn);
122			}
123			/* add "mandatory" rtrace options */
124	greg	1.1	for (i = 0; myrtopts[i] != NULL; i++)
125			rtargv[rtargc++] = myrtopts[i];
126	greg	2.20	/* finally, put back final argument */
127	greg	1.3	rtargv[rtargc++] = argv[gargc-1];
128	greg	1.1	rtargv[rtargc] = NULL;
129	greg	2.20	if (!nprocs) { /* just asking for defaults? */
130	greg	1.9	printopts(); fflush(stdout);
131	greg	1.1	rtpath = getpath(rtargv[0], getenv("PATH"), X_OK);
132			if (rtpath == NULL) {
133			eputs(rtargv[0]);
134			eputs(": command not found\n");
135			exit(1);
136			}
137			execv(rtpath, rtargv);
138	greg	2.20	perror(rtpath); /* execv() should not return */
139	greg	1.1	exit(1);
140			}
141	greg	2.3	if (gargc < 2 \|\| argv[gargc-1][0] == '-')
142			error(USER, "missing octree argument");
143	greg	1.1	/* else initialize and run our calculation */
144	greg	2.20	init(nprocs);
145	greg	1.3	if (gargc+1 < argc)
146			for (i = gargc+1; i < argc; i++) {
147	greg	1.5	if ((fp = fopen(argv[i], "r")) == NULL) {
148	greg	1.3	sprintf(errmsg,
149			"cannot open scene file \"%s\"", argv[i]);
150			error(SYSTEM, errmsg);
151			}
152			filter(fp, argv[i]);
153			fclose(fp);
154			}
155			else
156			filter(stdin, "standard input");
157	greg	2.20	quit(0);
158	greg	1.1	}
159
160	greg	2.20	static void
161	greg	2.21	killpersist(void) /* kill persistent rtrace process */
162	greg	2.20	{
163			FILE *fp = fopen(persistfn, "r");
164			int pid;
165
166			if (fp == NULL)
167			return;
168			if (fscanf(fp, "%*s %d", &pid) != 1 \|\| kill(pid, SIGALRM) < 0)
169			unlink(persistfn);
170			fclose(fp);
171			}
172	greg	1.1
173	greg	2.23	int
174			done_rprocs(struct rtproc *rtp)
175			{
176			int st0, st1 = 0;
177
178			if (rtp->next != NULL) { /* close last opened first! */
179			st1 = done_rprocs(rtp->next);
180			free((void *)rtp->next);
181			rtp->next = NULL;
182			}
183			st0 = close_process(&rtp->pd);
184			if (st0 < 0)
185			error(WARNING, "unknown return status from rtrace process");
186			else if (st0 > 0)
187			return(st0);
188			return(st1);
189			}
190
191	greg	2.11	void
192	greg	2.20	quit(int status) /* exit with status */
193	greg	1.3	{
194			int rtstat;
195
196	greg	2.20	if (rt0.next != NULL) /* terminate persistent rtrace */
197			killpersist();
198			/* clean up rtrace process(es) */
199	greg	2.23	rtstat = done_rprocs(&rt0);
200			if (status == 0)
201			status = rtstat;
202	greg	1.3	exit(status);
203			}
204
205	schorsch	2.19	void
206	greg	2.20	init(int np) /* start rtrace and set up buffers */
207	greg	1.1	{
208	greg	2.20	struct rtproc *rtp;
209			int i;
210	greg	1.1	int maxbytes;
211	greg	1.3	/* set up object functions */
212			ofun[OBJ_FACE].funp = o_face;
213	greg	1.4	ofun[OBJ_SPHERE].funp = o_sphere;
214	greg	1.3	ofun[OBJ_RING].funp = o_ring;
215			/* set up signal handling */
216	schorsch	2.14	#ifdef SIGPIPE /* not present on Windows */
217	greg	1.3	signal(SIGPIPE, quit);
218	schorsch	2.13	#endif
219	greg	2.20	rtp = &rt0; /* start rtrace process(es) */
220			for (i = 0; i++ < np; ) {
221			errno = 0;
222			maxbytes = open_process(&rtp->pd, rtargv);
223			if (maxbytes == 0) {
224			eputs(rtargv[0]);
225			eputs(": command not found\n");
226			exit(1);
227			}
228			if (maxbytes < 0)
229			error(SYSTEM, "cannot start rtrace process");
230			rtp->bsiz = maxbytes/(6*sizeof(float));
231			rtp->buf = (float )malloc(6sizeof(float)*rtp->bsiz--);
232			rtp->dest = (float *)calloc(rtp->bsiz, sizeof(float ));
233			if (rtp->buf == NULL \|\| rtp->dest == NULL)
234			error(SYSTEM, "out of memory in init");
235			rtp->nrays = 0;
236			if (i == np) /* last process? */
237			break;
238	greg	2.24	if (np > 1)
239			sleep(2); /* wait for persist file */
240	greg	2.20	rtp->next = (struct rtproc *)malloc(sizeof(struct rtproc));
241			if (rtp->next == NULL)
242			error(SYSTEM, "out of memory in init");
243			rtp = rtp->next;
244	greg	1.1	}
245	greg	2.20	rtp->next = NULL;
246	greg	1.6	/* set up urand */
247	greg	2.7	initurand(16384);
248	greg	1.1	}
249
250
251	greg	2.11	void
252	schorsch	2.19	eputs( /* put string to stderr */
253			register char *s
254			)
255	greg	1.1	{
256			static int midline = 0;
257
258			if (!*s) return;
259			if (!midline) {
260			fputs(progname, stderr);
261			fputs(": ", stderr);
262			}
263			fputs(s, stderr);
264			midline = s[strlen(s)-1] != '\n';
265			}
266
267
268	greg	2.11	void
269	greg	1.2	wputs(s) /* print warning if enabled */
270			char *s;
271			{
272			if (warnings)
273			eputs(s);
274			}
275
276
277	schorsch	2.19	void
278			filter( /* process stream */
279			register FILE *infp,
280			char *name
281			)
282	greg	1.1	{
283			char buf[512];
284			FILE *pfp;
285			register int c;
286
287			while ((c = getc(infp)) != EOF) {
288			if (isspace(c))
289			continue;
290			if (c == '#') { /* comment/options */
291			buf[0] = c;
292			fgets(buf+1, sizeof(buf)-1, infp);
293	greg	1.2	xoptions(buf, name);
294	greg	1.1	} else if (c == '!') { /* command */
295			buf[0] = c;
296			fgetline(buf+1, sizeof(buf)-1, infp);
297			if ((pfp = popen(buf+1, "r")) == NULL) {
298			sprintf(errmsg, "cannot execute \"%s\"", buf);
299			error(SYSTEM, errmsg);
300			}
301	greg	1.2	filter(pfp, buf);
302	greg	1.1	pclose(pfp);
303			} else { /* object */
304			ungetc(c, infp);
305			xobject(infp, name);
306			}
307			}
308			}
309
310
311	schorsch	2.19	void
312			xoptions( /* process options in string s */
313			char *s,
314			char *nm
315			)
316	greg	1.1	{
317			extern FILE *freopen();
318			char buf[64];
319			int nerrs = 0;
320			register char *cp;
321
322	greg	1.2	if (strncmp(s, "#@mkillum", 9) \|\| !isspace(s[9])) {
323			fputs(s, stdout); /* not for us */
324	greg	1.1	return;
325	greg	1.2	}
326	greg	1.1	cp = s+10;
327			while (*cp) {
328			switch (*cp) {
329			case ' ':
330			case '\t':
331			case '\n':
332	greg	2.11	case '\r':
333			case '\f':
334	greg	1.1	cp++;
335			continue;
336			case 'm': /* material name */
337	greg	1.2	if (*++cp != '=')
338	greg	1.1	break;
339	greg	2.8	if (!++cp \|\| isspace(cp))
340	greg	1.2	break;
341	greg	1.1	atos(thisillum.matname, MAXSTR, cp);
342			cp = sskip(cp);
343			if (!(thisillum.flags & IL_DATCLB)) {
344			strcpy(thisillum.datafile, thisillum.matname);
345			thisillum.dfnum = 0;
346			}
347			continue;
348			case 'f': /* data file name */
349	greg	1.2	if (*++cp != '=')
350	greg	1.1	break;
351	greg	2.8	if (!++cp \|\| isspace(cp)) {
352	greg	1.6	strcpy(thisillum.datafile,thisillum.matname);
353	greg	1.10	thisillum.dfnum = 0;
354	greg	1.1	thisillum.flags &= ~IL_DATCLB;
355			continue;
356			}
357			atos(thisillum.datafile, MAXSTR, cp);
358			cp = sskip(cp);
359			thisillum.dfnum = 0;
360			thisillum.flags \|= IL_DATCLB;
361			continue;
362			case 'i': /* include material */
363			case 'e': /* exclude material */
364	greg	1.6	if (cp[1] != '=')
365	greg	1.1	break;
366	greg	1.7	matselect = (*cp == 'i') ? S_ELEM : S_COMPL;
367			cp += 2;
368			atos(matcheck, MAXSTR, cp);
369	greg	1.1	cp = sskip(cp);
370			continue;
371			case 'a': /* use everything */
372			cp = sskip(cp);
373			matselect = S_ALL;
374			continue;
375			case 'n': /* use nothing (passive) */
376			cp = sskip(cp);
377			matselect = S_NONE;
378			continue;
379			case 'c': /* color calculation */
380	greg	1.2	if (*++cp != '=')
381	greg	1.1	break;
382	greg	1.2	switch (*++cp) {
383	greg	1.1	case 'a': /* average */
384	greg	1.6	thisillum.flags = (thisillum.flags\|IL_COLAVG)
385			& ~IL_COLDST;
386	greg	1.1	break;
387			case 'd': /* distribution */
388	greg	1.5	thisillum.flags \|= (IL_COLDST\|IL_COLAVG);
389	greg	1.1	break;
390			case 'n': /* none */
391			thisillum.flags &= ~(IL_COLAVG\|IL_COLDST);
392			break;
393			default:
394			goto opterr;
395			}
396			cp = sskip(cp);
397			continue;
398			case 'd': /* point sample density */
399	greg	1.2	if (*++cp != '=')
400	greg	1.1	break;
401	greg	2.11	if (!isintd(++cp, " \t\n\r"))
402	greg	1.1	break;
403	greg	1.3	thisillum.sampdens = atoi(cp);
404	greg	1.1	cp = sskip(cp);
405			continue;
406			case 's': /* point super-samples */
407	greg	1.2	if (*++cp != '=')
408	greg	1.1	break;
409	greg	2.11	if (!isintd(++cp, " \t\n\r"))
410	greg	1.1	break;
411			thisillum.nsamps = atoi(cp);
412			cp = sskip(cp);
413			continue;
414	greg	1.3	case 'l': /* light sources */
415	greg	1.4	cp++;
416			if (*cp == '+')
417	greg	1.3	thisillum.flags \|= IL_LIGHT;
418	greg	1.4	else if (*cp == '-')
419			thisillum.flags &= ~IL_LIGHT;
420	greg	1.3	else
421	greg	1.4	break;
422			cp++;
423	greg	1.3	continue;
424	greg	1.12	case 'b': /* brightness */
425			if (*++cp != '=')
426			break;
427	greg	2.11	if (!isfltd(++cp, " \t\n\r"))
428	greg	1.12	break;
429			thisillum.minbrt = atof(cp);
430	greg	1.13	if (thisillum.minbrt < 0.)
431			thisillum.minbrt = 0.;
432	greg	1.12	cp = sskip(cp);
433			continue;
434	greg	1.1	case 'o': /* output file */
435	greg	1.2	if (*++cp != '=')
436	greg	1.1	break;
437	greg	2.8	if (!++cp \|\| isspace(cp))
438	greg	1.2	break;
439	greg	1.1	atos(buf, sizeof(buf), cp);
440			cp = sskip(cp);
441			if (freopen(buf, "w", stdout) == NULL) {
442			sprintf(errmsg,
443			"cannot open output file \"%s\"", buf);
444			error(SYSTEM, errmsg);
445			}
446			doneheader = 0;
447			continue;
448	greg	1.6	case '!': /* processed file! */
449	greg	1.7	sprintf(errmsg, "(%s): already processed!", nm);
450	greg	1.6	error(WARNING, errmsg);
451	greg	1.7	matselect = S_NONE;
452	greg	1.6	return;
453	greg	1.1	}
454			opterr: /* skip faulty option */
455	greg	2.9	while (cp && !isspace(cp))
456			cp++;
457	greg	1.1	nerrs++;
458			}
459	greg	1.2	/* print header? */
460			checkhead();
461			/* issue warnings? */
462	greg	1.1	if (nerrs) {
463			sprintf(errmsg, "(%s): %d error(s) in option line:",
464			nm, nerrs);
465			error(WARNING, errmsg);
466	greg	1.2	wputs(s);
467			printf("# %s: the following option line has %d error(s):\n",
468			progname, nerrs);
469	greg	1.1	}
470	greg	1.2	/* print pure comment */
471	greg	1.4	printf("# %s", s+2);
472	greg	1.9	}
473
474	schorsch	2.19	void
475			printopts(void) /* print out option default values */
476	greg	1.9	{
477	greg	1.11	printf("m=%-15s\t\t# material name\n", thisillum.matname);
478			printf("f=%-15s\t\t# data file name\n", thisillum.datafile);
479	greg	2.2	if (thisillum.flags & IL_COLAVG)
480			if (thisillum.flags & IL_COLDST)
481			printf("c=d\t\t\t\t# color distribution\n");
482			else
483			printf("c=a\t\t\t\t# color average\n");
484			else
485			printf("c=n\t\t\t\t# color none\n");
486			if (thisillum.flags & IL_LIGHT)
487			printf("l+\t\t\t\t# light type on\n");
488			else
489			printf("l-\t\t\t\t# light type off\n");
490	greg	1.9	printf("d=%d\t\t\t\t# density of points\n", thisillum.sampdens);
491			printf("s=%d\t\t\t\t# samples per point\n", thisillum.nsamps);
492	greg	1.12	printf("b=%f\t\t\t# minimum average brightness\n", thisillum.minbrt);
493	greg	1.1	}
494	greg	1.2
495
496	schorsch	2.19	void
497			printhead( /* print out header */
498			register int ac,
499			register char **av
500			)
501	greg	1.2	{
502			putchar('#');
503			while (ac-- > 0) {
504			putchar(' ');
505			fputs(*av++, stdout);
506			}
507	greg	1.6	fputs("\n#@mkillum !\n", stdout);
508	greg	1.2	}
509
510
511	schorsch	2.19	void
512			xobject( /* translate an object from fp */
513			FILE *fp,
514			char *nm
515			)
516	greg	1.2	{
517			OBJREC thisobj;
518			char str[MAXSTR];
519			int doit;
520			/* read the object */
521			if (fgetword(thisillum.altmat, MAXSTR, fp) == NULL)
522			goto readerr;
523			if (fgetword(str, MAXSTR, fp) == NULL)
524			goto readerr;
525			/* is it an alias? */
526	greg	2.12	if (!strcmp(str, ALIASKEY)) {
527	greg	1.2	if (fgetword(str, MAXSTR, fp) == NULL)
528			goto readerr;
529	greg	2.12	printf("\n%s %s %s", thisillum.altmat, ALIASKEY, str);
530	greg	1.2	if (fgetword(str, MAXSTR, fp) == NULL)
531			goto readerr;
532	greg	1.4	printf("\t%s\n", str);
533	greg	1.2	return;
534			}
535	greg	1.4	thisobj.omod = OVOID; /* unused field */
536	greg	1.2	if ((thisobj.otype = otype(str)) < 0) {
537			sprintf(errmsg, "(%s): unknown type \"%s\"", nm, str);
538			error(USER, errmsg);
539			}
540			if (fgetword(str, MAXSTR, fp) == NULL)
541			goto readerr;
542			thisobj.oname = str;
543			if (readfargs(&thisobj.oargs, fp) != 1)
544			goto readerr;
545			thisobj.os = NULL;
546			/* check for translation */
547			switch (matselect) {
548			case S_NONE:
549			doit = 0;
550			break;
551			case S_ALL:
552	greg	1.4	doit = 1;
553	greg	1.2	break;
554			case S_ELEM:
555			doit = !strcmp(thisillum.altmat, matcheck);
556			break;
557			case S_COMPL:
558			doit = strcmp(thisillum.altmat, matcheck);
559			break;
560			}
561	greg	1.4	doit = doit && issurface(thisobj.otype);
562	greg	1.2	/* print header? */
563			checkhead();
564			/* process object */
565			if (doit)
566	greg	2.20	(*ofun[thisobj.otype].funp)(&thisobj, &thisillum, &rt0, nm);
567	greg	1.2	else
568			printobj(thisillum.altmat, &thisobj);
569			/* free arguments */
570			freefargs(&thisobj.oargs);
571			return;
572			readerr:
573			sprintf(errmsg, "(%s): error reading scene", nm);
574			error(USER, errmsg);
575			}