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.1 by greg, Thu Feb 2 11:16:29 1989 UTC vs.
Revision 2.7 by schorsch, Sun Nov 16 10:29:38 2003 UTC

# Line 1 | Line 1
1 – /*
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ LBL";
2 > static const char       RCSid[] = "$Id$";
3   #endif
4 + /*
5   *  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
# Line 12 | 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 > #include  "rterror.h"
20 > #include  "calcomp.h"
21  
22 + #define  ZNAME          "Z`SYS`"                /* z function name */
23 + #define  RNAME          "R`SYS`"                /* r function name */
24 +
25   #define  PI             3.14159265358979323846
26  
27   #define  FTINY          1e-9
# Line 26 | Line 31 | static char SCCSid[] = "$SunId$ LBL";
31   #define  UP             04
32   #define  DOWN           010
33  
29 – double  funvalue(), l_hermite(), argument();
34  
35 + void
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 + void
74 + printhead(ac, av)               /* print command header */
75 + register int  ac;
76 + register char  **av;
77 + {
78 +        putchar('#');
79 +        while (ac--) {
80 +                putchar(' ');
81 +                fputs(*av++, stdout);
82 +        }
83 +        putchar('\n');
84 + }
85 +
86 +
87 + double
88 + l_hermite(char *nm)
89 + {
90 +        double  t;
91 +        
92 +        t = argument(5);
93 +        return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
94 +                argument(2)*(-2.0*t+3.0)*t*t +
95 +                argument(3)*((t-2.0)*t+1.0)*t +
96 +                argument(4)*(t-1.0)*t*t );
97 + }
98 +
99 +
100 + double
101 + l_bezier(char *nm)
102 + {
103 +        double  t;
104 +
105 +        t = argument(5);
106 +        return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
107 +                argument(2) * 3.*t*(1.+t*(-2.+t)) +
108 +                argument(3) * 3.*t*t*(1.-t) +
109 +                argument(4) * t*t*t );
110 + }
111 +
112 +
113 + double
114 + l_bspline(char *nm)
115 + {
116 +        double  t;
117 +
118 +        t = argument(5);
119 +        return( argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
120 +                argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
121 +                argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
122 +                argument(4) * (1./6.*t*t*t) );
123 + }
124 +
125 +
126 + int
127   main(argc, argv)
128   int  argc;
129   char  *argv[];
# Line 41 | Line 136 | char  *argv[];
136          int  i, nseg;
137          int  orient;
138  
139 <        varset("PI", PI);
140 <        funset("hermite", 5, l_hermite);
139 >        varset("PI", ':', PI);
140 >        funset("hermite", 5, ':', l_hermite);
141 >        funset("bezier", 5, ':', l_bezier);
142 >        funset("bspline", 5, ':', l_bspline);
143  
144          if (argc < 6)
145                  goto userror;
146  
147          for (i = 6; i < argc; i++)
148                  if (!strcmp(argv[i], "-e"))
149 <                        scompile(NULL, argv[++i]);
149 >                        scompile(argv[++i], NULL, 0);
150                  else if (!strcmp(argv[i], "-f"))
151                          fcompile(argv[++i]);
152                  else if (!strcmp(argv[i], "-s"))
# Line 58 | Line 155 | char  *argv[];
155                          goto userror;
156  
157          sprintf(stmp, "%s(t)=%s;", ZNAME, argv[3]);
158 <        scompile(NULL, stmp);
158 >        scompile(stmp, NULL, 0);
159          sprintf(stmp, "%s(t)=%s;", RNAME, argv[4]);
160 <        scompile(NULL, stmp);
160 >        scompile(stmp, NULL, 0);
161          nseg = atoi(argv[5]);
162          if (nseg <= 0)
163                  goto userror;
164 <        modname = smooth ? "phong" : argv[1];
164 >        modname = smooth ? "Phong" : argv[1];
165  
166          printhead(argc, argv);
167 +        eclock = 0;
168  
169          lastnz = lastnr = 0.0;
170          t = 0.0;
# Line 96 | Line 194 | char  *argv[];
194                  if (!orient)
195                          goto endfor;
196                  if (smooth) {
197 <                        printf("\n%s texfunc phong\n", argv[1]);
197 >                        printf("\n%s texfunc Phong\n", argv[1]);
198                          printf("4 rev_dx rev_dy rev_dz rev.cal\n");
199                          printf("0\n4\n");
200                          if (orient&(UP|DOWN)) {
# Line 141 | Line 239 | char  *argv[];
239                  lastnz = nz; lastnr = nr;
240                  nz = nextnz; nr = nextnr;
241          }
242 <        quit(0);
242 >        return 0;
243  
244   userror:
245          fprintf(stderr,
246          "Usage: %s material name z(t) r(t) nseg [-e expr] [-f file] [-s]\n",
247                          argv[0]);
248 <        quit(1);
248 >        return 1;
249   }
250  
251  
154 – computen(nzp, nrp, z0, r0, z1, r1)              /* compute normal */
155 – double  *nzp, *nrp, z0, r0, z1, r1;
156 – {
157 –        extern double  sqrt();
158 –        double  dr, dz, len;
159 –
160 –        dz = r0 - r1;                           /* right angle vector */
161 –        dr = z1 - z0;
162 –        len = sqrt(dr*dr + dz*dz);
163 –        *nzp = dz/len;
164 –        *nrp = dr/len;
165 – }
166 –
167 –
168 – eputs(msg)
169 – char  *msg;
170 – {
171 –        fputs(msg, stderr);
172 – }
173 –
174 –
175 – wputs(msg)
176 – char  *msg;
177 – {
178 –        eputs(msg);
179 – }
180 –
181 –
182 – quit(code)
183 – {
184 –        exit(code);
185 – }
186 –
187 –
188 – printhead(ac, av)               /* print command header */
189 – register int  ac;
190 – register char  **av;
191 – {
192 –        putchar('#');
193 –        while (ac--) {
194 –                putchar(' ');
195 –                fputs(*av++, stdout);
196 –        }
197 –        putchar('\n');
198 – }
199 –
200 –
201 – double
202 – l_hermite()                    
203 – {
204 –        double  t;
205 –        
206 –        t = argument(5);
207 –        return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
208 –                argument(2)*(-2.0*t+3.0)*t*t +
209 –                argument(3)*((t-2.0)*t+1.0)*t +
210 –                argument(4)*(t-1.0)*t*t );
211 – }

Diff Legend

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