Interpolation function. Can be used for dense output after a step. It calls the low-level routine dintp.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
class(ddeabm_class), | intent(inout) | :: | me | |||
real(kind=wp), | intent(in) | :: | tc |
point at which solution is desired |
||
real(kind=wp), | intent(out), | dimension(me%neq) | :: | yc |
interpolated state at |
subroutine ddeabm_interp(me, tc, yc) implicit none class(ddeabm_class), intent(inout) :: me real(wp), intent(in) :: tc !! point at which solution is desired real(wp), dimension(me%neq), intent(out) :: yc !! interpolated state at `tc` ! interpolate to get the state at tc: call dintp(me%x, me%yy, tc, yc, & me%ypout, me%neq, me%kold, me%phi, & me%ivc, me%iv, me%kgi, me%gi, me%alpha, & me%g, me%w, me%xold, me%p) end subroutine ddeabm_interp