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

# Content
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 char curmat[80]; /* current material */
30 char curobj[128] = "Untitled"; /* current object name */
31
32 double unit_mult = 1.; /* units multiplier */
33
34 #define hasmult (unit_mult < .999 || unit_mult > 1.001)
35
36 /*
37 * Stuff for tracking and reusing vertices:
38 */
39
40 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 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 inp = "standard input";
112 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 if (!strcmp(id, VOIDID)) /* cannot set */
236 return;
237 printf("m %s\n", id);
238 strcpy(curmat, id);
239 }
240
241
242 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 init() /* initialize dispatch table and output */
266 {
267 lu_init(&vertab, NVERTS);
268 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 if (hasmult)
292 printf("xf -s %.4e\n", unit_mult);
293 printf("o %s\n", curobj);
294 }
295
296
297 uninit() /* mark end of MGF file */
298 {
299 puts("o");
300 if (hasmult)
301 puts("xf");
302 puts("# End of data converted from Radiance scene input");
303 lu_done(&rdispatch);
304 lu_done(&rmats);
305 lu_done(&vertab);
306 }
307
308
309 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 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 getvertid(vname, vp) /* get/set vertex ID for this point */
339 char *vname;
340 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 VCOPY(vert[vndx].p, vp); /* assign it */
367 printf("v v%d =\n\tp %.15g %.15g %.15g\n", /* print it */
368 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 vert[vndx].lused = clock; /* record this use */
373 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 register char *cp;
388 register int i;
389
390 if (fa->nfargs < 9 | fa->nfargs % 3)
391 return(-1);
392 setmat(mod);
393 setobj(id);
394 cp = entbuf;
395 *cp++ = 'f';
396 for (i = 0; i < fa->nfargs; i += 3) {
397 *cp++ = ' ';
398 getvertid(cp, fa->farg + i);
399 while (*cp)
400 cp++;
401 }
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 char v1[6], v2[6];
413
414 if (fa->nfargs != 8)
415 return(-1);
416 setmat(mod);
417 setobj(id);
418 getvertid(v1, fa->farg);
419 getvertid(v2, fa->farg + 3);
420 if (typ[1] == 'u') /* cup -> inverted cone */
421 printf("cone %s %.12g %s %.12g\n",
422 v1, -fa->farg[6], v2, -fa->farg[7]);
423 else
424 printf("cone %s %.12g %s %.12g\n",
425 v1, fa->farg[6], v2, fa->farg[7]);
426 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 char cent[6];
436
437 if (fa->nfargs != 4)
438 return(-1);
439 setmat(mod);
440 setobj(id);
441 printf("sph %s %.12g\n", getvertid(cent, fa->farg),
442 typ[0]=='b' ? -fa->farg[3] : fa->farg[3]);
443 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 char v1[6], v2[6];
453
454 if (fa->nfargs != 7)
455 return(-1);
456 setmat(mod);
457 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 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 setobj(id);
475 printf("v cent =\n\tp %.12g %.12g %.12g\n",
476 fa->farg[0], fa->farg[1], fa->farg[2]);
477 printf("\tn %.12g %.12g %.12g\n",
478 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 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 clrverts(); /* vertex id's no longer reliable */
526 return(0);
527 }
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 puts("\tts 1 0");
551 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 puts("\tc"); /* put diffuse component */
569 d = cxyz[0] + cxyz[1] + cxyz[2];
570 if (d > FTINY)
571 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 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 puts("\tc"); /* put diffuse component */
597 d = cxyz[0] + cxyz[1] + cxyz[2];
598 if (d > FTINY)
599 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
600 printf("\trd %.4f\n", cxyz[1]*(1. - fa->farg[3]));
601 /* put specular component */
602 printf("\trs %.4f %.4f\n", cxyz[1]*fa->farg[3],
603 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 trgb[i] = fa->farg[i] * (1. - F)*(1. - F) /
627 (1. - F*F*fa->farg[i]*fa->farg[i]);
628 rrgb[i] = F * (1. + (1. - 2.*F)*fa->farg[i]) /
629 (1. - F*F*fa->farg[i]*fa->farg[i]);
630 }
631 rgb_cie(cxyz, rrgb); /* put reflected component */
632 puts("\tc");
633 d = cxyz[0] + cxyz[1] + cxyz[2];
634 if (d > FTINY)
635 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
636 printf("\trs %.4f 0\n", cxyz[1]);
637 rgb_cie(cxyz, trgb); /* put transmitted component */
638 puts("\tc");
639 d = cxyz[0] + cxyz[1] + cxyz[2];
640 if (d > FTINY)
641 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
642 printf("\tts %.4f 0\n", cxyz[1]);
643 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 puts("\tc"); /* put specular component */
665 d = cxyz[0] + cxyz[1] + cxyz[2];
666 if (d > FTINY)
667 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
668 printf("\trs %.4f 0\n", cxyz[1]);
669 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 puts("\tc"); /* put transmitted diffuse */
698 d = cxyz[0] + cxyz[1] + cxyz[2];
699 if (d > FTINY)
700 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 /* put transmitted specular */
703 printf("\tts %.4f %.4f\n", cxyz[1]*trans*tspec*(1. - fa->farg[3]), rough);
704 /* put reflected diffuse */
705 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 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 puts("\tc");
727 if (d > FTINY)
728 printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d);
729 printf("\ted %.4g\n", cxyz[1]*WHTEFFICACY);
730 return(0);
731 }