ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cv/rad2mgf.c
Revision: 2.21
Committed: Mon Oct 27 10:26:45 2003 UTC (20 years, 5 months ago) by schorsch
Content type: text/plain
Branch: MAIN
Changes since 2.20: +2 -1 lines
Log Message:
Various compatibility fixes.

File Contents

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