Solving a linear system Ax=b with A being sparse matrix and b is a dense vector in Julia

Solution 1:

Define "doesn't work". Perhaps with a minimum (non)working example.

For instance, I just tried this and it seemed to work beautifully:

julia> using LinearAlgebra, SparseArrays
julia> m = sparse(repeat(1:1000, inner=10), rand(1:1000, 10000), rand(10000));

julia> rhs = rand(1000);

julia> m \ rhs
1000-element Vector{Float64}:
 -0.8976027048792632
  1.451352807805229
  3.6881731426288042
 -0.8042571078676164
 -2.7624771730498727

Without more information about your particular problem, there's nothing to be said.