ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cv/rad2mgf.c
Revision: 2.9
Committed: Tue Apr 11 13:33:10 1995 UTC (28 years, 11 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.8: +2 -2 lines
Log Message:
minor changes

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