From 8baa350ce7b716d0dfac30561393b3e108714e3a Mon Sep 17 00:00:00 2001 From: Chris Rackauckas Date: Mon, 13 Dec 2021 10:31:05 -0500 Subject: [PATCH] Avoid RecursiveFactorization except on simple number types --- src/linear_nonlinear.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/linear_nonlinear.jl b/src/linear_nonlinear.jl index bfe03e333..0b23d6b39 100644 --- a/src/linear_nonlinear.jl +++ b/src/linear_nonlinear.jl @@ -101,7 +101,7 @@ function (p::DefaultLinSolve)(x,A::Union{AbstractMatrix,AbstractDiffEqOperator}, p.openblas = isopenblas() end - if ArrayInterface.can_setindex(x) && (size(A,1) <= 100 || (p.openblas && size(A,1) <= 500)) + if eltype(A) <: Union{Float32,Float64,ComplexF32,ComplexF64} && ArrayInterface.can_setindex(x) && (size(A,1) <= 100 || (p.openblas && size(A,1) <= 500)) p.A = RecursiveFactorization.lu!(A) else p.A = lu!(A)