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.5 by greg, Wed Mar 7 11:14:40 1990 UTC vs.
Revision 2.17 by greg, Wed Apr 23 01:57:04 2025 UTC

# Line 1 | Line 1
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ LBL";
2 > static const char       RCSid[] = "$Id$";
3   #endif
4 –
5 – /* Copyright (c) 1989 Regents of the University of California */
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  <stdio.h>
14 > #include  <stdlib.h>
15 > #include  <math.h>
16  
17 < #define  ZNAME          "Z_"                    /* z function name */
18 < #define  RNAME          "R_"                    /* r function name */
17 > #include  "rtio.h"
18 > #include  "rterror.h"
19 > #include  "resolu.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 28 | Line 31 | static char SCCSid[] = "$SunId$ LBL";
31   #define  UP             04
32   #define  DOWN           010
33  
31 – double  funvalue(), l_hermite(), l_bezier(), l_bspline(), argument();
34  
35 + void
36 + computen(                       /* compute normal */
37 +        double  *nzp,
38 +        double  *nrp,
39 +        double z0,
40 +        double r0,
41 +        double z1,
42 +        double r1
43 + )
44 + {
45 +        double  dr, dz, len;
46  
47 < main(argc, argv)
48 < int  argc;
49 < char  *argv[];
47 >        dz = r0 - r1;                           /* right angle vector */
48 >        dr = z1 - z0;
49 >        len = sqrt(dr*dr + dz*dz);
50 >        *nzp = dz/len;
51 >        *nrp = dr/len;
52 > }
53 >
54 >
55 > double
56 > l_hermite(char *nm)
57   {
58 +        double  t;
59 +        
60 +        t = argument(5);
61 +        return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
62 +                argument(2)*(-2.0*t+3.0)*t*t +
63 +                argument(3)*((t-2.0)*t+1.0)*t +
64 +                argument(4)*(t-1.0)*t*t );
65 + }
66 +
67 +
68 + double
69 + l_bezier(char *nm)
70 + {
71 +        double  t;
72 +
73 +        t = argument(5);
74 +        return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
75 +                argument(2) * 3.*t*(1.+t*(-2.+t)) +
76 +                argument(3) * 3.*t*t*(1.-t) +
77 +                argument(4) * t*t*t );
78 + }
79 +
80 +
81 + double
82 + l_bspline(char *nm)
83 + {
84 +        double  t;
85 +
86 +        t = argument(5);
87 +        return( argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
88 +                argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
89 +                argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
90 +                argument(4) * (1./6.*t*t*t) );
91 + }
92 +
93 +
94 + int
95 + main(
96 +        int  argc,
97 +        char  *argv[]
98 + )
99 + {
100          char  stmp[256];
101          char  *modname;
102          int  smooth = 0;
# Line 43 | Line 105 | char  *argv[];
105          int  i, nseg;
106          int  orient;
107  
108 <        varset("PI", PI);
109 <        funset("hermite", 5, l_hermite);
110 <        funset("bezier", 5, l_bezier);
111 <        funset("bspline", 5, l_bspline);
108 >        esupport |= E_VARIABLE|E_FUNCTION|E_RCONST;
109 >        esupport &= ~(E_OUTCHAN|E_INCHAN);
110 >        varset("PI", ':', PI);
111 >        funset("hermite", 5, ':', l_hermite);
112 >        funset("bezier", 5, ':', l_bezier);
113 >        funset("bspline", 5, ':', l_bspline);
114  
115          if (argc < 6)
116                  goto userror;
117  
118          for (i = 6; i < argc; i++)
119                  if (!strcmp(argv[i], "-e"))
120 <                        scompile(NULL, argv[++i]);
121 <                else if (!strcmp(argv[i], "-f"))
122 <                        fcompile(argv[++i]);
123 <                else if (!strcmp(argv[i], "-s"))
120 >                        scompile(argv[++i], NULL, 0);
121 >                else if (!strcmp(argv[i], "-f")) {
122 >                        char  *fpath = getpath(argv[++i], getrlibpath(), 0);
123 >                        if (fpath == NULL) {
124 >                                fprintf(stderr, "%s: cannot find file '%s'\n",
125 >                                                argv[0], argv[i]);
126 >                                quit(1);
127 >                        }
128 >                        fcompile(fpath);
129 >                } else if (!strcmp(argv[i], "-s"))
130                          smooth = 1;
131                  else
132                          goto userror;
133  
134          sprintf(stmp, "%s(t)=%s;", ZNAME, argv[3]);
135 <        scompile(NULL, stmp);
135 >        scompile(stmp, NULL, 0);
136          sprintf(stmp, "%s(t)=%s;", RNAME, argv[4]);
137 <        scompile(NULL, stmp);
138 <        nseg = atoi(argv[5]);
137 >        scompile(stmp, NULL, 0);
138 >        nseg = eval(argv[5]) + .5;
139          if (nseg <= 0)
140                  goto userror;
141          modname = smooth ? "Phong" : argv[1];
142  
143 <        printhead(argc, argv);
143 >        fputs("# ", stdout);
144 >        printargs(argc, argv, stdout);
145 >        doptimize(1);
146 >        eclock++;
147  
148          lastnz = lastnr = 0.0;
149          t = 0.0;
# Line 105 | Line 178 | char  *argv[];
178                          printf("0\n4\n");
179                          if (orient&(UP|DOWN)) {
180                                  t = (nextnz - lastnz)/(z - lastz);
181 <                                printf("\t%15.9g\t%15.9g\n",
181 >                                printf("\t%18.12g\t%18.12g\n",
182                                                  t, lastnz - t*lastz);
183                          } else
184                                  printf("\t0\t%d\n", orient&IN ? 1 : -1);
185                          if (orient&(OUT|IN))  {
186                                  t = (nextnr - lastnr)/(r - lastr);
187 <                                printf("\t%15.9g\t%15.9g\n",
187 >                                printf("\t%18.12g\t%18.12g\n",
188                                                  t, lastnr - t*lastr);
189                          } else
190                                  printf("\t0\t%d\n", orient&UP ? 1 : -1);
# Line 121 | Line 194 | char  *argv[];
194                                          orient&DOWN ? "tube" : "cylinder",
195                                          argv[2], i);
196                          printf("0\n0\n7\n");
197 <                        printf("\t0\t0\t%15.9g\n", lastz);
198 <                        printf("\t0\t0\t%15.9g\n", z);
199 <                        printf("\t%15.9g\n", r);
197 >                        printf("\t0\t0\t%18.12g\n", lastz);
198 >                        printf("\t0\t0\t%18.12g\n", z);
199 >                        printf("\t%18.12g\n", r);
200                  } else if (!(orient&(UP|DOWN))) {
201                          printf("\n%s ring %s.%d\n", modname, argv[2], i);
202                          printf("0\n0\n8\n");
203 <                        printf("\t0\t0\t%15.9g\n", z);
204 <                        printf("\t0\t0\t%15.9g\n", orient&IN ? 1.0 : -1.0);
205 <                        printf("\t%15.9g\t%15.9g\n", lastr, r);
203 >                        printf("\t0\t0\t%18.12g\n", z);
204 >                        printf("\t0\t0\t%18.12g\n", orient&IN ? 1.0 : -1.0);
205 >                        printf("\t%18.12g\t%18.12g\n", lastr, r);
206                  } else {
207                          printf("\n%s %s %s.%d\n", modname,
208                                          orient&DOWN ? "cup" : "cone",
209                                          argv[2], i);
210                          printf("0\n0\n8\n");
211 <                        printf("\t0\t0\t%15.9g\n", lastz);
212 <                        printf("\t0\t0\t%15.9g\n", z);
213 <                        printf("\t%15.9g\t%15.9g\n", lastr, r);
211 >                        printf("\t0\t0\t%18.12g\n", lastz);
212 >                        printf("\t0\t0\t%18.12g\n", z);
213 >                        printf("\t%18.12g\t%18.12g\n", lastr, r);
214                  }
215          endfor:
216                  lastz = z; lastr = r;
# Line 145 | Line 218 | char  *argv[];
218                  lastnz = nz; lastnr = nr;
219                  nz = nextnz; nr = nextnr;
220          }
221 <        quit(0);
221 >        return 0;
222  
223   userror:
224          fprintf(stderr,
225          "Usage: %s material name z(t) r(t) nseg [-e expr] [-f file] [-s]\n",
226                          argv[0]);
227 <        quit(1);
155 < }
156 <
157 <
158 < computen(nzp, nrp, z0, r0, z1, r1)              /* compute normal */
159 < double  *nzp, *nrp, z0, r0, z1, r1;
160 < {
161 <        extern double  sqrt();
162 <        double  dr, dz, len;
163 <
164 <        dz = r0 - r1;                           /* right angle vector */
165 <        dr = z1 - z0;
166 <        len = sqrt(dr*dr + dz*dz);
167 <        *nzp = dz/len;
168 <        *nrp = dr/len;
169 < }
170 <
171 <
172 < eputs(msg)
173 < char  *msg;
174 < {
175 <        fputs(msg, stderr);
176 < }
177 <
178 <
179 < wputs(msg)
180 < char  *msg;
181 < {
182 <        eputs(msg);
183 < }
184 <
185 <
186 < quit(code)
187 < {
188 <        exit(code);
189 < }
190 <
191 <
192 < printhead(ac, av)               /* print command header */
193 < register int  ac;
194 < register char  **av;
195 < {
196 <        putchar('#');
197 <        while (ac--) {
198 <                putchar(' ');
199 <                fputs(*av++, stdout);
200 <        }
201 <        putchar('\n');
202 < }
203 <
204 <
205 < double
206 < l_hermite()                    
207 < {
208 <        double  t;
209 <        
210 <        t = argument(5);
211 <        return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
212 <                argument(2)*(-2.0*t+3.0)*t*t +
213 <                argument(3)*((t-2.0)*t+1.0)*t +
214 <                argument(4)*(t-1.0)*t*t );
215 < }
216 <
217 <
218 < double
219 < l_bezier()
220 < {
221 <        double  t;
222 <
223 <        t = argument(5);
224 <        return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
225 <                argument(2) * 3.*t*(1.+t*(-2.+t)) +
226 <                argument(3) * 3.*t*t*(1.-t) +
227 <                argument(4) * t*t*t );
228 < }
229 <
230 <
231 < double
232 < l_bspline()
233 < {
234 <        double  t;
235 <
236 <        t = argument(5);
237 <        return( argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
238 <                argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
239 <                argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
240 <                argument(4) * (1./6.*t*t*t) );
227 >        return 1;
228   }

Diff Legend

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