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

# Content
1 #ifndef lint
2 static const char RCSid[] = "$Id: rad2mgf.c,v 2.20 2003/10/22 02:06:34 greg Exp $";
3 #endif
4 /*
5 * Convert Radiance scene description to MGF
6 */
7
8 #include "platform.h"
9 #include "standard.h"
10 #include <ctype.h>
11 #include <string.h>
12 #include <stdio.h>
13
14 #include "rtprocess.h"
15 #include "object.h"
16 #include "color.h"
17 #include "lookup.h"
18
19 #define C_1SIDEDTHICK 0.005
20
21 int o_face(), o_cone(), o_sphere(), o_ring(), o_cylinder();
22 int o_instance(), o_illum();
23 int o_plastic(), o_metal(), o_glass(), o_dielectric(),
24 o_mirror(), o_trans(), o_light();
25
26 LUTAB rmats = LU_SINIT(free,NULL); /* defined material table */
27
28 LUTAB rdispatch = LU_SINIT(NULL,NULL); /* function dispatch table */
29
30 char curmat[80]; /* current material */
31 char curobj[128] = "Untitled"; /* current object name */
32
33 double unit_mult = 1.; /* units multiplier */
34
35 #define hasmult (unit_mult < .999 || unit_mult > 1.001)
36
37 /*
38 * Stuff for tracking and reusing vertices:
39 */
40
41 char VKFMT[] = "%+16.9e %+16.9e %+16.9e";
42 #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 long vclock; /* incremented at each vertex request */
49
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 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 default:
85 goto unkopt;
86 }
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 inp = "standard input";
115 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 } 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 }
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 if (!strcmp(id, VOIDID)) /* cannot set */
242 return;
243 printf("m %s\n", id);
244 strcpy(curmat, id);
245 }
246
247
248 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 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 if ((*cp < '!') | (*cp > '~')) /* limit to visible chars */
268 *cp = '?';
269 diff += *cp != *cp2;
270 }
271 if (!diff && !*cp2)
272 return;
273 *cp2 = '\0';
274 fputs("o\no ", stdout);
275 puts(curobj);
276 }
277
278
279 init() /* initialize dispatch table and output */
280 {
281 lu_init(&vertab, NVERTS);
282 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 add2dispatch("mesh", o_instance);
293 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 add2dispatch("dielectric", o_dielectric);
299 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 puts("# The following was converted from RADIANCE scene input");
307 if (hasmult)
308 printf("xf -s %.4e\n", unit_mult);
309 printf("o %s\n", curobj);
310 }
311
312
313 uninit() /* mark end of MGF file */
314 {
315 puts("o");
316 if (hasmult)
317 puts("xf");
318 puts("# End of data converted from RADIANCE scene input");
319 lu_done(&rdispatch);
320 lu_done(&rmats);
321 lu_done(&vertab);
322 }
323
324
325 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 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 getvertid(vname, vp) /* get/set vertex ID for this point */
355 char *vname;
356 FVECT vp;
357 {
358 static char vkey[VKLEN];
359 register LUENT *lp;
360 register int i, vndx;
361
362 vclock++; /* increment counter */
363 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 VCOPY(vert[vndx].p, vp); /* assign it */
383 printf("v v%d =\n\tp %.15g %.15g %.15g\n", /* print it */
384 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 vert[vndx].lused = vclock; /* record this use */
389 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 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 o_face(mod, typ, id, fa) /* print out a polygon */
426 char *mod, *typ, *id;
427 FUNARGS *fa;
428 {
429 char entbuf[2048], *linestart;
430 register char *cp;
431 register int i;
432
433 if ((fa->nfargs < 9) | (fa->nfargs % 3))
434 return(-1);
435 setmat(mod);
436 setobj(id);
437 cp = linestart = entbuf;
438 *cp++ = 'f';
439 for (i = 0; i < fa->nfargs; i += 3) {
440 *cp++ = ' ';
441 if (cp - linestart > 72) {
442 *cp++ = '\\'; *cp++ = '\n';
443 linestart = cp;
444 *cp++ = ' '; *cp++ = ' ';
445 }
446 getvertid(cp, fa->farg + i);
447 while (*cp)
448 cp++;
449 }
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 char v1[6], v2[6];
461
462 if (fa->nfargs != 8)
463 return(-1);
464 setmat(mod);
465 setobj(id);
466 getvertid(v1, fa->farg);
467 getvertid(v2, fa->farg + 3);
468 if (typ[1] == 'u') /* cup -> inverted cone */
469 printf("cone %s %.12g %s %.12g\n",
470 v1, -fa->farg[6], v2, -fa->farg[7]);
471 else
472 printf("cone %s %.12g %s %.12g\n",
473 v1, fa->farg[6], v2, fa->farg[7]);
474 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 char cent[6];
484
485 if (fa->nfargs != 4)
486 return(-1);
487 setmat(mod);
488 setobj(id);
489 printf("sph %s %.12g\n", getvertid(cent, fa->farg),
490 typ[0]=='b' ? -fa->farg[3] : fa->farg[3]);
491 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 char v1[6], v2[6];
501
502 if (fa->nfargs != 7)
503 return(-1);
504 setmat(mod);
505 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 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 setobj(id);
523 printf("v cent =\n\tp %.12g %.12g %.12g\n",
524 fa->farg[0], fa->farg[1], fa->farg[2]);
525 printf("\tn %.12g %.12g %.12g\n",
526 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 o_instance(mod, typ, id, fa) /* convert an instance (or mesh) */
539 char *mod, *typ, *id;
540 FUNARGS *fa;
541 {
542 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 clrverts(); /* vertex id's no longer reliable */
574 return(0);
575 }
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 puts("\tts 1 0");
590 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 puts("\tc"); /* put diffuse component */
608 d = cxyz[0] + cxyz[1] + cxyz[2];
609 if (d > FTINY)
610 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 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 puts("\tc"); /* put diffuse component */
636 d = cxyz[0] + cxyz[1] + cxyz[2];
637 if (d > FTINY)
638 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
639 printf("\trd %.4f\n", cxyz[1]*(1. - fa->farg[3]));
640 /* put specular component */
641 printf("\trs %.4f %.4f\n", cxyz[1]*fa->farg[3],
642 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 printf("\tir %f 0\n", nrfr);
663 F = (1. - nrfr)/(1. + nrfr); /* use normal incidence */
664 F *= F;
665 for (i = 0; i < 3; i++) {
666 trgb[i] = fa->farg[i] * (1. - F)*(1. - F) /
667 (1. - F*F*fa->farg[i]*fa->farg[i]);
668 rrgb[i] = F * (1. + (1. - 2.*F)*fa->farg[i]) /
669 (1. - F*F*fa->farg[i]*fa->farg[i]);
670 }
671 rgb_cie(cxyz, rrgb); /* put reflected component */
672 puts("\tc");
673 d = cxyz[0] + cxyz[1] + cxyz[2];
674 if (d > FTINY)
675 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
676 printf("\trs %.4f 0\n", cxyz[1]);
677 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 rgb_cie(cxyz, trgb); /* put transmitted component */
708 puts("\tc");
709 d = cxyz[0] + cxyz[1] + cxyz[2];
710 if (d > FTINY)
711 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
712 printf("\tts %.4f 0\n", cxyz[1]);
713 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 puts("\tc"); /* put specular component */
735 d = cxyz[0] + cxyz[1] + cxyz[2];
736 if (d > FTINY)
737 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
738 printf("\trs %.4f 0\n", cxyz[1]);
739 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 puts("\tc"); /* put transmitted diffuse */
768 d = cxyz[0] + cxyz[1] + cxyz[2];
769 if (d > FTINY)
770 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 /* put transmitted specular */
773 printf("\tts %.4f %.4f\n", cxyz[1]*trans*tspec*(1. - fa->farg[3]), rough);
774 /* put reflected diffuse */
775 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 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 puts("\tc");
797 if (d > FTINY)
798 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
799 printf("\ted %.4g\n", cxyz[1]*(PI*WHTEFFICACY));
800 return(0);
801 }