cubically interpolate p on x to q on xx
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer | :: | nx | ||||
real(kind=wp) | :: | p(nx) | ||||
integer | :: | mx | ||||
real(kind=wp) | :: | q(mx) | ||||
integer | :: | ix(mx) | ||||
real(kind=wp) | :: | dxm(mx) | ||||
real(kind=wp) | :: | dx(mx) | ||||
real(kind=wp) | :: | dxp(mx) | ||||
real(kind=wp) | :: | dxpp(mx) |
subroutine cubt1(nx,p,mx,q,ix,dxm,dx,dxp,dxpp) implicit none integer :: mx,ix(mx),nx,i,ii real(wp) :: p(nx),q(mx),dxm(mx),dx(mx),dxp(mx),dxpp(mx) do ii=1,mx i = ix(ii) q(ii) = dxm(ii)*p(i-1)+dx(ii)*p(i)+dxp(ii)*p(i+1)+dxpp(ii)*p(i+2) end do end subroutine cubt1