ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cv/rad2mgf.c
Revision: 2.5
Committed: Mon Aug 15 16:29:22 1994 UTC (29 years, 7 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.4: +12 -0 lines
Log Message:
fixed bug due to vertices changed in include file (instancing)

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