标签: python scipy petsc
我似乎找不到如何有效加载scipy稀疏矩阵的方法,例如: csr_matrix,加入petsc4py矩阵,例如: PETSc.Mat().createAIJ。我找到了this thread,但我无法应用它。
csr_matrix
PETSc.Mat().createAIJ
我也很欣赏这个东西实际记录的指针。 demo目录中的示例只解释了一个部分,我看不到任何文档字符串。
demo
答案 0 :(得分:8)
你的链接说要在PETSc中创建一个稀疏矩阵,你应该使用这样的命令:
PETSc.Mat().createAIJ(size=(nrows,ncols), csr=(ai,aj,aa))
根据this,ai,aj和aa在PETSc中发言:
ai
aj
aa
> i - row indices > j - column indices > a - matrix values
这些分别与.indptr的{{1}},.indices和.data属性相同,有关详细信息,请参阅docs。
.indptr
.indices
.data
因此,如果您的链接正确,则以下内容应该有效:
scypy.sparse.csr_matrix
不幸的是,我自己无法测试。