ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/gen/xform.c
Revision: 2.10
Committed: Fri Dec 8 18:29:06 1995 UTC (28 years, 4 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.9: +19 -0 lines
Log Message:
added M_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
412 if (readfargs(&fa, fin) != 1)
413 return(-1);
414 if (fa.nsargs != 0 || fa.nfargs != 4)
415 return(-1);
416 printf("0\n0\n4");
417 printf(" %18.12g %18.12g %18.12g", fa.farg[0]/tot.sca,
418 fa.farg[1]/tot.sca, fa.farg[2]/tot.sca);
419 printf(" %18.12g\n", fa.farg[3]);
420 freefargs(&fa);
421 return(0);
422 }
423
424
425 m_dielectric(fin) /* transform arguments for dielectric */
426 FILE *fin;
427 {
428 FUNARGS fa;
429
430 if (readfargs(&fa, fin) != 1)
431 return(-1);
432 if (fa.nsargs != 0 || fa.nfargs != 5)
433 return(-1);
434 printf("0\n0\n5");
435 printf(" %18.12g %18.12g %18.12g",
436 pow(fa.farg[0], 1.0/tot.sca),
437 pow(fa.farg[1], 1.0/tot.sca),
438 pow(fa.farg[2], 1.0/tot.sca));
439 printf(" %18.12g %18.12g\n", fa.farg[3], fa.farg[4]);
440 freefargs(&fa);
441 return(0);
442 }
443
444
445 m_interface(fin) /* transform arguments for interface */
446 FILE *fin;
447 {
448 FUNARGS fa;
449
450 if (readfargs(&fa, fin) != 1)
451 return(-1);
452 if (fa.nsargs != 0 || fa.nfargs != 8)
453 return(-1);
454 printf("0\n0\n8\n");
455 printf("%18.12g %18.12g %18.12g",
456 pow(fa.farg[0], 1.0/tot.sca),
457 pow(fa.farg[1], 1.0/tot.sca),
458 pow(fa.farg[2], 1.0/tot.sca));
459 printf(" %18.12g\n", fa.farg[3]);
460 printf("%18.12g %18.12g %18.12g",
461 pow(fa.farg[4], 1.0/tot.sca),
462 pow(fa.farg[5], 1.0/tot.sca),
463 pow(fa.farg[6], 1.0/tot.sca));
464 printf(" %18.12g\n", fa.farg[7]);
465 freefargs(&fa);
466 return(0);
467 }
468
469
470 text(fin) /* transform text arguments */
471 FILE *fin;
472 {
473 int i;
474 FVECT v;
475 FUNARGS fa;
476
477 if (readfargs(&fa, fin) != 1)
478 return(-1);
479 if (fa.nfargs < 9)
480 return(-1);
481 /* string arguments */
482 printf("%d", fa.nsargs);
483 for (i = 0; i < fa.nsargs; i++)
484 printf(" %s", fa.sarg[i]);
485 printf("\n0\n%d\n", fa.nfargs);
486 /* anchor point */
487 multp3(v, fa.farg, tot.xfm);
488 printf(" %18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
489 /* right vector */
490 multv3(v, fa.farg+3, tot.xfm);
491 printf(" %18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
492 /* down vector */
493 multv3(v, fa.farg+6, tot.xfm);
494 printf(" %18.12g %18.12g %18.12g", v[0], v[1], v[2]);
495 /* remaining arguments */
496 for (i = 9; i < fa.nfargs; i++) {
497 if (i%3 == 0)
498 putchar('\n');
499 printf(" %18.12g", fa.farg[i]);
500 }
501 putchar('\n');
502 freefargs(&fa);
503 return(0);
504 }
505
506
507 o_source(fin) /* transform source arguments */
508 FILE *fin;
509 {
510 FVECT dv;
511 FUNARGS fa;
512
513 if (readfargs(&fa, fin) != 1)
514 return(-1);
515 if (fa.nsargs != 0 || fa.nfargs != 4)
516 return(-1);
517 /* transform direction vector */
518 multv3(dv, fa.farg, tot.xfm);
519 /* output */
520 printf("0\n0\n4");
521 printf(" %18.12g %18.12g %18.12g %18.12g\n",
522 dv[0], dv[1], dv[2], fa.farg[3]);
523 freefargs(&fa);
524 return(0);
525 }
526
527
528 o_sphere(fin) /* transform sphere arguments */
529 FILE *fin;
530 {
531 FVECT cent;
532 double rad;
533 FUNARGS fa;
534
535 if (readfargs(&fa, fin) != 1)
536 return(-1);
537 if (fa.nsargs != 0 || fa.nfargs != 4)
538 return(-1);
539
540 multp3(cent, fa.farg, tot.xfm); /* transform center */
541
542 rad = fa.farg[3] * tot.sca; /* scale radius */
543
544 printf("0\n0\n4");
545 printf(" %18.12g %18.12g %18.12g %18.12g\n",
546 cent[0], cent[1], cent[2], rad);
547 freefargs(&fa);
548 return(0);
549 }
550
551
552 o_face(fin) /* transform face arguments */
553 FILE *fin;
554 {
555 FVECT p;
556 register int i;
557 FUNARGS fa;
558
559 if (readfargs(&fa, fin) != 1)
560 return(-1);
561 if (fa.nsargs != 0 || fa.nfargs % 3)
562 return(-1);
563
564 printf("0\n0\n%d\n", fa.nfargs);
565
566 for (i = 0; i < fa.nfargs; i += 3) {
567 if (reverse)
568 multp3(p, fa.farg+(fa.nfargs-i-3), tot.xfm);
569 else
570 multp3(p, fa.farg+i, tot.xfm);
571 printf(" %18.12g %18.12g %18.12g\n", p[0], p[1], p[2]);
572 }
573 freefargs(&fa);
574 return(0);
575 }
576
577
578 o_cone(fin) /* transform cone and cup arguments */
579 FILE *fin;
580 {
581 FVECT p0, p1;
582 double r0, r1;
583 FUNARGS fa;
584
585 if (readfargs(&fa, fin) != 1)
586 return(-1);
587 if (fa.nsargs != 0 || fa.nfargs != 8)
588 return(-1);
589
590 printf("0\n0\n8\n");
591
592 multp3(p0, fa.farg, tot.xfm);
593 multp3(p1, fa.farg+3, tot.xfm);
594 r0 = fa.farg[6] * tot.sca;
595 r1 = fa.farg[7] * tot.sca;
596 printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
597 printf(" %18.12g %18.12g %18.12g\n", p1[0], p1[1], p1[2]);
598 printf(" %18.12g %18.12g\n", r0, r1);
599
600 freefargs(&fa);
601 return(0);
602 }
603
604
605 o_cylinder(fin) /* transform cylinder and tube arguments */
606 FILE *fin;
607 {
608 FVECT p0, p1;
609 double rad;
610 FUNARGS fa;
611
612 if (readfargs(&fa, fin) != 1)
613 return(-1);
614 if (fa.nsargs != 0 || fa.nfargs != 7)
615 return(-1);
616
617 printf("0\n0\n7\n");
618
619 multp3(p0, fa.farg, tot.xfm);
620 multp3(p1, fa.farg+3, tot.xfm);
621 rad = fa.farg[6] * 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\n", rad);
625 freefargs(&fa);
626 return(0);
627 }
628
629
630 o_ring(fin) /* transform ring arguments */
631 FILE *fin;
632 {
633 FVECT p0, pd;
634 double r0, r1;
635 FUNARGS fa;
636
637 if (readfargs(&fa, fin) != 1)
638 return(-1);
639 if (fa.nsargs != 0 || fa.nfargs != 8)
640 return(-1);
641
642 printf("0\n0\n8\n");
643
644 multp3(p0, fa.farg, tot.xfm);
645 multv3(pd, fa.farg+3, tot.xfm);
646 if (invert) {
647 pd[0] = -pd[0];
648 pd[1] = -pd[1];
649 pd[2] = -pd[2];
650 }
651 r0 = fa.farg[6] * tot.sca;
652 r1 = fa.farg[7] * tot.sca;
653 printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
654 printf(" %18.12g %18.12g %18.12g\n", pd[0], pd[1], pd[2]);
655 printf(" %18.12g %18.12g\n", r0, r1);
656 freefargs(&fa);
657 return(0);
658 }
659
660
661 initotypes() /* initialize ofun[] array */
662 {
663 register int i;
664
665 if (ofun[OBJ_SOURCE].funp == o_source)
666 return; /* done already */
667 /* alias is additional */
668 ofun[ALIAS].funame = ALIASID;
669 ofun[ALIAS].flags = 0;
670 /* functions get new transform */
671 for (i = 0; i < NUMTYPES; i++)
672 if (hasfunc(i))
673 ofun[i].funp = addxform;
674 /* special cases */
675 ofun[OBJ_SOURCE].funp = o_source;
676 ofun[OBJ_SPHERE].funp =
677 ofun[OBJ_BUBBLE].funp = o_sphere;
678 ofun[OBJ_FACE].funp = o_face;
679 ofun[OBJ_CONE].funp =
680 ofun[OBJ_CUP].funp = o_cone;
681 ofun[OBJ_CYLINDER].funp =
682 ofun[OBJ_TUBE].funp = o_cylinder;
683 ofun[OBJ_RING].funp = o_ring;
684 ofun[OBJ_INSTANCE].funp = addxform;
685 ofun[MAT_GLOW].funp = m_glow;
686 ofun[MAT_SPOT].funp = m_spot;
687 ofun[MAT_DIELECTRIC].funp = m_dielectric;
688 ofun[MAT_INTERFACE].funp = m_interface;
689 ofun[MAT_MIST].funp = m_mist;
690 ofun[PAT_CTEXT].funp =
691 ofun[PAT_BTEXT].funp =
692 ofun[MIX_TEXT].funp = text;
693 ofun[ALIAS].funp = alias;
694 /* surface inverses */
695 tinvers[OBJ_FACE] = OBJ_FACE;
696 tinvers[OBJ_SOURCE] = OBJ_SOURCE;
697 tinvers[OBJ_CONE] = OBJ_CUP;
698 tinvers[OBJ_CUP] = OBJ_CONE;
699 tinvers[OBJ_SPHERE] = OBJ_BUBBLE;
700 tinvers[OBJ_BUBBLE] = OBJ_SPHERE;
701 tinvers[OBJ_RING] = OBJ_RING;
702 tinvers[OBJ_CYLINDER] = OBJ_TUBE;
703 tinvers[OBJ_TUBE] = OBJ_CYLINDER;
704 tinvers[OBJ_INSTANCE] = OBJ_INSTANCE; /* oh, well */
705 }