src/gen/genworm.c

#ifndef lint
static char SCCSid[] = "$SunId$ LBL";
#endif

/* Copyright (c) 1989 Regents of the University of California */

/*
 *  genworm.c - program to generate worms (strings with varying thickness).
 *
 *      The program takes as input the functions of t for x, y,
 *      z and r (the radius).  Cylinders and cones will be
 *      joined by spheres.  Negative radii are forbidden.
 *
 *      9/24/87
 */

#include  <stdio.h>

#include  "fvect.h"

#define  XNAME          "X`SYS`"                /* x function name */
#define  YNAME          "Y`SYS`"                /* y function name */
#define  ZNAME          "Z`SYS`"                /* z function name */
#define  RNAME          "R`SYS`"                /* r function name */

#define  PI             3.14159265358979323846

#define  FTINY          1e-7

#define  max(a,b)       ((a) > (b) ? (a) : (b))


double  funvalue(), l_hermite(), l_bezier(), l_bspline(), argument();


main(argc, argv)
int  argc;
char  *argv[];
{
        extern long     eclock;
        char  stmp[256];
        double  t, f, lastr, r;
        FVECT  lastp, p;
        int  i, nseg;

        varset("PI", ':', PI);
        funset("hermite", 5, ':', l_hermite);
        funset("bezier", 5, ':', l_bezier);
        funset("bspline", 5, ':', l_bspline);

        if (argc < 8)
                goto userror;

        for (i = 8; i < argc; i++)
                if (!strcmp(argv[i], "-e"))
                        scompile(argv[++i], NULL, 0);
                else if (!strcmp(argv[i], "-f"))
                        fcompile(argv[++i]);
                else
                        goto userror;

        sprintf(stmp, "%s(t)=%s;", XNAME, argv[3]);
        scompile(stmp, NULL, 0);
        sprintf(stmp, "%s(t)=%s;", YNAME, argv[4]);
        scompile(stmp, NULL, 0);
        sprintf(stmp, "%s(t)=%s;", ZNAME, argv[5]);
        scompile(stmp, NULL, 0);
        sprintf(stmp, "%s(t)=%s;", RNAME, argv[6]);
        scompile(stmp, NULL, 0);
        nseg = atoi(argv[7]);
        if (nseg <= 0)
                goto userror;

        printhead(argc, argv);
        eclock = 0;

        for (i = 0; i <= nseg; i++) {
                t = (double)i/nseg;
                p[0] = funvalue(XNAME, 1, &t);
                p[1] = funvalue(YNAME, 1, &t);
                p[2] = funvalue(ZNAME, 1, &t);
                r = funvalue(RNAME, 1, &t);
                if (i)
                        if (lastr <= r+FTINY && lastr >= r-FTINY) {
                                printf("\n%s cylinder %s.c%d\n",
                                                argv[1], argv[2], i);
                                printf("0\n0\n7\n");
                                printf("%18.12g %18.12g %18.12g\n",
                                                lastp[0], lastp[1], lastp[2]);
                                printf("%18.12g %18.12g %18.12g\n",
                                                p[0], p[1], p[2]);
                                printf("%18.12g\n", r);
                        } else {
                                printf("\n%s cone %s.c%d\n",
                                                argv[1], argv[2], i);
                                printf("0\n0\n8\n");
                                f = (lastr - r)/dist2(lastp,p);
                                printf("%18.12g %18.12g %18.12g\n",
                                        lastp[0] + f*lastr*(p[0] - lastp[0]),
                                        lastp[1] + f*lastr*(p[1] - lastp[1]),
                                        lastp[2] + f*lastr*(p[2] - lastp[2]));
                                printf("%18.12g %18.12g %18.12g\n",
                                        p[0] + f*r*(p[0] - lastp[0]),
                                        p[1] + f*r*(p[1] - lastp[1]),
                                        p[2] + f*r*(p[2] - lastp[2]));
                                f = 1.0 - (lastr-r)*f;
                                f = f <= 0.0 ? 0.0 : sqrt(f);
                                printf("%18.12g %18.12g\n", f*lastr, f*r);
                        }
                printf("\n%s sphere %s.s%d\n", argv[1], argv[2], i);
                printf("0\n0\n4 %18.12g %18.12g %18.12g %18.12g\n",
                                p[0], p[1], p[2], r);
                VCOPY(lastp, p);
                lastr = r;
        }
        quit(0);

userror:
        fprintf(stderr,
"Usage: %s material name x(t) y(t) z(t) r(t) nseg [-e expr] [-f file]\n",
                        argv[0]);
        quit(1);
}


eputs(msg)
char  *msg;
{
        fputs(msg, stderr);
}


wputs(msg)
char  *msg;
{
        eputs(msg);
}


quit(code)
{
        exit(code);
}


printhead(ac, av)               /* print command header */
register int  ac;
register char  **av;
{
        putchar('#');
        while (ac--) {
                putchar(' ');
                fputs(*av++, stdout);
        }
        putchar('\n');
}


double
l_hermite()                     
{
        double  t;
        
        t = argument(5);
        return( argument(1)*((2.0*t-3.0)*t*t+1.0) +
                argument(2)*(-2.0*t+3.0)*t*t +
                argument(3)*((t-2.0)*t+1.0)*t +
                argument(4)*(t-1.0)*t*t );
}


double
l_bezier()
{
        double  t;

        t = argument(5);
        return( argument(1) * (1.+t*(-3.+t*(3.-t))) +
                argument(2) * 3.*t*(1.+t*(-2.+t)) +
                argument(3) * 3.*t*t*(1.-t) +
                argument(4) * t*t*t );
}


double
l_bspline()
{
        double  t;

        t = argument(5);
        return( argument(1) * (1./6.+t*(-1./2.+t*(1./2.-1./6.*t))) +
                argument(2) * (2./3.+t*t*(-1.+1./2.*t)) +
                argument(3) * (1./6.+t*(1./2.+t*(1./2.-1./2.*t))) +
                argument(4) * (1./6.*t*t*t) );
}
Revision:	1.9
Committed:	Tue Aug 13 16:07:12 1991 UTC (32 years, 8 months ago) by greg
Content type:	text/plain
Branch:	MAIN
Changes since 1.8:	+4 -4 lines
Log Message:	changed private variable names with `SYS` context
#	User	Rev	Content
1	greg	1.1	#ifndef lint
2			static char SCCSid[] = "$SunId$ LBL";
3			#endif
4	greg	1.2
5	greg	1.4	/* Copyright (c) 1989 Regents of the University of California */
6
7	greg	1.2	/*
8	greg	1.1	* genworm.c - program to generate worms (strings with varying thickness).
9			*
10			* The program takes as input the functions of t for x, y,
11			* z and r (the radius). Cylinders and cones will be
12			* joined by spheres. Negative radii are forbidden.
13			*
14			* 9/24/87
15			*/
16
17			#include <stdio.h>
18
19			#include "fvect.h"
20
21	greg	1.9	#define XNAME "X`SYS`" /* x function name */
22			#define YNAME "Y`SYS`" /* y function name */
23			#define ZNAME "Z`SYS`" /* z function name */
24			#define RNAME "R`SYS`" /* r function name */
25	greg	1.1
26			#define PI 3.14159265358979323846
27
28			#define FTINY 1e-7
29
30			#define max(a,b) ((a) > (b) ? (a) : (b))
31
32
33	greg	1.4	double funvalue(), l_hermite(), l_bezier(), l_bspline(), argument();
34	greg	1.1
35
36			main(argc, argv)
37			int argc;
38			char *argv[];
39			{
40	greg	1.5	extern long eclock;
41	greg	1.1	char stmp[256];
42			double t, f, lastr, r;
43			FVECT lastp, p;
44			int i, nseg;
45
46	greg	1.7	varset("PI", ':', PI);
47	greg	1.8	funset("hermite", 5, ':', l_hermite);
48			funset("bezier", 5, ':', l_bezier);
49			funset("bspline", 5, ':', l_bspline);
50	greg	1.1
51			if (argc < 8)
52			goto userror;
53
54			for (i = 8; i < argc; i++)
55			if (!strcmp(argv[i], "-e"))
56	greg	1.6	scompile(argv[++i], NULL, 0);
57	greg	1.1	else if (!strcmp(argv[i], "-f"))
58			fcompile(argv[++i]);
59			else
60			goto userror;
61
62			sprintf(stmp, "%s(t)=%s;", XNAME, argv[3]);
63	greg	1.6	scompile(stmp, NULL, 0);
64	greg	1.1	sprintf(stmp, "%s(t)=%s;", YNAME, argv[4]);
65	greg	1.6	scompile(stmp, NULL, 0);
66	greg	1.1	sprintf(stmp, "%s(t)=%s;", ZNAME, argv[5]);
67	greg	1.6	scompile(stmp, NULL, 0);
68	greg	1.1	sprintf(stmp, "%s(t)=%s;", RNAME, argv[6]);
69	greg	1.6	scompile(stmp, NULL, 0);
70	greg	1.1	nseg = atoi(argv[7]);
71			if (nseg <= 0)
72			goto userror;
73
74			printhead(argc, argv);
75	greg	1.5	eclock = 0;
76	greg	1.1
77			for (i = 0; i <= nseg; i++) {
78			t = (double)i/nseg;
79			p[0] = funvalue(XNAME, 1, &t);
80			p[1] = funvalue(YNAME, 1, &t);
81			p[2] = funvalue(ZNAME, 1, &t);
82			r = funvalue(RNAME, 1, &t);
83			if (i)
84			if (lastr <= r+FTINY && lastr >= r-FTINY) {
85			printf("\n%s cylinder %s.c%d\n",
86			argv[1], argv[2], i);
87			printf("0\n0\n7\n");
88			printf("%18.12g %18.12g %18.12g\n",
89			lastp[0], lastp[1], lastp[2]);
90			printf("%18.12g %18.12g %18.12g\n",
91			p[0], p[1], p[2]);
92			printf("%18.12g\n", r);
93			} else {
94			printf("\n%s cone %s.c%d\n",
95			argv[1], argv[2], i);
96			printf("0\n0\n8\n");
97			f = (lastr - r)/dist2(lastp,p);
98			printf("%18.12g %18.12g %18.12g\n",
99			lastp[0] + flastr(p[0] - lastp[0]),
100			lastp[1] + flastr(p[1] - lastp[1]),
101			lastp[2] + flastr(p[2] - lastp[2]));
102			printf("%18.12g %18.12g %18.12g\n",
103			p[0] + fr(p[0] - lastp[0]),
104			p[1] + fr(p[1] - lastp[1]),
105			p[2] + fr(p[2] - lastp[2]));
106			f = 1.0 - (lastr-r)*f;
107			f = f <= 0.0 ? 0.0 : sqrt(f);
108			printf("%18.12g %18.12g\n", flastr, fr);
109			}
110			printf("\n%s sphere %s.s%d\n", argv[1], argv[2], i);
111			printf("0\n0\n4 %18.12g %18.12g %18.12g %18.12g\n",
112			p[0], p[1], p[2], r);
113			VCOPY(lastp, p);
114			lastr = r;
115			}
116			quit(0);
117
118			userror:
119			fprintf(stderr,
120			"Usage: %s material name x(t) y(t) z(t) r(t) nseg [-e expr] [-f file]\n",
121			argv[0]);
122			quit(1);
123			}
124
125
126			eputs(msg)
127			char *msg;
128			{
129			fputs(msg, stderr);
130			}
131
132
133			wputs(msg)
134			char *msg;
135			{
136			eputs(msg);
137			}
138
139
140			quit(code)
141			{
142			exit(code);
143			}
144
145
146			printhead(ac, av) /* print command header */
147			register int ac;
148			register char **av;
149			{
150			putchar('#');
151			while (ac--) {
152			putchar(' ');
153			fputs(*av++, stdout);
154			}
155			putchar('\n');
156			}
157
158
159			double
160			l_hermite()
161			{
162			double t;
163
164			t = argument(5);
165			return( argument(1)((2.0t-3.0)tt+1.0) +
166			argument(2)(-2.0t+3.0)tt +
167			argument(3)((t-2.0)t+1.0)*t +
168			argument(4)(t-1.0)t*t );
169	greg	1.3	}
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.tt*(1.-t) +
181			argument(4) * ttt );
182	greg	1.4	}
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.+tt(-1.+1./2.*t)) +
193			argument(3) * (1./6.+t(1./2.+t(1./2.-1./2.*t))) +
194			argument(4) * (1./6.tt*t) );
195	greg	1.1	}