[petsc-dev] PC ML Setup Error

[Matthew Knepley]

Is this the bad interaction Jed found between PCMG and ML?

   Matt

On Wed, Jun 29, 2011 at 12:45 PM, <domenico.borzacchiello at univ-st-etienne.fr
> wrote:

> Hi,
> I use ML as a preconditioner for the pressure poisson equation for Stokes
> problem. I solve for the pressure at each time step with ad DMMGSolve on
> the whole (Stokes) system. I randomly get the following error at some
> point.
>
> [0]PETSC ERROR: --------------------- Error Message
> ------------------------------------
> [0]PETSC ERROR: Operation done in wrong order!
> [0]PETSC ERROR: Number levels already set for MG
>  make sure that you call PCMGSetLevels() before KSPSetFromOptions()!
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: Petsc Development HG revision:
> d1471993ab06fb6a9072f156c420ef0a456a086c  HG Date: Sun May 01 09:59:41
> 2011 -0500
> [0]PETSC ERROR: See docs/changes/index.html for recent updates.
> [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
> [0]PETSC ERROR: See docs/index.html for manual pages.
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: ./2dbphlow on a arch-linu named
> domenico-Dell-System-XPS-L502X by domenico Wed Jun 29 14:32:00 2011
> [0]PETSC ERROR: Libraries linked from
> /opt/petsc-dev/arch-linux2-c-debug/lib
> [0]PETSC ERROR: Configure run at Thu May  5 03:28:56 2011
> [0]PETSC ERROR: Configure options --with-cc=icc --with-fc=ifort
> --download-f-blas-lapack=1 --download-mpich=1 --download-superlu_dist
> --download-parmetis --download-mumps --download-blacs --download-scalapack
> --download-hypre --download-ml --with-c2html=0 --download-prometheus
> --download-parms
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: PCMGSetLevels() line 151 in src/ksp/pc/impls/mg/mg.c
> [0]PETSC ERROR: PCSetUp_ML() line 610 in src/ksp/pc/impls/ml/ml.c
> [0]PETSC ERROR: PCSetUp() line 819 in src/ksp/pc/interface/precon.c
> [0]PETSC ERROR: KSPSetUp() line 261 in src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: KSPSolve() line 383 in src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: PCApply_LSC() line 85 in src/ksp/pc/impls/lsc/lsc.c
> [0]PETSC ERROR: PCApply() line 383 in src/ksp/pc/interface/precon.c
> [0]PETSC ERROR: FGMREScycle() line 174 in
> src/ksp/ksp/impls/gmres/fgmres/fgmres.c
> [0]PETSC ERROR: KSPSolve_FGMRES() line 300 in
> src/ksp/ksp/impls/gmres/fgmres/fgmres.c
> [0]PETSC ERROR: KSPSolve() line 426 in src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: PCApply_FieldSplit_Schur() line 572 in
> src/ksp/pc/impls/fieldsplit/fieldsplit.c
> [0]PETSC ERROR: PCApply() line 383 in src/ksp/pc/interface/precon.c
> [0]PETSC ERROR: FGMREScycle() line 174 in
> src/ksp/ksp/impls/gmres/fgmres/fgmres.c
> [0]PETSC ERROR: KSPSolve_FGMRES() line 300 in
> src/ksp/ksp/impls/gmres/fgmres/fgmres.c
> [0]PETSC ERROR: KSPSolve() line 426 in src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: SNES_KSPSolve() line 3124 in src/snes/interface/snes.c
> [0]PETSC ERROR: SNESSolve_LS() line 190 in src/snes/impls/ls/ls.c
> [0]PETSC ERROR: SNESSolve() line 2424 in src/snes/interface/snes.c
> [0]PETSC ERROR: DMMGSolveSNES() line 503 in src/snes/utils/damgsnes.c
> [0]PETSC ERROR: DMMGSolve() line 313 in src/snes/utils/damg.c
> [0]PETSC ERROR: main() line 726 in 2dbphlow.c
> application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
> [cli_0]: aborting job:
> application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
>
>
> =====================================================================================
> =   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
> =   EXIT CODE: 256
> =   CLEANING UP REMAINING PROCESSES
> =   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
>
> =====================================================================================
>
>
> I don't know the reason for it, it occurs at some random timestep if I
> change any of the input parameters like the viscosity, the geometry, the
> mesh size etc... Sometimes it just doesn't occur at all.
>
> ML is used with the following options:
>
>
> -stokes_fieldsplit_pressure_lsc_pc_type ml \
>  -stokes_fieldsplit_pressure_lsc_pc_ml_maxCoarseSize 1 \
>  -stokes_fieldsplit_pressure_lsc_pc_ml_Threshold 0.009 \
>  -stokes_fieldsplit_pressure_lsc_pc_ml_Symmetrize\
>  -stokes_fieldsplit_pressure_lsc_pc_ml_KeepAggInfo\
>  -stokes_fieldsplit_pressure_lsc_pc_ml_Reusable\
>  -stokes_fieldsplit_pressure_lsc_pc_ml_DampingFactor 1.3333 \
>
> Any idea?
>
>
>
>
>


-- 
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-dev/attachments/20110629/e09441c2/attachment.html>