[petsc-dev] [petsc-users] Bad memory scaling with PETSc 3.10

Mark Adams mfadams at lbl.gov
Thu Mar 21 18:37:18 CDT 2019

> Could you explain this more by adding some small examples?
Since you are considering implementing all-at-once (four nested loops,
right?) I'll give you my old code.

This code is hardwired for two AMG and for a geometric-AMG, where the
blocks of the R (and hence P) matrices are scaled identities and I only
store the scale. So you ignore those branches. This code also does
equivalent real form complex, so more branches to ignore.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20190321/1891005b/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: prom_mat_prod.C
Type: application/octet-stream
Size: 40192 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20190321/1891005b/attachment-0001.obj>

More information about the petsc-dev mailing list