Sparse solve with sparse rhs

Specific Domains Numerics
,
14

BasicLU.jl supports sparse right-hand side for information.

m = 1000
F = basiclu_object(m)
A = sprand(m, m, 5e-3) + I
err = factorize(F, A)
b = sparsevec([1], [1.0], m)
x = solve(F, b, 'N')  # sparse solution

Référence: L. Schork and J. Gondzio, Permuting Spiked Matrices to Triangular Form and its Application to the Forrest-Tomlin Update