ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/gen/genworm.c
Revision: 2.7
Committed: Sun Nov 16 10:29:38 2003 UTC (20 years, 4 months ago) by schorsch
Content type: text/plain
Branch: MAIN
Changes since 2.6: +17 -8 lines
Log Message:
Continued ANSIfication and reduced other compile warnings.

File Contents

# Content
1 #ifndef lint
2 static const char RCSid[] = "$Id: genworm.c,v 2.6 2003/08/04 22:37:53 greg Exp $";
3 #endif
4 /*
5 * genworm.c - program to generate worms (strings with varying thickness).
6 *
7 * The program takes as input the functions of t for x, y,
8 * z and r (the radius). Cylinders and cones will be
9 * joined by spheres. Negative radii are forbidden.
10 *
11 * 9/24/87
12 */
13
14 #include <stdlib.h>
15 #include <stdio.h>
16 #include <math.h>
17 #include <string.h>
18
19 #include "calcomp.h"
20 #include "fvect.h"
21
22 #define XNAME "X`SYS`" /* x function name */
23 #define YNAME "Y`SYS`" /* y function name */
24 #define ZNAME "Z`SYS`" /* z function name */
25 #define RNAME "R`SYS`" /* r function name */
26
27 #define PI 3.14159265358979323846
28
29 #define max(a,b) ((a) > (b) ? (a) : (b))
30
31
32 /* XXX redundant, move to library */
33 double l_hermite(char *), l_bezier(char *), l_bspline(char *);
34
35
36 static void printhead(int ac, char **av);
37
38
39 int
40 main(argc, argv)
41 int argc;
42 char *argv[];
43 {
44 char stmp[256];
45 double t, f, lastr = 0, r;
46 FVECT lastp, p;
47 int i, nseg;
48
49 varset("PI", ':', PI);
50 funset("hermite", 5, ':', l_hermite);
51 funset("bezier", 5, ':', l_bezier);
52 funset("bspline", 5, ':', l_bspline);
53
54 if (argc < 8)
55 goto userror;
56
57 for (i = 8; i < argc; i++)
58 if (!strcmp(argv[i], "-e"))
59 scompile(argv[++i], NULL, 0);
60 else if (!strcmp(argv[i], "-f"))
61 fcompile(argv[++i]);
62 else
63 goto userror;
64
65 sprintf(stmp, "%s(t)=%s;", XNAME, argv[3]);
66 scompile(stmp, NULL, 0);
67 sprintf(stmp, "%s(t)=%s;", YNAME, argv[4]);
68 scompile(stmp, NULL, 0);
69 sprintf(stmp, "%s(t)=%s;", ZNAME, argv[5]);
70 scompile(stmp, NULL, 0);
71 sprintf(stmp, "%s(t)=%s;", RNAME, argv[6]);
72 scompile(stmp, NULL, 0);
73 nseg = atoi(argv[7]);
74 if (nseg <= 0)
75 goto userror;
76
77 printhead(argc, argv);
78 eclock = 0;
79
80 for (i = 0; i <= nseg; i++) {
81 t = (double)i/nseg;
82 p[0] = funvalue(XNAME, 1, &t);
83 p[1] = funvalue(YNAME, 1, &t);
84 p[2] = funvalue(ZNAME, 1, &t);
85 r = funvalue(RNAME, 1, &t);
86 if (i) {
87 if (lastr <= r+FTINY && lastr >= r-FTINY) {
88 printf("\n%s cylinder %s.c%d\n",
89 argv[1], argv[2], i);
90 printf("0\n0\n7\n");
91 printf("%18.12g %18.12g %18.12g\n",
92 lastp[0], lastp[1], lastp[2]);
93 printf("%18.12g %18.12g %18.12g\n",
94 p[0], p[1], p[2]);
95 printf("%18.12g\n", r);
96 } else {
97 printf("\n%s cone %s.c%d\n",
98 argv[1], argv[2], i);
99 printf("0\n0\n8\n");
100 f = (lastr - r)/dist2(lastp,p);
101 printf("%18.12g %18.12g %18.12g\n",
102 lastp[0] + f*lastr*(p[0] - lastp[0]),
103 lastp[1] + f*lastr*(p[1] - lastp[1]),
104 lastp[2] + f*lastr*(p[2] - lastp[2]));
105 printf("%18.12g %18.12g %18.12g\n",
106 p[0] + f*r*(p[0] - lastp[0]),
107 p[1] + f*r*(p[1] - lastp[1]),
108 p[2] + f*r*(p[2] - lastp[2]));
109 f = 1.0 - (lastr-r)*f;
110 f = f <= 0.0 ? 0.0 : sqrt(f);
111 printf("%18.12g %18.12g\n", f*lastr, f*r);
112 }
113 }
114 printf("\n%s sphere %s.s%d\n", argv[1], argv[2], i);
115 printf("0\n0\n4 %18.12g %18.12g %18.12g %18.12g\n",
116 p[0], p[1], p[2], r);
117 VCOPY(lastp, p);
118 lastr = r;
119 }
120 return 0;
121
122 userror:
123 fprintf(stderr,
124 "Usage: %s material name x(t) y(t) z(t) r(t) nseg [-e expr] [-f file]\n",
125 argv[0]);
126 return 1;
127 }
128
129
130 void
131 eputs(msg)
132 char *msg;
133 {
134 fputs(msg, stderr);
135 }
136
137
138 void
139 wputs(msg)
140 char *msg;
141 {
142 eputs(msg);
143 }
144
145
146 void
147 quit(code)
148 int code;
149 {
150 exit(code);
151 }
152
153
154 void
155 printhead( /* print command header */
156 register int ac,
157 register char **av
158 )
159 {
160 putchar('#');
161 while (ac--) {
162 putchar(' ');
163 fputs(*av++, stdout);
164 }
165 putchar('\n');
166 }
167
168
169 double
170 l_hermite(char *nm)
171 {
172 double t;
173
174 t = argument(5);
175 return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
176 argument(2)*(-2.0*t+3.0)*t*t +
177 argument(3)*((t-2.0)*t+1.0)*t +
178 argument(4)*(t-1.0)*t*t );
179 }
180
181
182 double
183 l_bezier(char *nm)
184 {
185 double t;
186
187 t = argument(5);
188 return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
189 argument(2) * 3.*t*(1.+t*(-2.+t)) +
190 argument(3) * 3.*t*t*(1.-t) +
191 argument(4) * t*t*t );
192 }
193
194
195 double
196 l_bspline(char *nm)
197 {
198 double t;
199
200 t = argument(5);
201 return( argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
202 argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
203 argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
204 argument(4) * (1./6.*t*t*t) );
205 }