ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cv/rad2mgf.c
Revision: 2.26
Committed: Tue Oct 25 23:23:17 2011 UTC (12 years, 5 months ago) by greg
Content type: text/plain
Branch: MAIN
CVS Tags: rad4R2P2, rad5R0, rad5R1, rad4R2, rad4R1, rad4R2P1
Changes since 2.25: +2 -2 lines
Log Message:
Bug fix in handling of trans2 material

File Contents

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