From 9527c041d20a89218ebcf8f064b42b2bd7b40e43 Mon Sep 17 00:00:00 2001 From: Kanav Gupta Date: Sun, 14 Jul 2019 21:12:43 +0530 Subject: [PATCH] fix bdf and RadauIIA5 --- src/integrators/integrator_interface.jl | 1 + src/perform_step/bdf_perform_step.jl | 3 +-- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/src/integrators/integrator_interface.jl b/src/integrators/integrator_interface.jl index b49e18b710..1ca8034f40 100644 --- a/src/integrators/integrator_interface.jl +++ b/src/integrators/integrator_interface.jl @@ -75,6 +75,7 @@ end # the ordering of the cache arrays is important!!! @inline DiffEqBase.get_tmp_cache(integrator,alg::OrdinaryDiffEqAlgorithm,cache) = (cache.tmp,) +@inline DiffEqBase.get_tmp_cache(integrator,alg::RadauIIA5,cache) = (cache.tmp,cache.atmp) @inline DiffEqBase.get_tmp_cache(integrator,alg::OrdinaryDiffEqNewtonAdaptiveAlgorithm,cache) = (cache.nlsolver.tmp,cache.atmp) @inline DiffEqBase.get_tmp_cache(integrator,alg::OrdinaryDiffEqRosenbrockAdaptiveAlgorithm,cache) = (cache.tmp,cache.linsolve_tmp) @inline DiffEqBase.get_tmp_cache(integrator,alg::CompositeAlgorithm, cache) = get_tmp_cache(integrator, integrator.alg.algs[1], cache.caches[1]) diff --git a/src/perform_step/bdf_perform_step.jl b/src/perform_step/bdf_perform_step.jl index 08a5f144ae..7d78caf8ad 100644 --- a/src/perform_step/bdf_perform_step.jl +++ b/src/perform_step/bdf_perform_step.jl @@ -95,8 +95,7 @@ end cache.dtₙ₋₁ = dtₙ perform_step!(integrator, cache.eulercache, repeat_step) cache.fsalfirstprev .= integrator.fsalfirst - # why? - # nlsolver.tmp = tmp + nlsolver.tmp = tmp return end