ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/gen/genrev.c
(Generate patch)

Comparing ray/src/gen/genrev.c (file contents):
Revision 1.4 by greg, Fri Mar 2 17:24:07 1990 UTC vs.
Revision 2.6 by greg, Mon Aug 4 22:37:53 2003 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1989 Regents of the University of California */
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ LBL";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6 –
4   /*
5   *  genrev.c - program to generate functions of rotation about z
6   *
# Line 14 | Line 11 | static char SCCSid[] = "$SunId$ LBL";
11   *      8/6/86
12   */
13  
14 + #include  <stdlib.h>
15   #include  <stdio.h>
16 + #include  <string.h>
17 + #include  <math.h>
18 + #include  "calcomp.h"
19  
20 < #define  ZNAME          "Z_"                    /* z function name */
21 < #define  RNAME          "R_"                    /* r function name */
20 > #define  ZNAME          "Z`SYS`"                /* z function name */
21 > #define  RNAME          "R`SYS`"                /* r function name */
22  
23   #define  PI             3.14159265358979323846
24  
# Line 28 | Line 29 | static char SCCSid[] = "$SunId$ LBL";
29   #define  UP             04
30   #define  DOWN           010
31  
32 < double  funvalue(), l_hermite(), l_bezier(), argument();
32 > void    quit(), eputs(), wputs();
33  
34  
35 + computen(nzp, nrp, z0, r0, z1, r1)              /* compute normal */
36 + double  *nzp, *nrp, z0, r0, z1, r1;
37 + {
38 +        double  dr, dz, len;
39 +
40 +        dz = r0 - r1;                           /* right angle vector */
41 +        dr = z1 - z0;
42 +        len = sqrt(dr*dr + dz*dz);
43 +        *nzp = dz/len;
44 +        *nrp = dr/len;
45 + }
46 +
47 +
48 + void
49 + eputs(msg)
50 + char  *msg;
51 + {
52 +        fputs(msg, stderr);
53 + }
54 +
55 +
56 + void
57 + wputs(msg)
58 + char  *msg;
59 + {
60 +        eputs(msg);
61 + }
62 +
63 +
64 + void
65 + quit(code)
66 + int  code;
67 + {
68 +        exit(code);
69 + }
70 +
71 +
72 + printhead(ac, av)               /* print command header */
73 + register int  ac;
74 + register char  **av;
75 + {
76 +        putchar('#');
77 +        while (ac--) {
78 +                putchar(' ');
79 +                fputs(*av++, stdout);
80 +        }
81 +        putchar('\n');
82 + }
83 +
84 +
85 + double
86 + l_hermite(char *nm)
87 + {
88 +        double  t;
89 +        
90 +        t = argument(5);
91 +        return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
92 +                argument(2)*(-2.0*t+3.0)*t*t +
93 +                argument(3)*((t-2.0)*t+1.0)*t +
94 +                argument(4)*(t-1.0)*t*t );
95 + }
96 +
97 +
98 + double
99 + l_bezier(char *nm)
100 + {
101 +        double  t;
102 +
103 +        t = argument(5);
104 +        return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
105 +                argument(2) * 3.*t*(1.+t*(-2.+t)) +
106 +                argument(3) * 3.*t*t*(1.-t) +
107 +                argument(4) * t*t*t );
108 + }
109 +
110 +
111 + double
112 + l_bspline(char *nm)
113 + {
114 +        double  t;
115 +
116 +        t = argument(5);
117 +        return( argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
118 +                argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
119 +                argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
120 +                argument(4) * (1./6.*t*t*t) );
121 + }
122 +
123 +
124   main(argc, argv)
125   int  argc;
126   char  *argv[];
# Line 43 | Line 133 | char  *argv[];
133          int  i, nseg;
134          int  orient;
135  
136 <        varset("PI", PI);
137 <        funset("hermite", 5, l_hermite);
138 <        funset("bezier", 5, l_bezier);
136 >        varset("PI", ':', PI);
137 >        funset("hermite", 5, ':', l_hermite);
138 >        funset("bezier", 5, ':', l_bezier);
139 >        funset("bspline", 5, ':', l_bspline);
140  
141          if (argc < 6)
142                  goto userror;
143  
144          for (i = 6; i < argc; i++)
145                  if (!strcmp(argv[i], "-e"))
146 <                        scompile(NULL, argv[++i]);
146 >                        scompile(argv[++i], NULL, 0);
147                  else if (!strcmp(argv[i], "-f"))
148                          fcompile(argv[++i]);
149                  else if (!strcmp(argv[i], "-s"))
# Line 61 | Line 152 | char  *argv[];
152                          goto userror;
153  
154          sprintf(stmp, "%s(t)=%s;", ZNAME, argv[3]);
155 <        scompile(NULL, stmp);
155 >        scompile(stmp, NULL, 0);
156          sprintf(stmp, "%s(t)=%s;", RNAME, argv[4]);
157 <        scompile(NULL, stmp);
157 >        scompile(stmp, NULL, 0);
158          nseg = atoi(argv[5]);
159          if (nseg <= 0)
160                  goto userror;
161          modname = smooth ? "Phong" : argv[1];
162  
163          printhead(argc, argv);
164 +        eclock = 0;
165  
166          lastnz = lastnr = 0.0;
167          t = 0.0;
# Line 154 | Line 246 | userror:
246   }
247  
248  
157 – computen(nzp, nrp, z0, r0, z1, r1)              /* compute normal */
158 – double  *nzp, *nrp, z0, r0, z1, r1;
159 – {
160 –        extern double  sqrt();
161 –        double  dr, dz, len;
162 –
163 –        dz = r0 - r1;                           /* right angle vector */
164 –        dr = z1 - z0;
165 –        len = sqrt(dr*dr + dz*dz);
166 –        *nzp = dz/len;
167 –        *nrp = dr/len;
168 – }
169 –
170 –
171 – eputs(msg)
172 – char  *msg;
173 – {
174 –        fputs(msg, stderr);
175 – }
176 –
177 –
178 – wputs(msg)
179 – char  *msg;
180 – {
181 –        eputs(msg);
182 – }
183 –
184 –
185 – quit(code)
186 – {
187 –        exit(code);
188 – }
189 –
190 –
191 – printhead(ac, av)               /* print command header */
192 – register int  ac;
193 – register char  **av;
194 – {
195 –        putchar('#');
196 –        while (ac--) {
197 –                putchar(' ');
198 –                fputs(*av++, stdout);
199 –        }
200 –        putchar('\n');
201 – }
202 –
203 –
204 – double
205 – l_hermite()                    
206 – {
207 –        double  t;
208 –        
209 –        t = argument(5);
210 –        return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
211 –                argument(2)*(-2.0*t+3.0)*t*t +
212 –                argument(3)*((t-2.0)*t+1.0)*t +
213 –                argument(4)*(t-1.0)*t*t );
214 – }
215 –
216 –
217 – double
218 – l_bezier()
219 – {
220 –        double  t;
221 –
222 –        t = argument(5);
223 –        return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
224 –                argument(2) * 3.*t*(1.+t*(-2.+t)) +
225 –                argument(3) * 3.*t*t*(1.-t) +
226 –                argument(4) * t*t*t );
227 – }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines