[petsc-users] Meaning of error message (gamg & fieldsplit related)

Wed Aug 3 09:50:11 CDT 2016

Are you saying it works for a while but fails when the problem is
large, or that it never works with fieldsplit_1?

And how many processors are you using?

On Mon, Aug 1, 2016 at 9:13 PM, Safin, Artur <aks084000 at utdallas.edu> wrote:
> Hello,
>
>
> I am running some code that employs gamg preconditioning within a
> fieldsplit, and for sufficiently large/refined meshes, I am getting the
> following error:
>
>
> ------------------------------------------------------------------------------------------------------------------------------------------------
>
>     Residual norms for fieldsplit_0_ solve.
>     0 KSP unpreconditioned resid norm 1.019675281087e-08 true resid norm
> 1.019675281087e-08 ||r(i)||/||b|| 1.000000000000e+00
>     1 KSP unpreconditioned resid norm 3.855246547147e-09 true resid norm
> 3.855246547147e-09 ||r(i)||/||b|| 3.780857120550e-01
>     2 KSP unpreconditioned resid norm 1.438241386184e-09 true resid norm
> 1.438241386184e-09 ||r(i)||/||b|| 1.410489606701e-01
>     3 KSP unpreconditioned resid norm 3.624902894294e-10 true resid norm
> 3.624902894294e-10 ||r(i)||/||b|| 3.554958094531e-02
>     4 KSP unpreconditioned resid norm 1.267419175485e-10 true resid norm
> 1.267419175485e-10 ||r(i)||/||b|| 1.242963518870e-02
>     5 KSP unpreconditioned resid norm 2.929693449291e-11 true resid norm
> 2.929693449291e-11 ||r(i)||/||b|| 2.873163156576e-03
>     6 KSP unpreconditioned resid norm 9.520263854387e-12 true resid norm
> 9.520263854423e-12 ||r(i)||/||b|| 9.336564326903e-04
>     7 KSP unpreconditioned resid norm 1.679490979841e-12 true resid norm
> 1.679490979825e-12 ||r(i)||/||b|| 1.647084136466e-04
>     8 KSP unpreconditioned resid norm 3.608932906029e-13 true resid norm
> 3.608932905928e-13 ||r(i)||/||b|| 3.539296257217e-05
>     9 KSP unpreconditioned resid norm 9.297426160279e-14 true resid norm
> 9.297426159708e-14 ||r(i)||/||b|| 9.118026426799e-06
>
> [0]PETSC ERROR: --------------------- Error Message
> --------------------------------------------------------------
>
> [0]PETSC ERROR: Petsc has generated inconsistent data
> [0]PETSC ERROR: !(matA_1 && !matA_1->compressedrow.use)
> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for
> trouble shooting.
> [0]PETSC ERROR: Petsc Release Version 3.7.3, Jul, 24, 2016
> [0]PETSC ERROR: ./main_3D on a x86_64 named artur-ubuntu by artur Mon Aug  1
> 22:08:45 2016
> [0]PETSC ERROR: Configure options --with-scalar-type=complex --with-mpi=1
> --with-clanguage=c++ --with-cc=mpicc --with-fc=gfortran --with-cxx=mpic++
> --with-fc=mpif90 --download-mumps --download-scalapack
> [0]PETSC ERROR: #1 smoothAggs() line 354 in
> /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/pc/impls/gamg/agg.c
> [0]PETSC ERROR: #2 PCGAMGCoarsen_AGG() line 998 in
> /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/pc/impls/gamg/agg.c
> [0]PETSC ERROR: #3 PCSetUp_GAMG() line 571 in
> /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/pc/impls/gamg/gamg.c
> [0]PETSC ERROR: #4 PCSetUp() line 968 in
> /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/pc/interface/precon.c
> [0]PETSC ERROR: #5 KSPSetUp() line 390 in
> /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: #6 KSPSolve() line 599 in
> /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: #7 PCApply_FieldSplit() line 1016 in
> /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/pc/impls/fieldsplit/fieldsplit.c
> [0]PETSC ERROR: #8 PCApply() line 482 in
> /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/pc/interface/precon.c
> [0]PETSC ERROR: #9 KSP_PCApply() line 244 in
> /home/artur/Rorsrach/Packages/petsc-3.7.3/include/petsc/private/kspimpl.h
> [0]PETSC ERROR: #10 KSPInitialResidual() line 69 in
> /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/ksp/interface/itres.c
> [0]PETSC ERROR: #11 KSPSolve_GMRES() line 239 in
> /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/ksp/impls/gmres/gmres.c
> [0]PETSC ERROR: #12 KSPSolve() line 656 in
> /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: #13 solve() line 765 in
> /home/artur/Desktop/Preconditioned/MI_3D/Cpp/MorseI_PML.cpp
>
> terminate called after throwing an instance of 'std::runtime_error'
>   what():  Error detected in C PETSc
> [artur-ubuntu:07250] *** Process received signal ***
> [artur-ubuntu:07250] Signal: Aborted (6)
> [artur-ubuntu:07250] Signal code:  (-6)
> [artur-ubuntu:07250] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x36cb0)
> [0x7fe413815cb0]
> [artur-ubuntu:07250] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x37)
> [0x7fe413815c37]
> [artur-ubuntu:07250] [ 2] /lib/x86_64-linux-gnu/libc.so.6(abort+0x148)
> [0x7fe413819028]
> [artur-ubuntu:07250] [ 3]
> /usr/lib/x86_64-linux-gnu/libstdc++.so.6(_ZN9__gnu_cxx27__verbose_terminate_handlerEv+0x155)
> [0x7fe413f0a535]
> [artur-ubuntu:07250] [ 4] /usr/lib/x86_64-linux-gnu/libstdc++.so.6(+0x5e6d6)
> [0x7fe413f086d6]
> [artur-ubuntu:07250] [ 5] /usr/lib/x86_64-linux-gnu/libstdc++.so.6(+0x5e703)
> [0x7fe413f08703]
> [artur-ubuntu:07250] [ 6] /usr/lib/x86_64-linux-gnu/libstdc++.so.6(+0x5e922)
> [0x7fe413f08922]
> [artur-ubuntu:07250] [ 7]
> /home/artur/Rorsrach/Packages/petsc-3.7.3/x86_64/lib/libpetsc.so.3.7(+0x18d9ec)
> [0x7fe414c419ec]
> [artur-ubuntu:07250] [ 8]
> /home/artur/Rorsrach/Packages/petsc-3.7.3/x86_64/lib/libpetsc.so.3.7(PetscError+0x45b)
> [0x7fe414c41e94]
> [artur-ubuntu:07250] [ 9] ./main_3D(_ZN10MorseI_PMLILi3EE5solveEv+0x1cf0)
> [0x430d00]
> [artur-ubuntu:07250] [10] ./main_3D(_ZN10MorseI_PMLILi3EE3runEv+0xd9)
> [0x435bd9]
> [artur-ubuntu:07250] [11] ./main_3D(main+0x6c) [0x41a19c]
> [artur-ubuntu:07250] [12]
> /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7fe413800f45]
> [artur-ubuntu:07250] [13] ./main_3D() [0x41a223]
> [artur-ubuntu:07250] *** End of error message ***
> ------------------------------------------------------------------------------------------------------------------------------------------------
>
>
>
> The problem specifically appears when I attempt to precondition fieldsplit_1
> with gamg (no problems with gamg in fieldsplit_0 though for some reason). I
> am curious if someone can explain what this error actually means; this comes
> from line 354 in
> http://www.mcs.anl.gov/petsc/petsc-3.7.3/src/ksp/pc/impls/gamg/agg.c.html
>
>
> Thanks,
>
>
> Artur