diff options
Diffstat (limited to 'math/deboor/setdat2.f')
| -rw-r--r-- | math/deboor/setdat2.f | 29 |
1 files changed, 29 insertions, 0 deletions
diff --git a/math/deboor/setdat2.f b/math/deboor/setdat2.f new file mode 100644 index 00000000..cf767a59 --- /dev/null +++ b/math/deboor/setdat2.f @@ -0,0 +1,29 @@ + subroutine setdt2(icount) +c from * a practical guide to splines * by c. de boor +c to be called in main program l 2 m a i n . +c this routine is set up to provide the specific data for example 3 +c in chapter xiv. +c + integer icount, i,k,l,ntau + real break,coef,gtau,step,tau,totalw,weight,x,round +c parameter lpkmax=100,ntmax=200,ltkmax=2000 +c common / data / ntau, tau(ntmax),gtau(ntmax),weight(ntmax),totalw +c common /approx/ break(lpkmax),coef(ltkmax),l,k + common / data / ntau, tau(200),gtau(200),weight(200),totalw + common /approx/ break(100),coef(2000),l,k + round(x) = float(ifix(x*100.))/100. + if (icount .gt. 0) stop + icount = icount + 1 + ntau = 65 + step = 3./float(ntau-1) + do 10 i=1,ntau + tau(i) = i*step + gtau(i) = round(exp(tau(i))) + 10 weight(i) = 1. + totalw = ntau + l = 1 + break(1) = tau(1) + break(2) = tau(ntau) + k = 3 + return + end |