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.14 by greg, Sat Dec 28 18:05:14 2019 UTC

# Line 11 | Line 11 | static const char      RCSid[] = "$Id$";
11   *      8/6/86
12   */
13  
14 < #include  <stdio.h>
14 > #include  <stdlib.h>
15   #include  <math.h>
16  
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  
# Line 26 | Line 31 | static const char      RCSid[] = "$Id$";
31   #define  UP             04
32   #define  DOWN           010
33  
29 – double  funvalue(), l_hermite(), l_bezier(), l_bspline(), argument();
34  
35 < void    quit(), eputs(), wputs();
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 + double
50 + l_hermite(char *nm)
51 + {
52 +        double  t;
53 +        
54 +        t = argument(5);
55 +        return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
56 +                argument(2)*(-2.0*t+3.0)*t*t +
57 +                argument(3)*((t-2.0)*t+1.0)*t +
58 +                argument(4)*(t-1.0)*t*t );
59 + }
60 +
61 +
62 + double
63 + l_bezier(char *nm)
64 + {
65 +        double  t;
66 +
67 +        t = argument(5);
68 +        return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
69 +                argument(2) * 3.*t*(1.+t*(-2.+t)) +
70 +                argument(3) * 3.*t*t*(1.-t) +
71 +                argument(4) * t*t*t );
72 + }
73 +
74 +
75 + double
76 + l_bspline(char *nm)
77 + {
78 +        double  t;
79 +
80 +        t = argument(5);
81 +        return( argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
82 +                argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
83 +                argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
84 +                argument(4) * (1./6.*t*t*t) );
85 + }
86 +
87 +
88 + int
89   main(argc, argv)
90   int  argc;
91   char  *argv[];
92   {
38 –        extern long     eclock;
93          char  stmp[256];
94          char  *modname;
95          int  smooth = 0;
# Line 44 | Line 98 | char  *argv[];
98          int  i, nseg;
99          int  orient;
100  
101 +        esupport |= E_VARIABLE|E_FUNCTION|E_RCONST;
102 +        esupport &= ~(E_OUTCHAN|E_INCHAN);
103          varset("PI", ':', PI);
104          funset("hermite", 5, ':', l_hermite);
105          funset("bezier", 5, ':', l_bezier);
# Line 55 | Line 111 | char  *argv[];
111          for (i = 6; i < argc; i++)
112                  if (!strcmp(argv[i], "-e"))
113                          scompile(argv[++i], NULL, 0);
114 <                else if (!strcmp(argv[i], "-f"))
115 <                        fcompile(argv[++i]);
116 <                else if (!strcmp(argv[i], "-s"))
114 >                else if (!strcmp(argv[i], "-f")) {
115 >                        char  *fpath = getpath(argv[++i], getrlibpath(), 0);
116 >                        if (fpath == NULL) {
117 >                                fprintf(stderr, "%s: cannot find file '%s'\n",
118 >                                                argv[0], argv[i]);
119 >                                quit(1);
120 >                        }
121 >                        fcompile(fpath);
122 >                } else if (!strcmp(argv[i], "-s"))
123                          smooth = 1;
124                  else
125                          goto userror;
# Line 66 | Line 128 | char  *argv[];
128          scompile(stmp, NULL, 0);
129          sprintf(stmp, "%s(t)=%s;", RNAME, argv[4]);
130          scompile(stmp, NULL, 0);
131 <        nseg = atoi(argv[5]);
131 >        nseg = eval(argv[5]) + .5;
132          if (nseg <= 0)
133                  goto userror;
134          modname = smooth ? "Phong" : argv[1];
135  
136 <        printhead(argc, argv);
136 >        fputs("# ", stdout);
137 >        printargs(argc, argv, stdout);
138          eclock = 0;
139  
140          lastnz = lastnr = 0.0;
# Line 107 | Line 170 | char  *argv[];
170                          printf("0\n4\n");
171                          if (orient&(UP|DOWN)) {
172                                  t = (nextnz - lastnz)/(z - lastz);
173 <                                printf("\t%15.9g\t%15.9g\n",
173 >                                printf("\t%18.12g\t%18.12g\n",
174                                                  t, lastnz - t*lastz);
175                          } else
176                                  printf("\t0\t%d\n", orient&IN ? 1 : -1);
177                          if (orient&(OUT|IN))  {
178                                  t = (nextnr - lastnr)/(r - lastr);
179 <                                printf("\t%15.9g\t%15.9g\n",
179 >                                printf("\t%18.12g\t%18.12g\n",
180                                                  t, lastnr - t*lastr);
181                          } else
182                                  printf("\t0\t%d\n", orient&UP ? 1 : -1);
# Line 123 | Line 186 | char  *argv[];
186                                          orient&DOWN ? "tube" : "cylinder",
187                                          argv[2], i);
188                          printf("0\n0\n7\n");
189 <                        printf("\t0\t0\t%15.9g\n", lastz);
190 <                        printf("\t0\t0\t%15.9g\n", z);
191 <                        printf("\t%15.9g\n", r);
189 >                        printf("\t0\t0\t%18.12g\n", lastz);
190 >                        printf("\t0\t0\t%18.12g\n", z);
191 >                        printf("\t%18.12g\n", r);
192                  } else if (!(orient&(UP|DOWN))) {
193                          printf("\n%s ring %s.%d\n", modname, argv[2], i);
194                          printf("0\n0\n8\n");
195 <                        printf("\t0\t0\t%15.9g\n", z);
196 <                        printf("\t0\t0\t%15.9g\n", orient&IN ? 1.0 : -1.0);
197 <                        printf("\t%15.9g\t%15.9g\n", lastr, r);
195 >                        printf("\t0\t0\t%18.12g\n", z);
196 >                        printf("\t0\t0\t%18.12g\n", orient&IN ? 1.0 : -1.0);
197 >                        printf("\t%18.12g\t%18.12g\n", lastr, r);
198                  } else {
199                          printf("\n%s %s %s.%d\n", modname,
200                                          orient&DOWN ? "cup" : "cone",
201                                          argv[2], i);
202                          printf("0\n0\n8\n");
203 <                        printf("\t0\t0\t%15.9g\n", lastz);
204 <                        printf("\t0\t0\t%15.9g\n", z);
205 <                        printf("\t%15.9g\t%15.9g\n", lastr, r);
203 >                        printf("\t0\t0\t%18.12g\n", lastz);
204 >                        printf("\t0\t0\t%18.12g\n", z);
205 >                        printf("\t%18.12g\t%18.12g\n", lastr, r);
206                  }
207          endfor:
208                  lastz = z; lastr = r;
# Line 147 | Line 210 | char  *argv[];
210                  lastnz = nz; lastnr = nr;
211                  nz = nextnz; nr = nextnr;
212          }
213 <        quit(0);
213 >        return 0;
214  
215   userror:
216          fprintf(stderr,
217          "Usage: %s material name z(t) r(t) nseg [-e expr] [-f file] [-s]\n",
218                          argv[0]);
219 <        quit(1);
219 >        return 1;
220   }
221  
222  
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