From 6c058417c408e871242877391ab81ca9116c6792 Mon Sep 17 00:00:00 2001 From: ChrisRackauckas-Claude Date: Thu, 9 Jul 2026 21:41:20 -0400 Subject: [PATCH] fix(tests): use BrownFullBasicInit for KdV single-soliton test MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Over-specified boundary derivatives leave an O(dx^2) residual that OrdinaryDiffEq CheckInit rejects (normresid ≈ 3.6e-5 > 1e-6). Initialize algebraic variables with BrownFullBasicInit, matching the Complex/schroedinger tests. Co-Authored-By: Chris Rackauckas --- test/Higher_Order/MOL_1D_HigherOrder.jl | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/test/Higher_Order/MOL_1D_HigherOrder.jl b/test/Higher_Order/MOL_1D_HigherOrder.jl index 5f5820cea..54b6d7f9e 100644 --- a/test/Higher_Order/MOL_1D_HigherOrder.jl +++ b/test/Higher_Order/MOL_1D_HigherOrder.jl @@ -4,6 +4,7 @@ using ModelingToolkit, MethodOfLines, LinearAlgebra, Test, OrdinaryDiffEq, DomainSets using SciMLBase using ModelingToolkit: Differential +using DiffEqBase: BrownFullBasicInit # Beam Equation #! Broken due to overlapping inner corner bc eqs - determine state sharing heuristic; sort by dorder and give precedence to higher @@ -131,7 +132,7 @@ end @named pdesys = PDESystem(eq, bcs, domains, [x, t], [u(x, t)]) prob = discretize(pdesys, discretization) - sol = solve(prob, FBDF(), saveat = 0.1) + sol = solve(prob, FBDF(), saveat = 0.1, initializealg = BrownFullBasicInit()) @test SciMLBase.successful_retcode(sol)