1 subroutine r8evspline(xget,x,nx,ilinx,f,ict,fval,ier)
15 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
80 call r8herm1x(xget,x,nx,ilinx,i(1),xparam(1),hx(1),hxi(1),ier)
83 call r8fvspline(ict,1,1,fval,i,xparam,hx,hxi,f,nx)
91 subroutine r8fvspline(ict,ivec,ivecd,
92 > fval,ii,xparam,hx,hxi,fin,nx)
99 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
103 real*8 xp,xpi,xp2,xpi2,cx,cxi,hx2,cxd,cxdi
142 data sixth/0.166666666666666667_r8/
167 cxi=xpi*(xpi2-1.0_r8)
172 sum=xpi*fin(0,i) +xp*fin(0,i+1)
173 sum=sum+sixth*hx2*(cxi*fin(1,i) +cx*fin(1,i+1))
192 cxd=3.0_r8*xp2-1.0_r8
193 cxdi=-3.0_r8*xpi2+1.0_r8
195 sum=hxi(v)*(fin(0,i+1)-fin(0,i))
196 sum=sum+sixth*hx(v)*(cxdi*fin(1,i) +cxd*fin(1,i+1))
213 sum=xpi*fin(1,i) +xp*fin(1,i+1)
225 fval(v,iadr)=hxi(v)*(fin(1,i+1)-fin(1,i))