OpenFOAM、PETSc或其他稀疏矩阵乘法源代码

Question

有人能告诉我，在哪里可以找到由 OpenFOAM、PETSc 或类似的东西实现的矩阵乘法的源代码吗？这不可能是一个简单的算法。 我找到了 OpenFOAM 和 PETSc 的主页，但在文档中我找不到乘法方法和源代码。

Answer 1

PETSc 实现了多种格式的矩阵乘法，请查看 MatMult_SeqAIJ 用于最基本的实现。对于以压缩稀疏行形式存储的稀疏矩阵，行开头为 ai，列索引为 aj，条目为 aa，乘法由以下简单内容组成核心。

for (i=0; i<m; i++) {
  y[i] = 0;
  for (j=ai[i]; j<ai[i+1]; j++)
    y[i] += aa[j] * x[aj[j]];
}