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.2 by greg, Sun Sep 10 16:45:14 1989 UTC vs.
Revision 2.5 by schorsch, Sun Jun 8 12:03:09 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  
19 < #define  ZNAME          "Z_"                    /* z function name */
20 < #define  RNAME          "R_"                    /* r function name */
19 > #define  ZNAME          "Z`SYS`"                /* z function name */
20 > #define  RNAME          "R`SYS`"                /* r function name */
21  
22   #define  PI             3.14159265358979323846
23  
# Line 28 | Line 28 | static char SCCSid[] = "$SunId$ LBL";
28   #define  UP             04
29   #define  DOWN           010
30  
31 < double  funvalue(), l_hermite(), argument();
31 > double  funvalue(), l_hermite(), l_bezier(), l_bspline(), argument();
32  
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[];
128   {
129 +        extern long     eclock;
130          char  stmp[256];
131          char  *modname;
132          int  smooth = 0;
# Line 43 | Line 135 | char  *argv[];
135          int  i, nseg;
136          int  orient;
137  
138 <        varset("PI", PI);
139 <        funset("hermite", 5, l_hermite);
138 >        varset("PI", ':', PI);
139 >        funset("hermite", 5, ':', l_hermite);
140 >        funset("bezier", 5, ':', l_bezier);
141 >        funset("bspline", 5, ':', l_bspline);
142  
143          if (argc < 6)
144                  goto userror;
145  
146          for (i = 6; i < argc; i++)
147                  if (!strcmp(argv[i], "-e"))
148 <                        scompile(NULL, argv[++i]);
148 >                        scompile(argv[++i], NULL, 0);
149                  else if (!strcmp(argv[i], "-f"))
150                          fcompile(argv[++i]);
151                  else if (!strcmp(argv[i], "-s"))
# Line 60 | Line 154 | char  *argv[];
154                          goto userror;
155  
156          sprintf(stmp, "%s(t)=%s;", ZNAME, argv[3]);
157 <        scompile(NULL, stmp);
157 >        scompile(stmp, NULL, 0);
158          sprintf(stmp, "%s(t)=%s;", RNAME, argv[4]);
159 <        scompile(NULL, stmp);
159 >        scompile(stmp, NULL, 0);
160          nseg = atoi(argv[5]);
161          if (nseg <= 0)
162                  goto userror;
163 <        modname = smooth ? "phong" : argv[1];
163 >        modname = smooth ? "Phong" : argv[1];
164  
165          printhead(argc, argv);
166 +        eclock = 0;
167  
168          lastnz = lastnr = 0.0;
169          t = 0.0;
# Line 98 | Line 193 | char  *argv[];
193                  if (!orient)
194                          goto endfor;
195                  if (smooth) {
196 <                        printf("\n%s texfunc phong\n", argv[1]);
196 >                        printf("\n%s texfunc Phong\n", argv[1]);
197                          printf("4 rev_dx rev_dy rev_dz rev.cal\n");
198                          printf("0\n4\n");
199                          if (orient&(UP|DOWN)) {
# Line 153 | Line 248 | userror:
248   }
249  
250  
156 – computen(nzp, nrp, z0, r0, z1, r1)              /* compute normal */
157 – double  *nzp, *nrp, z0, r0, z1, r1;
158 – {
159 –        extern double  sqrt();
160 –        double  dr, dz, len;
161 –
162 –        dz = r0 - r1;                           /* right angle vector */
163 –        dr = z1 - z0;
164 –        len = sqrt(dr*dr + dz*dz);
165 –        *nzp = dz/len;
166 –        *nrp = dr/len;
167 – }
168 –
169 –
170 – eputs(msg)
171 – char  *msg;
172 – {
173 –        fputs(msg, stderr);
174 – }
175 –
176 –
177 – wputs(msg)
178 – char  *msg;
179 – {
180 –        eputs(msg);
181 – }
182 –
183 –
184 – quit(code)
185 – {
186 –        exit(code);
187 – }
188 –
189 –
190 – printhead(ac, av)               /* print command header */
191 – register int  ac;
192 – register char  **av;
193 – {
194 –        putchar('#');
195 –        while (ac--) {
196 –                putchar(' ');
197 –                fputs(*av++, stdout);
198 –        }
199 –        putchar('\n');
200 – }
201 –
202 –
203 – double
204 – l_hermite()                    
205 – {
206 –        double  t;
207 –        
208 –        t = argument(5);
209 –        return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
210 –                argument(2)*(-2.0*t+3.0)*t*t +
211 –                argument(3)*((t-2.0)*t+1.0)*t +
212 –                argument(4)*(t-1.0)*t*t );
213 – }

Diff Legend

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