src/gen/mkillum.c

#ifndef lint
static const char RCSid[] = "$Id: mkillum.c,v 2.21 2004/09/19 08:42:22 greg Exp $";
#endif
/*
 * Make illum sources for optimizing rendering process
 */

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

#include  "platform.h"
#include  "rtprocess.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? */

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);
}

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

        if (rt0.next != NULL)           /* terminate persistent rtrace */
                killpersist();
                                        /* clean up rtrace process(es) */
        for (rtp = &rt0; rtp != NULL; rtp = rtp->next) {
                rtstat = close_process(&rtp->pd);
                if (status == 0) {
                        if (rtstat < 0)
                                error(WARNING,
                                "unknown return status from rtrace process");
                        else
                                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");
                if (!i && np > 1)
                        sleep(2);       /* wait for persist file */
                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;
                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.22
Committed:	Sat Oct 23 18:55:52 2004 UTC (19 years, 6 months ago) by schorsch
Content type:	text/plain
Branch:	MAIN
Changes since 2.21:	+2 -1 lines
Log Message:	Compatibility fixes as suggested by Siegbert Debatin.
#	User	Rev	Content
1	greg	1.1	#ifndef lint
2	schorsch	2.22	static const char RCSid[] = "$Id: mkillum.c,v 2.21 2004/09/19 08:42:22 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	schorsch	2.22	#include "rtprocess.h"
13	greg	2.15	#include "mkillum.h"
14	schorsch	2.19	#include "random.h"
15	greg	1.1
16			/* default parameters */
17	greg	1.8	#define SAMPDENS 48 /* points per projected steradian */
18	greg	1.1	#define NSAMPS 32 /* samples per point */
19			#define DFLMAT "illum_mat" /* material name */
20	greg	1.10	#define DFLDAT "illum" /* data file name */
21	greg	1.1	/* selection options */
22			#define S_NONE 0 /* select none */
23			#define S_ELEM 1 /* select specified element */
24			#define S_COMPL 2 /* select all but element */
25			#define S_ALL 3 /* select all */
26
27			/* rtrace command and defaults */
28	greg	2.5	char *rtargv[64] = { "rtrace", "-dj", ".25", "-dr", "3", "-dv-",
29	greg	2.17	"-ab", "2", "-ad", "1024", "-as", "512", "-aa", ".1", };
30	greg	1.4	int rtargc = 14;
31	greg	1.1	/* overriding rtrace options */
32	gregl	2.10	char *myrtopts[] = { "-I-", "-i-", "-ld-", "-ov", "-h-",
33			"-fff", "-y", "0", NULL };
34	greg	1.1
35	greg	2.20	struct rtproc rt0; /* head of rtrace process list */
36	greg	1.1
37	greg	1.2	struct illum_args thisillum = { /* our illum and default values */
38			0,
39			DFLMAT,
40	greg	1.10	DFLDAT,
41	greg	1.2	0,
42			VOIDID,
43	greg	1.3	SAMPDENS,
44	greg	1.2	NSAMPS,
45	greg	1.12	0.,
46	greg	1.2	};
47	greg	1.1
48			char matcheck[MAXSTR]; /* current material to include or exclude */
49			int matselect = S_ALL; /* selection criterion */
50
51	greg	1.2	FUN ofun[NUMOTYPE] = INIT_OTYPE; /* object types */
52	greg	1.1
53	greg	2.20	char persistfn[] = "pfXXXXXX"; /* persist file name */
54
55	greg	1.1	int gargc; /* global argc */
56			char *gargv; / global argv */
57			#define progname gargv[0]
58
59			int doneheader = 0; /* printed header yet? */
60	greg	1.2	#define checkhead() if (!doneheader++) printhead(gargc,gargv)
61	greg	1.1
62	greg	1.2	int warnings = 1; /* print warnings? */
63
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.11	void
174	greg	2.20	quit(int status) /* exit with status */
175	greg	1.3	{
176	greg	2.20	struct rtproc *rtp;
177	greg	1.3	int rtstat;
178
179	greg	2.20	if (rt0.next != NULL) /* terminate persistent rtrace */
180			killpersist();
181			/* clean up rtrace process(es) */
182			for (rtp = &rt0; rtp != NULL; rtp = rtp->next) {
183			rtstat = close_process(&rtp->pd);
184			if (status == 0) {
185			if (rtstat < 0)
186			error(WARNING,
187			"unknown return status from rtrace process");
188			else
189			status = rtstat;
190			}
191	schorsch	2.16	}
192	greg	1.3	exit(status);
193			}
194
195	schorsch	2.19	void
196	greg	2.20	init(int np) /* start rtrace and set up buffers */
197	greg	1.1	{
198	greg	2.20	struct rtproc *rtp;
199			int i;
200	greg	1.1	int maxbytes;
201	greg	1.3	/* set up object functions */
202			ofun[OBJ_FACE].funp = o_face;
203	greg	1.4	ofun[OBJ_SPHERE].funp = o_sphere;
204	greg	1.3	ofun[OBJ_RING].funp = o_ring;
205			/* set up signal handling */
206	schorsch	2.14	#ifdef SIGPIPE /* not present on Windows */
207	greg	1.3	signal(SIGPIPE, quit);
208	schorsch	2.13	#endif
209	greg	2.20	rtp = &rt0; /* start rtrace process(es) */
210			for (i = 0; i++ < np; ) {
211			errno = 0;
212			maxbytes = open_process(&rtp->pd, rtargv);
213			if (maxbytes == 0) {
214			eputs(rtargv[0]);
215			eputs(": command not found\n");
216			exit(1);
217			}
218			if (maxbytes < 0)
219			error(SYSTEM, "cannot start rtrace process");
220			if (!i && np > 1)
221			sleep(2); /* wait for persist file */
222			rtp->bsiz = maxbytes/(6*sizeof(float));
223			rtp->buf = (float )malloc(6sizeof(float)*rtp->bsiz--);
224			rtp->dest = (float *)calloc(rtp->bsiz, sizeof(float ));
225			if (rtp->buf == NULL \|\| rtp->dest == NULL)
226			error(SYSTEM, "out of memory in init");
227			rtp->nrays = 0;
228			if (i == np) /* last process? */
229			break;
230			rtp->next = (struct rtproc *)malloc(sizeof(struct rtproc));
231			if (rtp->next == NULL)
232			error(SYSTEM, "out of memory in init");
233			rtp = rtp->next;
234	greg	1.1	}
235	greg	2.20	rtp->next = NULL;
236	greg	1.6	/* set up urand */
237	greg	2.7	initurand(16384);
238	greg	1.1	}
239
240
241	greg	2.11	void
242	schorsch	2.19	eputs( /* put string to stderr */
243			register char *s
244			)
245	greg	1.1	{
246			static int midline = 0;
247
248			if (!*s) return;
249			if (!midline) {
250			fputs(progname, stderr);
251			fputs(": ", stderr);
252			}
253			fputs(s, stderr);
254			midline = s[strlen(s)-1] != '\n';
255			}
256
257
258	greg	2.11	void
259	greg	1.2	wputs(s) /* print warning if enabled */
260			char *s;
261			{
262			if (warnings)
263			eputs(s);
264			}
265
266
267	schorsch	2.19	void
268			filter( /* process stream */
269			register FILE *infp,
270			char *name
271			)
272	greg	1.1	{
273			char buf[512];
274			FILE *pfp;
275			register int c;
276
277			while ((c = getc(infp)) != EOF) {
278			if (isspace(c))
279			continue;
280			if (c == '#') { /* comment/options */
281			buf[0] = c;
282			fgets(buf+1, sizeof(buf)-1, infp);
283	greg	1.2	xoptions(buf, name);
284	greg	1.1	} else if (c == '!') { /* command */
285			buf[0] = c;
286			fgetline(buf+1, sizeof(buf)-1, infp);
287			if ((pfp = popen(buf+1, "r")) == NULL) {
288			sprintf(errmsg, "cannot execute \"%s\"", buf);
289			error(SYSTEM, errmsg);
290			}
291	greg	1.2	filter(pfp, buf);
292	greg	1.1	pclose(pfp);
293			} else { /* object */
294			ungetc(c, infp);
295			xobject(infp, name);
296			}
297			}
298			}
299
300
301	schorsch	2.19	void
302			xoptions( /* process options in string s */
303			char *s,
304			char *nm
305			)
306	greg	1.1	{
307			extern FILE *freopen();
308			char buf[64];
309			int nerrs = 0;
310			register char *cp;
311
312	greg	1.2	if (strncmp(s, "#@mkillum", 9) \|\| !isspace(s[9])) {
313			fputs(s, stdout); /* not for us */
314	greg	1.1	return;
315	greg	1.2	}
316	greg	1.1	cp = s+10;
317			while (*cp) {
318			switch (*cp) {
319			case ' ':
320			case '\t':
321			case '\n':
322	greg	2.11	case '\r':
323			case '\f':
324	greg	1.1	cp++;
325			continue;
326			case 'm': /* material name */
327	greg	1.2	if (*++cp != '=')
328	greg	1.1	break;
329	greg	2.8	if (!++cp \|\| isspace(cp))
330	greg	1.2	break;
331	greg	1.1	atos(thisillum.matname, MAXSTR, cp);
332			cp = sskip(cp);
333			if (!(thisillum.flags & IL_DATCLB)) {
334			strcpy(thisillum.datafile, thisillum.matname);
335			thisillum.dfnum = 0;
336			}
337			continue;
338			case 'f': /* data file name */
339	greg	1.2	if (*++cp != '=')
340	greg	1.1	break;
341	greg	2.8	if (!++cp \|\| isspace(cp)) {
342	greg	1.6	strcpy(thisillum.datafile,thisillum.matname);
343	greg	1.10	thisillum.dfnum = 0;
344	greg	1.1	thisillum.flags &= ~IL_DATCLB;
345			continue;
346			}
347			atos(thisillum.datafile, MAXSTR, cp);
348			cp = sskip(cp);
349			thisillum.dfnum = 0;
350			thisillum.flags \|= IL_DATCLB;
351			continue;
352			case 'i': /* include material */
353			case 'e': /* exclude material */
354	greg	1.6	if (cp[1] != '=')
355	greg	1.1	break;
356	greg	1.7	matselect = (*cp == 'i') ? S_ELEM : S_COMPL;
357			cp += 2;
358			atos(matcheck, MAXSTR, cp);
359	greg	1.1	cp = sskip(cp);
360			continue;
361			case 'a': /* use everything */
362			cp = sskip(cp);
363			matselect = S_ALL;
364			continue;
365			case 'n': /* use nothing (passive) */
366			cp = sskip(cp);
367			matselect = S_NONE;
368			continue;
369			case 'c': /* color calculation */
370	greg	1.2	if (*++cp != '=')
371	greg	1.1	break;
372	greg	1.2	switch (*++cp) {
373	greg	1.1	case 'a': /* average */
374	greg	1.6	thisillum.flags = (thisillum.flags\|IL_COLAVG)
375			& ~IL_COLDST;
376	greg	1.1	break;
377			case 'd': /* distribution */
378	greg	1.5	thisillum.flags \|= (IL_COLDST\|IL_COLAVG);
379	greg	1.1	break;
380			case 'n': /* none */
381			thisillum.flags &= ~(IL_COLAVG\|IL_COLDST);
382			break;
383			default:
384			goto opterr;
385			}
386			cp = sskip(cp);
387			continue;
388			case 'd': /* point sample density */
389	greg	1.2	if (*++cp != '=')
390	greg	1.1	break;
391	greg	2.11	if (!isintd(++cp, " \t\n\r"))
392	greg	1.1	break;
393	greg	1.3	thisillum.sampdens = atoi(cp);
394	greg	1.1	cp = sskip(cp);
395			continue;
396			case 's': /* point super-samples */
397	greg	1.2	if (*++cp != '=')
398	greg	1.1	break;
399	greg	2.11	if (!isintd(++cp, " \t\n\r"))
400	greg	1.1	break;
401			thisillum.nsamps = atoi(cp);
402			cp = sskip(cp);
403			continue;
404	greg	1.3	case 'l': /* light sources */
405	greg	1.4	cp++;
406			if (*cp == '+')
407	greg	1.3	thisillum.flags \|= IL_LIGHT;
408	greg	1.4	else if (*cp == '-')
409			thisillum.flags &= ~IL_LIGHT;
410	greg	1.3	else
411	greg	1.4	break;
412			cp++;
413	greg	1.3	continue;
414	greg	1.12	case 'b': /* brightness */
415			if (*++cp != '=')
416			break;
417	greg	2.11	if (!isfltd(++cp, " \t\n\r"))
418	greg	1.12	break;
419			thisillum.minbrt = atof(cp);
420	greg	1.13	if (thisillum.minbrt < 0.)
421			thisillum.minbrt = 0.;
422	greg	1.12	cp = sskip(cp);
423			continue;
424	greg	1.1	case 'o': /* output file */
425	greg	1.2	if (*++cp != '=')
426	greg	1.1	break;
427	greg	2.8	if (!++cp \|\| isspace(cp))
428	greg	1.2	break;
429	greg	1.1	atos(buf, sizeof(buf), cp);
430			cp = sskip(cp);
431			if (freopen(buf, "w", stdout) == NULL) {
432			sprintf(errmsg,
433			"cannot open output file \"%s\"", buf);
434			error(SYSTEM, errmsg);
435			}
436			doneheader = 0;
437			continue;
438	greg	1.6	case '!': /* processed file! */
439	greg	1.7	sprintf(errmsg, "(%s): already processed!", nm);
440	greg	1.6	error(WARNING, errmsg);
441	greg	1.7	matselect = S_NONE;
442	greg	1.6	return;
443	greg	1.1	}
444			opterr: /* skip faulty option */
445	greg	2.9	while (cp && !isspace(cp))
446			cp++;
447	greg	1.1	nerrs++;
448			}
449	greg	1.2	/* print header? */
450			checkhead();
451			/* issue warnings? */
452	greg	1.1	if (nerrs) {
453			sprintf(errmsg, "(%s): %d error(s) in option line:",
454			nm, nerrs);
455			error(WARNING, errmsg);
456	greg	1.2	wputs(s);
457			printf("# %s: the following option line has %d error(s):\n",
458			progname, nerrs);
459	greg	1.1	}
460	greg	1.2	/* print pure comment */
461	greg	1.4	printf("# %s", s+2);
462	greg	1.9	}
463
464	schorsch	2.19	void
465			printopts(void) /* print out option default values */
466	greg	1.9	{
467	greg	1.11	printf("m=%-15s\t\t# material name\n", thisillum.matname);
468			printf("f=%-15s\t\t# data file name\n", thisillum.datafile);
469	greg	2.2	if (thisillum.flags & IL_COLAVG)
470			if (thisillum.flags & IL_COLDST)
471			printf("c=d\t\t\t\t# color distribution\n");
472			else
473			printf("c=a\t\t\t\t# color average\n");
474			else
475			printf("c=n\t\t\t\t# color none\n");
476			if (thisillum.flags & IL_LIGHT)
477			printf("l+\t\t\t\t# light type on\n");
478			else
479			printf("l-\t\t\t\t# light type off\n");
480	greg	1.9	printf("d=%d\t\t\t\t# density of points\n", thisillum.sampdens);
481			printf("s=%d\t\t\t\t# samples per point\n", thisillum.nsamps);
482	greg	1.12	printf("b=%f\t\t\t# minimum average brightness\n", thisillum.minbrt);
483	greg	1.1	}
484	greg	1.2
485
486	schorsch	2.19	void
487			printhead( /* print out header */
488			register int ac,
489			register char **av
490			)
491	greg	1.2	{
492			putchar('#');
493			while (ac-- > 0) {
494			putchar(' ');
495			fputs(*av++, stdout);
496			}
497	greg	1.6	fputs("\n#@mkillum !\n", stdout);
498	greg	1.2	}
499
500
501	schorsch	2.19	void
502			xobject( /* translate an object from fp */
503			FILE *fp,
504			char *nm
505			)
506	greg	1.2	{
507			OBJREC thisobj;
508			char str[MAXSTR];
509			int doit;
510			/* read the object */
511			if (fgetword(thisillum.altmat, MAXSTR, fp) == NULL)
512			goto readerr;
513			if (fgetword(str, MAXSTR, fp) == NULL)
514			goto readerr;
515			/* is it an alias? */
516	greg	2.12	if (!strcmp(str, ALIASKEY)) {
517	greg	1.2	if (fgetword(str, MAXSTR, fp) == NULL)
518			goto readerr;
519	greg	2.12	printf("\n%s %s %s", thisillum.altmat, ALIASKEY, str);
520	greg	1.2	if (fgetword(str, MAXSTR, fp) == NULL)
521			goto readerr;
522	greg	1.4	printf("\t%s\n", str);
523	greg	1.2	return;
524			}
525	greg	1.4	thisobj.omod = OVOID; /* unused field */
526	greg	1.2	if ((thisobj.otype = otype(str)) < 0) {
527			sprintf(errmsg, "(%s): unknown type \"%s\"", nm, str);
528			error(USER, errmsg);
529			}
530			if (fgetword(str, MAXSTR, fp) == NULL)
531			goto readerr;
532			thisobj.oname = str;
533			if (readfargs(&thisobj.oargs, fp) != 1)
534			goto readerr;
535			thisobj.os = NULL;
536			/* check for translation */
537			switch (matselect) {
538			case S_NONE:
539			doit = 0;
540			break;
541			case S_ALL:
542	greg	1.4	doit = 1;
543	greg	1.2	break;
544			case S_ELEM:
545			doit = !strcmp(thisillum.altmat, matcheck);
546			break;
547			case S_COMPL:
548			doit = strcmp(thisillum.altmat, matcheck);
549			break;
550			}
551	greg	1.4	doit = doit && issurface(thisobj.otype);
552	greg	1.2	/* print header? */
553			checkhead();
554			/* process object */
555			if (doit)
556	greg	2.20	(*ofun[thisobj.otype].funp)(&thisobj, &thisillum, &rt0, nm);
557	greg	1.2	else
558			printobj(thisillum.altmat, &thisobj);
559			/* free arguments */
560			freefargs(&thisobj.oargs);
561			return;
562			readerr:
563			sprintf(errmsg, "(%s): error reading scene", nm);
564			error(USER, errmsg);
565			}