SciML / LinearSolve.jl

LinearSolve.jl: High-Performance Unified Interface for Linear Solvers in Julia. Easily switch between factorization and Krylov methods, add preconditioners, and all in one interface.
https://docs.sciml.ai/LinearSolve/stable/
Other
245 stars 52 forks source link

Can't use Pivoted QR Factorization #335

Closed samuelsonric closed 1 year ago

samuelsonric commented 1 year ago

LinearSolve v2.1.13

When I try to solve a linear problem using the pivoted QR factorization, I see the following exception:

TypeError: in setfield!, expected LinearAlgebra.QRCompactWY{Int64, Matrix{Int64}, Matrix{Int64}}, got a value of type QRPivoted{Int64, Matrix{Int64}, Vector{Int64}, Vector{Int64}}

MWE:

using LinearAlgebra
using LinearSolve

lp = LinearProblem([1;;], [1])
ls = init(lp, QRFactorization(ColumnNorm(), 16, true))

solve!(ls)