朱莉娅:将CHOLMOD因子转换为稀疏矩阵，然后再转换回来

Question

我有一个稀疏矩阵H的CHOLMOD分解，我想编辑上、下和块对角线因子的稀疏表示。我该怎么做？当我运行以下内容时，最后一行不起作用。

H = sprand(10,10,0.5)
fac = ldltfact(H; shift=0.0)
fD = fac[:D]
D = Base.SparseArrays.CHOLMOD.Sparse(fD)

有没有办法从稀疏矩阵到CHOLMOD.factor相反的方向？

Answer 1

提取ldltfact的相关分解矩阵可能有点乏味。以下示例显示了一个类似于问题中的示例，其中包含提取的矩阵恢复原始分解矩阵的最终测试：

srand(1)
pre = sprand(10,10,0.5)
H = pre + pre' + speye(10,10)
fac = ldltfact(H; shift=0.0)
P = sparse(1:size(H,1),fac[:p],ones(size(H,1)))
LD = sparse(fac[:LD]) # this matrix contains both D and L embedded in it
L = copy(LD)
for i=1:size(L,1)
  L[i,i] = 1.0
end
D = sparse(1:size(L,1),1:size(L,1),diag(LD))
PHP = P*H*P'
LDL = L*D*L'
using Base.Test
@test PHP ≈ LDL

预期输出(以及 Julia v0.6.3 上的实际输出(：

julia> @test PHP ≈ LDL
Test Passed

希望这有帮助。