[petsc-users] Question on matrix assembly

Barry Smith bsmith at petsc.dev
Sat Aug 15 21:38:19 CDT 2020


    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.


> 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>

More information about the petsc-users mailing list