ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/gen/xform.c
Revision: 1.17
Committed: Wed Oct 23 15:37:47 1991 UTC (32 years, 6 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 1.16: +12 -8 lines
Log Message:
changed double var[3] to FVECT var

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 (issurface(fn))
252 printf("\n%s %s ", newmod != NULL ? newmod : nam,
253 invert ? ofun[tinvers[fn]].funame : typ);
254 else
255 printf("\n%s %s ", nam, typ);
256 /* object name */
257 fgetword(nam, sizeof(nam), fin);
258 if (idprefix != NULL && issurface(fn))
259 printf("%s.%s\n", idprefix, nam);
260 else
261 printf("%s\n", 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 FUNARGS fa;
308
309 if (readfargs(&fa, fin) != 1)
310 return(-1);
311 /* string arguments */
312 printf("%d", fa.nsargs + xac-xfa);
313 for (i = 0; i < fa.nsargs; i++)
314 printf(" %s", fa.sarg[i]);
315 for (i = xfa; i < xac; i++) /* add xf arguments */
316 printf(" %s", xav[i]);
317 printf("\n");
318 #ifdef IARGS
319 /* integer arguments */
320 printf("%d", fa.niargs);
321 for (i = 0; i < fa.niargs; i++)
322 printf(" %d", fa.iarg[i]);
323 printf("\n");
324 #else
325 printf("0\n");
326 #endif
327 /* float arguments */
328 printf("%d", fa.nfargs);
329 for (i = 0; i < fa.nfargs; i++)
330 printf(" %18.12g", fa.farg[i]);
331 printf("\n");
332 freefargs(&fa);
333 return(0);
334 }
335
336
337 int
338 otype(ofname) /* get object function number from its name */
339 register char *ofname;
340 {
341 register int i;
342
343 for (i = 0; i < NUMTYPES; i++)
344 if (!strcmp(ofun[i].funame, ofname))
345 return(i);
346
347 return(-1); /* not found */
348 }
349
350
351 alias(fin) /* transfer alias */
352 FILE *fin;
353 {
354 char aliasnm[MAXSTR];
355
356 if (fgetword(aliasnm, MAXSTR, fin) == NULL)
357 return(-1);
358 printf("\t%s\n", aliasnm);
359 return(0);
360 }
361
362
363 m_glow(fin) /* transform arguments for proximity light */
364 FILE *fin;
365 {
366 FUNARGS fa;
367
368 if (readfargs(&fa, fin) != 1)
369 return(-1);
370 if (fa.nsargs != 0 || fa.nfargs != 4)
371 return(-1);
372 printf("0\n0\n4");
373 printf(" %18.12g %18.12g %18.12g",
374 fa.farg[0], fa.farg[1], fa.farg[2]);
375 printf(" %18.12g\n", fa.farg[3] * tot.sca);
376 freefargs(&fa);
377 return(0);
378 }
379
380
381 m_spot(fin) /* transform arguments for spotlight */
382 FILE *fin;
383 {
384 FVECT v;
385 FUNARGS fa;
386
387 if (readfargs(&fa, fin) != 1)
388 return(-1);
389 if (fa.nsargs != 0 || fa.nfargs != 7)
390 return(-1);
391 printf("0\n0\n7");
392 printf(" %18.12g %18.12g %18.12g %18.12g\n",
393 fa.farg[0], fa.farg[1], fa.farg[2], fa.farg[3]);
394 multv3(v, fa.farg+4, tot.xfm);
395 printf("\t%18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
396 freefargs(&fa);
397 return(0);
398 }
399
400
401 m_dielectric(fin) /* transform arguments for dielectric */
402 FILE *fin;
403 {
404 double pow();
405 FUNARGS fa;
406
407 if (readfargs(&fa, fin) != 1)
408 return(-1);
409 if (fa.nsargs != 0 || fa.nfargs != 5)
410 return(-1);
411 printf("0\n0\n5");
412 printf(" %18.12g %18.12g %18.12g",
413 pow(fa.farg[0], 1.0/tot.sca),
414 pow(fa.farg[1], 1.0/tot.sca),
415 pow(fa.farg[2], 1.0/tot.sca));
416 printf(" %18.12g %18.12g\n", fa.farg[3], fa.farg[4]);
417 freefargs(&fa);
418 return(0);
419 }
420
421
422 m_interface(fin) /* transform arguments for interface */
423 FILE *fin;
424 {
425 double pow();
426 FUNARGS fa;
427
428 if (readfargs(&fa, fin) != 1)
429 return(-1);
430 if (fa.nsargs != 0 || fa.nfargs != 8)
431 return(-1);
432 printf("0\n0\n8\n");
433 printf("%18.12g %18.12g %18.12g",
434 pow(fa.farg[0], 1.0/tot.sca),
435 pow(fa.farg[1], 1.0/tot.sca),
436 pow(fa.farg[2], 1.0/tot.sca));
437 printf(" %18.12g\n", fa.farg[3]);
438 printf("%18.12g %18.12g %18.12g",
439 pow(fa.farg[4], 1.0/tot.sca),
440 pow(fa.farg[5], 1.0/tot.sca),
441 pow(fa.farg[6], 1.0/tot.sca));
442 printf(" %18.12g\n", fa.farg[7]);
443 freefargs(&fa);
444 return(0);
445 }
446
447
448 text(fin) /* transform text arguments */
449 FILE *fin;
450 {
451 int i;
452 FVECT v;
453 FUNARGS fa;
454
455 if (readfargs(&fa, fin) != 1)
456 return(-1);
457 if (fa.nfargs != 9 && fa.nfargs != 11 && fa.nfargs != 15)
458 return(-1);
459 /* string arguments */
460 printf("%d", fa.nsargs);
461 for (i = 0; i < fa.nsargs; i++)
462 printf(" %s", fa.sarg[i]);
463 printf("\n0\n%d\n", fa.nfargs);
464 /* anchor point */
465 multp3(v, fa.farg, tot.xfm);
466 printf(" %18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
467 /* right vector */
468 multv3(v, fa.farg+3, tot.xfm);
469 printf(" %18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
470 /* down vector */
471 multv3(v, fa.farg+6, tot.xfm);
472 printf(" %18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
473 /* forground and background */
474 if (fa.nfargs == 11)
475 printf(" %18.12g %18.12g\n", fa.farg[9], fa.farg[10]);
476 else if (fa.nfargs == 15) {
477 printf(" %18.12g %18.12g %18.12g\n",
478 fa.farg[9], fa.farg[10], fa.farg[11]);
479 printf(" %18.12g %18.12g %18.12g\n",
480 fa.farg[12], fa.farg[13], fa.farg[14]);
481 }
482 freefargs(&fa);
483 return(0);
484 }
485
486
487 o_source(fin) /* transform source arguments */
488 FILE *fin;
489 {
490 FVECT dv;
491 FUNARGS fa;
492
493 if (readfargs(&fa, fin) != 1)
494 return(-1);
495 if (fa.nsargs != 0 || fa.nfargs != 4)
496 return(-1);
497 /* transform direction vector */
498 multv3(dv, fa.farg, tot.xfm);
499 /* output */
500 printf("0\n0\n4");
501 printf(" %18.12g %18.12g %18.12g %18.12g\n",
502 dv[0], dv[1], dv[2], fa.farg[3]);
503 freefargs(&fa);
504 return(0);
505 }
506
507
508 o_sphere(fin) /* transform sphere arguments */
509 FILE *fin;
510 {
511 FVECT cent;
512 double rad;
513 FUNARGS fa;
514
515 if (readfargs(&fa, fin) != 1)
516 return(-1);
517 if (fa.nsargs != 0 || fa.nfargs != 4)
518 return(-1);
519
520 multp3(cent, fa.farg, tot.xfm); /* transform center */
521
522 rad = fa.farg[3] * tot.sca; /* scale radius */
523
524 printf("0\n0\n4");
525 printf(" %18.12g %18.12g %18.12g %18.12g\n",
526 cent[0], cent[1], cent[2], rad);
527 freefargs(&fa);
528 return(0);
529 }
530
531
532 o_face(fin) /* transform face arguments */
533 FILE *fin;
534 {
535 FVECT p;
536 register int i;
537 FUNARGS fa;
538
539 if (readfargs(&fa, fin) != 1)
540 return(-1);
541 if (fa.nsargs != 0 || fa.nfargs % 3)
542 return(-1);
543
544 printf("0\n0\n%d\n", fa.nfargs);
545
546 for (i = 0; i < fa.nfargs; i += 3) {
547 if (reverse)
548 multp3(p, fa.farg+(fa.nfargs-i-3), tot.xfm);
549 else
550 multp3(p, fa.farg+i, tot.xfm);
551 printf(" %18.12g %18.12g %18.12g\n", p[0], p[1], p[2]);
552 }
553 freefargs(&fa);
554 return(0);
555 }
556
557
558 o_cone(fin) /* transform cone and cup arguments */
559 FILE *fin;
560 {
561 FVECT p0, p1;
562 double r0, r1;
563 FUNARGS fa;
564
565 if (readfargs(&fa, fin) != 1)
566 return(-1);
567 if (fa.nsargs != 0 || fa.nfargs != 8)
568 return(-1);
569
570 printf("0\n0\n8\n");
571
572 multp3(p0, fa.farg, tot.xfm);
573 multp3(p1, fa.farg+3, tot.xfm);
574 r0 = fa.farg[6] * tot.sca;
575 r1 = fa.farg[7] * tot.sca;
576 printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
577 printf(" %18.12g %18.12g %18.12g\n", p1[0], p1[1], p1[2]);
578 printf(" %18.12g %18.12g\n", r0, r1);
579
580 freefargs(&fa);
581 return(0);
582 }
583
584
585 o_cylinder(fin) /* transform cylinder and tube arguments */
586 FILE *fin;
587 {
588 FVECT p0, p1;
589 double rad;
590 FUNARGS fa;
591
592 if (readfargs(&fa, fin) != 1)
593 return(-1);
594 if (fa.nsargs != 0 || fa.nfargs != 7)
595 return(-1);
596
597 printf("0\n0\n7\n");
598
599 multp3(p0, fa.farg, tot.xfm);
600 multp3(p1, fa.farg+3, tot.xfm);
601 rad = fa.farg[6] * tot.sca;
602 printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
603 printf(" %18.12g %18.12g %18.12g\n", p1[0], p1[1], p1[2]);
604 printf(" %18.12g\n", rad);
605 freefargs(&fa);
606 return(0);
607 }
608
609
610 o_ring(fin) /* transform ring arguments */
611 FILE *fin;
612 {
613 FVECT p0, pd;
614 double r0, r1;
615 FUNARGS fa;
616
617 if (readfargs(&fa, fin) != 1)
618 return(-1);
619 if (fa.nsargs != 0 || fa.nfargs != 8)
620 return(-1);
621
622 printf("0\n0\n8\n");
623
624 multp3(p0, fa.farg, tot.xfm);
625 multv3(pd, fa.farg+3, tot.xfm);
626 if (invert) {
627 pd[0] = -pd[0];
628 pd[1] = -pd[1];
629 pd[2] = -pd[2];
630 }
631 r0 = fa.farg[6] * tot.sca;
632 r1 = fa.farg[7] * tot.sca;
633 printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
634 printf(" %18.12g %18.12g %18.12g\n", pd[0], pd[1], pd[2]);
635 printf(" %18.12g %18.12g\n", r0, r1);
636 freefargs(&fa);
637 return(0);
638 }
639
640
641 initotypes() /* initialize ofun[] array */
642 {
643 extern int o_source();
644 extern int o_sphere();
645 extern int o_face();
646 extern int o_cone();
647 extern int o_cylinder();
648 extern int o_ring();
649 extern int m_glow();
650 extern int m_spot();
651 extern int m_dielectric();
652 extern int m_interface();
653 extern int text();
654 extern int alias();
655 extern int passargs();
656 extern int addxform();
657 register int i;
658
659 if (ofun[OBJ_SOURCE].funp == o_source)
660 return; /* done already */
661 /* alias is additional */
662 ofun[ALIAS].funame = ALIASID;
663 ofun[ALIAS].flags = 0;
664 /* functions get new transform */
665 for (i = 0; i < NUMTYPES; i++)
666 if (hasfunc(i))
667 ofun[i].funp = addxform;
668 /* special cases */
669 ofun[OBJ_SOURCE].funp = o_source;
670 ofun[OBJ_SPHERE].funp =
671 ofun[OBJ_BUBBLE].funp = o_sphere;
672 ofun[OBJ_FACE].funp = o_face;
673 ofun[OBJ_CONE].funp =
674 ofun[OBJ_CUP].funp = o_cone;
675 ofun[OBJ_CYLINDER].funp =
676 ofun[OBJ_TUBE].funp = o_cylinder;
677 ofun[OBJ_RING].funp = o_ring;
678 ofun[OBJ_INSTANCE].funp = addxform;
679 ofun[MAT_GLOW].funp = m_glow;
680 ofun[MAT_SPOT].funp = m_spot;
681 ofun[MAT_DIELECTRIC].funp = m_dielectric;
682 ofun[MAT_INTERFACE].funp = m_interface;
683 ofun[PAT_CTEXT].funp =
684 ofun[PAT_BTEXT].funp =
685 ofun[MIX_TEXT].funp = text;
686 ofun[ALIAS].funp = alias;
687 /* surface inverses */
688 tinvers[OBJ_SOURCE] = OBJ_SOURCE;
689 tinvers[OBJ_CONE] = OBJ_CUP;
690 tinvers[OBJ_CUP] = OBJ_CONE;
691 tinvers[OBJ_SPHERE] = OBJ_BUBBLE;
692 tinvers[OBJ_BUBBLE] = OBJ_SPHERE;
693 tinvers[OBJ_RING] = OBJ_RING;
694 tinvers[OBJ_CYLINDER] = OBJ_TUBE;
695 tinvers[OBJ_TUBE] = OBJ_CYLINDER;
696 tinvers[OBJ_INSTANCE] = OBJ_INSTANCE; /* oh, well */
697 }