[petsc-users] Question on matrix assembly

[Barry Smith]

   Sajid,

    In the branch barry/2020-08-13/matsetvalues-nopreallocation-seqaij I have provided a new function
MatSeqAIJSetTotalPreallocation() that requires only a single estimate for the entire number of nonzeros
in the matrix (not row by row). If you provide the rows in order with sorted column indices for each 
row the assembly process will be super fast and you won't need to compute anything twice.

   Barry


> On Aug 14, 2020, at 2:59 PM, Sajid Ali <sajidsyed2021 at u.northwestern.edu> wrote:
> 
> Hi Barry, 
> 
> All entries for a row are available together, but there is no requirement to compute them in order.
> 
> Thank You, 
> Sajid Ali | PhD Candidate
> Applied Physics
> Northwestern University
> s-sajid-ali.github.io <http://s-sajid-ali.github.io/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20200815/15a42364/attachment.html>