fix benchmarks

Author: ChrisRackauckas

src/benchmark.jl

@@ -196,7 +196,7 @@ function WorkPrecision(prob,alg,abstols,reltols,dts=nothing;
     t = t/numruns
 
     if appxsol != nothing
-      errsol = calculate_errsol(prob,sol,approxsol)
+      errsol = calculate_errsol(prob,sol,appxsol)
       errors[i] = mean(errsol.errors[error_estimate])
     else
       errors[i] = mean(sol.errors[error_estimate])
@@ -254,7 +254,7 @@ function WorkPrecision(prob::Union{AbstractRODEProblem,AbstractSDEProblem},
                                   dense_errors = dense_errors)
       end
       if appxsol != nothing
-        errsol = calculate_errsol(prob,sol,approxsol)
+        errsol = calculate_errsol(prob,sol,appxsol)
         local_errors[j] = errsol.errors[error_estimate]
       else
         local_errors[j] = sol.errors[error_estimate]
@@ -268,22 +268,22 @@ function WorkPrecision(prob::Union{AbstractRODEProblem,AbstractSDEProblem},
   return WorkPrecision(prob,abstols,reltols,errors,times,name,N)
 end
 
-function calculate_errsol(prob,sol::AbstractODESolution,approxsol_setup::Dict)
-  true_sol = solve(prob,approxsol_setup[i][:alg];approxsol_setup[i]...)
+function calculate_errsol(prob,sol::AbstractODESolution,appxsol_setup::Dict)
+  true_sol = solve(prob,appxsol_setup[i][:alg];appxsol_setup[i]...)
   appxtrue(sol,true_sol)
 end
 
-function calculate_errsol(prob::AbstractSDEProblem,sol::AbstractRODESolution,approxsol_setup::Dict)
+function calculate_errsol(prob::AbstractSDEProblem,sol::AbstractRODESolution,appxsol_setup::Dict)
   prob2 = SDEProblem(prob.f,prob.g,prob.u0,prob.tspan,noise=NoiseWrapper(sol.W))
-  true_sol = solve(prob2,approxsol_setup[i][:alg];approxsol_setup[i]...)
+  true_sol = solve(prob2,appxsol_setup[i][:alg];appxsol_setup[i]...)
   appxtrue(sol,true_sol)
 end
 
-function calculate_errsol(prob,sol::AbstractODESolution,approxsol_setup::AbstractTimeseriesSolution)
+function calculate_errsol(prob,sol::AbstractODESolution,true_sol::AbstractTimeseriesSolution)
   appxtrue(sol,true_sol)
 end
 
-function calculate_errsol(prob::MonteCarloProblem,sol,approxsol_setup)
+function calculate_errsol(prob::MonteCarloProblem,sol,true_sol)
   appxtrue(sol,true_sol)
 end
 

src/ode_tableaus.jl

@@ -81,7 +81,8 @@ end
 Euler's method.
 """
 function constructEuler(T::Type = Float64)
-  A = [0]
+  A = Matrix{T}(1,1)
+  A[1] = 0
   c = [0]
   α = [1]
   A = map(T,A)
@@ -141,7 +142,8 @@ end
 Implicit Euler Method
 """
 function constructImplicitEuler(T::Type = Float64)
-  A = [1]
+  A = Matrix{T}(1,1)
+  A[1] = 1
   c = [1]
   α = [1]
   A = map(T,A)
@@ -154,7 +156,8 @@ end
 Order 2 Midpoint Method
 """
 function constructMidpointRule(T::Type = Float64)
-  A = [1//2]
+  A = Matrix{T}(1,1)
+  A[1] = 1//2
   c = [1//2]
   α = [1]
   A = map(T,A)

src/test_solution.jl

@@ -7,19 +7,20 @@ type TestSolution{T,N,hasinterp,tType,uType,iType} <: AbstractTimeseriesSolution
   u::uType
   interp::iType
   dense::Bool
+  retcode::Symbol
 end
 (T::TestSolution)(t) = T.interp(t)
 function TestSolution(t,u)
   T = eltype(eltype(u))
   N = length((size(u[1])..., length(u)))
-  TestSolution{T,N,false,typeof(t),typeof(u),Void}(t,u,nothing,false)
+  TestSolution{T,N,false,typeof(t),typeof(u),Void}(t,u,nothing,false,:Success)
 end
 function TestSolution(t,u,interp)
   T = eltype(eltype(u))
   N = length((size(u[1])..., length(u)))
-  TestSolution{T,N,true,typeof(t),typeof(u),typeof(interp)}(t,u,interp,true)
+  TestSolution{T,N,true,typeof(t),typeof(u),typeof(interp)}(t,u,interp,true,:Success)
 end
-TestSolution(interp::DESolution) = TestSolution{Void,0,true,Void,Void,typeof(interp)}(nothing,nothing,interp,true)
+TestSolution(interp::DESolution) = TestSolution{Void,0,true,Void,Void,typeof(interp)}(nothing,nothing,interp,true,:Success)
 hasinterp{T,N,hi,tType,uType,iType}(::TestSolution{T,N,hi,tType,uType,iType}) = hi
 """
 `appxtrue(sol::AbstractODESolution,sol2::TestSolution)`