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 2.2 by greg, Tue Jan 5 10:06:59 1993 UTC vs.
Revision 2.13 by greg, Tue Dec 10 17:22:39 2024 UTC

# Line 1 | Line 1
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ LBL";
2 > static const char       RCSid[] = "$Id$";
3   #endif
4 –
5 – /* Copyright (c) 1989 Regents of the University of California */
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  <stdio.h>
14 > #include  <stdlib.h>
15   #include  <math.h>
16 +
17 + #include  "calcomp.h"
18 + #include  "rtio.h"
19 + #include  "resolu.h"
20 + #include  "rterror.h"
21   #include  "fvect.h"
22  
23   #define  XNAME          "X`SYS`"                /* x function name */
# 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(), l_bezier(), l_bspline(), argument();
33 > /* XXX redundant, move to library */
34 > double  l_hermite(char *), l_bezier(char *), l_bspline(char *);
35  
36  
37 + int
38   main(argc, argv)
39   int  argc;
40   char  *argv[];
41   {
38 –        extern long     eclock;
42          char  stmp[256];
43 <        double  t, f, lastr, r;
43 >        double  t, f, lastr = 0, r;
44          FVECT  lastp, p;
45          int  i, nseg;
46  
47 +        esupport |= E_VARIABLE|E_FUNCTION|E_RCONST;
48 +        esupport &= ~(E_OUTCHAN|E_INCHAN);
49          varset("PI", ':', PI);
50          funset("hermite", 5, ':', l_hermite);
51          funset("bezier", 5, ':', l_bezier);
# Line 52 | Line 57 | char  *argv[];
57          for (i = 8; i < argc; i++)
58                  if (!strcmp(argv[i], "-e"))
59                          scompile(argv[++i], NULL, 0);
60 <                else if (!strcmp(argv[i], "-f"))
61 <                        fcompile(argv[++i]);
62 <                else
60 >                else if (!strcmp(argv[i], "-f")) {
61 >                        char  *fpath = getpath(argv[++i], getrlibpath(), 0);
62 >                        if (fpath == NULL) {
63 >                                fprintf(stderr, "%s: cannot find file '%s'\n",
64 >                                                argv[0], argv[i]);
65 >                                quit(1);
66 >                        }
67 >                        fcompile(fpath);
68 >                } else
69                          goto userror;
70  
71          sprintf(stmp, "%s(t)=%s;", XNAME, argv[3]);
# Line 65 | Line 76 | char  *argv[];
76          scompile(stmp, NULL, 0);
77          sprintf(stmp, "%s(t)=%s;", RNAME, argv[6]);
78          scompile(stmp, NULL, 0);
79 <        nseg = atoi(argv[7]);
79 >        nseg = eval(argv[7]) + .5;
80          if (nseg <= 0)
81                  goto userror;
82  
83 <        printhead(argc, argv);
84 <        eclock = 0;
83 >        fputs("# ", stdout);
84 >        printargs(argc, argv, stdout);
85 >        eclock = 1;
86  
87          for (i = 0; i <= nseg; i++) {
88                  t = (double)i/nseg;
# Line 78 | Line 90 | char  *argv[];
90                  p[1] = funvalue(YNAME, 1, &t);
91                  p[2] = funvalue(ZNAME, 1, &t);
92                  r = funvalue(RNAME, 1, &t);
93 <                if (i)
93 >                if (i) {
94                          if (lastr <= r+FTINY && lastr >= r-FTINY) {
95                                  printf("\n%s cylinder %s.c%d\n",
96                                                  argv[1], argv[2], i);
# Line 105 | Line 117 | char  *argv[];
117                                  f = f <= 0.0 ? 0.0 : sqrt(f);
118                                  printf("%18.12g %18.12g\n", f*lastr, f*r);
119                          }
120 +                }
121                  printf("\n%s sphere %s.s%d\n", argv[1], argv[2], i);
122                  printf("0\n0\n4 %18.12g %18.12g %18.12g %18.12g\n",
123                                  p[0], p[1], p[2], r);
124                  VCOPY(lastp, p);
125                  lastr = r;
126          }
127 <        quit(0);
127 >        return 0;
128  
129   userror:
130          fprintf(stderr,
131   "Usage: %s material name x(t) y(t) z(t) r(t) nseg [-e expr] [-f file]\n",
132                          argv[0]);
133 <        quit(1);
133 >        return 1;
134   }
135  
136  
124 – eputs(msg)
125 – char  *msg;
126 – {
127 –        fputs(msg, stderr);
128 – }
129 –
130 –
131 – wputs(msg)
132 – char  *msg;
133 – {
134 –        eputs(msg);
135 – }
136 –
137 –
138 – quit(code)
139 – {
140 –        exit(code);
141 – }
142 –
143 –
144 – printhead(ac, av)               /* print command header */
145 – register int  ac;
146 – register char  **av;
147 – {
148 –        putchar('#');
149 –        while (ac--) {
150 –                putchar(' ');
151 –                fputs(*av++, stdout);
152 –        }
153 –        putchar('\n');
154 – }
155 –
156 –
137   double
138 < l_hermite()                    
138 > l_hermite(char *nm)
139   {
140          double  t;
141          
# Line 168 | Line 148 | l_hermite()                    
148  
149  
150   double
151 < l_bezier()
151 > l_bezier(char *nm)
152   {
153          double  t;
154  
# Line 181 | Line 161 | l_bezier()
161  
162  
163   double
164 < l_bspline()
164 > l_bspline(char *nm)
165   {
166          double  t;
167  

Diff Legend

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