ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/gen/xform.c
Revision: 2.42
Committed: Sun Jul 26 18:36:17 2009 UTC (14 years, 9 months ago) by greg
Content type: text/plain
Branch: MAIN
CVS Tags: rad4R0
Changes since 2.41: +5 -5 lines
Log Message:
Simplification to bug fix

File Contents

# Content
1 #ifndef lint
2 static const char RCSid[] = "$Id: xform.c,v 2.41 2009/07/26 18:34:59 greg Exp $";
3 #endif
4 /*
5 * xform.c - program to transform object files.
6 * Transformations must preserve aspect ratio.
7 *
8 * 10/19/85
9 * 11/6/86 Finally added error checking!
10 */
11
12 #include <ctype.h>
13
14 #include "platform.h"
15 #include "paths.h"
16 #include "rtprocess.h" /* win_popen() */
17 #include "rtio.h"
18 #include "rtmath.h"
19 #include "object.h"
20 #include "otypes.h"
21
22 int xac; /* global xform argument count */
23 char **xav; /* global xform argument pointer */
24 int xfa; /* start of xf arguments */
25
26 XF tot; /* total transformation */
27 int reverse; /* boolean true if scene mirrored */
28
29 int invert; /* boolean true to invert surfaces */
30
31 int expand; /* boolean true to expand commands */
32
33 char *newmod; /* new modifier for surfaces */
34
35 char *idprefix; /* prefix for object identifiers */
36
37 FUN ofun[NUMOTYPE] = INIT_OTYPE; /* default types and actions */
38
39 short tinvers[NUMOTYPE]; /* inverse types for surfaces */
40
41 int nrept = 1; /* number of array repetitions */
42
43 int stdinused = 0; /* stdin has been used by -f option? */
44
45 char mainfn[PATH_MAX]; /* main file name */
46 FILE *mainfp = NULL; /* main file pointer */
47
48 #define progname (xav[0])
49
50 static int doargf(int ac, char **av, int fi);
51 static int doarray(int ac, char **av, int ai);
52 static void xform(char *name, FILE *fin);
53 static void xfcomm(char *fname, FILE *fin);
54 static void xfobject(char *fname, FILE *fin);
55 static int addxform(FILE *fin);
56 static int alias(FILE *fin);
57 void initotypes(void); /* XXX conflict with otypes.h */
58 static void openmain(char *iname);
59
60
61 int
62 main( /* get transform options and transform file */
63 int argc,
64 char *argv[]
65 )
66 {
67 int mal_prefix = 0;
68 int a;
69 /* check for argument list file */
70 for (a = 1; a < argc; a++)
71 if (!strcmp(argv[a], "-f"))
72 return(doargf(argc, argv, a));
73 /* check for regular array */
74 for (a = 1; a < argc; a++)
75 if (!strcmp(argv[a], "-a"))
76 return(doarray(argc, argv, a));
77
78 initotypes(); /* initialize */
79 invert = 0;
80 expand = 1;
81 newmod = NULL;
82 idprefix = NULL;
83
84 for (a = 1; a < argc; a++) {
85 if (argv[a][0] == '-')
86 switch (argv[a][1]) {
87 case 'm':
88 if (argv[a][2] | (a+1 >= argc))
89 break;
90 a++;
91 if (newmod == NULL)
92 newmod = argv[a];
93 continue;
94 case 'n':
95 if (argv[a][2] | (a+1 >= argc))
96 break;
97 a++;
98 if (idprefix == NULL)
99 idprefix = argv[a];
100 else {
101 register char *newp;
102 newp = (char *)malloc(strlen(idprefix)+
103 strlen(argv[a])+2);
104 if (newp == NULL)
105 exit(2);
106 sprintf(newp, "%s.%s",
107 idprefix, argv[a]);
108 if (mal_prefix++)
109 free((void *)idprefix);
110 idprefix = newp;
111 }
112 continue;
113 case 'c':
114 if (argv[a][2])
115 break;
116 expand = 0;
117 continue;
118 case 'e':
119 if (argv[a][2])
120 break;
121 expand = 1;
122 continue;
123 case 'I':
124 if (argv[a][2])
125 break;
126 invert = !invert;
127 continue;
128 }
129 break;
130 }
131
132 xav = argv;
133 xfa = a;
134
135 a += xf(&tot, argc-a, argv+a);
136
137 if ( (reverse = tot.sca < 0.0) )
138 tot.sca = -tot.sca;
139 if (invert)
140 reverse = !reverse;
141
142 if (a < argc && argv[a][0] == '-') {
143 fprintf(stderr, "%s: command line error at '%s'\n",
144 argv[0], argv[a]);
145 exit(1);
146 }
147
148 xac = a;
149 /* simple header */
150 putchar('#');
151 for (a = 0; a < xac; a++)
152 printf(" %s", xav[a]);
153 putchar('\n');
154 /* transform input */
155 if (xac == argc) {
156 if (stdinused) {
157 fprintf(stderr, "%s: cannot use stdin more than once\n",
158 argv[0]);
159 exit(1);
160 }
161 openmain(NULL);
162 xform(mainfn, mainfp);
163 } else
164 for (a = xac; a < argc; a++) {
165 openmain(argv[a]);
166 xform(mainfn, mainfp);
167 }
168
169 if (mal_prefix)
170 free((void *)idprefix);
171 return(0);
172 }
173
174
175 int
176 doargf( /* take argument list from file */
177 int ac,
178 char **av,
179 int fi
180 )
181 {
182 int inquote;
183 char *newav[256], **avp;
184 char argbuf[2048];
185 char *newid, newidbuf[128];
186 char *oldid;
187 register char *cp;
188 FILE *argfp;
189 int n, i, k, newac, err;
190
191 if (fi >= ac-1 || (av[fi+1][0] == '-' && av[fi+1][1] != '\0')) {
192 fprintf(stderr, "%s: missing file for -f option\n", av[0]);
193 exit(1);
194 }
195 if (av[fi+1][0] == '-' && av[fi+1][1] == '\0') {
196 if (stdinused++) {
197 fprintf(stderr,
198 "%s: cannot use stdin more than once\n",
199 av[0]);
200 exit(1);
201 }
202 argfp = stdin;
203 n = 100; /* we just don't know! */
204 } else {
205 if ((argfp = fopen(av[fi+1], "r")) == NULL) {
206 fprintf(stderr,
207 "%s: cannot open argument file \"%s\"\n",
208 av[0], av[fi+1]);
209 exit(1);
210 }
211 n = 0; /* count number of lines in file */
212 while (fgetline(argbuf,sizeof(argbuf),argfp) != NULL)
213 n += argbuf[0] && argbuf[0] != '#';
214 if (!n) {
215 fprintf(stderr, "%s: empty argument file \"%s\"\n",
216 av[0], av[fi+1]);
217 exit(1);
218 }
219 nrept *= n;
220 rewind(argfp);
221 }
222 err = 0; k = 0; /* read each arg list and call main */
223 while (fgetline(argbuf,sizeof(argbuf),argfp) != NULL) {
224 if (!argbuf[0] || argbuf[0] == '#')
225 continue;
226 avp = newav+2;
227 avp[0] = av[0];
228 for (i = 1; i < fi; i++)
229 avp[i] = av[i];
230 newac = i;
231 cp = argbuf; /* parse new words */
232 if (*cp == '!') cp++;
233 if (!strncmp(cp, "xform ", 6)) cp += 6;
234 inquote = 0;
235 for ( ; ; ) {
236 skipspaces:
237 while (isspace(*cp)) /* nullify spaces */
238 *cp++ = '\0';
239 if ((*cp == '"') | (*cp == '\''))
240 inquote = *cp++;
241 if (!*cp) /* all done? */
242 break;
243 if (cp[0] == '\\' && cp[1])
244 if (*++cp == '\n')
245 goto skipspaces;
246 avp[newac++] = cp; /* add argument to list */
247 if (inquote) {
248 while (*++cp)
249 if (*cp == inquote) {
250 *cp++ = '\0';
251 break;
252 }
253 } else {
254 while (*++cp && !isspace(*cp))
255 ;
256 }
257 }
258 for (i = fi+2; i < ac; i++)
259 avp[newac++] = av[i];
260 avp[newac] = NULL;
261 newid = newidbuf;
262 oldid = NULL;
263 for (i = 2; i < newac; i++)
264 if (!strcmp(avp[i-1], "-n")) {
265 oldid = avp[i];
266 if (strlen(oldid)+32 > sizeof(newidbuf)) {
267 newid = (char *)malloc(strlen(oldid)+32);
268 if (newid == NULL)
269 exit(2);
270 }
271 avp[i] = newid;
272 break;
273 }
274 if (oldid == NULL) {
275 newav[0] = av[0];
276 newav[1] = "-n";
277 newav[2] = newid;
278 avp = newav;
279 newac += 2;
280 }
281 if (oldid == NULL)
282 sprintf(newid, "i%d", k);
283 else
284 sprintf(newid, "%s.%d", oldid, k);
285 err |= main(newac, avp);
286 if (newid != newidbuf)
287 free((void *)newid);
288 k++;
289 }
290 fclose(argfp);
291 return(err);
292 }
293
294
295 int
296 doarray( /* make array */
297 int ac,
298 char **av,
299 int ai
300 )
301 {
302 char *newav[256], **avp;
303 char newid[128], repts[32];
304 char *oldid = NULL;
305 int n, i, err;
306
307 if (ai >= ac-1 || (n = atoi(av[ai+1])) <= 0) {
308 fprintf(stderr, "%s: missing count for -a option\n", av[0]);
309 exit(1);
310 }
311 nrept *= n;
312 avp = newav+2;
313 avp[0] = av[0];
314 for (i = 1; i < ac; i++)
315 if (!strcmp(av[i-1], "-n")) {
316 oldid = av[i];
317 avp[i] = newid;
318 } else
319 avp[i] = av[i];
320 avp[ai] = "-i";
321 avp[ai+1] = repts;
322 avp[i] = NULL;
323 if (oldid == NULL) {
324 newav[0] = av[0];
325 newav[1] = "-n";
326 newav[2] = newid;
327 avp = newav;
328 ac += 2;
329 }
330 err = 0;
331 for (i = 0; i < n; i++) {
332 if (oldid == NULL)
333 sprintf(newid, "a%d", i);
334 else
335 sprintf(newid, "%s.%d", oldid, i);
336 sprintf(repts, "%d", i);
337 err |= main(ac, avp);
338 }
339 return(err);
340 }
341
342
343 void
344 xform( /* transform stream by tot.xfm */
345 char *name,
346 register FILE *fin
347 )
348 {
349 int nobjs = 0;
350 register int c;
351
352 while ((c = getc(fin)) != EOF) {
353 if (isspace(c)) /* blank */
354 continue;
355 if (c == '#') { /* comment */
356 putchar(c);
357 do {
358 if ((c = getc(fin)) == EOF)
359 return;
360 putchar(c);
361 } while (c != '\n');
362 } else if (c == '!') { /* command */
363 ungetc(c, fin);
364 xfcomm(name, fin);
365 nobjs++;
366 } else { /* object */
367 ungetc(c, fin);
368 xfobject(name, fin);
369 nobjs++;
370 }
371 }
372 if (nobjs == 0)
373 fprintf(stderr, "%s: (%s): warning - empty file\n",
374 progname, name);
375 }
376
377
378 void
379 xfcomm( /* transform a command */
380 char *fname,
381 FILE *fin
382 )
383 {
384 FILE *pin;
385 char buf[2048];
386 int i;
387
388 fgetline(buf, sizeof(buf), fin);
389 if (expand) {
390 if ((pin = popen(buf+1, "r")) == NULL) {
391 fprintf(stderr, "%s: (%s): cannot execute \"%s\"\n",
392 progname, fname, buf);
393 exit(1);
394 }
395 xform(buf, pin);
396 pclose(pin);
397 } else {
398 printf("\n%s", buf);
399 if (xac > 1) {
400 printf(" | %s -e", xav[0]); /* expand next time */
401 for (i = 1; i < xac; i++)
402 if (i >= xfa || strcmp(xav[i], "-c"))
403 printf(" %s", xav[i]);
404 }
405 putchar('\n');
406 }
407 }
408
409
410 void
411 xfobject( /* transform an object */
412 char *fname,
413 FILE *fin
414 )
415 {
416 char typ[16], nam[MAXSTR];
417 int fn;
418 /* modifier and type */
419 strcpy(typ, "EOF");
420 fgetword(nam, sizeof(nam), fin);
421 fgetword(typ, sizeof(typ), fin);
422 if ((fn = otype(typ)) < 0) {
423 fprintf(stderr, "%s: (%s): unknown object type \"%s\"\n",
424 progname, fname, typ);
425 exit(1);
426 }
427 putchar('\n');
428 if (ismodifier(fn)) {
429 fputword(nam, stdout);
430 printf(" %s ", typ);
431 } else {
432 fputword(newmod != NULL ? newmod : nam, stdout);
433 printf(" %s ", invert ? ofun[tinvers[fn]].funame : typ);
434 }
435 /* object name */
436 fgetword(nam, sizeof(nam), fin);
437 if (idprefix == NULL || ismodifier(fn))
438 fputword(nam, stdout);
439 else {
440 char nnam[MAXSTR];
441 sprintf(nnam, "%s.%s", idprefix, nam);
442 fputword(nnam, stdout);
443 }
444 putchar('\n');
445 /* transform arguments */
446 if ((*ofun[fn].funp)(fin) < 0) {
447 fprintf(stderr, "%s: (%s): bad %s \"%s\"\n",
448 progname, fname, ofun[fn].funame, nam);
449 exit(1);
450 }
451 }
452
453
454 int
455 o_default( /* pass on arguments unchanged */
456 FILE *fin
457 )
458 {
459 register int i;
460 FUNARGS fa;
461
462 if (readfargs(&fa, fin) != 1)
463 return(-1);
464 /* string arguments */
465 printf("%d", fa.nsargs);
466 for (i = 0; i < fa.nsargs; i++) {
467 fputc(' ', stdout);
468 fputword(fa.sarg[i], stdout);
469 }
470 printf("\n");
471 #ifdef IARGS
472 /* integer arguments */
473 printf("%d", fa.niargs);
474 for (i = 0; i < fa.niargs; i++)
475 printf(" %d", fa.iarg[i]);
476 printf("\n");
477 #else
478 printf("0\n");
479 #endif
480 /* float arguments */
481 printf("%d", fa.nfargs);
482 for (i = 0; i < fa.nfargs; i++)
483 printf(" %18.12g", fa.farg[i]);
484 printf("\n");
485 freefargs(&fa);
486 return(0);
487 }
488
489
490 int
491 addxform( /* add xf arguments to strings */
492 FILE *fin
493 )
494 {
495 register int i;
496 int resetarr = 0;
497 FUNARGS fa;
498
499 if (readfargs(&fa, fin) != 1)
500 return(-1);
501 /* string arguments */
502 if (xac > xfa && strcmp(xav[xfa], "-i"))
503 resetarr = 2;
504 printf("%d", fa.nsargs + resetarr + xac-xfa);
505 for (i = 0; i < fa.nsargs; i++) {
506 fputc(' ', stdout);
507 fputword(fa.sarg[i], stdout);
508 }
509 if (resetarr)
510 printf(" -i 1");
511 for (i = xfa; i < xac; i++) /* add xf arguments */
512 printf(" %s", xav[i]);
513 printf("\n");
514 #ifdef IARGS
515 /* integer arguments */
516 printf("%d", fa.niargs);
517 for (i = 0; i < fa.niargs; i++)
518 printf(" %d", fa.iarg[i]);
519 printf("\n");
520 #else
521 printf("0\n");
522 #endif
523 /* float arguments */
524 printf("%d", fa.nfargs);
525 for (i = 0; i < fa.nfargs; i++)
526 printf(" %18.12g", fa.farg[i]);
527 printf("\n");
528 freefargs(&fa);
529 return(0);
530 }
531
532
533 int
534 alias( /* transfer alias */
535 FILE *fin
536 )
537 {
538 char aliasnm[MAXSTR];
539
540 if (fgetword(aliasnm, MAXSTR, fin) == NULL)
541 return(-1);
542 printf("\t%s\n", aliasnm);
543 return(0);
544 }
545
546
547 int
548 m_glow( /* transform arguments for proximity light */
549 FILE *fin
550 )
551 {
552 FUNARGS fa;
553
554 if (readfargs(&fa, fin) != 1)
555 return(-1);
556 if (fa.nsargs != 0 || fa.nfargs != 4)
557 return(-1);
558 printf("0\n0\n4");
559 printf(" %18.12g %18.12g %18.12g",
560 fa.farg[0], fa.farg[1], fa.farg[2]);
561 printf(" %18.12g\n", fa.farg[3] * tot.sca);
562 freefargs(&fa);
563 return(0);
564 }
565
566
567 int
568 m_spot( /* transform arguments for spotlight */
569 FILE *fin
570 )
571 {
572 FVECT v;
573 FUNARGS fa;
574
575 if (readfargs(&fa, fin) != 1)
576 return(-1);
577 if (fa.nsargs != 0 || fa.nfargs != 7)
578 return(-1);
579 printf("0\n0\n7");
580 printf(" %18.12g %18.12g %18.12g %18.12g\n",
581 fa.farg[0], fa.farg[1], fa.farg[2], fa.farg[3]);
582 multv3(v, fa.farg+4, tot.xfm);
583 printf("\t%18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
584 freefargs(&fa);
585 return(0);
586 }
587
588
589 int
590 m_mist( /* transform arguments for mist */
591 FILE *fin
592 )
593 {
594 FUNARGS fa;
595 int i;
596
597 if (readfargs(&fa, fin) != 1)
598 return(-1);
599 if (fa.nfargs > 7)
600 return(-1);
601 printf("%d", fa.nsargs);
602 if (idprefix == NULL)
603 for (i = 0; i < fa.nsargs; i++)
604 printf(" %s", fa.sarg[i]);
605 else
606 for (i = 0; i < fa.nsargs; i++) {
607 char sname[256], *sp;
608 register char *cp1, *cp2 = sname;
609 /* add idprefix */
610 for (sp = fa.sarg[i]; *sp; sp = cp1) {
611 for (cp1 = idprefix; *cp1; )
612 *cp2++ = *cp1++;
613 *cp2++ = '.';
614 for (cp1 = sp; *cp1 &&
615 (*cp2++ = *cp1++) != '>'; )
616 ;
617 }
618 *cp2 = '\0';
619 printf(" %s", sname);
620 }
621 printf("\n0\n%d", fa.nfargs);
622 if (fa.nfargs > 2)
623 printf(" %12.6g %12.6g %12.6g", fa.farg[0]/tot.sca,
624 fa.farg[1]/tot.sca, fa.farg[2]/tot.sca);
625 for (i = 3; i < fa.nfargs; i++)
626 printf(" %12.6g", fa.farg[i]);
627 printf("\n");
628 freefargs(&fa);
629 return(0);
630 }
631
632
633 int
634 m_dielectric( /* transform arguments for dielectric */
635 FILE *fin
636 )
637 {
638 FUNARGS fa;
639
640 if (readfargs(&fa, fin) != 1)
641 return(-1);
642 if (fa.nsargs != 0 || fa.nfargs != 5)
643 return(-1);
644 printf("0\n0\n5");
645 printf(" %12.6g %12.6g %12.6g",
646 pow(fa.farg[0], 1.0/tot.sca),
647 pow(fa.farg[1], 1.0/tot.sca),
648 pow(fa.farg[2], 1.0/tot.sca));
649 printf(" %12.6g %12.6g\n", fa.farg[3], fa.farg[4]);
650 freefargs(&fa);
651 return(0);
652 }
653
654
655 int
656 m_interface( /* transform arguments for interface */
657 FILE *fin
658 )
659 {
660 FUNARGS fa;
661
662 if (readfargs(&fa, fin) != 1)
663 return(-1);
664 if (fa.nsargs != 0 || fa.nfargs != 8)
665 return(-1);
666 printf("0\n0\n8\n");
667 printf("%12.6g %12.6g %12.6g",
668 pow(fa.farg[0], 1.0/tot.sca),
669 pow(fa.farg[1], 1.0/tot.sca),
670 pow(fa.farg[2], 1.0/tot.sca));
671 printf(" %12.6g\n", fa.farg[3]);
672 printf("%12.6g %12.6g %12.6g",
673 pow(fa.farg[4], 1.0/tot.sca),
674 pow(fa.farg[5], 1.0/tot.sca),
675 pow(fa.farg[6], 1.0/tot.sca));
676 printf(" %12.6g\n", fa.farg[7]);
677 freefargs(&fa);
678 return(0);
679 }
680
681
682 int
683 text( /* transform text arguments */
684 FILE *fin
685 )
686 {
687 int i;
688 FVECT v;
689 FUNARGS fa;
690
691 if (readfargs(&fa, fin) != 1)
692 return(-1);
693 if (fa.nfargs < 9)
694 return(-1);
695 /* string arguments */
696 printf("%d", fa.nsargs);
697 for (i = 0; i < fa.nsargs; i++) {
698 fputc(' ', stdout);
699 fputword(fa.sarg[i], stdout);
700 }
701 printf("\n0\n%d\n", fa.nfargs);
702 /* anchor point */
703 multp3(v, fa.farg, tot.xfm);
704 printf(" %18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
705 /* right vector */
706 multv3(v, fa.farg+3, tot.xfm);
707 printf(" %18.12g %18.12g %18.12g\n", v[0], v[1], v[2]);
708 /* down vector */
709 multv3(v, fa.farg+6, tot.xfm);
710 printf(" %18.12g %18.12g %18.12g", v[0], v[1], v[2]);
711 /* remaining arguments */
712 for (i = 9; i < fa.nfargs; i++) {
713 if (i%3 == 0)
714 putchar('\n');
715 printf(" %18.12g", fa.farg[i]);
716 }
717 putchar('\n');
718 freefargs(&fa);
719 return(0);
720 }
721
722
723 int
724 o_source( /* transform source arguments */
725 FILE *fin
726 )
727 {
728 FVECT dv;
729 FUNARGS fa;
730
731 if (readfargs(&fa, fin) != 1)
732 return(-1);
733 if (fa.nsargs != 0 || fa.nfargs != 4)
734 return(-1);
735 /* transform direction vector */
736 multv3(dv, fa.farg, tot.xfm);
737 /* output */
738 printf("0\n0\n4");
739 printf(" %18.12g %18.12g %18.12g %18.12g\n",
740 dv[0], dv[1], dv[2], fa.farg[3]);
741 freefargs(&fa);
742 return(0);
743 }
744
745
746 int
747 o_sphere( /* transform sphere arguments */
748 FILE *fin
749 )
750 {
751 FVECT cent;
752 double rad;
753 FUNARGS fa;
754
755 if (readfargs(&fa, fin) != 1)
756 return(-1);
757 if (fa.nsargs != 0 || fa.nfargs != 4)
758 return(-1);
759
760 multp3(cent, fa.farg, tot.xfm); /* transform center */
761
762 rad = fa.farg[3] * tot.sca; /* scale radius */
763
764 printf("0\n0\n4");
765 printf(" %18.12g %18.12g %18.12g %18.12g\n",
766 cent[0], cent[1], cent[2], rad);
767 freefargs(&fa);
768 return(0);
769 }
770
771
772 int
773 o_face( /* transform face arguments */
774 FILE *fin
775 )
776 {
777 FVECT p;
778 register int i;
779 FUNARGS fa;
780
781 if (readfargs(&fa, fin) != 1)
782 return(-1);
783 if (fa.nsargs != 0 || fa.nfargs % 3)
784 return(-1);
785
786 printf("0\n0\n%d\n", fa.nfargs);
787
788 for (i = 0; i < fa.nfargs; i += 3) {
789 if (reverse)
790 multp3(p, fa.farg+(fa.nfargs-i-3), tot.xfm);
791 else
792 multp3(p, fa.farg+i, tot.xfm);
793 printf(" %18.12g %18.12g %18.12g\n", p[0], p[1], p[2]);
794 }
795 freefargs(&fa);
796 return(0);
797 }
798
799
800 int
801 o_cone( /* transform cone and cup arguments */
802 FILE *fin
803 )
804 {
805 FVECT p0, p1;
806 double r0, r1;
807 FUNARGS fa;
808
809 if (readfargs(&fa, fin) != 1)
810 return(-1);
811 if (fa.nsargs != 0 || fa.nfargs != 8)
812 return(-1);
813
814 printf("0\n0\n8\n");
815
816 multp3(p0, fa.farg, tot.xfm);
817 multp3(p1, fa.farg+3, tot.xfm);
818 r0 = fa.farg[6] * tot.sca;
819 r1 = fa.farg[7] * tot.sca;
820 printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
821 printf(" %18.12g %18.12g %18.12g\n", p1[0], p1[1], p1[2]);
822 printf(" %18.12g %18.12g\n", r0, r1);
823
824 freefargs(&fa);
825 return(0);
826 }
827
828
829 int
830 o_cylinder( /* transform cylinder and tube arguments */
831 FILE *fin
832 )
833 {
834 FVECT p0, p1;
835 double rad;
836 FUNARGS fa;
837
838 if (readfargs(&fa, fin) != 1)
839 return(-1);
840 if (fa.nsargs != 0 || fa.nfargs != 7)
841 return(-1);
842
843 printf("0\n0\n7\n");
844
845 multp3(p0, fa.farg, tot.xfm);
846 multp3(p1, fa.farg+3, tot.xfm);
847 rad = fa.farg[6] * tot.sca;
848 printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
849 printf(" %18.12g %18.12g %18.12g\n", p1[0], p1[1], p1[2]);
850 printf(" %18.12g\n", rad);
851 freefargs(&fa);
852 return(0);
853 }
854
855
856 int
857 o_ring( /* transform ring arguments */
858 FILE *fin
859 )
860 {
861 FVECT p0, pd;
862 double r0, r1;
863 FUNARGS fa;
864
865 if (readfargs(&fa, fin) != 1)
866 return(-1);
867 if (fa.nsargs != 0 || fa.nfargs != 8)
868 return(-1);
869
870 printf("0\n0\n8\n");
871
872 multp3(p0, fa.farg, tot.xfm);
873 multv3(pd, fa.farg+3, tot.xfm);
874 if (invert) {
875 pd[0] = -pd[0];
876 pd[1] = -pd[1];
877 pd[2] = -pd[2];
878 }
879 r0 = fa.farg[6] * tot.sca;
880 r1 = fa.farg[7] * tot.sca;
881 printf(" %18.12g %18.12g %18.12g\n", p0[0], p0[1], p0[2]);
882 printf(" %18.12g %18.12g %18.12g\n", pd[0], pd[1], pd[2]);
883 printf(" %18.12g %18.12g\n", r0, r1);
884 freefargs(&fa);
885 return(0);
886 }
887
888
889 void
890 initotypes(void) /* initialize ofun[] array */
891 {
892 register int i;
893
894 if (ofun[OBJ_SOURCE].funp == o_source)
895 return; /* done already */
896 /* functions get new transform */
897 for (i = 0; i < NUMOTYPE; i++)
898 if (hasfunc(i))
899 ofun[i].funp = addxform;
900 /* special cases */
901 ofun[OBJ_SOURCE].funp = o_source;
902 ofun[OBJ_SPHERE].funp =
903 ofun[OBJ_BUBBLE].funp = o_sphere;
904 ofun[OBJ_FACE].funp = o_face;
905 ofun[OBJ_CONE].funp =
906 ofun[OBJ_CUP].funp = o_cone;
907 ofun[OBJ_CYLINDER].funp =
908 ofun[OBJ_TUBE].funp = o_cylinder;
909 ofun[OBJ_RING].funp = o_ring;
910 ofun[OBJ_INSTANCE].funp =
911 ofun[OBJ_MESH].funp = addxform;
912 ofun[MAT_GLOW].funp = m_glow;
913 ofun[MAT_SPOT].funp = m_spot;
914 ofun[MAT_DIELECTRIC].funp = m_dielectric;
915 ofun[MAT_INTERFACE].funp = m_interface;
916 ofun[MAT_MIST].funp = m_mist;
917 ofun[PAT_CTEXT].funp =
918 ofun[PAT_BTEXT].funp =
919 ofun[MIX_TEXT].funp = text;
920 ofun[MOD_ALIAS].funp = alias;
921 /* surface inverses */
922 tinvers[OBJ_FACE] = OBJ_FACE;
923 tinvers[OBJ_SOURCE] = OBJ_SOURCE;
924 tinvers[OBJ_CONE] = OBJ_CUP;
925 tinvers[OBJ_CUP] = OBJ_CONE;
926 tinvers[OBJ_SPHERE] = OBJ_BUBBLE;
927 tinvers[OBJ_BUBBLE] = OBJ_SPHERE;
928 tinvers[OBJ_RING] = OBJ_RING;
929 tinvers[OBJ_CYLINDER] = OBJ_TUBE;
930 tinvers[OBJ_TUBE] = OBJ_CYLINDER;
931 tinvers[OBJ_INSTANCE] = OBJ_INSTANCE; /* oh, well */
932 }
933
934
935 #ifdef OLDXFORM
936 void
937 openmain(
938 char *fname
939 )
940 {
941 if (fname == NULL) {
942 strcpy(mainfn, "standard input");
943 mainfp = stdin;
944 return;
945 }
946 if (mainfp != NULL) {
947 if (!strcmp(fname, mainfn)) {
948 rewind(mainfp);
949 return;
950 }
951 fclose(mainfp);
952 }
953 if ((mainfp = fopen(fname, "r")) == NULL) {
954 fprintf(stderr, "%s: cannot open file \"%s\"\n",
955 progname, fname);
956 exit(1);
957 }
958 strcpy(mainfn, fname);
959 }
960 #else
961 void
962 openmain( /* open input, changing directory for file */
963 char *iname
964 )
965 {
966 static char origdir[PATH_MAX];
967 static char curfn[PATH_MAX];
968 static int diffdir;
969 register char *fpath;
970
971 if (iname == NULL) { /* standard input */
972 if (mainfp == NULL) {
973 register int c;
974 strcpy(mainfn, "standard input");
975 if (nrept <= 1) {
976 mainfp = stdin;
977 return; /* just read once */
978 }
979 /* else copy */
980 if ((mainfp = tmpfile()) == NULL) {
981 fprintf(stderr,
982 "%s: cannot create temporary file\n",
983 progname);
984 exit(1);
985 }
986 while ((c = getc(stdin)) != EOF)
987 putc(c, mainfp);
988 }
989 rewind(mainfp); /* rewind copy */
990 return;
991 }
992 if (mainfp == NULL) { /* first call, initialize */
993 getcwd(origdir, sizeof(origdir));
994 } else if (!strcmp(iname, curfn)) { /* just need to rewind? */
995 rewind(mainfp);
996 return;
997 } else { /* else close old stream */
998 fclose(mainfp);
999 mainfp = NULL;
1000 if (diffdir) { /* return to our directory */
1001 chdir(origdir);
1002 diffdir = 0;
1003 }
1004 }
1005 strcpy(curfn, iname); /* remember input name */
1006 /* get full path for file */
1007 if ((fpath = getpath(iname, getrlibpath(), R_OK)) == NULL) {
1008 fprintf(stderr, "%s: cannot find file \"%s\"\n",
1009 progname, iname);
1010 exit(1);
1011 }
1012 if (fpath[0] == '.' && ISDIRSEP(fpath[1])) /* remove leading ./ */
1013 fpath += 2;
1014 /* record path name */
1015 strcpy(mainfn, fpath);
1016 if (expand) { /* change to local directory */
1017 register char *cp = fpath + strlen(fpath); /* get dir. */
1018 while (cp > fpath) {
1019 cp--;
1020 if (ISDIRSEP(*cp)) {
1021 if (cp == fpath)
1022 cp++; /* root special case */
1023 break;
1024 }
1025 }
1026 *cp = '\0';
1027 if (fpath[0]) { /* change to new directory? */
1028 if (chdir(fpath) < 0) {
1029 fprintf(stderr,
1030 "%s: cannot change directory to \"%s\"\n",
1031 progname, fpath);
1032 exit(1);
1033 }
1034 diffdir++;
1035 }
1036 /* get final path component */
1037 for (fpath = iname+strlen(iname);
1038 fpath > iname && !ISDIRSEP(fpath[-1]); fpath--)
1039 ;
1040 }
1041 /* finally, open the file */
1042 if ((mainfp = fopen(fpath, "r")) == NULL) {
1043 fprintf(stderr, "%s: cannot open file \"%s\"\n",
1044 progname, mainfn);
1045 exit(1);
1046 }
1047 }
1048 #endif