Documentation

BSeries.Numerics.Reindex

Stage relabelling and the antisymmetric two-stage family #

Relabelling the stages of a Runge–Kutta tableau leaves all elementary weights unchanged (P-equivalence; Butcher, Numerical Methods for ODEs, Theorem 381H). As an application we construct the antisymmetric two-stage family Ω_λ of arXiv:2507.21006, Section 7 — the tableau

A = [λ λ; -λ -λ], b = (λ, -λ)

whose negation is a stage swap of itself, so its character is even.

def BSeries.RungeKutta.reindex {ι : Type u} {κ : Type v} {R : Type w} (e : κ ι) (rk : RungeKutta ι R) :

Relabel the stages of a tableau along an equivalence.

Equations
Instances For
    theorem BSeries.RungeKutta.stageWeight_reindex {ι : Type u} {κ : Type v} {R : Type w} [Fintype ι] [Fintype κ] [CommSemiring R] (e : κ ι) (rk : RungeKutta ι R) (t : HopfAlgebras.PTree) (i : κ) :
    (reindex e rk).stageWeight t i = rk.stageWeight t (e i)
    theorem BSeries.RungeKutta.stageWeightList_reindex {ι : Type u} {κ : Type v} {R : Type w} [Fintype ι] [Fintype κ] [CommSemiring R] (e : κ ι) (rk : RungeKutta ι R) (ts : List HopfAlgebras.PTree) (i : κ) :
    (reindex e rk).stageWeightList ts i = rk.stageWeightList ts (e i)
    theorem BSeries.RungeKutta.weight_reindex {ι : Type u} {κ : Type v} {R : Type w} [Fintype ι] [Fintype κ] [CommSemiring R] (e : κ ι) (rk : RungeKutta ι R) (t : HopfAlgebras.PTree) :
    (reindex e rk).weight t = rk.weight t

    Stage relabelling preserves the elementary weights (P-equivalence; Butcher, Theorem 381H).

    theorem BSeries.RungeKutta.treeWeight_reindex {ι : Type u} {κ : Type v} {R : Type w} [Fintype ι] [Fintype κ] [CommSemiring R] (e : κ ι) (rk : RungeKutta ι R) (τ : HopfAlgebras.RootedTree) :
    (reindex e rk).treeWeight τ = rk.treeWeight τ

    The stage swap of a two-stage scheme.

    Equations
    Instances For

      The antisymmetric two-stage tableau Ω_λ of arXiv:2507.21006, Section 7: A = [λ λ; -λ -λ], b = (λ, -λ).

      Equations
      Instances For

        Negating Ω_λ is exactly the swap of its two stages.

        The antisymmetric family is even: the character of Ω_λ is fixed by the canonical involution (arXiv:2507.21006, Section 7).

        The quadrature weights of Ω_λ sum to zero: Ω_λ is inconsistent (a "do-nothing to first order" scheme).

        The elementary weight of the two-chain detects λ: ω_λ(chain₂) = 4λ².