ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/gen/xform.c
Revision: 2.11
Committed: Mon Dec 11 09:52:11 1995 UTC (28 years, 4 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.10: +37 -11 lines
Log Message:
bug fixes associated with new mist type

File Contents

# Content
1 /* Copyright (c) 1995 Regents of the University of California */
2
3 #ifndef lint
4 static char SCCSid[] = "$SunId$ LBL";
5 #endif
6
7 /*
8 * xform.c - program to transform object files.
9 * Transformations must preserve aspect ratio.
10 *
11 * 10/19/85
12 * 11/6/86 Finally added error checking!
13 */
14
15 #include "standard.h"
16
17 #include <ctype.h>
18
19 #include "object.h"
20
21 #include "otypes.h"
22
23 int xac; /* global xform argument count */
24 char **xav; /* global xform argument pointer */
25 int xfa; /* start of xf arguments */
26
27 XF tot; /* total transformation */
28 int reverse; /* boolean true if scene mirrored */
29
30 int invert = 0; /* boolean true to invert surfaces */
31
32 int expand = 0; /* boolean true to expand commands */
33
34 char *newmod = NULL; /* new modifier for surfaces */
35
36 char *idprefix = NULL; /* prefix for object identifiers */
37
38 #define ALIAS NUMOTYPE /* put alias at end of array */
39
40 #define NUMTYPES (NUMOTYPE+1) /* total number of object types */
41
42 FUN ofun[NUMTYPES] = INIT_OTYPE; /* default types and actions */
43
44 short tinvers[NUMOTYPE]; /* inverse types for surfaces */
45
46 extern char *malloc(), *fgetword();
47
48 #define progname (xav[0])
49
50
51 main(argc, argv) /* get transform options and transform file */
52 int argc;
53 char *argv[];
54 {
55 FILE *fopen();
56 FILE *fp;
57 int a;
58 /* check for array */
59 for (a = 1; a < argc; a++)
60 if (!strcmp(argv[a], "-a"))
61 return(doarray(argc, argv, a));
62
63 initotypes();
64
65 for (a = 1; a < argc; a++) {
66 if (argv[a][0] == '-')
67 switch (argv[a][1]) {
68 case 'm':
69 if (argv[a][2] || a+1 >= argc)
70 break;
71 newmod = argv[++a];
72 continue;
73 case 'n':
74 if (argv[a][2] || a+1 >= argc)
75 break;
76 idprefix = argv[++a];
77 continue;
78 case 'e':
79 if (argv[a][2])
80 break;
81 expand = 1;
82 continue;
83 case 'I':
84 if (argv[a][2])
85 break;
86 invert = 1;
87 continue;
88 }
89 break;
90 }
91
92 xav = argv;
93 xfa = a;
94
95 a += xf(&tot, argc-a, argv+a);
96
97 if (reverse = tot.sca < 0.0)
98 tot.sca = -tot.sca;
99 if (invert)
100 reverse = !reverse;
101
102 if (a < argc && argv[a][0] == '-') {
103 fprintf(stderr, "%s: command line error at '%s'\n",
104 argv[0], argv[a]);
105 exit(1);
106 }
107
108 xac = a;
109 /* simple header */
110 putchar('#');
111 for (a = 0; a < xac; a++)
112 printf(" %s", xav[a]);
113 putchar('\n');
114 /* transform input */
115 if (xac == argc)
116 xform("standard input", stdin);
117 else
118 for (a = xac; a < argc; a++) {
119 if ((fp = fopen(argv[a], "r")) == NULL) {
120 fprintf(stderr, "%s: cannot open \"%s\"\n",
121 progname, argv[a]);
122 exit(1);
123 }
124 xform(argv[a], fp);
125 fclose(fp);
126 }
127
128 return(0);
129 }
130
131
132 doarray(ac, av, ai) /* make array */
133 char **av;
134 int ac, ai;
135 {
136 char *newav[256], **avp;
137 char newid[128], repts[32];
138 char *oldid = NULL;
139 int i, err;
140
141 avp = newav+2;
142 avp[0] = av[0];
143 for (i = 1; i < ac; i++)
144 if (!strcmp(av[i-1], "-n")) {
145 oldid = av[i];
146 avp[i] = newid;
147 } else
148 avp[i] = av[i];
149 avp[ai] = "-i";
150 avp[ai+1] = repts;
151 avp[i] = NULL;
152 if (oldid == NULL) {
153 newav[0] = av[0];
154 newav[1] = "-n";
155 newav[2] = newid;
156 avp = newav;
157 ac += 2;
158 }
159 err = 0;
160 for (i = 0; i < atoi(av[ai+1]); i++) {
161 if (oldid == NULL)
162 sprintf(newid, "a%d", i);
163 else
164 sprintf(newid, "%s.%d", oldid, i);
165 sprintf(repts, "%d", i);
166 err |= main(ac, avp);
167 }
168 return(err);
169 }
170
171
172 xform(name, fin) /* transform stream by tot.xfm */
173 char *name;
174 register FILE *fin;
175 {
176 int nobjs = 0;
177 register int c;
178
179 while ((c = getc(fin)) != EOF) {
180 if (isspace(c)) /* blank */
181 continue;
182 if (c == '#') { /* comment */
183 putchar(c);
184 do {
185 if ((c = getc(fin)) == EOF)
186 return;
187 putchar(c);
188 } while (c != '\n');
189 } else if (c == '!') { /* command */
190 ungetc(c, fin);
191 xfcomm(name, fin);
192 nobjs++;
193 } else { /* object */
194 ungetc(c, fin);
195 xfobject(name, fin);
196 nobjs++;
197 }
198 }
199 if (nobjs == 0)
200 fprintf(stderr, "%s: (%s): warning - empty file\n",
201 progname, name);
202 }
203
204
205 xfcomm(fname, fin) /* transform a command */
206 char *fname;
207 FILE *fin;
208 {
209 FILE *popen();
210 char *fgetline();
211 FILE *pin;
212 char buf[512];
213 int i;
214
215 fgetline(buf, sizeof(buf), fin);
216 if (expand) {
217 if ((pin = popen(buf+1, "r")) == NULL) {
218 fprintf(stderr, "%s: (%s): cannot execute \"%s\"\n",
219 progname, fname, buf);
220 exit(1);
221 }
222 xform(buf, pin);
223 pclose(pin);
224 } else {
225 printf("\n%s", buf);
226 if (xac > 1) {
227 printf(" | %s -e", xav[0]);
228 for (i = 1; i < xac; i++)
229 printf(" %s", xav[i]);
230 }
231 putchar('\n');
232 }
233 }
234
235
236 xfobject(fname, fin) /* transform an object */
237 char *fname;
238 FILE *fin;
239 {
240 extern char *strcpy();
241 char typ[16], nam[MAXSTR];
242 int fn;
243 /* modifier and type */
244 strcpy(typ, "EOF");
245 fgetword(nam, sizeof(nam), fin);
246 fgetword(typ, sizeof(typ), fin);
247 if ((fn = otype(typ)) < 0) {
248 fprintf(stderr, "%s: (%s): unknown object type \"%s\"\n",
249 progname, fname, typ);
250 exit(1);
251 }
252 if (ismodifier(fn))
253 printf("\n%s %s ", nam, typ);
254 else
255 printf("\n%s %s ", newmod != NULL ? newmod : nam,
256 invert ? ofun[tinvers[fn]].funame : typ);
257 /* object name */
258 fgetword(nam, sizeof(nam), fin);
259 if (idprefix == NULL || ismodifier(fn))
260 printf("%s\n", nam);
261 else
262 printf("%s.%s\n", idprefix, nam);
263 /* transform arguments */
264 if ((*ofun[fn].funp)(fin) < 0) {
265 fprintf(stderr, "%s: (%s): bad %s \"%s\"\n",
266 progname, fname, ofun[fn].funame, nam);
267 exit(1);
268 }
269 }
270
271
272 o_default(fin) /* pass on arguments unchanged */
273 FILE *fin;
274 {
275 register int i;
276 FUNARGS fa;
277
278 if (readfargs(&fa, fin) != 1)
279 return(-1);
280 /* string arguments */
281 printf("%d", fa.nsargs);
282 for (i = 0; i < fa.nsargs; i++)
283 printf(" %s", fa.sarg[i]);
284 printf("\n");
285 #ifdef IARGS
286 /* integer arguments */
287 printf("%d", fa.niargs);
288 for (i = 0; i < fa.niargs; i++)
289 printf(" %d", fa.iarg[i]);
290 printf("\n");
291 #else
292 printf("0\n");
293 #endif
294 /* float arguments */
295 printf("%d", fa.nfargs);
296 for (i = 0; i < fa.nfargs; i++)
297 printf(" %18.12g", fa.farg[i]);
298 printf("\n");
299 freefargs(&fa);
300 return(0);
301 }
302
303
304 addxform(fin) /* add xf arguments to strings */
305 FILE *fin;
306 {
307 register int i;
308 int resetarr = 0;
309 FUNARGS fa;
310
311 if (readfargs(&fa, fin) != 1)
312 return(-1);
313 /* string arguments */
314 if (xac > xfa && strcmp(xav[xfa], "-i"))
315 resetarr = 2;
316 printf("%d", fa.nsargs + resetarr + xac-xfa);
317 for (i = 0; i < fa.nsargs; i++)
318 printf(" %s", fa.sarg[i]);
319 if (resetarr)
320 printf(" -i 1");
321 for (i = xfa; i < xac; i++) /* add xf arguments */
322 printf(" %s", xav[i]);
323 printf("\n");
324 #ifdef IARGS
325 /* integer arguments */
326 printf("%d", fa.niargs);
327 for (i = 0; i < fa.niargs; i++)
328 printf(" %d", fa.iarg[i]);
329 printf("\n");
330 #else
331 printf("0\n");
332 #endif
333 /* float arguments */
334 printf("%d", fa.nfargs);
335 for (i = 0; i < fa.nfargs; i++)
336 printf(" %18.12g", fa.farg[i]);
337 printf("\n");
338 freefargs(&fa);
339 return(0);
340 }
341
342
343 int
344 otype(ofname) /* get object function number from its name */
345 register char *ofname;
346 {
347 register int i;
348
349 for (i = 0; i < NUMTYPES; i++)
350 if (!strcmp(ofun[i].funame, ofname))
351 return(i);
352
353 return(-1); /* not found */
354 }
355
356
357 alias(fin) /* transfer alias */
358 FILE *fin;
359 {
360 char aliasnm[MAXSTR];
361
362 if (fgetword(aliasnm, MAXSTR, fin) == NULL)
363 return(-1);
364 printf("\t%s\n", aliasnm);
365 return(0);
366 }
367
368
369 m_glow(fin) /* transform arguments for proximity light */
370 FILE *fin;
371 {
372 FUNARGS fa;
373
374 if (readfargs(&fa, fin) != 1)
375 return(-1);
376 if (fa.nsargs != 0 || fa.nfargs != 4)
377 return(-1);
378 printf("0\n0\n4");
379 printf(" %18.12g %18.12g %18.12g",
380 fa.farg[0], fa.farg[1], fa.farg[2]);
381 printf(" %18.12g\n", fa.farg[3] * tot.sca);
382 freefargs(&fa);
383 return(0);
384 }
385
386
387 m_spot(fin) /* transform arguments for spotlight */
388 FILE *fin;
389 {
390 FVECT v;
391 FUNARGS fa;
392
393 if (readfargs(&fa, fin) != 1)
394 return(-1);
395 if (fa.nsargs != 0 || fa.nfargs != 7)
396 return(-1);
397 printf("0\n0\n7");
398 printf(" %18.12g %18.12g %18.12g %18.12g\n",
399 fa.farg[0], fa.farg[1], fa.farg[2], fa.farg[3]);
400 multv3(v, fa.farg+4, tot.xfm);
401 printf("\t%18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
402 freefargs(&fa);
403 return(0);
404 }
405
406
407 m_mist(fin) /* transform arguments for mist */
408 FILE *fin;
409 {
410 FUNARGS fa;
411 int i;
412
413 if (readfargs(&fa, fin) != 1)
414 return(-1);
415 if (fa.nfargs > 5)
416 return(-1);
417 printf("%d", fa.nsargs);
418 if (idprefix == NULL)
419 for (i = 0; i < fa.nsargs; i++)
420 printf(" %s", fa.sarg[i]);
421 else
422 for (i = 0; i < fa.nsargs; i++) {
423 char sname[256], *sp;
424 register char *cp1, *cp2 = sname;
425 /* add idprefix */
426 for (sp = fa.sarg[i]; *sp; sp = cp1) {
427 for (cp1 = idprefix; *cp1; )
428 *cp2++ = *cp1++;
429 *cp2++ = '.';
430 for (cp1 = sp; *cp1 &&
431 (*cp2++ = *cp1++) != '>'; )
432 ;
433 }
434 *cp2 = '\0';
435 printf(" %s", sname);
436 }
437 printf("\n0\n%d", fa.nfargs);
438 if (fa.nfargs > 2)
439 printf(" %12.6g %12.6g %12.6g", fa.farg[0]/tot.sca,
440 fa.farg[1]/tot.sca, fa.farg[2]/tot.sca);
441 if (fa.nfargs > 3)
442 printf(" %12.6g", fa.farg[3]);
443 if (fa.nfargs > 4)
444 printf(" %12.6g", fa.farg[4]);
445 printf("\n");
446 freefargs(&fa);
447 return(0);
448 }
449
450
451 m_dielectric(fin) /* transform arguments for dielectric */
452 FILE *fin;
453 {
454 FUNARGS fa;
455
456 if (readfargs(&fa, fin) != 1)
457 return(-1);
458 if (fa.nsargs != 0 || fa.nfargs != 5)
459 return(-1);
460 printf("0\n0\n5");
461 printf(" %12.6g %12.6g %12.6g",
462 pow(fa.farg[0], 1.0/tot.sca),
463 pow(fa.farg[1], 1.0/tot.sca),
464 pow(fa.farg[2], 1.0/tot.sca));
465 printf(" %12.6g %12.6g\n", fa.farg[3], fa.farg[4]);
466 freefargs(&fa);
467 return(0);
468 }
469
470
471 m_interface(fin) /* transform arguments for interface */
472 FILE *fin;
473 {
474 FUNARGS fa;
475
476 if (readfargs(&fa, fin) != 1)
477 return(-1);
478 if (fa.nsargs != 0 || fa.nfargs != 8)
479 return(-1);
480 printf("0\n0\n8\n");
481 printf("%12.6g %12.6g %12.6g",
482 pow(fa.farg[0], 1.0/tot.sca),
483 pow(fa.farg[1], 1.0/tot.sca),
484 pow(fa.farg[2], 1.0/tot.sca));
485 printf(" %12.6g\n", fa.farg[3]);
486 printf("%12.6g %12.6g %12.6g",
487 pow(fa.farg[4], 1.0/tot.sca),
488 pow(fa.farg[5], 1.0/tot.sca),
489 pow(fa.farg[6], 1.0/tot.sca));
490 printf(" %12.6g\n", fa.farg[7]);
491 freefargs(&fa);
492 return(0);
493 }
494
495
496 text(fin) /* transform text arguments */
497 FILE *fin;
498 {
499 int i;
500 FVECT v;
501 FUNARGS fa;
502
503 if (readfargs(&fa, fin) != 1)
504 return(-1);
505 if (fa.nfargs < 9)
506 return(-1);
507 /* string arguments */
508 printf("%d", fa.nsargs);
509 for (i = 0; i < fa.nsargs; i++)
510 printf(" %s", fa.sarg[i]);
511 printf("\n0\n%d\n", fa.nfargs);
512 /* anchor point */
513 multp3(v, fa.farg, tot.xfm);
514 printf(" %18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
515 /* right vector */
516 multv3(v, fa.farg+3, tot.xfm);
517 printf(" %18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
518 /* down vector */
519 multv3(v, fa.farg+6, tot.xfm);
520 printf(" %18.12g %18.12g %18.12g", v[0], v[1], v[2]);
521 /* remaining arguments */
522 for (i = 9; i < fa.nfargs; i++) {
523 if (i%3 == 0)
524 putchar('\n');
525 printf(" %18.12g", fa.farg[i]);
526 }
527 putchar('\n');
528 freefargs(&fa);
529 return(0);
530 }
531
532
533 o_source(fin) /* transform source arguments */
534 FILE *fin;
535 {
536 FVECT dv;
537 FUNARGS fa;
538
539 if (readfargs(&fa, fin) != 1)
540 return(-1);
541 if (fa.nsargs != 0 || fa.nfargs != 4)
542 return(-1);
543 /* transform direction vector */
544 multv3(dv, fa.farg, tot.xfm);
545 /* output */
546 printf("0\n0\n4");
547 printf(" %18.12g %18.12g %18.12g %18.12g\n",
548 dv[0], dv[1], dv[2], fa.farg[3]);
549 freefargs(&fa);
550 return(0);
551 }
552
553
554 o_sphere(fin) /* transform sphere arguments */
555 FILE *fin;
556 {
557 FVECT cent;
558 double rad;
559 FUNARGS fa;
560
561 if (readfargs(&fa, fin) != 1)
562 return(-1);
563 if (fa.nsargs != 0 || fa.nfargs != 4)
564 return(-1);
565
566 multp3(cent, fa.farg, tot.xfm); /* transform center */
567
568 rad = fa.farg[3] * tot.sca; /* scale radius */
569
570 printf("0\n0\n4");
571 printf(" %18.12g %18.12g %18.12g %18.12g\n",
572 cent[0], cent[1], cent[2], rad);
573 freefargs(&fa);
574 return(0);
575 }
576
577
578 o_face(fin) /* transform face arguments */
579 FILE *fin;
580 {
581 FVECT p;
582 register int i;
583 FUNARGS fa;
584
585 if (readfargs(&fa, fin) != 1)
586 return(-1);
587 if (fa.nsargs != 0 || fa.nfargs % 3)
588 return(-1);
589
590 printf("0\n0\n%d\n", fa.nfargs);
591
592 for (i = 0; i < fa.nfargs; i += 3) {
593 if (reverse)
594 multp3(p, fa.farg+(fa.nfargs-i-3), tot.xfm);
595 else
596 multp3(p, fa.farg+i, tot.xfm);
597 printf(" %18.12g %18.12g %18.12g\n", p[0], p[1], p[2]);
598 }
599 freefargs(&fa);
600 return(0);
601 }
602
603
604 o_cone(fin) /* transform cone and cup arguments */
605 FILE *fin;
606 {
607 FVECT p0, p1;
608 double r0, r1;
609 FUNARGS fa;
610
611 if (readfargs(&fa, fin) != 1)
612 return(-1);
613 if (fa.nsargs != 0 || fa.nfargs != 8)
614 return(-1);
615
616 printf("0\n0\n8\n");
617
618 multp3(p0, fa.farg, tot.xfm);
619 multp3(p1, fa.farg+3, tot.xfm);
620 r0 = fa.farg[6] * tot.sca;
621 r1 = fa.farg[7] * tot.sca;
622 printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
623 printf(" %18.12g %18.12g %18.12g\n", p1[0], p1[1], p1[2]);
624 printf(" %18.12g %18.12g\n", r0, r1);
625
626 freefargs(&fa);
627 return(0);
628 }
629
630
631 o_cylinder(fin) /* transform cylinder and tube arguments */
632 FILE *fin;
633 {
634 FVECT p0, p1;
635 double rad;
636 FUNARGS fa;
637
638 if (readfargs(&fa, fin) != 1)
639 return(-1);
640 if (fa.nsargs != 0 || fa.nfargs != 7)
641 return(-1);
642
643 printf("0\n0\n7\n");
644
645 multp3(p0, fa.farg, tot.xfm);
646 multp3(p1, fa.farg+3, tot.xfm);
647 rad = fa.farg[6] * tot.sca;
648 printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
649 printf(" %18.12g %18.12g %18.12g\n", p1[0], p1[1], p1[2]);
650 printf(" %18.12g\n", rad);
651 freefargs(&fa);
652 return(0);
653 }
654
655
656 o_ring(fin) /* transform ring arguments */
657 FILE *fin;
658 {
659 FVECT p0, pd;
660 double r0, r1;
661 FUNARGS fa;
662
663 if (readfargs(&fa, fin) != 1)
664 return(-1);
665 if (fa.nsargs != 0 || fa.nfargs != 8)
666 return(-1);
667
668 printf("0\n0\n8\n");
669
670 multp3(p0, fa.farg, tot.xfm);
671 multv3(pd, fa.farg+3, tot.xfm);
672 if (invert) {
673 pd[0] = -pd[0];
674 pd[1] = -pd[1];
675 pd[2] = -pd[2];
676 }
677 r0 = fa.farg[6] * tot.sca;
678 r1 = fa.farg[7] * tot.sca;
679 printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
680 printf(" %18.12g %18.12g %18.12g\n", pd[0], pd[1], pd[2]);
681 printf(" %18.12g %18.12g\n", r0, r1);
682 freefargs(&fa);
683 return(0);
684 }
685
686
687 initotypes() /* initialize ofun[] array */
688 {
689 register int i;
690
691 if (ofun[OBJ_SOURCE].funp == o_source)
692 return; /* done already */
693 /* alias is additional */
694 ofun[ALIAS].funame = ALIASID;
695 ofun[ALIAS].flags = 0;
696 /* functions get new transform */
697 for (i = 0; i < NUMTYPES; i++)
698 if (hasfunc(i))
699 ofun[i].funp = addxform;
700 /* special cases */
701 ofun[OBJ_SOURCE].funp = o_source;
702 ofun[OBJ_SPHERE].funp =
703 ofun[OBJ_BUBBLE].funp = o_sphere;
704 ofun[OBJ_FACE].funp = o_face;
705 ofun[OBJ_CONE].funp =
706 ofun[OBJ_CUP].funp = o_cone;
707 ofun[OBJ_CYLINDER].funp =
708 ofun[OBJ_TUBE].funp = o_cylinder;
709 ofun[OBJ_RING].funp = o_ring;
710 ofun[OBJ_INSTANCE].funp = addxform;
711 ofun[MAT_GLOW].funp = m_glow;
712 ofun[MAT_SPOT].funp = m_spot;
713 ofun[MAT_DIELECTRIC].funp = m_dielectric;
714 ofun[MAT_INTERFACE].funp = m_interface;
715 ofun[MAT_MIST].funp = m_mist;
716 ofun[PAT_CTEXT].funp =
717 ofun[PAT_BTEXT].funp =
718 ofun[MIX_TEXT].funp = text;
719 ofun[ALIAS].funp = alias;
720 /* surface inverses */
721 tinvers[OBJ_FACE] = OBJ_FACE;
722 tinvers[OBJ_SOURCE] = OBJ_SOURCE;
723 tinvers[OBJ_CONE] = OBJ_CUP;
724 tinvers[OBJ_CUP] = OBJ_CONE;
725 tinvers[OBJ_SPHERE] = OBJ_BUBBLE;
726 tinvers[OBJ_BUBBLE] = OBJ_SPHERE;
727 tinvers[OBJ_RING] = OBJ_RING;
728 tinvers[OBJ_CYLINDER] = OBJ_TUBE;
729 tinvers[OBJ_TUBE] = OBJ_CYLINDER;
730 tinvers[OBJ_INSTANCE] = OBJ_INSTANCE; /* oh, well */
731 }