scipy稀疏矩阵作为petsc4py的输入

时间:2013-03-15 21:47:45

标签: python scipy petsc

我似乎找不到如何有效加载scipy稀疏矩阵的方法,例如: csr_matrix,加入petsc4py矩阵,例如: PETSc.Mat().createAIJ。我找到了this thread,但我无法应用它。

我也很欣赏这个东西实际记录的指针。 demo目录中的示例只解释了一个部分,我看不到任何文档字符串。

1 个答案:

答案 0 :(得分:8)

你的链接说要在PETSc中创建一个稀疏矩阵,你应该使用这样的命令:

PETSc.Mat().createAIJ(size=(nrows,ncols), csr=(ai,aj,aa))

根据thisaiajaa在PETSc中发言:

> i - row indices
> j - column indices
> a - matrix values

这些分别与.indptr的{​​{1}},.indices.data属性相同,有关详细信息,请参阅docs

因此,如果您的链接正确,则以下内容应该有效:

scypy.sparse.csr_matrix

不幸的是,我自己无法测试。