LinearSolve.jl

LinearSolve.jl copied to clipboard

in default alg branch on ```julia if SciMLBase.has_ldiv!(A) # do ldiv!(u, A, b) via solve_function elseif SciMLBase.has_ldiv(A) # do u .= A \ b via solve_function end ```

vpuri3

WIP adjoint

2

should close https://github.com/SciML/LinearSolve.jl/issues/92

vpuri3

linearsolve using the transpose of the factorization

10

For defining rrules, I need to make use of the factorization in the forward pass. The rrules is also a linearsolve with the transpose of the linear system, how to...

MKAbdElrahman

AMGCL wrapper?

5

https://github.com/ddemidov/amgcl Has OpenMP/MPI/CUDA implementations, similar to Hypre but much more lightweight. Already has a [Python interface](https://github.com/ddemidov/amgcl/tree/1.4.2/pyamgcl) (see [test](https://github.com/ddemidov/amgcl/blob/1.4.2/tests/test_pyamgcl.py))

learning-chip

https://github.com/JuliaParallel/Elemental.jl

ChrisRackauckas

http://icl.cs.utk.edu/projectsfiles/magma/doxygen/_m_a_g_m_a-sparse.html#sparsesolvers

ChrisRackauckas

Some exist as jlls but not in Trilinos.jl. The newest ones do not.

ChrisRackauckas

http://icl.eecs.utk.edu/slate/

ChrisRackauckas

Better initial guesses for reinitialized iterative solvers

1

https://scicomp.stackexchange.com/questions/23323/role-of-initial-guess-for-iterative-linear-solver https://krypy.readthedocs.io/en/latest/krypy.utils.html#krypy.utils.hegedus http://math.aalto.fi/~kuiduhut/spkry.pdf https://fenicsproject.org/qa/7253/convergence-of-iterative-solver-with-initial-guess/

ChrisRackauckas

LU factorizations use `lu_instance`, but the other factorizations do not init a cache. This means that blank `A`, such as a just initialized Jacobian, can call the factorization which costs...

ChrisRackauckas