Hello,

I'm trying to solve a linear elasticity problem based on a tetraheadral
mesh with over a million degrees of freedom.

I've been using the LU solver in MUMPS but I'm now forced to switch to an
iterative solver because of the memory requirements of that solver.

CG and BoomerAMG in hypre seem to be a good combination, but I've not been
able to achieve convergence. The residual never goes below 1 even after
several thousand iterations and shows high oscillations.

The parameters I'm using are
-ksp_type cg
-ksp_monitor
-ksp_converged_reason
-pc_type hypre
-pc_hypre_type boomeramg
-pc_hypre_boomeramg_strong_threshold 0.7
-log_summary
-eps_monitor
-st_type sinvert

Does anyone have recommendations / prior experience on successfully using
BoomerAMG in libMesh?

Thanks,
Dana
------------------------------------------------------------------------------
WatchGuard Dimension instantly turns raw network data into actionable 
security intelligence. It gives you real-time visual feedback on key
security issues and trends.  Skip the complicated setup - simply import
a virtual appliance and go from zero to informed in seconds.
http://pubads.g.doubleclick.net/gampad/clk?id=123612991&iu=/4140/ostg.clktrk
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to