求解线性系统Ax=b，其中a是稀疏矩阵，b是Julia中的稠密向量

Question

我有一个矩阵A，它是以CSC格式存储的稀疏矩阵，还有一个密集向量b。并且我想为x求解系统Ax = b。我该怎么做？

使用运算符和inv(A)运算符都不起作用。我还尝试了链接中建议的factorize功能。

有人建议另一种方法吗？使用Julia v1.7.1。此外，矩阵A不是三对角的。请查看链接中的问题以获取错误信息。

Answer 1

定义"不起作用"；。也许有一个最低限度的(非(工作示例。

例如，我刚刚尝试了这个，它似乎运行得很好：

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

如果没有更多关于你特定问题的信息，就没有什么好说的了。