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

Comparing src/gen/genworm.c (file contents):
Revision 1.1 by greg, Thu Feb 2 11:16:31 1989 UTC vs.
Revision 2.6 by greg, Mon Aug 4 22:37:53 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   *  genworm.c - program to generate worms (strings with varying thickness).
6   *
7   *      The program takes as input the functions of t for x, y,
# Line 13 | Line 12 | static char SCCSid[] = "$SunId$ LBL";
12   */
13  
14   #include  <stdio.h>
15 <
15 > #include  <math.h>
16 > #include  "calcomp.h"
17   #include  "fvect.h"
18  
19 < #define  XNAME          "X_"                    /* x function name */
20 < #define  YNAME          "Y_"                    /* y function name */
21 < #define  ZNAME          "Z_"                    /* z function name */
22 < #define  RNAME          "R_"                    /* r function name */
19 > #define  XNAME          "X`SYS`"                /* x function name */
20 > #define  YNAME          "Y`SYS`"                /* y function name */
21 > #define  ZNAME          "Z`SYS`"                /* z function name */
22 > #define  RNAME          "R`SYS`"                /* r function name */
23  
24   #define  PI             3.14159265358979323846
25  
26 – #define  FTINY          1e-7
27 –
26   #define  max(a,b)       ((a) > (b) ? (a) : (b))
27  
28  
29 < double  funvalue(), l_hermite(), argument();
29 > double  l_hermite(char *), l_bezier(char *), l_bspline(char *);
30 > void  quit();
31  
32  
33   main(argc, argv)
# Line 40 | Line 39 | char  *argv[];
39          FVECT  lastp, p;
40          int  i, nseg;
41  
42 <        varset("PI", PI, NULL);
43 <        funset("hermite", 5, l_hermite);
42 >        varset("PI", ':', PI);
43 >        funset("hermite", 5, ':', l_hermite);
44 >        funset("bezier", 5, ':', l_bezier);
45 >        funset("bspline", 5, ':', l_bspline);
46  
47          if (argc < 8)
48                  goto userror;
49  
50          for (i = 8; i < argc; i++)
51                  if (!strcmp(argv[i], "-e"))
52 <                        scompile(NULL, argv[++i]);
52 >                        scompile(argv[++i], NULL, 0);
53                  else if (!strcmp(argv[i], "-f"))
54                          fcompile(argv[++i]);
55                  else
56                          goto userror;
57  
58          sprintf(stmp, "%s(t)=%s;", XNAME, argv[3]);
59 <        scompile(NULL, stmp);
59 >        scompile(stmp, NULL, 0);
60          sprintf(stmp, "%s(t)=%s;", YNAME, argv[4]);
61 <        scompile(NULL, stmp);
61 >        scompile(stmp, NULL, 0);
62          sprintf(stmp, "%s(t)=%s;", ZNAME, argv[5]);
63 <        scompile(NULL, stmp);
63 >        scompile(stmp, NULL, 0);
64          sprintf(stmp, "%s(t)=%s;", RNAME, argv[6]);
65 <        scompile(NULL, stmp);
65 >        scompile(stmp, NULL, 0);
66          nseg = atoi(argv[7]);
67          if (nseg <= 0)
68                  goto userror;
69  
70          printhead(argc, argv);
71 +        eclock = 0;
72  
73          for (i = 0; i <= nseg; i++) {
74                  t = (double)i/nseg;
# Line 74 | Line 76 | char  *argv[];
76                  p[1] = funvalue(YNAME, 1, &t);
77                  p[2] = funvalue(ZNAME, 1, &t);
78                  r = funvalue(RNAME, 1, &t);
79 <                if (i)
79 >                if (i) {
80                          if (lastr <= r+FTINY && lastr >= r-FTINY) {
81                                  printf("\n%s cylinder %s.c%d\n",
82                                                  argv[1], argv[2], i);
# Line 101 | Line 103 | char  *argv[];
103                                  f = f <= 0.0 ? 0.0 : sqrt(f);
104                                  printf("%18.12g %18.12g\n", f*lastr, f*r);
105                          }
106 +                }
107                  printf("\n%s sphere %s.s%d\n", argv[1], argv[2], i);
108                  printf("0\n0\n4 %18.12g %18.12g %18.12g %18.12g\n",
109                                  p[0], p[1], p[2], r);
# Line 117 | Line 120 | userror:
120   }
121  
122  
123 + void
124   eputs(msg)
125   char  *msg;
126   {
# Line 124 | Line 128 | char  *msg;
128   }
129  
130  
131 + void
132   wputs(msg)
133   char  *msg;
134   {
# Line 131 | Line 136 | char  *msg;
136   }
137  
138  
139 + void
140   quit(code)
141 + int  code;
142   {
143          exit(code);
144   }
# Line 151 | Line 158 | register char  **av;
158  
159  
160   double
161 < l_hermite()                    
161 > l_hermite(char *nm)
162   {
163          double  t;
164          
# Line 160 | Line 167 | l_hermite()                    
167                  argument(2)*(-2.0*t+3.0)*t*t +
168                  argument(3)*((t-2.0)*t+1.0)*t +
169                  argument(4)*(t-1.0)*t*t );
170 + }
171 +
172 +
173 + double
174 + l_bezier(char *nm)
175 + {
176 +        double  t;
177 +
178 +        t = argument(5);
179 +        return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
180 +                argument(2) * 3.*t*(1.+t*(-2.+t)) +
181 +                argument(3) * 3.*t*t*(1.-t) +
182 +                argument(4) * t*t*t );
183 + }
184 +
185 +
186 + double
187 + l_bspline(char *nm)
188 + {
189 +        double  t;
190 +
191 +        t = argument(5);
192 +        return( argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
193 +                argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
194 +                argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
195 +                argument(4) * (1./6.*t*t*t) );
196   }

Diff Legend

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