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

Wed Aug 3 12:12:20 CDT 2016

   This should be fixed in the master branch 

> On Aug 1, 2016, at 11: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