我想将一个小的cholmod_sparse
矩阵插入一个更大的矩阵。是否有一种优雅的方法可以做到这一点,或者我必须通过循环以及直接对内部列压缩表示形式进行操作来解决这个问题?
我想在CHOLMOD中做些什么,我可以用Matlab
伪代码轻松地说明一下:
n = 24;
L = sparse(n,n);
blocksize = 4;
for i = 1:blocksize:n
... % compute some sparse matrix S
L(i:end, i:i+blocksize) = S; % insert S into L
end
这个Matlab示例创建的稀疏矩阵S在每次迭代中变得越来越小,并插入到L的下部三角形部分。我想为此创建一个CHOLMOD(C ++)版本。