--- ray/src/gen/genworm.c	1989/09/10 16:45:18	1.2
+++ ray/src/gen/genworm.c	2003/07/21 22:30:18	2.5
@@ -1,9 +1,6 @@
-/* Copyright (c) 1989 Regents of the University of California */
-
 #ifndef lint
-static char SCCSid[] = "$SunId$ LBL";
+static const char	RCSid[] = "$Id: genworm.c,v 2.5 2003/07/21 22:30:18 schorsch Exp $";
 #endif
-
 /*
  *  genworm.c - program to generate worms (strings with varying thickness).
  *
@@ -15,60 +12,63 @@ static char SCCSid[] = "$SunId$ LBL";
  */
 
 #include  <stdio.h>
-
+#include  <math.h>
 #include  "fvect.h"
 
-#define  XNAME		"X_"			/* x function name */
-#define  YNAME		"Y_"			/* y function name */
-#define  ZNAME		"Z_"			/* z function name */
-#define  RNAME		"R_"			/* r function name */
+#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(), argument();
+double  funvalue(), l_hermite(), l_bezier(), l_bspline(), argument();
+void  quit();
 
 
 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, NULL);
-	funset("hermite", 5, l_hermite);
+	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(NULL, argv[++i]);
+			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(NULL, stmp);
+	scompile(stmp, NULL, 0);
 	sprintf(stmp, "%s(t)=%s;", YNAME, argv[4]);
-	scompile(NULL, stmp);
+	scompile(stmp, NULL, 0);
 	sprintf(stmp, "%s(t)=%s;", ZNAME, argv[5]);
-	scompile(NULL, stmp);
+	scompile(stmp, NULL, 0);
 	sprintf(stmp, "%s(t)=%s;", RNAME, argv[6]);
-	scompile(NULL, stmp);
+	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;
@@ -76,7 +76,7 @@ char  *argv[];
 		p[1] = funvalue(YNAME, 1, &t);
 		p[2] = funvalue(ZNAME, 1, &t);
 		r = funvalue(RNAME, 1, &t);
-		if (i)
+		if (i) {
 			if (lastr <= r+FTINY && lastr >= r-FTINY) {
 				printf("\n%s cylinder %s.c%d\n",
 						argv[1], argv[2], i);
@@ -103,6 +103,7 @@ char  *argv[];
 				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);
@@ -119,6 +120,7 @@ userror:
 }
 
 
+void
 eputs(msg)
 char  *msg;
 {
@@ -126,6 +128,7 @@ char  *msg;
 }
 
 
+void
 wputs(msg)
 char  *msg;
 {
@@ -133,7 +136,9 @@ char  *msg;
 }
 
 
+void
 quit(code)
+int  code;
 {
 	exit(code);
 }
@@ -162,4 +167,30 @@ l_hermite()			
 		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) );
 }