rkdp87_dense_coefficients Module Procedure

module procedure rkdp87_dense_coefficients pure module subroutine rkdp87_dense_coefficients(me, b, form)

Dense output for the Prince-Dormand 8(7) (RK8(7)13M) method rkdp87. This is an interpolant of order 8 constructed by P. Stone, which requires 8 additional function evaluations (see rkdp87_dense_extra_stages).

Note

The interpolant is published in the polynomial form. Here it has been converted (exactly, using high precision arithmetic) to the form used in DOPRI5 and DOP853 (see dense_coefficients), which has much smaller roundoff error for these coefficients. This code is generated by scripts/generate_dense.py.

References

Arguments

Type IntentOptional Attributes Name
class(rkdp87_class), intent(in) :: me
real(kind=wp), intent(out), dimension(:,:), allocatable :: b

dense output coefficients

integer, intent(out) :: form

the form of b