Re: [petsc-dev] [petsc-users] performance regression with GAMG

2023-10-06 Thread Pierre Jolivet
Hello Stephan, Could you share the Amat/Pmat in binary format of the specific fieldsplit block, as well as all inputs needed to generate the same grid hierarchy (block size, options used, near kernel)? Alternatively, have you been able to generate the same error in a plain PETSc example? I’m

Re: [petsc-dev] [petsc-users] performance regression with GAMG

2023-10-06 Thread Pierre Jolivet
On Oct 6, 2023, at 3:48 AM, Mark Adams wrote:Pierre, (moved to dev)It looks like there is a subtle bug in the new MatFilter.My guess is that after the compression/filter the communication buffers and lists need to be recomputed because the graph has changed.Maybe an issue with

Re: [petsc-dev] [petsc-users] performance regression with GAMG

2023-10-06 Thread Mark Adams
Stephan Kramer has a reproducer (see this thread) with 1500 processors. They have lots of tests and this seemed to be the only, or the smallest, one that failed. I added a switch in GAMG to use your "low memory filter" to use this. Stephan is accommodating and would test a fix. Stephan has stack