# Benchmarks for our own new differential equation solver

**URL:** <https://discourse.julialang.org/t/benchmarks-for-our-own-new-differential-equation-solver/81109>\
**Category:** Modelling & Simulations\
**Created:** [May 15, 2022, 6:00pm UTC](https://discourse.julialang.org/t/benchmarks-for-our-own-new-differential-equation-solver/81109 "2022-05-15T18:00:20Z")\
**Posts on this page:** 3\
**Page:** 1

<div class="post-metadata">

**Author:** ![Mikel\_Antonana](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/mikel_antonana/32/3095_2.png) [@Mikel\_Antonana](https://discourse.julialang.org/u/Mikel_Antonana)\
**Post date:** [May 15, 2022, 6:00pm UTC](https://discourse.julialang.org/t/benchmarks-for-our-own-new-differential-equation-solver/81109/1 "2022-05-15T18:00:20Z")

</div>

Hi everyone,

We are developing a new differential equation solver,

````julia
function DiffEqBase.__solve(prob::DiffEqBase.AbstractODEProblem{uType,tType,isinplace},
    alg::IRKAlgorithm, udim::Array{floatType,dim},args...;
    s=8,
    dt,
    m=1,
    initial_interp=true,
    save_everystep=true,
    partitioned=false) where {floatType<: Union{Float32,Float64}, uType, tType, isinplace, dim}```

but when I try to execute a benchmark

````

abstols=reltols=[1e-4, 1e-6]  
dts=[10.0, 20.0]

setups = [  
Dict(:alg=\>Vern9())  
Dict(:alg=\>IRKGL\_simd(),:udim=\>udim,:s=\>8,:dts=\>dts)  
]

solnames =[“Vern9”, “IRKGL\_simd”]  
wp = WorkPrecisionSet(prob,abstols,reltols,setups;appxsol=test\_sol,save\_everystep=false,numruns=1,maxiters=10^9, names=solnames);

```julia

I get the next error,

```

MethodError: no method matching alg\_cache(::IRKGL\_simd, ::Array{Float64, 3}, ::Array{Float64, 3}, ::Type{Float64}, …```

Do I something wrong?

Thanks  
M.

---

<div class="post-metadata">

**Author:** ![ChrisRackauckas](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/chrisrackauckas/32/77_2.png) [@ChrisRackauckas](https://discourse.julialang.org/u/ChrisRackauckas)\
**Post date:** [May 15, 2022, 8:11pm UTC](https://discourse.julialang.org/t/benchmarks-for-our-own-new-differential-equation-solver/81109/2 "2022-05-15T20:11:25Z")

</div>

> [@Mikel\_Antonana](#):
>
> `udim::Array{floatType,dim}`

`udim` isn’t a supported argument, this should go into `IRKAlgorithm`. `m`,`s`, `partitioned`, and `initial_interp` are also not common args, so they should be in the alg.

Other than that, I’d need more details (like the error message) to help more.

---

<div class="post-metadata">

**Author:** ![Mikel\_Antonana](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/mikel_antonana/32/3095_2.png) [@Mikel\_Antonana](https://discourse.julialang.org/u/Mikel_Antonana)\
**Post date:** [May 16, 2022, 5:31pm UTC](https://discourse.julialang.org/t/benchmarks-for-our-own-new-differential-equation-solver/81109/3 "2022-05-16T17:31:11Z")

</div>

You can download the notebook (Benchmark error.html), where you can check the details of the experiment (I do not know if it is enough to help us) from: [IRKGL\_SIMD.jl/Benchmarks at master · mikelehu/IRKGL\_SIMD.jl · GitHub](https://github.com/mikelehu/IRKGL_SIMD.jl/tree/master/Benchmarks)

If we cannot use the BenchmarkTools package, the alternative will be to perform the comparisons manually.

Thanks