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

Safin, Artur aks084000 at utdallas.edu
Mon Aug 1 23:13:32 CDT 2016


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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160802/66af7fb8/attachment.html>


More information about the petsc-users mailing list