[petsc-users] The PETSC ERROR: VecMAXPY() when using GMRES

Song Gao song.gao2 at mail.mcgill.ca
Wed Mar 11 14:33:24 CDT 2015


Hello,

I'm solving Navier-Stokes equations by finite element method.  I use KSP as
the linear system solver. I run with 2 cpu. The code runs fine in non-debug
version. But when I switch to the debug version, the code gives the
following error.

I output the matrix and rhs before calling KSPSolve to make sure no NAN or
INF in them. The condition number of matrix is ~2e4.  Seems okay.

I also run the code with valgrind, but didn't find any other errors. The
valgrind output is attached. Any ideas of what I can do next? Thanks in
advance.


[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: Invalid argument!
[0]PETSC ERROR: Scalar value must be same on all processes, argument # 3!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 7, Sat May 11 22:15:24
CDT 2013
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR:
/home/cfd/sgao/mycodes/fensap_new_edge_coefficient/fensapng-mf-newmuscl-overledg_org/bin/fensapMPI_LINUX64_DEBUG
on a linux named anakin by sgao Wed Mar 11 15:07:53 2015
[0]PETSC ERROR: Libraries linked from /tmp/PETSC33/petsc-3.3-p7/linux/lib
[0]PETSC ERROR: Configure run at Wed Jan 15 12:04:54 2014
[0]PETSC ERROR: Configure options
--with-mpi-dir=/usr/local.linux64/lib64/MPI-openmpi-1.4.5/
--with-shared-libraries=0 --COPTFLAGS=-g --FOPTFLAGS=-g --with-debugging=yes
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: VecMAXPY() line 1186 in src/vec/vec/interface/rvector.c
[0]PETSC ERROR: KSPGMRESBuildSoln() line 345 in
src/ksp/ksp/impls/gmres/gmres.c
[0]PETSC ERROR: KSPGMRESCycle() line 206 in src/ksp/ksp/impls/gmres/gmres.c
[0]PETSC ERROR: KSPSolve_GMRES() line 231 in src/ksp/ksp/impls/gmres/gmres.c
[0]PETSC ERROR: KSPSolve() line 446 in src/ksp/ksp/interface/itfunc.c
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150311/eead4b97/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: valgrindout
Type: application/octet-stream
Size: 9775 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150311/eead4b97/attachment-0001.obj>


More information about the petsc-users mailing list