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

[Matthew Knepley]

On Wed, Mar 11, 2015 at 2:33 PM, Song Gao <song.gao2 at mail.mcgill.ca> wrote:

> 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?
>

Is there any chance you could spawn the debugger, -start_in_debugger, and
when you get the error,
print out the value of 'alpha' on each process?

Otherwise the best thing to do is output your Mat and RHS in binary and
send them so we can try to reproduce.

    Matt

  Thanks,

     Matt


> 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
>
>


-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150311/474ff375/attachment.html>