Merge pull request #2525 from SciML/precision

Fix precision mixing in RodasTableau and better test precision mix

Author: ChrisRackauckas

lib/OrdinaryDiffEqRosenbrock/src/rosenbrock_tableaus.jl

@@ -225,7 +225,7 @@ end
 struct RodasTableau{T, T2}
     A::Matrix{T}
     C::Matrix{T}
-    gamma::T
+    gamma::T2
     c::Vector{T2}
     d::Vector{T}
     H::Matrix{T}

test/interface/precision_mixing.jl

@@ -0,0 +1,22 @@
+using OrdinaryDiffEq
+function ODE(du, u, t, R, K)
+    du .= u
+end
+params = BigFloat[1. 0.91758707304098; 1.48439909482661 1.]
+u0 = BigFloat[0.1, 0.1]
+tspan = (1.0, 31.0)
+R = BigFloat[0.443280390004304303, 1.172917082211452]
+K = BigFloat[13.470600276901400604, 12.52980757005]
+ODE_ = (du, u, params, t) -> ODE(du, u, t, R, K)
+odeProblem = ODEProblem(ODE_, u0, tspan, params)
+for alg in [AutoVern8(Rodas5(), nonstifftol = 11 / 10)
+            FBDF()
+            QNDF()
+            Tsit5()
+            Rodas5P()
+            TRBDF2()
+            KenCarp4()
+            RadauIIA5()
+    ]
+    Solution = solve(odeProblem, alg, saveat = 1, abstol = 1.e-12, reltol = 1.e-6)
+end

test/runtests.jl

@@ -85,6 +85,7 @@ end
     if !is_APPVEYOR && (GROUP == "All" || GROUP == "InterfaceIV" || GROUP == "Interface")
         @time @safetestset "Autodiff Error Tests" include("interface/autodiff_error_tests.jl")
         @time @safetestset "Ambiguity Tests" include("interface/ambiguity_tests.jl")
+        @time @safetestset "Precision Mixing Tests" include("interface/precision_mixing.jl")
         @time @safetestset "Sized Matrix Tests" include("interface/sized_matrix_tests.jl")
         @time @safetestset "Second Order with First Order Solver Tests" include("interface/second_order_with_first_order_solvers.jl")
     end