fixing type stability, running benchmarks.

QuEraComputing · Dec 22, 2023 · 04b834d · 04b834d
1 parent 3df38f3
commit 04b834d
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 6 deletions.
diff --git a/benchmarks/rabi-diffeq.jl b/benchmarks/rabi-diffeq.jl
@@ -13,7 +13,7 @@ u0 = ComplexF64[1.0, 0.0]
 tspan = (0.0, 2π)
 prob = ODEProblem(f, u0, tspan)
 # get precompilation out of the way
-sol = solve(prob, DP5(), reltol=1e-10, abstol=1e-10)
+sol = solve(prob, DP8(), reltol=1e-10, abstol=1e-10)
 
 # terminate benchmark after maximum of 5 minutes
-@benchmark solve(prob, DP5(), reltol=1e-10, abstol=1e-10) samples=10000 evals=5 seconds=300
+@benchmark solve(prob, DP8(), reltol=1e-10, abstol=1e-10) samples=10000 evals=5 seconds=300
diff --git a/benchmarks/rabi-dormand-prince.jl b/benchmarks/rabi-dormand-prince.jl
@@ -1,5 +1,5 @@
 using BenchmarkTools
-using DormandPrince:DP5Solver,  integrate
+using DormandPrince:DP5Solver, DP8Solver,  integrate
 
 function fcn(x, y, f)
     g(x) = 2.2*2π*sin(2π*x)
@@ -8,12 +8,12 @@ function fcn(x, y, f)
     f[2] = -1im * g(x)/2 * y[1]
 end
 
-solver = DP5Solver(
+solver = DP8Solver(
     fcn,
     0.0,
     ComplexF64[1.0, 0.0]
 )
 
 integrate(solver, 2π)
 
-@benchmark  integrate(clean_solver, 2π) setup=(clean_solver = DP5Solver(fcn, 0.0, ComplexF64[1.0, 0.0])) samples=10000 evals=5 seconds=500
+@benchmark  integrate(clean_solver, 2π) setup=(clean_solver = DP8Solver(fcn, 0.0, ComplexF64[1.0, 0.0])) samples=10000 evals=5 seconds=500
diff --git a/src/types.jl b/src/types.jl
@@ -113,7 +113,7 @@ struct DP5Solver{T, StateType , F, OptionsType, ConstsType, VarsType} <: Abstrac
         consts = Consts(expo1, options)
         vars = Vars{T}(;x=x, h=options.initial_step_size)
 
-        new{T, StateType, F, typeof(options), typeof(consts), typeof(varS)}(
+        new{T, StateType, F, typeof(options), typeof(consts), typeof(vars)}(
             f, y, k1, k2, k3, k4, k5, k6, y1, ysti, options, consts, vars
         )
     end