H A D | spline.f90 | 2 subroutine spline(f,x,n,t,d2f) argument 8 double precision x(n),f(n),d2f(n),t local 17 d2f(1)=(f(2)-f(1))/(x(2)-x(1)) & 23 d2f(jj)=(f(jj+1)-f(jj))/(x(jj+1)-x(jj)) & 29 d2f(n)=(f(1)-f(n))/(x(1)-x(n)+t) & 37 call baksublu(cmat,d2f,n,n,iperm) 44 subroutine dspline(f,x,d2f,n,df) argument 49 double precision f(n),x(n),d2f(n),df(n) local 67 subroutine dsplv(f,x,t,d2f,n,xv,dfv) argument 107 subroutine evalspline(f,x,t,d2f,n,xv,fv) argument [all …]
|