在Python中加载matlab稀疏矩阵（matlab v 7.3）

Question

我有一个存储大型matlab稀疏矩阵（Matlab7.3）的数据文件，需要在python程序中使用。我使用h5py加载这个稀疏矩阵，发现有3个数据结构与稀疏矩阵相关。

假设稀疏矩阵的名称为M，3个数据结构为M['data']、M['ir']、M['jc']。最初我认为M['ir']和M['jc']存储非零项的行索引和列索引，但我刚刚发现M['jc']中存在一些大于稀疏矩阵行数的值。谁能解释一下3数据结构中存储了哪些类型的信息？

Answer 1

ir正如您所猜到的，是非空行的row-index。对于列索引来说，事情有些复杂，但在Mathworks mex-Function留档中有充分的记录。

粘贴自http://www.mathworks.de/de/help/matlab/apiref/mxsetir.html:

If the jth column of the sparse mxArray has any nonzero elements:

jc[j] is the index in ir, pr, and pi (if it exists) of the first nonzero element in the jth column.
jc[j+1]-1 is the index of the last nonzero element in the jth column.
For the jth column of the sparse matrix, jc[j] is the total number of nonzero elements in all preceding columns.
The number of nonzero elements in the jth column of the sparse mxArray is:

jc[j+1] - jc[j];

还要检查mxSetIr上的文档。假设您也可以访问matlab，您可能应该检查从文档链接的mex示例。