From 920fade3ae01daacd343669e86590d7196ffdd18 Mon Sep 17 00:00:00 2001 From: Christopher Rackauckas Date: Fri, 6 Oct 2023 08:43:30 +0200 Subject: [PATCH 1/3] HOTFIX: Don't use MKL with numbers other that Float64 --- src/default.jl | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/default.jl b/src/default.jl index 94ec4d658..36dadb21b 100644 --- a/src/default.jl +++ b/src/default.jl @@ -170,7 +170,7 @@ function defaultalg(A, b, assump::OperatorAssumptions) DefaultAlgorithmChoice.RFLUFactorization #elseif A === nothing || A isa Matrix # alg = FastLUFactorization() - elseif usemkl + elseif usemkl && eltype(A) <: Float64 DefaultAlgorithmChoice.MKLLUFactorization else DefaultAlgorithmChoice.LUFactorization @@ -179,7 +179,7 @@ function defaultalg(A, b, assump::OperatorAssumptions) DefaultAlgorithmChoice.QRFactorization elseif __conditioning(assump) === OperatorCondition.SuperIllConditioned DefaultAlgorithmChoice.SVDFactorization - elseif usemkl + elseif usemkl && eltype(A) <: Float64 DefaultAlgorithmChoice.MKLLUFactorization else DefaultAlgorithmChoice.LUFactorization From ed9d147992aafa0329985e703339b1d81e0b6d88 Mon Sep 17 00:00:00 2001 From: Christopher Rackauckas Date: Fri, 6 Oct 2023 08:50:53 +0200 Subject: [PATCH 2/3] Update src/default.jl --- src/default.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/default.jl b/src/default.jl index 36dadb21b..66a73e373 100644 --- a/src/default.jl +++ b/src/default.jl @@ -179,7 +179,7 @@ function defaultalg(A, b, assump::OperatorAssumptions) DefaultAlgorithmChoice.QRFactorization elseif __conditioning(assump) === OperatorCondition.SuperIllConditioned DefaultAlgorithmChoice.SVDFactorization - elseif usemkl && eltype(A) <: Float64 + elseif usemkl && eltype(A) <: Union{Float64, Float32} DefaultAlgorithmChoice.MKLLUFactorization else DefaultAlgorithmChoice.LUFactorization From 473cb4707b23d65e35e1749a4de2f47d4d54b7e9 Mon Sep 17 00:00:00 2001 From: Christopher Rackauckas Date: Fri, 6 Oct 2023 08:50:58 +0200 Subject: [PATCH 3/3] Update src/default.jl --- src/default.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/default.jl b/src/default.jl index 66a73e373..60cbfcc48 100644 --- a/src/default.jl +++ b/src/default.jl @@ -170,7 +170,7 @@ function defaultalg(A, b, assump::OperatorAssumptions) DefaultAlgorithmChoice.RFLUFactorization #elseif A === nothing || A isa Matrix # alg = FastLUFactorization() - elseif usemkl && eltype(A) <: Float64 + elseif usemkl && eltype(A) <: Union{Float64, Float32} DefaultAlgorithmChoice.MKLLUFactorization else DefaultAlgorithmChoice.LUFactorization