diff --git a/ext/LinearSolveBandedMatricesExt.jl b/ext/LinearSolveBandedMatricesExt.jl index 7b4d6a054..f402e06fa 100644 --- a/ext/LinearSolveBandedMatricesExt.jl +++ b/ext/LinearSolveBandedMatricesExt.jl @@ -6,7 +6,7 @@ import LinearSolve: defaultalg, # Defaults for BandedMatrices function defaultalg(A::BandedMatrix, b, ::OperatorAssumptions) - return DefaultLinearSolver(DefaultAlgorithmChoice.QRFactorization) + return DefaultLinearSolver(DefaultAlgorithmChoice.DirectLdiv!) end function defaultalg(A::Symmetric{<:Number, <:BandedMatrix}, b, ::OperatorAssumptions)