ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cv/rad2mgf.c
Revision: 2.19
Committed: Sat Aug 16 15:54:30 2003 UTC (20 years, 8 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.18: +3 -2 lines
Log Message:
Added basic mesh primitive support

File Contents

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