petsc4py

Question

için bir girdi olarak scipy seyrek matrisleri Scipy seyrek matrislerini verimli bir şekilde nasıl yükleyeceğinizin bir yolunu bulmak gibi görünmüyor, ör. csr_matrix, bir petsc4py matrisine, örn. PETSc.Mat().createAIJ. this thread'u buldum, ancak uygulayamıyorum.

Ayrıca bu belgelerin belgelendiği bir işaretçiyi de takdir ediyorum. demo dizindeki örnekler yalnızca bir bölümü açıklar ve herhangi bir belge dizesi göremiyorum. PETSc-speak

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

göre this için, ai, aj ve aa vardır:

Answer 1

Bağlantınız bir böyle komutu kullanmalısınız, PETSc bir seyrek matris oluşturmak için söylüyor

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

Bunlar scypy.sparse.csr_matrix ait .indptr, .indices ve .data özelliklerine sırasıyla eşdeğerdir,bakınız Detaylar için. Bağlantı doğru ise

Yani, şu çalışmalıdır:

>>> from petsc4py import PETSc 
>>> import scipy.sparse 
>>> csr_mat = scipy.sparse.rand(1000, 1000, density=0.001, format='csr') 
>>> petsc_mat = PETSc.Mat().createAIJ(size=csr_mat.shape, 
...         csr=(csr_mat.indptr, csr_mat.indices, 
...          csr_mat.data)) 

Maalesef kendim test edemezsiniz.