1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
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
|
