<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Hello,</p>
<p><br>
</p>
<p>I am running some code that employs gamg preconditioning within a fieldsplit, and for sufficiently large/refined meshes, I am getting the following error:</p>
<p><br>
</p>
<p>------------------------------------------------------------------------------------------------------------------------------------------------</p>
<p></p>
<div> Residual norms for fieldsplit_0_ solve.</div>
<div> 0 KSP unpreconditioned resid norm 1.019675281087e-08 true resid norm 1.019675281087e-08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div> 1 KSP unpreconditioned resid norm 3.855246547147e-09 true resid norm 3.855246547147e-09 ||r(i)||/||b|| 3.780857120550e-01</div>
<div> 2 KSP unpreconditioned resid norm 1.438241386184e-09 true resid norm 1.438241386184e-09 ||r(i)||/||b|| 1.410489606701e-01</div>
<div> 3 KSP unpreconditioned resid norm 3.624902894294e-10 true resid norm 3.624902894294e-10 ||r(i)||/||b|| 3.554958094531e-02</div>
<div> 4 KSP unpreconditioned resid norm 1.267419175485e-10 true resid norm 1.267419175485e-10 ||r(i)||/||b|| 1.242963518870e-02</div>
<div> 5 KSP unpreconditioned resid norm 2.929693449291e-11 true resid norm 2.929693449291e-11 ||r(i)||/||b|| 2.873163156576e-03</div>
<div> 6 KSP unpreconditioned resid norm 9.520263854387e-12 true resid norm 9.520263854423e-12 ||r(i)||/||b|| 9.336564326903e-04</div>
<div> 7 KSP unpreconditioned resid norm 1.679490979841e-12 true resid norm 1.679490979825e-12 ||r(i)||/||b|| 1.647084136466e-04</div>
<div> 8 KSP unpreconditioned resid norm 3.608932906029e-13 true resid norm 3.608932905928e-13 ||r(i)||/||b|| 3.539296257217e-05</div>
<div> 9 KSP unpreconditioned resid norm 9.297426160279e-14 true resid norm 9.297426159708e-14 ||r(i)||/||b|| 9.118026426799e-06</div>
<p></p>
<p><span style="font-size: 12pt;">[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</span><br>
</p>
<p></p>
<div>[0]PETSC ERROR: Petsc has generated inconsistent data</div>
<div>[0]PETSC ERROR: <b>!(matA_1 && !matA_1->compressedrow.use)</b></div>
<div>[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.</div>
<div>[0]PETSC ERROR: Petsc Release Version 3.7.3, Jul, 24, 2016 </div>
<div>[0]PETSC ERROR: ./main_3D on a x86_64 named artur-ubuntu by artur Mon Aug 1 22:08:45 2016</div>
<div>[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</div>
<div>[0]PETSC ERROR: #1 smoothAggs() line 354 in /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/pc/impls/gamg/agg.c</div>
<div>[0]PETSC ERROR: #2 PCGAMGCoarsen_AGG() line 998 in /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/pc/impls/gamg/agg.c</div>
<div>[0]PETSC ERROR: #3 PCSetUp_GAMG() line 571 in /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/pc/impls/gamg/gamg.c</div>
<div>[0]PETSC ERROR: #4 PCSetUp() line 968 in /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/pc/interface/precon.c</div>
<div>[0]PETSC ERROR: #5 KSPSetUp() line 390 in /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/ksp/interface/itfunc.c</div>
<div>[0]PETSC ERROR: #6 KSPSolve() line 599 in /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/ksp/interface/itfunc.c</div>
<div>[0]PETSC ERROR: #7 PCApply_FieldSplit() line 1016 in /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/pc/impls/fieldsplit/fieldsplit.c</div>
<div>[0]PETSC ERROR: #8 PCApply() line 482 in /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/pc/interface/precon.c</div>
<div>[0]PETSC ERROR: #9 KSP_PCApply() line 244 in /home/artur/Rorsrach/Packages/petsc-3.7.3/include/petsc/private/kspimpl.h</div>
<div>[0]PETSC ERROR: #10 KSPInitialResidual() line 69 in /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/ksp/interface/itres.c</div>
<div>[0]PETSC ERROR: #11 KSPSolve_GMRES() line 239 in /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/ksp/impls/gmres/gmres.c</div>
<div>[0]PETSC ERROR: #12 KSPSolve() line 656 in /home/artur/Rorsrach/Packages/petsc-3.7.3/src/ksp/ksp/interface/itfunc.c</div>
<div>[0]PETSC ERROR: #13 solve() line 765 in /home/artur/Desktop/Preconditioned/MI_3D/Cpp/MorseI_PML.cpp</div>
<div>
<p style="font-family: Calibri, Arial, Helvetica, sans-serif, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 16px;">
</p>
<div>terminate called after throwing an instance of 'std::runtime_error'</div>
<div> what(): Error detected in C PETSc</div>
<div>[artur-ubuntu:07250] *** Process received signal ***</div>
<div>[artur-ubuntu:07250] Signal: Aborted (6)</div>
<div>[artur-ubuntu:07250] Signal code: (-6)</div>
<div>[artur-ubuntu:07250] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x36cb0) [0x7fe413815cb0]</div>
<div>[artur-ubuntu:07250] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x37) [0x7fe413815c37]</div>
<div>[artur-ubuntu:07250] [ 2] /lib/x86_64-linux-gnu/libc.so.6(abort+0x148) [0x7fe413819028]</div>
<div>[artur-ubuntu:07250] [ 3] /usr/lib/x86_64-linux-gnu/libstdc++.so.6(_ZN9__gnu_cxx27__verbose_terminate_handlerEv+0x155) [0x7fe413f0a535]</div>
<div>[artur-ubuntu:07250] [ 4] /usr/lib/x86_64-linux-gnu/libstdc++.so.6(+0x5e6d6) [0x7fe413f086d6]</div>
<div>[artur-ubuntu:07250] [ 5] /usr/lib/x86_64-linux-gnu/libstdc++.so.6(+0x5e703) [0x7fe413f08703]</div>
<div>[artur-ubuntu:07250] [ 6] /usr/lib/x86_64-linux-gnu/libstdc++.so.6(+0x5e922) [0x7fe413f08922]</div>
<div>[artur-ubuntu:07250] [ 7] /home/artur/Rorsrach/Packages/petsc-3.7.3/x86_64/lib/libpetsc.so.3.7(+0x18d9ec) [0x7fe414c419ec]</div>
<div>[artur-ubuntu:07250] [ 8] /home/artur/Rorsrach/Packages/petsc-3.7.3/x86_64/lib/libpetsc.so.3.7(PetscError+0x45b) [0x7fe414c41e94]</div>
<div>[artur-ubuntu:07250] [ 9] ./main_3D(_ZN10MorseI_PMLILi3EE5solveEv+0x1cf0) [0x430d00]</div>
<div>[artur-ubuntu:07250] [10] ./main_3D(_ZN10MorseI_PMLILi3EE3runEv+0xd9) [0x435bd9]</div>
<div>[artur-ubuntu:07250] [11] ./main_3D(main+0x6c) [0x41a19c]</div>
<div>[artur-ubuntu:07250] [12] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7fe413800f45]</div>
<div>[artur-ubuntu:07250] [13] ./main_3D() [0x41a223]</div>
<div>[artur-ubuntu:07250] *** End of error message ***</div>
------------------------------------------------------------------------------------------------------------------------------------------------
<p></p>
<div><br>
</div>
<p style="font-family: Calibri, Arial, Helvetica, sans-serif, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 16px;">
</p>
<br>
</div>
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 <a href="http://www.mcs.anl.gov/petsc/petsc-3.7.3/src/ksp/pc/impls/gamg/agg.c.html" id="LPlnk343556" title="http://www.mcs.anl.gov/petsc/petsc-3.7.3/src/ksp/pc/impls/gamg/agg.c.html
Ctrl+Click or tap to follow the link">http://www.mcs.anl.gov/petsc/petsc-3.7.3/src/ksp/pc/impls/gamg/agg.c.html</a>
<p></p>
<p><br>
</p>
<p>Thanks,</p>
<p><br>
</p>
<p>Artur</p>
</div>
</body>
</html>