1 SUBROUTINE initspline(amat, splnot, h, weight, nots)
8 REAL(rprec),
DIMENSION(nots,nots) :: amat
9 REAL(rprec),
DIMENSION(*) :: splnot, h, weight
17 IF (nots .lt. 3) stop
'nots<3'
18 eps = 1.0/(splnot(nots)-splnot(1))
26 h(i) = splnot(i+1) - splnot(i)
27 IF (eps*h(i) .le. 1.e-8_dp) stop
'h(i)<1.e-8'
30 END SUBROUTINE initspline