ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/gen/genrev.c
Revision: 2.3
Committed: Thu Nov 18 09:33:01 1993 UTC (30 years, 4 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.2: +1 -0 lines
Log Message:
minor compiler warning fixes

File Contents

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