ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/gen/genrev.c
Revision: 2.5
Committed: Sun Jun 8 12:03:09 2003 UTC (20 years, 9 months ago) by schorsch
Content type: text/plain
Branch: MAIN
Changes since 2.4: +92 -88 lines
Log Message:
Reduced compile warnings/errors on Windows.

File Contents

# User Rev Content
1 greg 1.1 #ifndef lint
2 schorsch 2.5 static const char RCSid[] = "$Id: genrev.c,v 2.4 2003/02/22 02:07:23 greg Exp $";
3 greg 1.1 #endif
4 greg 1.2 /*
5 greg 1.1 * genrev.c - program to generate functions of rotation about z
6     *
7     * The program takes as input the functions of t for z and r
8     * (the radius). If z is increasing, the normal points
9     * outward, inward for decreasing z. Negative radii are forbidden.
10     *
11     * 8/6/86
12     */
13    
14 schorsch 2.5 #include <stdlib.h>
15 greg 1.1 #include <stdio.h>
16 schorsch 2.5 #include <string.h>
17 greg 2.2 #include <math.h>
18 greg 1.1
19 greg 1.10 #define ZNAME "Z`SYS`" /* z function name */
20     #define RNAME "R`SYS`" /* r function name */
21 greg 1.1
22     #define PI 3.14159265358979323846
23    
24     #define FTINY 1e-9
25     /* orientation */
26     #define OUT 01
27     #define IN 02
28     #define UP 04
29     #define DOWN 010
30    
31 greg 1.5 double funvalue(), l_hermite(), l_bezier(), l_bspline(), argument();
32 greg 1.1
33 greg 2.4 void quit(), eputs(), wputs();
34    
35 greg 1.1
36 schorsch 2.5 computen(nzp, nrp, z0, r0, z1, r1) /* compute normal */
37     double *nzp, *nrp, z0, r0, z1, r1;
38     {
39     double dr, dz, len;
40    
41     dz = r0 - r1; /* right angle vector */
42     dr = z1 - z0;
43     len = sqrt(dr*dr + dz*dz);
44     *nzp = dz/len;
45     *nrp = dr/len;
46     }
47    
48    
49     void
50     eputs(msg)
51     char *msg;
52     {
53     fputs(msg, stderr);
54     }
55    
56    
57     void
58     wputs(msg)
59     char *msg;
60     {
61     eputs(msg);
62     }
63    
64    
65     void
66     quit(code)
67     int code;
68     {
69     exit(code);
70     }
71    
72    
73     printhead(ac, av) /* print command header */
74     register int ac;
75     register char **av;
76     {
77     putchar('#');
78     while (ac--) {
79     putchar(' ');
80     fputs(*av++, stdout);
81     }
82     putchar('\n');
83     }
84    
85    
86     double
87     l_hermite()
88     {
89     double t;
90    
91     t = argument(5);
92     return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
93     argument(2)*(-2.0*t+3.0)*t*t +
94     argument(3)*((t-2.0)*t+1.0)*t +
95     argument(4)*(t-1.0)*t*t );
96     }
97    
98    
99     double
100     l_bezier()
101     {
102     double t;
103    
104     t = argument(5);
105     return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
106     argument(2) * 3.*t*(1.+t*(-2.+t)) +
107     argument(3) * 3.*t*t*(1.-t) +
108     argument(4) * t*t*t );
109     }
110    
111    
112     double
113     l_bspline()
114     {
115     double t;
116    
117     t = argument(5);
118     return( argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
119     argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
120     argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
121     argument(4) * (1./6.*t*t*t) );
122     }
123    
124    
125 greg 1.1 main(argc, argv)
126     int argc;
127     char *argv[];
128     {
129 greg 1.6 extern long eclock;
130 greg 1.1 char stmp[256];
131     char *modname;
132     int smooth = 0;
133     double t, lastz, z, nextz, lastr, r, nextr;
134     double lastnz, lastnr, nz, nr, nextnz, nextnr;
135     int i, nseg;
136     int orient;
137    
138 greg 1.8 varset("PI", ':', PI);
139 greg 1.9 funset("hermite", 5, ':', l_hermite);
140     funset("bezier", 5, ':', l_bezier);
141     funset("bspline", 5, ':', l_bspline);
142 greg 1.1
143     if (argc < 6)
144     goto userror;
145    
146     for (i = 6; i < argc; i++)
147     if (!strcmp(argv[i], "-e"))
148 greg 1.7 scompile(argv[++i], NULL, 0);
149 greg 1.1 else if (!strcmp(argv[i], "-f"))
150     fcompile(argv[++i]);
151     else if (!strcmp(argv[i], "-s"))
152     smooth = 1;
153     else
154     goto userror;
155    
156     sprintf(stmp, "%s(t)=%s;", ZNAME, argv[3]);
157 greg 1.7 scompile(stmp, NULL, 0);
158 greg 1.1 sprintf(stmp, "%s(t)=%s;", RNAME, argv[4]);
159 greg 1.7 scompile(stmp, NULL, 0);
160 greg 1.1 nseg = atoi(argv[5]);
161     if (nseg <= 0)
162     goto userror;
163 greg 1.3 modname = smooth ? "Phong" : argv[1];
164 greg 1.1
165     printhead(argc, argv);
166 greg 1.6 eclock = 0;
167 greg 1.1
168     lastnz = lastnr = 0.0;
169     t = 0.0;
170     lastz = funvalue(ZNAME, 1, &t);
171     lastr = funvalue(RNAME, 1, &t);
172     t = 1.0/nseg;
173     z = funvalue(ZNAME, 1, &t);
174     r = funvalue(RNAME, 1, &t);
175     computen(&nz, &nr, lastz, lastr, z, r);
176     for (i = 1; i <= nseg; i++) {
177     if (i < nseg) {
178     t = (double)(i+1)/nseg;
179     nextz = funvalue(ZNAME, 1, &t);
180     nextr = funvalue(RNAME, 1, &t);
181     computen(&nextnz, &nextnr, z, r, nextz, nextr);
182     } else
183     nextnz = nextnr = 0.0;
184     orient = 0;
185     if (z < lastz-FTINY)
186     orient |= DOWN;
187     else if (z > lastz+FTINY)
188     orient |= UP;
189     if (r < lastr-FTINY)
190     orient |= IN;
191     else if (r > lastr+FTINY)
192     orient |= OUT;
193     if (!orient)
194     goto endfor;
195     if (smooth) {
196 greg 1.3 printf("\n%s texfunc Phong\n", argv[1]);
197 greg 1.1 printf("4 rev_dx rev_dy rev_dz rev.cal\n");
198     printf("0\n4\n");
199     if (orient&(UP|DOWN)) {
200     t = (nextnz - lastnz)/(z - lastz);
201     printf("\t%15.9g\t%15.9g\n",
202     t, lastnz - t*lastz);
203     } else
204     printf("\t0\t%d\n", orient&IN ? 1 : -1);
205     if (orient&(OUT|IN)) {
206     t = (nextnr - lastnr)/(r - lastr);
207     printf("\t%15.9g\t%15.9g\n",
208     t, lastnr - t*lastr);
209     } else
210     printf("\t0\t%d\n", orient&UP ? 1 : -1);
211     }
212     if (!(orient&(IN|OUT))) {
213     printf("\n%s %s %s.%d\n", modname,
214     orient&DOWN ? "tube" : "cylinder",
215     argv[2], i);
216     printf("0\n0\n7\n");
217     printf("\t0\t0\t%15.9g\n", lastz);
218     printf("\t0\t0\t%15.9g\n", z);
219     printf("\t%15.9g\n", r);
220     } else if (!(orient&(UP|DOWN))) {
221     printf("\n%s ring %s.%d\n", modname, argv[2], i);
222     printf("0\n0\n8\n");
223     printf("\t0\t0\t%15.9g\n", z);
224     printf("\t0\t0\t%15.9g\n", orient&IN ? 1.0 : -1.0);
225     printf("\t%15.9g\t%15.9g\n", lastr, r);
226     } else {
227     printf("\n%s %s %s.%d\n", modname,
228     orient&DOWN ? "cup" : "cone",
229     argv[2], i);
230     printf("0\n0\n8\n");
231     printf("\t0\t0\t%15.9g\n", lastz);
232     printf("\t0\t0\t%15.9g\n", z);
233     printf("\t%15.9g\t%15.9g\n", lastr, r);
234     }
235     endfor:
236     lastz = z; lastr = r;
237     z = nextz; r = nextr;
238     lastnz = nz; lastnr = nr;
239     nz = nextnz; nr = nextnr;
240     }
241     quit(0);
242    
243     userror:
244     fprintf(stderr,
245     "Usage: %s material name z(t) r(t) nseg [-e expr] [-f file] [-s]\n",
246     argv[0]);
247     quit(1);
248     }
249    
250