try QR

Author: ChrisRackauckas

test/Hydraulic/isothermal_compressible.jl

@@ -151,7 +151,7 @@ end
     @mtkbuild sys = System()
     prob = ODEProblem(sys, [], (0, 5))
     @show prob.u0
-    sol = solve(prob, Rodas5P(); abstol = 1e-6, reltol = 1e-9)
+    sol = solve(prob, Rodas5P(linsolve = OrdinaryDiffEq.OrdinaryDiffEqNonlinearSolve.LinearSolve.QRFactorization(OrdinaryDiffEq.OrdinaryDiffEqNonlinearSolve.LinearSolve.ColumnNorm())); abstol = 1e-6, reltol = 1e-9)
     @show sol[1]
     # begin
     #     fig = Figure()