--- ray/src/gen/genrev.c 1990/03/02 17:24:07 1.4 +++ ray/src/gen/genrev.c 1990/07/20 14:33:10 1.6 @@ -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(), l_bezier(), 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; @@ -46,6 +47,7 @@ char *argv[]; varset("PI", PI); funset("hermite", 5, l_hermite); funset("bezier", 5, l_bezier); + funset("bspline", 5, l_bspline); if (argc < 6) goto userror; @@ -70,6 +72,7 @@ char *argv[]; modname = smooth ? "Phong" : argv[1]; printhead(argc, argv); + eclock = 0; lastnz = lastnr = 0.0; t = 0.0; @@ -224,4 +227,17 @@ l_bezier() 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) ); }