ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cv/rad2mgf.c
Revision: 2.11
Committed: Thu Apr 13 15:33:48 1995 UTC (29 years ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.10: +36 -15 lines
Log Message:
added commenting of unsupported Radiance primitives

File Contents

# User Rev Content
1 greg 2.1 /* Copyright (c) 1994 Regents of the University of California */
2    
3     #ifndef lint
4     static char SCCSid[] = "$SunId$ LBL";
5     #endif
6    
7     /*
8     * Convert Radiance scene description to MGF
9     */
10    
11     #include <stdio.h>
12 greg 2.10 #include <math.h>
13 greg 2.1 #include <string.h>
14     #include "fvect.h"
15     #include "object.h"
16     #include "color.h"
17     #include "lookup.h"
18    
19 greg 2.7 #define PI 3.14159265358979323846
20    
21 greg 2.10 #define C_1SIDEDTHICK 0.005
22    
23 greg 2.1 int o_face(), o_cone(), o_sphere(), o_ring(), o_cylinder();
24 greg 2.11 int o_instance(), o_illum();
25 greg 2.10 int o_plastic(), o_metal(), o_glass(), o_dielectric(),
26 greg 2.11 o_mirror(), o_trans(), o_light();
27 greg 2.1
28     extern void free();
29     extern char *malloc();
30    
31     LUTAB rmats = LU_SINIT(free,NULL); /* defined material table */
32    
33     LUTAB rdispatch = LU_SINIT(NULL,NULL); /* function dispatch table */
34    
35 greg 2.2 char curmat[80]; /* current material */
36     char curobj[128] = "Untitled"; /* current object name */
37 greg 2.1
38 greg 2.2 double unit_mult = 1.; /* units multiplier */
39 greg 2.1
40 greg 2.2 #define hasmult (unit_mult < .999 || unit_mult > 1.001)
41 greg 2.1
42 greg 2.4 /*
43     * Stuff for tracking and reusing vertices:
44     */
45 greg 2.2
46 greg 2.9 char VKFMT[] = "%+16.9e %+16.9e %+16.9e";
47 greg 2.4 #define VKLEN 64
48    
49     #define mkvkey(k,v) sprintf(k, VKFMT, (v)[0], (v)[1], (v)[2])
50    
51     #define NVERTS 256
52    
53 greg 2.8 long vclock; /* incremented at each vertex request */
54 greg 2.4
55     struct vert {
56     long lused; /* when last used (0 if unassigned) */
57     FVECT p; /* track point position only */
58     } vert[NVERTS]; /* our vertex cache */
59    
60     LUTAB vertab = LU_SINIT(free,NULL); /* our vertex lookup table */
61    
62    
63 greg 2.1 main(argc, argv)
64     int argc;
65     char **argv;
66     {
67     int i;
68    
69     for (i = 1; i < argc && argv[i][0] == '-'; i++)
70     switch (argv[i][1]) {
71     case 'd': /* units */
72     switch (argv[i][2]) {
73     case 'm': /* meters */
74     unit_mult = 1.;
75     break;
76     case 'c': /* centimeters */
77     unit_mult = .01;
78     break;
79     case 'f': /* feet */
80     unit_mult = 12.*.0254;
81     break;
82     case 'i': /* inches */
83     unit_mult = .0254;
84     break;
85     default:
86     goto unkopt;
87     }
88     break;
89 greg 2.6 default:
90     goto unkopt;
91 greg 2.1 }
92     init();
93     if (i >= argc)
94     rad2mgf(NULL);
95     else
96     for ( ; i < argc; i++)
97     rad2mgf(argv[i]);
98     uninit();
99     exit(0);
100     unkopt:
101     fprintf(stderr, "Usage: %s [-d{m|c|f|i}] file ..\n", argv[0]);
102     exit(1);
103     }
104    
105    
106     rad2mgf(inp) /* convert a Radiance file to MGF */
107     char *inp;
108     {
109     #define mod buf
110     #define typ (buf+128)
111     #define id (buf+256)
112     #define alias (buf+384)
113     char buf[512];
114     FUNARGS fa;
115     register FILE *fp;
116     register int c;
117    
118     if (inp == NULL) {
119 greg 2.2 inp = "standard input";
120 greg 2.1 fp = stdin;
121     } else if (inp[0] == '!') {
122     if ((fp = popen(inp+1, "r")) == NULL) {
123     fputs(inp, stderr);
124     fputs(": cannot execute\n", stderr);
125     exit(1);
126     }
127     } else if ((fp = fopen(inp, "r")) == NULL) {
128     fputs(inp, stderr);
129     fputs(": cannot open\n", stderr);
130     exit(1);
131     }
132     printf("# Begin conversion from: %s\n", inp);
133     while ((c = getc(fp)) != EOF)
134     switch (c) {
135     case ' ': /* white space */
136     case '\t':
137     case '\n':
138     case '\r':
139     case '\f':
140     break;
141     case '#': /* comment */
142     if (fgets(buf, sizeof(buf), fp) != NULL)
143     printf("# %s", buf);
144     break;
145     case '!': /* inline command */
146     ungetc(c, fp);
147     fgetline(buf, sizeof(buf), fp);
148     rad2mgf(buf);
149     break;
150     default: /* Radiance primitive */
151     ungetc(c, fp);
152     if (fscanf(fp, "%s %s %s", mod, typ, id) != 3) {
153     fputs(inp, stderr);
154     fputs(": unexpected EOF\n", stderr);
155     exit(1);
156     }
157     if (!strcmp(typ, "alias")) {
158     strcpy(alias, "EOF");
159     fscanf(fp, "%s", alias);
160     newmat(id, alias);
161     } else {
162     if (!readfargs(&fa, fp)) {
163     fprintf(stderr,
164     "%s: bad argument syntax for %s \"%s\"\n",
165     inp, typ, id);
166     exit(1);
167     }
168     cvtprim(inp, mod, typ, id, &fa);
169     freefargs(&fa);
170     }
171     break;
172     }
173     printf("# End conversion from: %s\n", inp);
174     if (inp[0] == '!')
175     pclose(fp);
176     else
177     fclose(fp);
178     #undef mod
179     #undef typ
180     #undef id
181     #undef alias
182     }
183    
184    
185     cvtprim(inp, mod, typ, id, fa) /* process Radiance primitive */
186     char *inp, *mod, *typ, *id;
187     FUNARGS *fa;
188     {
189     int (*df)();
190    
191     df = (int (*)())lu_find(&rdispatch, typ)->data;
192     if (df != NULL) { /* convert */
193     if ((*df)(mod, typ, id, fa) < 0) {
194     fprintf(stderr, "%s: bad %s \"%s\"\n", typ, id);
195     exit(1);
196     }
197 greg 2.11 } else { /* unsupported */
198     o_unsupported(mod, typ, id, fa);
199     if (lu_find(&rmats, mod)->data != NULL) /* make alias */
200     newmat(id, mod);
201     }
202 greg 2.1 }
203    
204    
205     newmat(id, alias) /* add a modifier to the alias list */
206     char *id;
207     char *alias;
208     {
209     register LUENT *lp, *lpa;
210    
211     if (alias != NULL) { /* look up alias */
212     if ((lpa = lu_find(&rmats, alias)) == NULL)
213     goto memerr;
214     if (lpa->data == NULL)
215     alias = NULL; /* doesn't exist! */
216     }
217     if ((lp = lu_find(&rmats, id)) == NULL) /* look up material */
218     goto memerr;
219     if (alias != NULL && lp->data == lpa->key)
220     return; /* alias set already */
221     if (lp->data == NULL) { /* allocate material */
222     if ((lp->key = (char *)malloc(strlen(id)+1)) == NULL)
223     goto memerr;
224     strcpy(lp->key, id);
225     }
226     if (alias == NULL) { /* set this material */
227     lp->data = lp->key;
228     printf("m %s =\n", id);
229     } else { /* set this alias */
230     lp->data = lpa->key;
231     printf("m %s = %s\n", id, alias);
232     }
233     strcpy(curmat, id);
234     return;
235     memerr:
236     fputs("Out of memory in newmat!\n", stderr);
237     exit(1);
238     }
239    
240    
241     setmat(id) /* set material to this one */
242     char *id;
243     {
244     if (!strcmp(id, curmat)) /* already set? */
245     return;
246 greg 2.2 if (!strcmp(id, VOIDID)) /* cannot set */
247     return;
248 greg 2.1 printf("m %s\n", id);
249     strcpy(curmat, id);
250     }
251    
252    
253 greg 2.2 setobj(id) /* set object name to this one */
254     char *id;
255     {
256     register char *cp, *cp2;
257     char *end = NULL;
258     int diff = 0;
259     /* use all but final suffix */
260     for (cp = id; *cp; cp++)
261     if (*cp == '.')
262     end = cp;
263     if (end == NULL)
264     end = cp;
265     /* copy to current object */
266     for (cp = id, cp2 = curobj; cp < end; *cp2++ = *cp++)
267     diff += *cp != *cp2;
268     if (!diff && !*cp2)
269     return;
270     *cp2 = '\0';
271     fputs("o\no ", stdout);
272     puts(curobj);
273     }
274    
275    
276 greg 2.1 init() /* initialize dispatch table and output */
277     {
278 greg 2.4 lu_init(&vertab, NVERTS);
279 greg 2.1 lu_init(&rdispatch, 22);
280     add2dispatch("polygon", o_face);
281     add2dispatch("cone", o_cone);
282     add2dispatch("cup", o_cone);
283     add2dispatch("sphere", o_sphere);
284     add2dispatch("bubble", o_sphere);
285     add2dispatch("cylinder", o_cylinder);
286     add2dispatch("tube", o_cylinder);
287     add2dispatch("ring", o_ring);
288     add2dispatch("instance", o_instance);
289     add2dispatch("plastic", o_plastic);
290     add2dispatch("plastic2", o_plastic);
291     add2dispatch("metal", o_metal);
292     add2dispatch("metal2", o_metal);
293     add2dispatch("glass", o_glass);
294 greg 2.10 add2dispatch("dielectric", o_dielectric);
295 greg 2.1 add2dispatch("trans", o_trans);
296     add2dispatch("trans2", o_trans);
297     add2dispatch("mirror", o_mirror);
298     add2dispatch("light", o_light);
299     add2dispatch("spotlight", o_light);
300     add2dispatch("glow", o_light);
301     add2dispatch("illum", o_illum);
302 greg 2.11 puts("# The following was converted from RADIANCE scene input");
303 greg 2.2 if (hasmult)
304 greg 2.1 printf("xf -s %.4e\n", unit_mult);
305 greg 2.2 printf("o %s\n", curobj);
306 greg 2.1 }
307    
308    
309     uninit() /* mark end of MGF file */
310     {
311 greg 2.2 puts("o");
312     if (hasmult)
313 greg 2.1 puts("xf");
314 greg 2.11 puts("# End of data converted from RADIANCE scene input");
315 greg 2.1 lu_done(&rdispatch);
316     lu_done(&rmats);
317 greg 2.4 lu_done(&vertab);
318 greg 2.1 }
319    
320    
321 greg 2.5 clrverts() /* clear vertex table */
322     {
323     register int i;
324    
325     lu_done(&vertab);
326     for (i = 0; i < NVERTS; i++)
327     vert[i].lused = 0;
328     lu_init(&vertab, NVERTS);
329     }
330    
331    
332 greg 2.1 add2dispatch(name, func) /* add function to dispatch table */
333     char *name;
334     int (*func)();
335     {
336     register LUENT *lp;
337    
338     lp = lu_find(&rdispatch, name);
339     if (lp->key != NULL) {
340     fputs(name, stderr);
341     fputs(": duplicate dispatch entry!\n", stderr);
342     exit(1);
343     }
344     lp->key = name;
345     lp->data = (char *)func;
346     }
347    
348    
349     char *
350 greg 2.2 getvertid(vname, vp) /* get/set vertex ID for this point */
351     char *vname;
352 greg 2.1 FVECT vp;
353     {
354 greg 2.9 static char vkey[VKLEN];
355 greg 2.1 register LUENT *lp;
356     register int i, vndx;
357    
358 greg 2.8 vclock++; /* increment counter */
359 greg 2.1 mkvkey(vkey, vp);
360     if ((lp = lu_find(&vertab, vkey)) == NULL)
361     goto memerr;
362     if (lp->data == NULL) { /* allocate new vertex entry */
363     if (lp->key != NULL) /* reclaim deleted entry */
364     vertab.ndel--;
365     else {
366     if ((lp->key = (char *)malloc(VKLEN)) == NULL)
367     goto memerr;
368     strcpy(lp->key, vkey);
369     }
370     vndx = 0; /* find oldest vertex */
371     for (i = 1; i < NVERTS; i++)
372     if (vert[i].lused < vert[vndx].lused)
373     vndx = i;
374     if (vert[vndx].lused) { /* free old entry first */
375     mkvkey(vkey, vert[vndx].p);
376     lu_delete(&vertab, vkey);
377     }
378 greg 2.3 VCOPY(vert[vndx].p, vp); /* assign it */
379 greg 2.2 printf("v v%d =\n\tp %.15g %.15g %.15g\n", /* print it */
380 greg 2.1 vndx, vp[0], vp[1], vp[2]);
381     lp->data = (char *)&vert[vndx]; /* set it */
382     } else
383     vndx = (struct vert *)lp->data - vert;
384 greg 2.8 vert[vndx].lused = vclock; /* record this use */
385 greg 2.1 sprintf(vname, "v%d", vndx);
386     return(vname);
387     memerr:
388     fputs("Out of memory in getvertid!\n", stderr);
389     exit(1);
390     }
391    
392    
393     int
394 greg 2.11 o_unsupported(mod, typ, id, fa) /* mark unsupported primitive */
395     char *mod, *typ, *id;
396     FUNARGS *fa;
397     {
398     register int i;
399    
400     fputs("\n# Unsupported RADIANCE primitive:\n", stdout);
401     printf("# %s %s %s", mod, typ, id);
402     printf("\n# %d", fa->nsargs);
403     for (i = 0; i < fa->nsargs; i++)
404     printf(" %s", fa->sarg[i]);
405     #ifdef IARGS
406     printf("\n# %d", fa->niargs);
407     for (i = 0; i < fa->niargs; i++)
408     printf(" %ld", fa->iarg[i]);
409     #else
410     fputs("\n# 0", stdout);
411     #endif
412     printf("\n# %d", fa->nfargs);
413     for (i = 0; i < fa->nfargs; i++)
414     printf(" %g", fa->farg[i]);
415     fputs("\n\n", stdout);
416     return(0);
417     }
418    
419    
420     int
421 greg 2.1 o_face(mod, typ, id, fa) /* print out a polygon */
422     char *mod, *typ, *id;
423     FUNARGS *fa;
424     {
425     char entbuf[512];
426 greg 2.2 register char *cp;
427 greg 2.1 register int i;
428    
429     if (fa->nfargs < 9 | fa->nfargs % 3)
430     return(-1);
431     setmat(mod);
432 greg 2.2 setobj(id);
433     cp = entbuf;
434     *cp++ = 'f';
435 greg 2.1 for (i = 0; i < fa->nfargs; i += 3) {
436 greg 2.2 *cp++ = ' ';
437     getvertid(cp, fa->farg + i);
438     while (*cp)
439     cp++;
440 greg 2.1 }
441     puts(entbuf);
442     return(0);
443     }
444    
445    
446     int
447     o_cone(mod, typ, id, fa) /* print out a cone */
448     char *mod, *typ, *id;
449     register FUNARGS *fa;
450     {
451 greg 2.2 char v1[6], v2[6];
452    
453 greg 2.1 if (fa->nfargs != 8)
454     return(-1);
455     setmat(mod);
456 greg 2.2 setobj(id);
457     getvertid(v1, fa->farg);
458     getvertid(v2, fa->farg + 3);
459 greg 2.1 if (typ[1] == 'u') /* cup -> inverted cone */
460 greg 2.2 printf("cone %s %.12g %s %.12g\n",
461     v1, -fa->farg[6], v2, -fa->farg[7]);
462 greg 2.1 else
463 greg 2.2 printf("cone %s %.12g %s %.12g\n",
464     v1, fa->farg[6], v2, fa->farg[7]);
465 greg 2.1 return(0);
466     }
467    
468    
469     int
470     o_sphere(mod, typ, id, fa) /* print out a sphere */
471     char *mod, *typ, *id;
472     register FUNARGS *fa;
473     {
474 greg 2.2 char cent[6];
475    
476 greg 2.1 if (fa->nfargs != 4)
477     return(-1);
478     setmat(mod);
479 greg 2.2 setobj(id);
480     printf("sph %s %.12g\n", getvertid(cent, fa->farg),
481     typ[0]=='b' ? -fa->farg[3] : fa->farg[3]);
482 greg 2.1 return(0);
483     }
484    
485    
486     int
487     o_cylinder(mod, typ, id, fa) /* print out a cylinder */
488     char *mod, *typ, *id;
489     register FUNARGS *fa;
490     {
491 greg 2.2 char v1[6], v2[6];
492    
493 greg 2.1 if (fa->nfargs != 7)
494     return(-1);
495     setmat(mod);
496 greg 2.2 setobj(id);
497     getvertid(v1, fa->farg);
498     getvertid(v2, fa->farg + 3);
499     printf("cyl %s %.12g %s\n", v1,
500     typ[0]=='t' ? -fa->farg[6] : fa->farg[6], v2);
501 greg 2.1 return(0);
502     }
503    
504    
505     int
506     o_ring(mod, typ, id, fa) /* print out a ring */
507     char *mod, *typ, *id;
508     register FUNARGS *fa;
509     {
510     if (fa->nfargs != 8)
511     return(-1);
512     setmat(mod);
513 greg 2.2 setobj(id);
514     printf("v cent =\n\tp %.12g %.12g %.12g\n",
515 greg 2.1 fa->farg[0], fa->farg[1], fa->farg[2]);
516 greg 2.2 printf("\tn %.12g %.12g %.12g\n",
517 greg 2.1 fa->farg[3], fa->farg[4], fa->farg[5]);
518     if (fa->farg[6] < fa->farg[7])
519     printf("ring cent %.12g %.12g\n",
520     fa->farg[6], fa->farg[7]);
521     else
522     printf("ring cent %.12g %.12g\n",
523     fa->farg[7], fa->farg[6]);
524     return(0);
525     }
526    
527    
528     int
529     o_instance(mod, typ, id, fa) /* convert an instance */
530     char *mod, *typ, *id;
531     FUNARGS *fa;
532     {
533 greg 2.2 register int i;
534     register char *cp;
535     char *start = NULL, *end = NULL;
536     /*
537     * We don't really know how to do this, so we just create
538     * a reference to an undefined MGF file and it's the user's
539     * responsibility to create this file and put the appropriate
540     * stuff into it.
541     */
542     if (fa->nsargs < 1)
543     return(-1);
544     setmat(mod); /* only works if surfaces are void */
545     setobj(id);
546     for (cp = fa->sarg[0]; *cp; cp++) /* construct MGF file name */
547     if (*cp == '/')
548     start = cp+1;
549     else if (*cp == '.')
550     end = cp;
551     if (start == NULL)
552     start = fa->sarg[0];
553     if (end == NULL || start >= end)
554     end = cp;
555     fputs("i ", stdout); /* print include entity */
556     for (cp = start; cp < end; cp++)
557     putchar(*cp);
558     fputs(".mgf", stdout); /* add MGF suffix */
559     for (i = 1; i < fa->nsargs; i++) { /* add transform */
560     putchar(' ');
561     fputs(fa->sarg[i], stdout);
562     }
563     putchar('\n');
564 greg 2.5 clrverts(); /* vertex id's no longer reliable */
565 greg 2.2 return(0);
566 greg 2.1 }
567    
568    
569     int
570     o_illum(mod, typ, id, fa) /* convert an illum material */
571     char *mod, *typ, *id;
572     FUNARGS *fa;
573     {
574     if (fa->nsargs == 1 && strcmp(fa->sarg[0], VOIDID)) {
575     newmat(id, fa->sarg[0]); /* just create alias */
576     return(0);
577     }
578     /* else create invisible material */
579     newmat(id, NULL);
580 greg 2.2 puts("\tts 1 0");
581 greg 2.1 return(0);
582     }
583    
584    
585     int
586     o_plastic(mod, typ, id, fa) /* convert a plastic material */
587     char *mod, *typ, *id;
588     register FUNARGS *fa;
589     {
590     COLOR cxyz, rrgb;
591     double d;
592    
593     if (fa->nfargs != (typ[7]=='2' ? 6 : 5))
594     return(-1);
595     newmat(id, NULL);
596     rrgb[0] = fa->farg[0]; rrgb[1] = fa->farg[1]; rrgb[2] = fa->farg[2];
597     rgb_cie(cxyz, rrgb);
598 greg 2.2 puts("\tc"); /* put diffuse component */
599 greg 2.1 d = cxyz[0] + cxyz[1] + cxyz[2];
600     if (d > FTINY)
601 greg 2.2 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
602     printf("\trd %.4f\n", cxyz[1]*(1. - fa->farg[3]));
603     if (fa->farg[3] > FTINY) { /* put specular component */
604     puts("\tc");
605     printf("\trs %.4f %.4f\n", fa->farg[3],
606     typ[7]=='2' ? .5*(fa->farg[4] + fa->farg[5]) :
607     fa->farg[4]);
608     }
609 greg 2.1 return(0);
610     }
611    
612    
613     int
614     o_metal(mod, typ, id, fa) /* convert a metal material */
615     char *mod, *typ, *id;
616     register FUNARGS *fa;
617     {
618     COLOR cxyz, rrgb;
619     double d;
620    
621     if (fa->nfargs != (typ[5]=='2' ? 6 : 5))
622     return(-1);
623     newmat(id, NULL);
624     rrgb[0] = fa->farg[0]; rrgb[1] = fa->farg[1]; rrgb[2] = fa->farg[2];
625     rgb_cie(cxyz, rrgb);
626 greg 2.2 puts("\tc"); /* put diffuse component */
627 greg 2.1 d = cxyz[0] + cxyz[1] + cxyz[2];
628     if (d > FTINY)
629 greg 2.2 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
630     printf("\trd %.4f\n", cxyz[1]*(1. - fa->farg[3]));
631 greg 2.1 /* put specular component */
632 greg 2.2 printf("\trs %.4f %.4f\n", cxyz[1]*fa->farg[3],
633 greg 2.1 typ[5]=='2' ? .5*(fa->farg[4] + fa->farg[5]) :
634     fa->farg[4]);
635     return(0);
636     }
637    
638    
639     int
640     o_glass(mod, typ, id, fa) /* convert a glass material */
641     char *mod, *typ, *id;
642     register FUNARGS *fa;
643     {
644     COLOR cxyz, rrgb, trgb;
645     double nrfr = 1.52, F, d;
646     register int i;
647    
648     if (fa->nfargs != 3 && fa->nfargs != 4)
649     return(-1);
650     newmat(id, NULL);
651     if (fa->nfargs == 4)
652     nrfr = fa->farg[3];
653 greg 2.10 printf("\tir %f 0\n", nrfr);
654 greg 2.1 F = (1. - nrfr)/(1. + nrfr); /* use normal incidence */
655     F *= F;
656     for (i = 0; i < 3; i++) {
657 greg 2.2 trgb[i] = fa->farg[i] * (1. - F)*(1. - F) /
658 greg 2.1 (1. - F*F*fa->farg[i]*fa->farg[i]);
659 greg 2.2 rrgb[i] = F * (1. + (1. - 2.*F)*fa->farg[i]) /
660     (1. - F*F*fa->farg[i]*fa->farg[i]);
661 greg 2.1 }
662     rgb_cie(cxyz, rrgb); /* put reflected component */
663 greg 2.2 puts("\tc");
664 greg 2.1 d = cxyz[0] + cxyz[1] + cxyz[2];
665     if (d > FTINY)
666 greg 2.2 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
667     printf("\trs %.4f 0\n", cxyz[1]);
668 greg 2.10 rgb_cie(cxyz, trgb); /* put transmitted component */
669     puts("\tc");
670     d = cxyz[0] + cxyz[1] + cxyz[2];
671     if (d > FTINY)
672     printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
673     printf("\tts %.4f 0\n", cxyz[1]);
674     return(0);
675     }
676    
677    
678     int
679     o_dielectric(mod, typ, id, fa) /* convert a dielectric material */
680     char *mod, *typ, *id;
681     register FUNARGS *fa;
682     {
683     COLOR cxyz, trgb;
684     double F, d;
685     register int i;
686    
687     if (fa->nfargs != 5)
688     return(-1);
689     newmat(id, NULL);
690     F = (1. - fa->farg[3])/(1. + fa->farg[3]); /* normal incidence */
691     F *= F;
692     for (i = 0; i < 3; i++)
693     trgb[i] = (1. - F)*pow(fa->farg[i], C_1SIDEDTHICK/unit_mult);
694     printf("\tir %f 0\n", fa->farg[3]); /* put index of refraction */
695     printf("\tsides 1\n");
696     puts("\tc"); /* put reflected component */
697     printf("\trs %.4f 0\n", F);
698 greg 2.1 rgb_cie(cxyz, trgb); /* put transmitted component */
699 greg 2.2 puts("\tc");
700 greg 2.1 d = cxyz[0] + cxyz[1] + cxyz[2];
701     if (d > FTINY)
702 greg 2.2 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
703     printf("\tts %.4f 0\n", cxyz[1]);
704 greg 2.1 return(0);
705     }
706    
707    
708     int
709     o_mirror(mod, typ, id, fa) /* convert a mirror material */
710     char *mod, *typ, *id;
711     register FUNARGS *fa;
712     {
713     COLOR cxyz, rrgb;
714     double d;
715    
716     if (fa->nsargs == 1) { /* use alternate material */
717     newmat(id, fa->sarg[0]);
718     return(0);
719     }
720     if (fa->nfargs != 3)
721     return(-1);
722     newmat(id, NULL);
723     rrgb[0] = fa->farg[0]; rrgb[1] = fa->farg[1]; rrgb[2] = fa->farg[2];
724     rgb_cie(cxyz, rrgb);
725 greg 2.2 puts("\tc"); /* put specular component */
726 greg 2.1 d = cxyz[0] + cxyz[1] + cxyz[2];
727     if (d > FTINY)
728 greg 2.2 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
729     printf("\trs %.4f 0\n", cxyz[1]);
730 greg 2.1 return(0);
731     }
732    
733    
734     int
735     o_trans(mod, typ, id, fa) /* convert a trans material */
736     char *mod, *typ, *id;
737     register FUNARGS *fa;
738     {
739     COLOR cxyz, rrgb;
740     double rough, trans, tspec, d;
741    
742     if (typ[4] == '2') { /* trans2 */
743     if (fa->nfargs != 8)
744     return(-1);
745     rough = .5*(fa->farg[4] + fa->farg[5]);
746     trans = fa->farg[6];
747     tspec = fa->farg[7];
748     } else { /* trans */
749     if (fa->nfargs != 7)
750     return(-1);
751     rough = fa->farg[4];
752     trans = fa->farg[5];
753     tspec = fa->farg[6];
754     }
755     newmat(id, NULL);
756     rrgb[0] = fa->farg[0]; rrgb[1] = fa->farg[1]; rrgb[2] = fa->farg[2];
757     rgb_cie(cxyz, rrgb);
758 greg 2.2 puts("\tc"); /* put transmitted diffuse */
759 greg 2.1 d = cxyz[0] + cxyz[1] + cxyz[2];
760     if (d > FTINY)
761 greg 2.2 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
762     printf("\ttd %.4f\n", cxyz[1]*trans*(1. - fa->farg[3])*(1. - tspec));
763 greg 2.1 /* put transmitted specular */
764 greg 2.2 printf("\tts %.4f %.4f\n", cxyz[1]*trans*tspec*(1. - fa->farg[3]), rough);
765 greg 2.1 /* put reflected diffuse */
766 greg 2.2 printf("\trd %.4f\n", cxyz[1]*(1. - fa->farg[3])*(1. - trans));
767     puts("\tc"); /* put reflected specular */
768     printf("\trs %.4f %.4f\n", fa->farg[3], rough);
769 greg 2.1 return(0);
770     }
771    
772    
773     int
774     o_light(mod, typ, id, fa) /* convert a light type */
775     char *mod, *typ, *id;
776     register FUNARGS *fa;
777     {
778     COLOR cxyz, rrgb;
779     double d;
780    
781     if (fa->nfargs < 3)
782     return(-1);
783     newmat(id, NULL);
784     rrgb[0] = fa->farg[0]; rrgb[1] = fa->farg[1]; rrgb[2] = fa->farg[2];
785     rgb_cie(cxyz, rrgb);
786     d = cxyz[0] + cxyz[1] + cxyz[2];
787 greg 2.2 puts("\tc");
788 greg 2.1 if (d > FTINY)
789 greg 2.2 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
790 greg 2.7 printf("\ted %.4g\n", cxyz[1]*(PI*WHTEFFICACY));
791 greg 2.1 return(0);
792     }