--- ray/src/gen/genrev.c 1989/09/10 16:45:14 1.2 +++ ray/src/gen/genrev.c 1991/04/23 15:51:17 1.9 @@ -1,9 +1,9 @@ -/* Copyright (c) 1989 Regents of the University of California */ - #ifndef lint static char SCCSid[] = "$SunId$ LBL"; #endif +/* Copyright (c) 1989 Regents of the University of California */ + /* * genrev.c - program to generate functions of rotation about z * @@ -28,13 +28,14 @@ static char SCCSid[] = "$SunId$ LBL"; #define UP 04 #define DOWN 010 -double funvalue(), l_hermite(), argument(); +double funvalue(), l_hermite(), l_bezier(), l_bspline(), argument(); main(argc, argv) int argc; char *argv[]; { + extern long eclock; char stmp[256]; char *modname; int smooth = 0; @@ -43,15 +44,17 @@ char *argv[]; int i, nseg; int orient; - varset("PI", PI); - 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 < 6) goto userror; for (i = 6; 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 if (!strcmp(argv[i], "-s")) @@ -60,15 +63,16 @@ char *argv[]; goto userror; sprintf(stmp, "%s(t)=%s;", ZNAME, argv[3]); - scompile(NULL, stmp); + scompile(stmp, NULL, 0); sprintf(stmp, "%s(t)=%s;", RNAME, argv[4]); - scompile(NULL, stmp); + scompile(stmp, NULL, 0); nseg = atoi(argv[5]); if (nseg <= 0) goto userror; - modname = smooth ? "phong" : argv[1]; + modname = smooth ? "Phong" : argv[1]; printhead(argc, argv); + eclock = 0; lastnz = lastnr = 0.0; t = 0.0; @@ -98,7 +102,7 @@ char *argv[]; if (!orient) goto endfor; if (smooth) { - printf("\n%s texfunc phong\n", argv[1]); + printf("\n%s texfunc Phong\n", argv[1]); printf("4 rev_dx rev_dy rev_dz rev.cal\n"); printf("0\n4\n"); if (orient&(UP|DOWN)) { @@ -210,4 +214,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) ); }