[petsc-users] GMRES error

Peter Brune prbrune at gmail.com
Fri Mar 28 08:10:47 CDT 2014


That's what happens when you try to VecAXPY/MAXPY/whatever with a NaN or
inf.  We should really have a better error message for this.

- Peter


On Fri, Mar 28, 2014 at 7:50 AM, Mark Adams <mfadams at lbl.gov> wrote:

> I get this error in the middle of a GMRES iteration.  Any ideas?
>
>   0 SNES Function norm 6.067786990073e+17
>
>
>     Residual norms for fsa_ solve.
>
>
>     0 KSP Residual norm 6.067786990073e+17
>
>
>     1 KSP Residual norm 6.037192129591e+17
>
>
>     2 KSP Residual norm 5.852390318109e+17
>
>
>     3 KSP Residual norm 5.562063100790e+17
>
>
>     4 KSP Residual norm 5.284787094948e+17
>
>
>     5 KSP Residual norm 5.183753281867e+17
>
>
>     6 KSP Residual norm 4.813052731918e+17
>
>
>     7 KSP Residual norm 4.571658388437e+17
>
>
>     8 KSP Residual norm 4.326765203018e+17
>
>
>     9 KSP Residual norm 4.044786063878e+17
>
>
>    10 KSP Residual norm 3.852338418686e+17
>
>
>    11 KSP Residual norm 3.533696733608e+17
>
>
>    12 KSP Residual norm 3.242503119640e+17
>
>
>    13 KSP Residual norm 2.946010215248e+17
>
>
>    14 KSP Residual norm 2.642699328492e+17
>
>
>    15 KSP Residual norm 2.365537803535e+17
>
>
>    16 KSP Residual norm 1.841438363878e+17
>
>
>    17 KSP Residual norm 1.333683638359e+17
>
>
>    18 KSP Residual norm 7.325375319129e+16
>
>
>    19 KSP Residual norm 2.964254455676e+16
>
>
>    20 KSP Residual norm 1.216650547542e+16                    [3]PETSC
> ERROR: --------------------- Error Message
> ------------------------------------
>
> [3]PETSC ERROR: Invalid argument!
>
>
> [3]PETSC ERROR: Scalar value must be same on all processes, argument # 3!
>
>
> [3]PETSC ERROR:
> ------------------------------------------------------------------------
>
>
> [3]PETSC ERROR: Petsc Development GIT revision:
> a9e422f38a085385bc6ffdad704c9d6a22342ac9  GIT Date: 2013-12-31 01:26:24
> -0700
> [3]PETSC ERROR: See docs/changes/index.html for recent updates.
>
>
> [3]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
>
>
> [3]PETSC ERROR: See docs/index.html for manual pages.
>
>
> [3]PETSC ERROR:
> ------------------------------------------------------------------------
>
>
> [3]PETSC ERROR: ../../XGC1_3/xgc2 on a arch-xc30-dbg named nid04306 by
> madams Fri Mar 28 05:34:54 2014
>
> [3]PETSC ERROR: Libraries linked from
> /global/homes/m/madams/petsc_private/arch-xc30-dbg/lib
>
> [3]PETSC ERROR: Configure run at Tue Dec 31 14:53:42 2013
>
>
> [3]PETSC ERROR: Configure options --COPTFLAGS="-g -no-ipo"
> --CXXOPTFLAGS="-g -no-ipo" --FOPTFLAGS="-g -no-ipo" --download-hypre
> --download-metis --download-parmetis --with-cc=cc --with-clib-autodetect=0
> --with-cxx=CC --with-cxxlib-autodetect=0 --with-debugging=1 --with-fc=ftn
> --with-fortranlib-autodetect=0
> --with-hdf5-dir=/opt/cray/hdf5-parallel/1.8.9/intel/120/
> --with-shared-libraries=0 --with-x=0 LIBS=-lstdc++ PETSC_ARCH=arch-xc30-dbg
>
> [3]PETSC ERROR:
> ------------------------------------------------------------------------
>
>
> [3]PETSC ERROR: VecMAXPY() line 1257 in
> /global/u2/m/madams/petsc_private/src/vec/vec/interface/rvector.c
>
> [3]PETSC ERROR: KSPGMRESBuildSoln() line 353 in
> /global/u2/m/madams/petsc_private/src/ksp/ksp/impls/gmres/gmres.c
>
> [3]PETSC ERROR: KSPGMRESCycle() line 211 in
> /global/u2/m/madams/petsc_private/src/ksp/ksp/impls/gmres/gmres.c
>
> [3]PETSC ERROR: KSPSolve_GMRES() line 235 in
> /global/u2/m/madams/petsc_private/src/ksp/ksp/impls/gmres/gmres.c
>
> [3]PETSC ERROR: KSPSolve() line 432 in
> /global/u2/m/madams/petsc_private/src/ksp/ksp/interface/itfunc.c
>
> [3]PETSC ERROR: SNESSolve_KSPONLY() line 44 in
> /global/u2/m/madams/petsc_private/src/snes/impls/ksponly/ksponly.c
>
> [3]PETSC ERROR: SNESSolve() line 3812 in
> /global/u2/m/madams/petsc_private/src/snes/interface/snes.c
>
> Rank 3 [Fri Mar 28 05:35:04 2014] [c6-2c1s4n2] application called
> MPI_Abort(MPI_COMM_WORLD, 62) - process 3
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140328/4d9de35f/attachment.html>


More information about the petsc-users mailing list