MRKISS 2025-09-10
A tiny library with zero dependencies that aims to make it easy to use & experiment with explicit Runge-Kutta methods.
Loading...
Searching...
No Matches
mrkiss_erk_ralston_2 Module Reference

Butcher tableau for ralston's 2 stage order (4) Runge-Kutta method. More...

Variables

integer, parameter s = 2
 The order of the overall method.
integer, parameter m = 1
 Number of methods.
real(kind=rk), dimension(s, s), parameter a
 The \(\mathbf{a}\) matrix for the Butcher Tableau.
real(kind=rk), dimension(s, m), parameter b
 The \(\mathbf{b}\) matrix for the Butcher Tableau.
real(kind=rk), dimension(s), parameter c
 The \(\mathbf{c}\) matrix for the Butcher Tableau.
integer, dimension(m), parameter p = [2]
 The method orders.
integer, dimension(m), parameter se = [2]
 Number of stages for each method.

Detailed Description

Butcher tableau for ralston's 2 stage order (4) Runge-Kutta method.

IMO
This method was designed to minimize truncation error over every other consideration includeing stability, roundoff error, performance, and storage.
Stability Image Links
References:
  • Ralston (1962); Runge-Kutta Methods with Minimum Error Bounds; Math. of Comp. 16 (80); p431-437; zotero://select/items/0_2FK7F4J3

Variable Documentation

◆ s

integer, parameter mrkiss_erk_ralston_2::s = 2

The order of the overall method.

Definition at line 60 of file mrkiss_erk_ralston_2.f90.

◆ m

integer, parameter mrkiss_erk_ralston_2::m = 1

Number of methods.

Definition at line 62 of file mrkiss_erk_ralston_2.f90.

◆ a

real(kind=rk), dimension(s,s), parameter mrkiss_erk_ralston_2::a

The \(\mathbf{a}\) matrix for the Butcher Tableau.

Definition at line 64 of file mrkiss_erk_ralston_2.f90.

64 real(kind=rk), parameter :: a(s,s) = reshape([0.0_rk, 0.0_rk, &
65 2.0_rk, 0.0_rk], [s, s]) / 3.0_rk

◆ b

real(kind=rk), dimension(s,m), parameter mrkiss_erk_ralston_2::b

The \(\mathbf{b}\) matrix for the Butcher Tableau.

Definition at line 67 of file mrkiss_erk_ralston_2.f90.

67 real(kind=rk), parameter :: b(s,m) = reshape([1.0_rk, 3.0_rk], [s, m]) / 4.0_rk

◆ c

real(kind=rk), dimension(s), parameter mrkiss_erk_ralston_2::c

The \(\mathbf{c}\) matrix for the Butcher Tableau.

Definition at line 69 of file mrkiss_erk_ralston_2.f90.

69 real(kind=rk), parameter :: c(s) = [ 0.0_rk, 2.0_rk] / 3.0_rk

◆ p

integer, dimension(m), parameter mrkiss_erk_ralston_2::p = [2]

The method orders.

Definition at line 71 of file mrkiss_erk_ralston_2.f90.

◆ se

integer, dimension(m), parameter mrkiss_erk_ralston_2::se = [2]

Number of stages for each method.

Definition at line 73 of file mrkiss_erk_ralston_2.f90.