LinearSolve.jl
LinearSolve.jl copied to clipboard
SciML
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.
If I want to solve a 12 parameter model: ``` A = rand(100,12) b = rand(100,1) x = A\b ``` works just fine in Julia or matlab. But the following...
It would be nice to be able to give solver parameters as arguments when defining the solver type. Currently I think one has to give the solver parameters as keywordargument...
Is this expected to work? (related to https://github.com/PainterQubits/Unitful.jl/issues/46): ```julia julia> using LinearSolve, Unitful julia> A = rand(3,3)u"nm"; x = rand(3)u"nm"; julia> solve(init(LinearProblem(A, x))) ERROR: DimensionError: nm and false are not...
Strumpack (https://github.com/pghysels/STRUMPACK) contains dense and sparse linear solvers, with OpenMP/MPI/CUDA backends. It is a relatively new package written by SuperLU developers. It implements a variety of hierarchical low-rank compression algorithms...
As discussed [here](https://discourse.julialang.org/t/discussion-on-why-i-no-longer-recommend-julia-by-yuri-vishnevsky/81151/240?u=bruno_amorim) it would be nice to support eigenvalues problems in the SciML ecosystem. Some of the things that should be included: - dense eigenproblems via the Julia's built...
We should put together a benchmark script and have a bunch of people run it. It should just run LUFactorization, RFLUFactorization, and FastLUFactorization (and MKLFactorization when that exists). It would...
According to the documentation, `RFLUFactorization` is supposed to be faster than backslash for matrices on the order of a few hundred. I did some quick benchmarking and I'm not finding...
https://github.com/SciML/LinearSolve.jl/issues/159#issuecomment-1188398811