Preallocating Matrix

Andreas Grassl Andreas.Grassl at
Thu Apr 23 10:38:59 CDT 2009


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?



 /"\                               Grassl Andreas
 \ /    ASCII Ribbon Campaign      Uni Innsbruck Institut f. Mathematik
  X      against HTML email        Technikerstr. 13 Zi 709
 / \                               +43 (0)512 507 6091

More information about the petsc-users mailing list