From 2713004355a92c7778ade6a65a66acc0baff2776 Mon Sep 17 00:00:00 2001 From: Xingjian Guo Date: Mon, 30 Jul 2018 12:45:23 -0400 Subject: [PATCH] Bugfix for `calc_W!` --- src/derivative_utils.jl | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/derivative_utils.jl b/src/derivative_utils.jl index b4aeae7873..e05443c368 100644 --- a/src/derivative_utils.jl +++ b/src/derivative_utils.jl @@ -282,7 +282,7 @@ function calc_W!(integrator, cache::OrdinaryDiffEqConstantCache, dtgamma, repeat if !W_transform if DiffEqBase.has_jac(f) J = f.jac(uprev, p, t) - if !isa(f, DiffEqBase.AbstractDiffEqLinearOperator) + if !isa(J, DiffEqBase.AbstractDiffEqLinearOperator) J = DiffEqArrayOperator(J) end W = WOperator(mass_matrix, dtgamma, J; transform=false) @@ -297,7 +297,7 @@ function calc_W!(integrator, cache::OrdinaryDiffEqConstantCache, dtgamma, repeat else if DiffEqBase.has_jac(f) J = f.jac(uprev, p, t) - if !isa(f, DiffEqBase.AbstractDiffEqLinearOperator) + if !isa(J, DiffEqBase.AbstractDiffEqLinearOperator) J = DiffEqArrayOperator(J) end W = WOperator(mass_matrix, dtgamma, J; transform=true)