ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/gen/xform.c
Revision: 2.16
Committed: Thu Jun 3 14:15:01 1999 UTC (24 years, 11 months ago) by gwlarson
Content type: text/plain
Branch: MAIN
Changes since 2.15: +88 -3 lines
Log Message:
added -f option to read transform arguments from a file

File Contents

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