Hello, I'm assembling large matrices giving just the numbers of zero per row and wondering if it is possible to extract the nonzero-structure in array-format it can be fed again into
MatSeqAIJSetPreallocation(Mat B,PetscInt nz,const PetscInt nnz[]) to detect the bottleneck? cheers ando -- /"\ Grassl Andreas \ / ASCII Ribbon Campaign Uni Innsbruck Institut f. Mathematik X against HTML email Technikerstr. 13 Zi 709 / \ +43 (0)512 507 6091
