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 2.4 by greg, Sat Feb 22 02:07:23 2003 UTC vs.
Revision 2.5 by schorsch, Sun Jun 8 12:03:09 2003 UTC

# Line 11 | Line 11 | static const char      RCSid[] = "$Id$";
11   *      8/6/86
12   */
13  
14 + #include  <stdlib.h>
15   #include  <stdio.h>
16 + #include  <string.h>
17   #include  <math.h>
18  
19   #define  ZNAME          "Z`SYS`"                /* z function name */
# Line 31 | Line 33 | double  funvalue(), l_hermite(), l_bezier(), l_bspline
33   void    quit(), eputs(), wputs();
34  
35  
36 + 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   main(argc, argv)
126   int  argc;
127   char  *argv[];
# Line 157 | Line 248 | userror:
248   }
249  
250  
160 – computen(nzp, nrp, z0, r0, z1, r1)              /* compute normal */
161 – double  *nzp, *nrp, z0, r0, z1, r1;
162 – {
163 –        double  dr, dz, len;
164 –
165 –        dz = r0 - r1;                           /* right angle vector */
166 –        dr = z1 - z0;
167 –        len = sqrt(dr*dr + dz*dz);
168 –        *nzp = dz/len;
169 –        *nrp = dr/len;
170 – }
171 –
172 –
173 – void
174 – eputs(msg)
175 – char  *msg;
176 – {
177 –        fputs(msg, stderr);
178 – }
179 –
180 –
181 – void
182 – wputs(msg)
183 – char  *msg;
184 – {
185 –        eputs(msg);
186 – }
187 –
188 –
189 – void
190 – quit(code)
191 – int  code;
192 – {
193 –        exit(code);
194 – }
195 –
196 –
197 – printhead(ac, av)               /* print command header */
198 – register int  ac;
199 – register char  **av;
200 – {
201 –        putchar('#');
202 –        while (ac--) {
203 –                putchar(' ');
204 –                fputs(*av++, stdout);
205 –        }
206 –        putchar('\n');
207 – }
208 –
209 –
210 – double
211 – l_hermite()                    
212 – {
213 –        double  t;
214 –        
215 –        t = argument(5);
216 –        return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
217 –                argument(2)*(-2.0*t+3.0)*t*t +
218 –                argument(3)*((t-2.0)*t+1.0)*t +
219 –                argument(4)*(t-1.0)*t*t );
220 – }
221 –
222 –
223 – double
224 – l_bezier()
225 – {
226 –        double  t;
227 –
228 –        t = argument(5);
229 –        return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
230 –                argument(2) * 3.*t*(1.+t*(-2.+t)) +
231 –                argument(3) * 3.*t*t*(1.-t) +
232 –                argument(4) * t*t*t );
233 – }
234 –
235 –
236 – double
237 – l_bspline()
238 – {
239 –        double  t;
240 –
241 –        t = argument(5);
242 –        return( argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
243 –                argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
244 –                argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
245 –                argument(4) * (1./6.*t*t*t) );
246 – }

Diff Legend

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