ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/gen/xform.c
Revision: 2.3
Committed: Thu Mar 12 11:43:58 1992 UTC (32 years, 1 month ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.2: +6 -6 lines
Log Message:
fixed handling of surfaces

File Contents

# Content
1 /* Copyright (c) 1991 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 register int c;
177
178 while ((c = getc(fin)) != EOF) {
179 if (isspace(c)) /* blank */
180 continue;
181 if (c == '#') { /* comment */
182 putchar(c);
183 do {
184 if ((c = getc(fin)) == EOF)
185 return;
186 putchar(c);
187 } while (c != '\n');
188 } else if (c == '!') { /* command */
189 ungetc(c, fin);
190 xfcomm(name, fin);
191 } else { /* object */
192 ungetc(c, fin);
193 xfobject(name, fin);
194 }
195 }
196 }
197
198
199 xfcomm(fname, fin) /* transform a command */
200 FILE *fin;
201 {
202 FILE *popen();
203 char *fgetline();
204 FILE *pin;
205 char buf[512];
206 int i;
207
208 fgetline(buf, sizeof(buf), fin);
209 if (expand) {
210 if (xac > 2) {
211 if ((pin = popen(buf+1, "r")) == NULL) {
212 fprintf(stderr,
213 "%s: (%s): cannot execute \"%s\"\n",
214 progname, fname, buf);
215 exit(1);
216 }
217 xform(buf, pin);
218 pclose(pin);
219 } else {
220 fflush(stdout);
221 system(buf+1);
222 }
223 } else {
224 printf("\n%s", buf);
225 if (xac > 1) {
226 printf(" | %s -e", xav[0]);
227 for (i = 1; i < xac; i++)
228 printf(" %s", xav[i]);
229 }
230 putchar('\n');
231 }
232 }
233
234
235 xfobject(fname, fin) /* transform an object */
236 char *fname;
237 FILE *fin;
238 {
239 extern char *strcpy();
240 char typ[16], nam[MAXSTR];
241 int fn;
242 /* modifier and type */
243 strcpy(typ, "EOF");
244 fgetword(nam, sizeof(nam), fin);
245 fgetword(typ, sizeof(typ), fin);
246 if ((fn = otype(typ)) < 0) {
247 fprintf(stderr, "%s: (%s): unknown object type \"%s\"\n",
248 progname, fname, typ);
249 exit(1);
250 }
251 if (ismodifier(fn))
252 printf("\n%s %s ", nam, typ);
253 else
254 printf("\n%s %s ", newmod != NULL ? newmod : nam,
255 invert ? ofun[tinvers[fn]].funame : typ);
256 /* object name */
257 fgetword(nam, sizeof(nam), fin);
258 if (idprefix == NULL || ismodifier(fn))
259 printf("%s\n", nam);
260 else
261 printf("%s.%s\n", idprefix, nam);
262 /* transform arguments */
263 if ((*ofun[fn].funp)(fin) < 0) {
264 fprintf(stderr, "%s: (%s): bad %s \"%s\"\n",
265 progname, fname, ofun[fn].funame, nam);
266 exit(1);
267 }
268 }
269
270
271 o_default(fin) /* pass on arguments unchanged */
272 FILE *fin;
273 {
274 register int i;
275 FUNARGS fa;
276
277 if (readfargs(&fa, fin) != 1)
278 return(-1);
279 /* string arguments */
280 printf("%d", fa.nsargs);
281 for (i = 0; i < fa.nsargs; i++)
282 printf(" %s", fa.sarg[i]);
283 printf("\n");
284 #ifdef IARGS
285 /* integer arguments */
286 printf("%d", fa.niargs);
287 for (i = 0; i < fa.niargs; i++)
288 printf(" %d", fa.iarg[i]);
289 printf("\n");
290 #else
291 printf("0\n");
292 #endif
293 /* float arguments */
294 printf("%d", fa.nfargs);
295 for (i = 0; i < fa.nfargs; i++)
296 printf(" %18.12g", fa.farg[i]);
297 printf("\n");
298 freefargs(&fa);
299 return(0);
300 }
301
302
303 addxform(fin) /* add xf arguments to strings */
304 FILE *fin;
305 {
306 register int i;
307 int resetarr = 0;
308 FUNARGS fa;
309
310 if (readfargs(&fa, fin) != 1)
311 return(-1);
312 /* string arguments */
313 if (xac > xfa && strcmp(xav[xfa], "-i"))
314 resetarr = 2;
315 printf("%d", fa.nsargs + resetarr + xac-xfa);
316 for (i = 0; i < fa.nsargs; i++)
317 printf(" %s", fa.sarg[i]);
318 if (resetarr)
319 printf(" -i 1");
320 for (i = xfa; i < xac; i++) /* add xf arguments */
321 printf(" %s", xav[i]);
322 printf("\n");
323 #ifdef IARGS
324 /* integer arguments */
325 printf("%d", fa.niargs);
326 for (i = 0; i < fa.niargs; i++)
327 printf(" %d", fa.iarg[i]);
328 printf("\n");
329 #else
330 printf("0\n");
331 #endif
332 /* float arguments */
333 printf("%d", fa.nfargs);
334 for (i = 0; i < fa.nfargs; i++)
335 printf(" %18.12g", fa.farg[i]);
336 printf("\n");
337 freefargs(&fa);
338 return(0);
339 }
340
341
342 int
343 otype(ofname) /* get object function number from its name */
344 register char *ofname;
345 {
346 register int i;
347
348 for (i = 0; i < NUMTYPES; i++)
349 if (!strcmp(ofun[i].funame, ofname))
350 return(i);
351
352 return(-1); /* not found */
353 }
354
355
356 alias(fin) /* transfer alias */
357 FILE *fin;
358 {
359 char aliasnm[MAXSTR];
360
361 if (fgetword(aliasnm, MAXSTR, fin) == NULL)
362 return(-1);
363 printf("\t%s\n", aliasnm);
364 return(0);
365 }
366
367
368 m_glow(fin) /* transform arguments for proximity light */
369 FILE *fin;
370 {
371 FUNARGS fa;
372
373 if (readfargs(&fa, fin) != 1)
374 return(-1);
375 if (fa.nsargs != 0 || fa.nfargs != 4)
376 return(-1);
377 printf("0\n0\n4");
378 printf(" %18.12g %18.12g %18.12g",
379 fa.farg[0], fa.farg[1], fa.farg[2]);
380 printf(" %18.12g\n", fa.farg[3] * tot.sca);
381 freefargs(&fa);
382 return(0);
383 }
384
385
386 m_spot(fin) /* transform arguments for spotlight */
387 FILE *fin;
388 {
389 FVECT v;
390 FUNARGS fa;
391
392 if (readfargs(&fa, fin) != 1)
393 return(-1);
394 if (fa.nsargs != 0 || fa.nfargs != 7)
395 return(-1);
396 printf("0\n0\n7");
397 printf(" %18.12g %18.12g %18.12g %18.12g\n",
398 fa.farg[0], fa.farg[1], fa.farg[2], fa.farg[3]);
399 multv3(v, fa.farg+4, tot.xfm);
400 printf("\t%18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
401 freefargs(&fa);
402 return(0);
403 }
404
405
406 m_dielectric(fin) /* transform arguments for dielectric */
407 FILE *fin;
408 {
409 double pow();
410 FUNARGS fa;
411
412 if (readfargs(&fa, fin) != 1)
413 return(-1);
414 if (fa.nsargs != 0 || fa.nfargs != 5)
415 return(-1);
416 printf("0\n0\n5");
417 printf(" %18.12g %18.12g %18.12g",
418 pow(fa.farg[0], 1.0/tot.sca),
419 pow(fa.farg[1], 1.0/tot.sca),
420 pow(fa.farg[2], 1.0/tot.sca));
421 printf(" %18.12g %18.12g\n", fa.farg[3], fa.farg[4]);
422 freefargs(&fa);
423 return(0);
424 }
425
426
427 m_interface(fin) /* transform arguments for interface */
428 FILE *fin;
429 {
430 double pow();
431 FUNARGS fa;
432
433 if (readfargs(&fa, fin) != 1)
434 return(-1);
435 if (fa.nsargs != 0 || fa.nfargs != 8)
436 return(-1);
437 printf("0\n0\n8\n");
438 printf("%18.12g %18.12g %18.12g",
439 pow(fa.farg[0], 1.0/tot.sca),
440 pow(fa.farg[1], 1.0/tot.sca),
441 pow(fa.farg[2], 1.0/tot.sca));
442 printf(" %18.12g\n", fa.farg[3]);
443 printf("%18.12g %18.12g %18.12g",
444 pow(fa.farg[4], 1.0/tot.sca),
445 pow(fa.farg[5], 1.0/tot.sca),
446 pow(fa.farg[6], 1.0/tot.sca));
447 printf(" %18.12g\n", fa.farg[7]);
448 freefargs(&fa);
449 return(0);
450 }
451
452
453 text(fin) /* transform text arguments */
454 FILE *fin;
455 {
456 int i;
457 FVECT v;
458 FUNARGS fa;
459
460 if (readfargs(&fa, fin) != 1)
461 return(-1);
462 if (fa.nfargs != 9 && fa.nfargs != 11 && fa.nfargs != 15)
463 return(-1);
464 /* string arguments */
465 printf("%d", fa.nsargs);
466 for (i = 0; i < fa.nsargs; i++)
467 printf(" %s", fa.sarg[i]);
468 printf("\n0\n%d\n", fa.nfargs);
469 /* anchor point */
470 multp3(v, fa.farg, tot.xfm);
471 printf(" %18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
472 /* right vector */
473 multv3(v, fa.farg+3, tot.xfm);
474 printf(" %18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
475 /* down vector */
476 multv3(v, fa.farg+6, tot.xfm);
477 printf(" %18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
478 /* forground and background */
479 if (fa.nfargs == 11)
480 printf(" %18.12g %18.12g\n", fa.farg[9], fa.farg[10]);
481 else if (fa.nfargs == 15) {
482 printf(" %18.12g %18.12g %18.12g\n",
483 fa.farg[9], fa.farg[10], fa.farg[11]);
484 printf(" %18.12g %18.12g %18.12g\n",
485 fa.farg[12], fa.farg[13], fa.farg[14]);
486 }
487 freefargs(&fa);
488 return(0);
489 }
490
491
492 o_source(fin) /* transform source arguments */
493 FILE *fin;
494 {
495 FVECT dv;
496 FUNARGS fa;
497
498 if (readfargs(&fa, fin) != 1)
499 return(-1);
500 if (fa.nsargs != 0 || fa.nfargs != 4)
501 return(-1);
502 /* transform direction vector */
503 multv3(dv, fa.farg, tot.xfm);
504 /* output */
505 printf("0\n0\n4");
506 printf(" %18.12g %18.12g %18.12g %18.12g\n",
507 dv[0], dv[1], dv[2], fa.farg[3]);
508 freefargs(&fa);
509 return(0);
510 }
511
512
513 o_sphere(fin) /* transform sphere arguments */
514 FILE *fin;
515 {
516 FVECT cent;
517 double rad;
518 FUNARGS fa;
519
520 if (readfargs(&fa, fin) != 1)
521 return(-1);
522 if (fa.nsargs != 0 || fa.nfargs != 4)
523 return(-1);
524
525 multp3(cent, fa.farg, tot.xfm); /* transform center */
526
527 rad = fa.farg[3] * tot.sca; /* scale radius */
528
529 printf("0\n0\n4");
530 printf(" %18.12g %18.12g %18.12g %18.12g\n",
531 cent[0], cent[1], cent[2], rad);
532 freefargs(&fa);
533 return(0);
534 }
535
536
537 o_face(fin) /* transform face arguments */
538 FILE *fin;
539 {
540 FVECT p;
541 register int i;
542 FUNARGS fa;
543
544 if (readfargs(&fa, fin) != 1)
545 return(-1);
546 if (fa.nsargs != 0 || fa.nfargs % 3)
547 return(-1);
548
549 printf("0\n0\n%d\n", fa.nfargs);
550
551 for (i = 0; i < fa.nfargs; i += 3) {
552 if (reverse)
553 multp3(p, fa.farg+(fa.nfargs-i-3), tot.xfm);
554 else
555 multp3(p, fa.farg+i, tot.xfm);
556 printf(" %18.12g %18.12g %18.12g\n", p[0], p[1], p[2]);
557 }
558 freefargs(&fa);
559 return(0);
560 }
561
562
563 o_cone(fin) /* transform cone and cup arguments */
564 FILE *fin;
565 {
566 FVECT p0, p1;
567 double r0, r1;
568 FUNARGS fa;
569
570 if (readfargs(&fa, fin) != 1)
571 return(-1);
572 if (fa.nsargs != 0 || fa.nfargs != 8)
573 return(-1);
574
575 printf("0\n0\n8\n");
576
577 multp3(p0, fa.farg, tot.xfm);
578 multp3(p1, fa.farg+3, tot.xfm);
579 r0 = fa.farg[6] * tot.sca;
580 r1 = fa.farg[7] * tot.sca;
581 printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
582 printf(" %18.12g %18.12g %18.12g\n", p1[0], p1[1], p1[2]);
583 printf(" %18.12g %18.12g\n", r0, r1);
584
585 freefargs(&fa);
586 return(0);
587 }
588
589
590 o_cylinder(fin) /* transform cylinder and tube arguments */
591 FILE *fin;
592 {
593 FVECT p0, p1;
594 double rad;
595 FUNARGS fa;
596
597 if (readfargs(&fa, fin) != 1)
598 return(-1);
599 if (fa.nsargs != 0 || fa.nfargs != 7)
600 return(-1);
601
602 printf("0\n0\n7\n");
603
604 multp3(p0, fa.farg, tot.xfm);
605 multp3(p1, fa.farg+3, tot.xfm);
606 rad = fa.farg[6] * tot.sca;
607 printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
608 printf(" %18.12g %18.12g %18.12g\n", p1[0], p1[1], p1[2]);
609 printf(" %18.12g\n", rad);
610 freefargs(&fa);
611 return(0);
612 }
613
614
615 o_ring(fin) /* transform ring arguments */
616 FILE *fin;
617 {
618 FVECT p0, pd;
619 double r0, r1;
620 FUNARGS fa;
621
622 if (readfargs(&fa, fin) != 1)
623 return(-1);
624 if (fa.nsargs != 0 || fa.nfargs != 8)
625 return(-1);
626
627 printf("0\n0\n8\n");
628
629 multp3(p0, fa.farg, tot.xfm);
630 multv3(pd, fa.farg+3, tot.xfm);
631 if (invert) {
632 pd[0] = -pd[0];
633 pd[1] = -pd[1];
634 pd[2] = -pd[2];
635 }
636 r0 = fa.farg[6] * tot.sca;
637 r1 = fa.farg[7] * tot.sca;
638 printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
639 printf(" %18.12g %18.12g %18.12g\n", pd[0], pd[1], pd[2]);
640 printf(" %18.12g %18.12g\n", r0, r1);
641 freefargs(&fa);
642 return(0);
643 }
644
645
646 initotypes() /* initialize ofun[] array */
647 {
648 extern int o_source();
649 extern int o_sphere();
650 extern int o_face();
651 extern int o_cone();
652 extern int o_cylinder();
653 extern int o_ring();
654 extern int m_glow();
655 extern int m_spot();
656 extern int m_dielectric();
657 extern int m_interface();
658 extern int text();
659 extern int alias();
660 extern int passargs();
661 extern int addxform();
662 register int i;
663
664 if (ofun[OBJ_SOURCE].funp == o_source)
665 return; /* done already */
666 /* alias is additional */
667 ofun[ALIAS].funame = ALIASID;
668 ofun[ALIAS].flags = 0;
669 /* functions get new transform */
670 for (i = 0; i < NUMTYPES; i++)
671 if (hasfunc(i))
672 ofun[i].funp = addxform;
673 /* special cases */
674 ofun[OBJ_SOURCE].funp = o_source;
675 ofun[OBJ_SPHERE].funp =
676 ofun[OBJ_BUBBLE].funp = o_sphere;
677 ofun[OBJ_FACE].funp = o_face;
678 ofun[OBJ_CONE].funp =
679 ofun[OBJ_CUP].funp = o_cone;
680 ofun[OBJ_CYLINDER].funp =
681 ofun[OBJ_TUBE].funp = o_cylinder;
682 ofun[OBJ_RING].funp = o_ring;
683 ofun[OBJ_INSTANCE].funp = addxform;
684 ofun[MAT_GLOW].funp = m_glow;
685 ofun[MAT_SPOT].funp = m_spot;
686 ofun[MAT_DIELECTRIC].funp = m_dielectric;
687 ofun[MAT_INTERFACE].funp = m_interface;
688 ofun[PAT_CTEXT].funp =
689 ofun[PAT_BTEXT].funp =
690 ofun[MIX_TEXT].funp = text;
691 ofun[ALIAS].funp = alias;
692 /* surface inverses */
693 tinvers[OBJ_SOURCE] = OBJ_SOURCE;
694 tinvers[OBJ_CONE] = OBJ_CUP;
695 tinvers[OBJ_CUP] = OBJ_CONE;
696 tinvers[OBJ_SPHERE] = OBJ_BUBBLE;
697 tinvers[OBJ_BUBBLE] = OBJ_SPHERE;
698 tinvers[OBJ_RING] = OBJ_RING;
699 tinvers[OBJ_CYLINDER] = OBJ_TUBE;
700 tinvers[OBJ_TUBE] = OBJ_CYLINDER;
701 tinvers[OBJ_INSTANCE] = OBJ_INSTANCE; /* oh, well */
702 }