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 <[email protected]> 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
