[petsc-users] Bad memory scaling with PETSc 3.10

Zhang, Hong hzhang at mcs.anl.gov
Fri Apr 12 11:18:50 CDT 2019

Myriam :
Thanks for your effort. It will help us improve PETSc.

Hi all,

I used the wrong script, that's why it diverged... Sorry about that.
I tried again with the right script applied on a tiny problem (~200
elements). I can see a small difference in memory usage (gain ~ 1mB).
when adding the -mat_freeintermediatestructures option. I still have to
execute larger cases to plot the scaling. The supercomputer I am used to
run my jobs on is really busy at the moment so it takes a while. I hope
I'll send you the results on Monday.

Thanks everyone,


Le 04/11/19 à 06:01, Jed Brown a écrit :
> "Zhang, Hong" <hzhang at mcs.anl.gov<mailto:hzhang at mcs.anl.gov>> writes:
>> Jed:
>>>> Myriam,
>>>> Thanks for the plot. '-mat_freeintermediatedatastructures' should not affect solution. It releases almost half of memory in C=PtAP if C is not reused.
>>> And yet if turning it on causes divergence, that would imply a bug.
>>> Hong, are you able to reproduce the experiment to see the memory
>>> scaling?
>> I like to test his code using an alcf machine, but my hands are full now. I'll try it as soon as I find time, hopefully next week.
> I have now compiled and run her code locally.
> Myriam, thanks for your last mail adding configuration and removing the
> MemManager.h dependency.  I ran with and without
> -mat_freeintermediatedatastructures and don't see a difference in
> convergence.  What commands did you run to observe that difference?

Myriam Peyrounette

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190412/5d524ad4/attachment.html>

More information about the petsc-users mailing list