FIXME : Add documentation
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real, | intent(in), | dimension(:) | :: | r | ||
real, | intent(in), | dimension(:) | :: | data | ||
real, | intent(in), | dimension(:) | :: | x | ||
real, | intent(out), | dimension(:) | :: | dint | ||
real, | intent(out), | dimension(:) | :: | ddint |
subroutine inter_d_cspl(r,data,x,dint,ddint)
implicit none
real, dimension(:), intent(in) :: r, data
real, dimension(:), intent(in) :: x
real, dimension(:), intent(out) :: dint, ddint
integer :: i, m
type(spline) :: spl
real, parameter :: tension = 1.0
spl = new_spline(r, data, tension)
m = size(x)
do i = 1, m
dint(i) = spl%interpolate(x(i))
ddint(i)= spl%derivative(x(i))
end do
end subroutine inter_d_cspl