ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/gen/xform.c
Revision: 2.54
Committed: Sat Jun 7 05:09:45 2025 UTC (2 weeks, 2 days ago) by greg
Content type: text/plain
Branch: MAIN
CVS Tags: HEAD
Changes since 2.53: +1 -2 lines
Log Message:
refactor: Put some declarations into "paths.h" and included in "platform.h"

File Contents

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