ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/gen/xform.c
Revision: 2.49
Committed: Tue Feb 2 18:02:32 2016 UTC (8 years, 3 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.48: +1 -2 lines
Log Message:
Moved declaration of popen to paths.h and put convert_command() into module

File Contents

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