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

Wed Aug 17 16:04:44 CDT 2016

> On Aug 17, 2016, at 11:02 AM, Safin, Artur <aks084000 at utdallas.edu> wrote:
> 
> Hi Barry,
> 
> Well, this is a bit strange: my run produces a completely different output.

    I am a complete idiot. I did not notice you were running with complex. I switched to complex and get the same error you got. 

    I've attached a patch that will remove the problem; it was some unneeded checks that did not belong there. This is also fixed in the master branch. I can't fix it in maint because of code changes since the maint release.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: gamg.patch
Type: application/octet-stream
Size: 918 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160817/7b4b0605/attachment.obj>
-------------- next part --------------


    Barry


> What I have is
> 
> --------------------------------------------------------------------------------------------------------------------------------------------------------
> artur at artur-ubuntu:~/Downloads/temp$ ./ex -options_file options.txt
> 
>    Residual norms for fieldsplit_PA_ solve.
>    0 KSP unpreconditioned resid norm 3.037294001981e-01 true resid norm 3.037294001981e-01 ||r(i)||/||b|| 1.000000000000e+00
>    1 KSP unpreconditioned resid norm 2.205402771720e-02 true resid norm 2.205402771720e-02 ||r(i)||/||b|| 7.261077690475e-02
>    2 KSP unpreconditioned resid norm 2.342871234277e-03 true resid norm 2.342871234279e-03 ||r(i)||/||b|| 7.713679455299e-03
>    3 KSP unpreconditioned resid norm 2.847545158015e-04 true resid norm 2.847545157994e-04 ||r(i)||/||b|| 9.375270079671e-04
>    4 KSP unpreconditioned resid norm 8.896004812140e-05 true resid norm 8.896004811976e-05 ||r(i)||/||b|| 2.928924498640e-04
>    5 KSP unpreconditioned resid norm 4.813451098741e-06 true resid norm 4.813451098816e-06 ||r(i)||/||b|| 1.584782736105e-05
>    6 KSP unpreconditioned resid norm 1.029168681786e-06 true resid norm 1.029168680690e-06 ||r(i)||/||b|| 3.388439446491e-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 Development GIT revision: v3.7.3-1165-gfeaa1dd  GIT Date: 2016-08-16 11:58:28 -0500
> [0]PETSC ERROR: ./ex on a x86_64 named artur-ubuntu by artur Wed Aug 17 10:47:27 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 336 in /home/artur/Rorsrach/Packages/petsc/src/ksp/pc/impls/gamg/agg.c
> [0]PETSC ERROR: #2 PCGAMGCoarsen_AGG() line 1004 in /home/artur/Rorsrach/Packages/petsc/src/ksp/pc/impls/gamg/agg.c
> [0]PETSC ERROR: #3 PCSetUp_GAMG() line 526 in /home/artur/Rorsrach/Packages/petsc/src/ksp/pc/impls/gamg/gamg.c
> [0]PETSC ERROR: #4 PCSetUp() line 968 in /home/artur/Rorsrach/Packages/petsc/src/ksp/pc/interface/precon.c
> [0]PETSC ERROR: #5 KSPSetUp() line 393 in /home/artur/Rorsrach/Packages/petsc/src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: #6 KSPSolve() line 602 in /home/artur/Rorsrach/Packages/petsc/src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: #7 PCApply_FieldSplit() line 1040 in /home/artur/Rorsrach/Packages/petsc/src/ksp/pc/impls/fieldsplit/fieldsplit.c
> [0]PETSC ERROR: #8 PCApply() line 482 in /home/artur/Rorsrach/Packages/petsc/src/ksp/pc/interface/precon.c
> [0]PETSC ERROR: #9 KSP_PCApply() line 245 in /home/artur/Rorsrach/Packages/petsc/include/petsc/private/kspimpl.h
> [0]PETSC ERROR: #10 KSPInitialResidual() line 69 in /home/artur/Rorsrach/Packages/petsc/src/ksp/ksp/interface/itres.c
> [0]PETSC ERROR: #11 KSPSolve_GMRES() line 239 in /home/artur/Rorsrach/Packages/petsc/src/ksp/ksp/impls/gmres/gmres.c
> [0]PETSC ERROR: #12 KSPSolve() line 659 in /home/artur/Rorsrach/Packages/petsc/src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: #13 main() line 67 in /home/artur/Downloads/temp/ex.c
> [0]PETSC ERROR: PETSc Option Table entries:
> [0]PETSC ERROR: -fieldsplit_PA_gamg_est_ksp_pmat_view
> [0]PETSC ERROR: -fieldsplit_PA_ksp_monitor_true_residual
> [0]PETSC ERROR: -fieldsplit_PA_ksp_type fgmres
> [0]PETSC ERROR: -fieldsplit_PA_pc_type gamg
> [0]PETSC ERROR: -fieldsplit_PB_ksp_monitor_true_residual
> [0]PETSC ERROR: -fieldsplit_PB_ksp_type fgmres
> [0]PETSC ERROR: -fieldsplit_PB_pc_type gamg
> [0]PETSC ERROR: -fieldsplit_TA_ksp_monitor_true_residual
> [0]PETSC ERROR: -fieldsplit_TA_ksp_type fgmres
> [0]PETSC ERROR: -fieldsplit_TA_pc_type gamg
> [0]PETSC ERROR: -fieldsplit_TB_ksp_monitor_true_residual
> [0]PETSC ERROR: -fieldsplit_TB_ksp_type fgmres
> [0]PETSC ERROR: -fieldsplit_TB_pc_type gamg
> [0]PETSC ERROR: -ksp_converged_reason
> [0]PETSC ERROR: -pc_fieldsplit_type multiplicative
> [0]PETSC ERROR: -pc_type fieldsplit
> [0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint at mcs.anl.gov----------
> --------------------------------------------------------------------------------------------------------------------------------------------------------
> 
> Do you have any idea as to why our results are completely different?
> 
> Artur