89 real(kind=
rp),
intent(out) :: coeffs(4)
90 real(kind=
rp),
intent(in) :: dt(10)
91 integer,
intent(in) :: order
92 real(kind=
rp) dta, dtb, dtc, dtd, dte, dts
100 coeffs(2) = -0.5_rp * dt(1) / dt(2)
101 coeffs(1) = 1.0_rp - coeffs(2)
109 coeffs(3) = dte*( 0.5_rp*dtb + dtc/3.0_rp )
110 coeffs(2) = -0.5_rp * dta - coeffs(3) * dtd
111 coeffs(1) = 1.0_rp - coeffs(2) - coeffs(3)
113 call neko_error(
"The order of the AB time scheme must be 1 to 3.")
Interface for setting the scheme coefficients.
Adam-Bashforth scheme for time integration.
subroutine ab_time_scheme_compute_coeffs(coeffs, dt, order)
Compute the scheme coefficients.
subroutine, public rzero(a, n)
Zero a real vector.
integer, parameter, public rp
Global precision used in computations.
Base class for time integration schemes.
Explicit Adam-Bashforth scheme for time integration.
Base abstract class for time integration schemes.