我想为稀疏矩阵重新创建一个求解函数（求解 Ax = b for x）。在 Julia 文档中，它说当我们将稀疏矩阵应用于 lufact() 时，它返回以下内容：
L, U, p, q, Rs = F[:(:)]
使用 Julia doc 中给定的公式：LU = Rs.*A[p,q]，我做了一些代数并得到以下公式：
x = U \ ( L \ (Rs.*b[p]) )
ipermute!(x,q)
当矩阵密集时，此公式与 Julia 中的默认 F\b 求解器匹配，但当矩阵稀疏时，结果关闭。有谁知道为什么？