[petsc-users] How to speed up geometric multigrid

[Matthew Knepley]

On Mon, Sep 23, 2013 at 8:55 AM, Michele Rosso <mrosso at uci.edu> wrote:

>  Hi,
>
> I am successfully using PETSc to solve a 3D Poisson's equation with CG +
> MG .  Such equation arises from a projection algorithm for a multiphase
> incompressible flow simulation.
> I set up the solver  as I was suggested to do in a previous thread(title: "GAMG speed") and run a test case (liquid droplet with surface
> tension falling under the effect of gravity in a quiescent fluid).
> The solution of the Poisson Equation via multigrid is correct but it
> becomes progressively slower and slower as the simulation progresses (I am
> performing successive solves) due to an increase in the number of
> iterations.
> Since the solution of the Poisson equation is mission-critical, I need to
> speed it up as much as I can.
> Could you please help me out with this?
>

First, what does the coefficient look like?

Second, what are the boundary conditions?

   Matt


> I run the test case with the following options:
>
> -pc_type mg  -pc_mg_galerkin  -pc_mg_levels 5   -mg_levels_ksp_type
> richardson -mg_levels_ksp_max_it 1
> -mg_coarse_pc_type lu   -mg_coarse_pc_factor_mat_solver_package
> superlu_dist
> -log_summary -ksp_view  -ksp_monitor_true_residual  -options_left
>
> Please find the diagnostic for the final solve in the attached file
> "final.txt'.
> Thank you,
>
> Michele
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130923/87fbc91a/attachment.html>