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

[Mark Adams]

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