On 06/18/2013 03:05 PM, John Peterson wrote: > > Do you still run out of memory if you run without GAMG? > > > > There could be some GAMG options that control memory consumption, I > don't know too much about it. > > > > I am able to to solve the problem with -pc_type bjacobi and > -sub_pc_type icc, but that still uses a lot of memory, around 60Gb. > And also this required more than 5000 CG iterations which is why I > moved to AMG. > > > OK, switching from bjacobi -> GAMG caused memory consumption to more > than double?! I'd definitely look into the GAMG options... > > Another possibility is that you could build PETSc with Hypre and run > with: > > -pc_type hypre -pc_hypre_type boomeramg > -pc_hypre_boomeramg_strong_threshold 0.7
Thanks. This works, but the number of CG iterations still don't seem quite right (more than 1000 to get to rtol 1e-4). Do you know if I somehow have to tell hypre that I am solving a vector problem (I'm guessing otherwise it would not be able to exploit the elliptic structure)? Jens > > -- > John ------------------------------------------------------------------------------ This SF.net email is sponsored by Windows: Build for Windows Store. http://p.sf.net/sfu/windows-dev2dev _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
