在 MEX 中创建稀疏矩阵

Question

如何在用 C 编写的 MEX 文件中创建二维稀疏矩阵。创建矩阵后如何像 C 中那样单独访问元素，比如 mat[i][j]？

我厌倦了使用 mxCreateNumericArray 函数，但我不是无法访问元素并将其设置为稀疏矩阵。

请帮忙

Answer 1

查看 mxCreateSparse 上的此页面。 然后您将需要查看 < a href="http://www.mathworks.com/help/matlab/apiref/mxsetpr.html" rel="noreferrer">mxSetPr, mxSetIr 和 mxSetJc 和相应的“get”版本。

这是一个示例如何分配稀疏矩阵。我意识到这是一个旧链接，但据我所知，它没有改变。

基本上，它的工作原理是ir 数据包含行索引。 jr 数据包含 ir 数组的索引列表。例如，在如何分配稀疏矩阵的链接中，代码：

...
static double  static_pr_data[NZMAX] = {5.8, 6.2, 5.9, 6.1};
static int     static_ir_data[NZMAX] = {0, 2, 1, 3};
static int     static_jc_data[COLS+1] = {0, 2, 4};
...

数组 static_jc_data 告诉您索引 static_jc_data[c] 到 static_jc_data[c+ static_pr_data和static_ir_data的1]-1对应于矩阵的c列。在该范围内（static_jc_data[c] 到 static_jc_data[c+1]-1），static_pr_data 的条目为您提供矩阵中的值static_ir_data 为您提供正确的行。

例如，这里的矩阵为：

A = [ 5.8  0
      0    5.9
      6.2  0
      0    6.1];

要回答有关如何单独访问元素的问题，您必须搜索第 i,j 个元素是否存在，如果存在则返回该元素，否则返回 0为此，您需要从 static_ir_data[static_jc_data[j]] 到 static_ir_data[static_jc_data[j+1]-1] 进行搜索，以查看您的 我存在。如果是，则 static_pr_data 中的相应条目将包含您的条目。如果不存在，则返回 0。

但是，通常在使用稀疏矩阵时，如果您要对矩阵进行大量搜索以查看某个元素是否存在，您可能需要考虑如何使用它。通常，执行您正在执行的任何操作时，最好只遍历一次非零元素，而不是搜索每个第 i,j 个条目。

哦，还有最后一件事。请记住，在 MEX 代码中，所有索引都是基于 0 的，但在 MATLAB 中它们都是基于 1 的。这应该会增加乐趣。