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