ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/Development/ray/src/gen/genworm.c
(Generate patch)

Comparing ray/src/gen/genworm.c (file contents):
Revision 1.2 by greg, Sun Sep 10 16:45:18 1989 UTC vs.
Revision 2.7 by schorsch, Sun Nov 16 10:29:38 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   *  genworm.c - program to generate worms (strings with varying thickness).
6   *
# Line 14 | Line 11 | static char SCCSid[] = "$SunId$ LBL";
11   *      9/24/87
12   */
13  
14 + #include  <stdlib.h>
15   #include  <stdio.h>
16 + #include  <math.h>
17 + #include  <string.h>
18  
19 + #include  "calcomp.h"
20   #include  "fvect.h"
21  
22 < #define  XNAME          "X_"                    /* x function name */
23 < #define  YNAME          "Y_"                    /* y function name */
24 < #define  ZNAME          "Z_"                    /* z function name */
25 < #define  RNAME          "R_"                    /* r function name */
22 > #define  XNAME          "X`SYS`"                /* x function name */
23 > #define  YNAME          "Y`SYS`"                /* y function name */
24 > #define  ZNAME          "Z`SYS`"                /* z function name */
25 > #define  RNAME          "R`SYS`"                /* r function name */
26  
27   #define  PI             3.14159265358979323846
28  
28 – #define  FTINY          1e-7
29 –
29   #define  max(a,b)       ((a) > (b) ? (a) : (b))
30  
31  
32 < double  funvalue(), l_hermite(), argument();
32 > /* XXX redundant, move to library */
33 > double  l_hermite(char *), l_bezier(char *), l_bspline(char *);
34  
35  
36 + static void printhead(int  ac, char  **av);
37 +
38 +
39 + int
40   main(argc, argv)
41   int  argc;
42   char  *argv[];
43   {
44          char  stmp[256];
45 <        double  t, f, lastr, r;
45 >        double  t, f, lastr = 0, r;
46          FVECT  lastp, p;
47          int  i, nseg;
48  
49 <        varset("PI", PI, NULL);
50 <        funset("hermite", 5, l_hermite);
49 >        varset("PI", ':', PI);
50 >        funset("hermite", 5, ':', l_hermite);
51 >        funset("bezier", 5, ':', l_bezier);
52 >        funset("bspline", 5, ':', l_bspline);
53  
54          if (argc < 8)
55                  goto userror;
56  
57          for (i = 8; i < argc; i++)
58                  if (!strcmp(argv[i], "-e"))
59 <                        scompile(NULL, argv[++i]);
59 >                        scompile(argv[++i], NULL, 0);
60                  else if (!strcmp(argv[i], "-f"))
61                          fcompile(argv[++i]);
62                  else
63                          goto userror;
64  
65          sprintf(stmp, "%s(t)=%s;", XNAME, argv[3]);
66 <        scompile(NULL, stmp);
66 >        scompile(stmp, NULL, 0);
67          sprintf(stmp, "%s(t)=%s;", YNAME, argv[4]);
68 <        scompile(NULL, stmp);
68 >        scompile(stmp, NULL, 0);
69          sprintf(stmp, "%s(t)=%s;", ZNAME, argv[5]);
70 <        scompile(NULL, stmp);
70 >        scompile(stmp, NULL, 0);
71          sprintf(stmp, "%s(t)=%s;", RNAME, argv[6]);
72 <        scompile(NULL, stmp);
72 >        scompile(stmp, NULL, 0);
73          nseg = atoi(argv[7]);
74          if (nseg <= 0)
75                  goto userror;
76  
77          printhead(argc, argv);
78 +        eclock = 0;
79  
80          for (i = 0; i <= nseg; i++) {
81                  t = (double)i/nseg;
# Line 76 | Line 83 | char  *argv[];
83                  p[1] = funvalue(YNAME, 1, &t);
84                  p[2] = funvalue(ZNAME, 1, &t);
85                  r = funvalue(RNAME, 1, &t);
86 <                if (i)
86 >                if (i) {
87                          if (lastr <= r+FTINY && lastr >= r-FTINY) {
88                                  printf("\n%s cylinder %s.c%d\n",
89                                                  argv[1], argv[2], i);
# Line 103 | Line 110 | char  *argv[];
110                                  f = f <= 0.0 ? 0.0 : sqrt(f);
111                                  printf("%18.12g %18.12g\n", f*lastr, f*r);
112                          }
113 +                }
114                  printf("\n%s sphere %s.s%d\n", argv[1], argv[2], i);
115                  printf("0\n0\n4 %18.12g %18.12g %18.12g %18.12g\n",
116                                  p[0], p[1], p[2], r);
117                  VCOPY(lastp, p);
118                  lastr = r;
119          }
120 <        quit(0);
120 >        return 0;
121  
122   userror:
123          fprintf(stderr,
124   "Usage: %s material name x(t) y(t) z(t) r(t) nseg [-e expr] [-f file]\n",
125                          argv[0]);
126 <        quit(1);
126 >        return 1;
127   }
128  
129  
130 + void
131   eputs(msg)
132   char  *msg;
133   {
# Line 126 | Line 135 | char  *msg;
135   }
136  
137  
138 + void
139   wputs(msg)
140   char  *msg;
141   {
# Line 133 | Line 143 | char  *msg;
143   }
144  
145  
146 + void
147   quit(code)
148 + int  code;
149   {
150          exit(code);
151   }
152  
153  
154 < printhead(ac, av)               /* print command header */
155 < register int  ac;
156 < register char  **av;
154 > void
155 > printhead(              /* print command header */
156 >        register int  ac,
157 >        register char  **av
158 > )
159   {
160          putchar('#');
161          while (ac--) {
# Line 153 | Line 167 | register char  **av;
167  
168  
169   double
170 < l_hermite()                    
170 > l_hermite(char *nm)
171   {
172          double  t;
173          
# Line 162 | Line 176 | l_hermite()                    
176                  argument(2)*(-2.0*t+3.0)*t*t +
177                  argument(3)*((t-2.0)*t+1.0)*t +
178                  argument(4)*(t-1.0)*t*t );
179 + }
180 +
181 +
182 + double
183 + l_bezier(char *nm)
184 + {
185 +        double  t;
186 +
187 +        t = argument(5);
188 +        return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
189 +                argument(2) * 3.*t*(1.+t*(-2.+t)) +
190 +                argument(3) * 3.*t*t*(1.-t) +
191 +                argument(4) * t*t*t );
192 + }
193 +
194 +
195 + double
196 + l_bspline(char *nm)
197 + {
198 +        double  t;
199 +
200 +        t = argument(5);
201 +        return( argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
202 +                argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
203 +                argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
204 +                argument(4) * (1./6.*t*t*t) );
205   }

Diff Legend

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