从 MATLAB 稀疏矩阵中提取压缩列存储向量

Question

我在 MATLAB 中有以下矩阵。

A =

     1     3     0     2     0
     0     1     0     0     2
     0     0     4     0     0
     1     0     0     0     4
     0     2     0     2     0

我可以使用sparse(A) 命令获得三元组形式，但是是否有内置方法可以从 MATLAB 中获取 column_pointer、row_index 和值，或者我应该编写代码吗？ （我在 C++ 中有一个工作代码，可以从三元组转换为压缩列存储，我想我必须将其移植到 .m 文件中）

MATLAB 的 find() 仍然给出三元组形式。但我想要压缩列存储形式，如下图所示。 （MATLAB 是基于 1 个索引的。所以我认为 col_ptr 和 row_index 都可能增加 1）

col_ptr =   [0 2 5 6 8 10];
row_index = [ 0 3 0 1 4 2 0 4 1 3]
values  =   [ 1 1 3 1 2 4 2 2 2 4]

Answer 1

我认为这样做：

col_ptr = [0 cumsum(sum(A~=0,1))];
[row_index, ~, values] = find(A);
row_index = row_index.'-1;
values = values.';