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.4 by greg, Fri Mar 2 17:24:07 1990 UTC vs.
Revision 2.10 by greg, Sat Oct 11 04:39:34 2008 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 > #include  "rterror.h"
20 > #include  "resolu.h"
21 > #include  "calcomp.h"
22  
23 + #define  ZNAME          "Z`SYS`"                /* z function name */
24 + #define  RNAME          "R`SYS`"                /* r function name */
25 +
26   #define  PI             3.14159265358979323846
27  
28   #define  FTINY          1e-9
# Line 28 | Line 32 | static char SCCSid[] = "$SunId$ LBL";
32   #define  UP             04
33   #define  DOWN           010
34  
31 – double  funvalue(), l_hermite(), l_bezier(), argument();
35  
36 + void
37 + computen(nzp, nrp, z0, r0, z1, r1)              /* compute normal */
38 + double  *nzp, *nrp, z0, r0, z1, r1;
39 + {
40 +        double  dr, dz, len;
41  
42 +        dz = r0 - r1;                           /* right angle vector */
43 +        dr = z1 - z0;
44 +        len = sqrt(dr*dr + dz*dz);
45 +        *nzp = dz/len;
46 +        *nrp = dr/len;
47 + }
48 +
49 +
50 + double
51 + l_hermite(char *nm)
52 + {
53 +        double  t;
54 +        
55 +        t = argument(5);
56 +        return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
57 +                argument(2)*(-2.0*t+3.0)*t*t +
58 +                argument(3)*((t-2.0)*t+1.0)*t +
59 +                argument(4)*(t-1.0)*t*t );
60 + }
61 +
62 +
63 + double
64 + l_bezier(char *nm)
65 + {
66 +        double  t;
67 +
68 +        t = argument(5);
69 +        return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
70 +                argument(2) * 3.*t*(1.+t*(-2.+t)) +
71 +                argument(3) * 3.*t*t*(1.-t) +
72 +                argument(4) * t*t*t );
73 + }
74 +
75 +
76 + double
77 + l_bspline(char *nm)
78 + {
79 +        double  t;
80 +
81 +        t = argument(5);
82 +        return( argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
83 +                argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
84 +                argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
85 +                argument(4) * (1./6.*t*t*t) );
86 + }
87 +
88 +
89 + int
90   main(argc, argv)
91   int  argc;
92   char  *argv[];
# Line 43 | Line 99 | char  *argv[];
99          int  i, nseg;
100          int  orient;
101  
102 <        varset("PI", PI);
103 <        funset("hermite", 5, l_hermite);
104 <        funset("bezier", 5, l_bezier);
102 >        varset("PI", ':', PI);
103 >        funset("hermite", 5, ':', l_hermite);
104 >        funset("bezier", 5, ':', l_bezier);
105 >        funset("bspline", 5, ':', l_bspline);
106  
107          if (argc < 6)
108                  goto userror;
109  
110          for (i = 6; i < argc; i++)
111                  if (!strcmp(argv[i], "-e"))
112 <                        scompile(NULL, argv[++i]);
112 >                        scompile(argv[++i], NULL, 0);
113                  else if (!strcmp(argv[i], "-f"))
114                          fcompile(argv[++i]);
115                  else if (!strcmp(argv[i], "-s"))
# Line 61 | Line 118 | char  *argv[];
118                          goto userror;
119  
120          sprintf(stmp, "%s(t)=%s;", ZNAME, argv[3]);
121 <        scompile(NULL, stmp);
121 >        scompile(stmp, NULL, 0);
122          sprintf(stmp, "%s(t)=%s;", RNAME, argv[4]);
123 <        scompile(NULL, stmp);
123 >        scompile(stmp, NULL, 0);
124          nseg = atoi(argv[5]);
125          if (nseg <= 0)
126                  goto userror;
127          modname = smooth ? "Phong" : argv[1];
128  
129 <        printhead(argc, argv);
129 >        fputs("# ", stdout);
130 >        printargs(argc, argv, stdout);
131 >        eclock = 0;
132  
133          lastnz = lastnr = 0.0;
134          t = 0.0;
# Line 104 | Line 163 | char  *argv[];
163                          printf("0\n4\n");
164                          if (orient&(UP|DOWN)) {
165                                  t = (nextnz - lastnz)/(z - lastz);
166 <                                printf("\t%15.9g\t%15.9g\n",
166 >                                printf("\t%18.12g\t%18.12g\n",
167                                                  t, lastnz - t*lastz);
168                          } else
169                                  printf("\t0\t%d\n", orient&IN ? 1 : -1);
170                          if (orient&(OUT|IN))  {
171                                  t = (nextnr - lastnr)/(r - lastr);
172 <                                printf("\t%15.9g\t%15.9g\n",
172 >                                printf("\t%18.12g\t%18.12g\n",
173                                                  t, lastnr - t*lastr);
174                          } else
175                                  printf("\t0\t%d\n", orient&UP ? 1 : -1);
# Line 120 | Line 179 | char  *argv[];
179                                          orient&DOWN ? "tube" : "cylinder",
180                                          argv[2], i);
181                          printf("0\n0\n7\n");
182 <                        printf("\t0\t0\t%15.9g\n", lastz);
183 <                        printf("\t0\t0\t%15.9g\n", z);
184 <                        printf("\t%15.9g\n", r);
182 >                        printf("\t0\t0\t%18.12g\n", lastz);
183 >                        printf("\t0\t0\t%18.12g\n", z);
184 >                        printf("\t%18.12g\n", r);
185                  } else if (!(orient&(UP|DOWN))) {
186                          printf("\n%s ring %s.%d\n", modname, argv[2], i);
187                          printf("0\n0\n8\n");
188 <                        printf("\t0\t0\t%15.9g\n", z);
189 <                        printf("\t0\t0\t%15.9g\n", orient&IN ? 1.0 : -1.0);
190 <                        printf("\t%15.9g\t%15.9g\n", lastr, r);
188 >                        printf("\t0\t0\t%18.12g\n", z);
189 >                        printf("\t0\t0\t%18.12g\n", orient&IN ? 1.0 : -1.0);
190 >                        printf("\t%18.12g\t%18.12g\n", lastr, r);
191                  } else {
192                          printf("\n%s %s %s.%d\n", modname,
193                                          orient&DOWN ? "cup" : "cone",
194                                          argv[2], i);
195                          printf("0\n0\n8\n");
196 <                        printf("\t0\t0\t%15.9g\n", lastz);
197 <                        printf("\t0\t0\t%15.9g\n", z);
198 <                        printf("\t%15.9g\t%15.9g\n", lastr, r);
196 >                        printf("\t0\t0\t%18.12g\n", lastz);
197 >                        printf("\t0\t0\t%18.12g\n", z);
198 >                        printf("\t%18.12g\t%18.12g\n", lastr, r);
199                  }
200          endfor:
201                  lastz = z; lastr = r;
# Line 144 | Line 203 | char  *argv[];
203                  lastnz = nz; lastnr = nr;
204                  nz = nextnz; nr = nextnr;
205          }
206 <        quit(0);
206 >        return 0;
207  
208   userror:
209          fprintf(stderr,
210          "Usage: %s material name z(t) r(t) nseg [-e expr] [-f file] [-s]\n",
211                          argv[0]);
212 <        quit(1);
212 >        return 1;
213   }
214  
215  
157 – computen(nzp, nrp, z0, r0, z1, r1)              /* compute normal */
158 – double  *nzp, *nrp, z0, r0, z1, r1;
159 – {
160 –        extern double  sqrt();
161 –        double  dr, dz, len;
162 –
163 –        dz = r0 - r1;                           /* right angle vector */
164 –        dr = z1 - z0;
165 –        len = sqrt(dr*dr + dz*dz);
166 –        *nzp = dz/len;
167 –        *nrp = dr/len;
168 – }
169 –
170 –
171 – eputs(msg)
172 – char  *msg;
173 – {
174 –        fputs(msg, stderr);
175 – }
176 –
177 –
178 – wputs(msg)
179 – char  *msg;
180 – {
181 –        eputs(msg);
182 – }
183 –
184 –
185 – quit(code)
186 – {
187 –        exit(code);
188 – }
189 –
190 –
191 – printhead(ac, av)               /* print command header */
192 – register int  ac;
193 – register char  **av;
194 – {
195 –        putchar('#');
196 –        while (ac--) {
197 –                putchar(' ');
198 –                fputs(*av++, stdout);
199 –        }
200 –        putchar('\n');
201 – }
202 –
203 –
204 – double
205 – l_hermite()                    
206 – {
207 –        double  t;
208 –        
209 –        t = argument(5);
210 –        return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
211 –                argument(2)*(-2.0*t+3.0)*t*t +
212 –                argument(3)*((t-2.0)*t+1.0)*t +
213 –                argument(4)*(t-1.0)*t*t );
214 – }
215 –
216 –
217 – double
218 – l_bezier()
219 – {
220 –        double  t;
221 –
222 –        t = argument(5);
223 –        return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
224 –                argument(2) * 3.*t*(1.+t*(-2.+t)) +
225 –                argument(3) * 3.*t*t*(1.-t) +
226 –                argument(4) * t*t*t );
227 – }

Diff Legend

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