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

Comparing ray/src/gen/genworm.c (file contents):
Revision 1.1 by greg, Thu Feb 2 11:16:31 1989 UTC vs.
Revision 1.5 by greg, Mon Jul 9 09:55:50 1990 UTC

# Line 1 | Line 1
1 – /*
2 –
1   #ifndef lint
2   static char SCCSid[] = "$SunId$ LBL";
3   #endif
4 +
5 + /* Copyright (c) 1989 Regents of the University of California */
6 +
7 + /*
8   *  genworm.c - program to generate worms (strings with varying thickness).
9   *
10   *      The program takes as input the functions of t for x, y,
# Line 28 | Line 30 | static char SCCSid[] = "$SunId$ LBL";
30   #define  max(a,b)       ((a) > (b) ? (a) : (b))
31  
32  
33 < double  funvalue(), l_hermite(), argument();
33 > double  funvalue(), l_hermite(), l_bezier(), l_bspline(), argument();
34  
35  
36   main(argc, argv)
37   int  argc;
38   char  *argv[];
39   {
40 +        extern long     eclock;
41          char  stmp[256];
42          double  t, f, lastr, r;
43          FVECT  lastp, p;
# Line 42 | Line 45 | char  *argv[];
45  
46          varset("PI", PI, NULL);
47          funset("hermite", 5, l_hermite);
48 +        funset("bezier", 5, l_bezier);
49 +        funset("bspline", 5, l_bspline);
50  
51          if (argc < 8)
52                  goto userror;
# Line 67 | Line 72 | char  *argv[];
72                  goto userror;
73  
74          printhead(argc, argv);
75 +        eclock = 0;
76  
77          for (i = 0; i <= nseg; i++) {
78                  t = (double)i/nseg;
# Line 160 | Line 166 | l_hermite()                    
166                  argument(2)*(-2.0*t+3.0)*t*t +
167                  argument(3)*((t-2.0)*t+1.0)*t +
168                  argument(4)*(t-1.0)*t*t );
169 + }
170 +
171 +
172 + double
173 + l_bezier()
174 + {
175 +        double  t;
176 +
177 +        t = argument(5);
178 +        return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
179 +                argument(2) * 3.*t*(1.+t*(-2.+t)) +
180 +                argument(3) * 3.*t*t*(1.-t) +
181 +                argument(4) * t*t*t );
182 + }
183 +
184 +
185 + double
186 + l_bspline()
187 + {
188 +        double  t;
189 +
190 +        t = argument(5);
191 +        return( argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
192 +                argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
193 +                argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
194 +                argument(4) * (1./6.*t*t*t) );
195   }

Diff Legend

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