dintdy Subroutine

public subroutine dintdy(T, K, Yh, Nyh, Dky, Iflag)

DESCRIPTION

DINTDY computes interpolated values of the K-th derivative of the dependent variable vector y, and stores it in DKY. This routine is called within the package with K = 0 and T = TOUT, but may also be called by the user for any K up to the current order. (See detailed instructions in the usage documentation.)

The computed values in DKY are gotten by interpolation using the Nordsieck history array YH. This array corresponds uniquely to a vector-valued polynomial of degree NQCUR or less, and DKY is set to the K-th derivative of this polynomial at T. The formula for DKY is:

               q
   DKY(i)  =  sum  c(j,K) * (T - tn)**(j-K) * h**(-j) * YH(i,j+1)
              j=K

where c(j,K) = j(j-1)…*(j-K+1), q = NQCUR, tn = TCUR, h = HCUR.

The quantities nq = NQCUR, l = nq+1, N = NEQ, tn, and h are communicated by COMMON. The above sum is done in reverse order. IFLAG is returned negative if either K or T is out of bounds.

Arguments

Type IntentOptional Attributes Name
real(kind=dp) :: T
integer :: K
real(kind=dp), intent(in) :: Yh(Nyh,*)
integer, intent(in) :: Nyh
real(kind=dp), intent(inout) :: Dky(*)
integer, intent(out) :: Iflag

Calls

proc~~dintdy~2~~CallsGraph proc~dintdy~2 M_odepack::dintdy proc~xerrwd~2 M_odepack::xerrwd proc~dintdy~2->proc~xerrwd~2

Called by

proc~~dintdy~2~~CalledByGraph proc~dintdy~2 M_odepack::dintdy proc~dlsodar~2 M_odepack::dlsodar proc~dlsodar~2->proc~dintdy~2 proc~drchek M_odepack::drchek proc~dlsodar~2->proc~drchek proc~dlsoda~2 M_odepack::dlsoda proc~dlsoda~2->proc~dintdy~2 proc~dlsodes~2 M_odepack::dlsodes proc~dlsodes~2->proc~dintdy~2 proc~dlsode~2 M_odepack::dlsode proc~dlsode~2->proc~dintdy~2 proc~dlsodis~2 M_odepack::dlsodis proc~dlsodis~2->proc~dintdy~2 proc~dlsodi~2 M_odepack::dlsodi proc~dlsodi~2->proc~dintdy~2 proc~dlsodkr~2 M_odepack::dlsodkr proc~dlsodkr~2->proc~dintdy~2 proc~dlsodkr~2->proc~drchek proc~dlsodpk~2 M_odepack::dlsodpk proc~dlsodpk~2->proc~dintdy~2 proc~dlsoibt~2 M_odepack::dlsoibt proc~dlsoibt~2->proc~dintdy~2 proc~drchek->proc~dintdy~2