[petsc-users] KSP no longer returning NaNs?

Kirk, Benjamin (JSC-EG311) benjamin.kirk-1 at nasa.gov
Wed Jul 25 17:59:58 CDT 2012

On 7/25/12 5:04 PM, "Jed Brown" <jedbrown at mcs.anl.gov> wrote:

> This is probably the main complaint:
> http://petsc.cs.iit.edu/petsc/petsc-dev/rev/7d6f5cbe67bc

Thanks, Jed.  The trace is below.  Indeed, a call to VecNorm exits with an error, where previously it would have returned the NaN.   So this may be in user code before handing b off to the KSP object.

Is there an alternative, recommended way to ask a PETSc vector if it has NaNs without causing an abort?

MPI: (gdb) #0  0x000000310240effe in waitpid () from /lib64/libpthread.so.0
MPI: #1  0x00002aaaaafaa3fc in mpi_sgi_system (header=<value optimized out>) at sig.c:89
MPI: #2  MPI_SGI_stacktraceback (header=<value optimized out>) at sig.c:272
MPI: #3  0x00002aaaaafaae8e in first_arriver_handler (signo=6, stack_trace_sem=0x2aaab2fc0500) at sig.c:415
MPI: #4  0x00002aaaaafab000 in slave_sig_handler (signo=6, siginfo=<value optimized out>, extra=<value optimized out>) at sig.c:494
MPI: #5  <signal handler called>
MPI: #6  0x0000003101c32885 in raise () from /lib64/libc.so.6
MPI: #7  0x0000003101c34065 in abort () from /lib64/libc.so.6
MPI: #8  0x00002aaaadc1b244 in PetscTraceBackErrorHandler () from /software/x86_64/petsc/3.2-p5/aerolab_hpc_l1-mpt-2.05-intel-12.1/lib/libpetsc.so
MPI: #9  0x00002aaaadc1ac83 in PetscError () from /software/x86_64/petsc/3.2-p5/aerolab_hpc_l1-mpt-2.05-intel-12.1/lib/libpetsc.so
MPI: #10 0x00002aaaadc7ff3f in VecNorm () from /software/x86_64/petsc/3.2-p5/aerolab_hpc_l1-mpt-2.05-intel-12.1/lib/libpetsc.so
MPI: #11 0x00002aaaad0afb57 in libMesh::PetscVector<double>::l2_norm() const () from /lustre/work/benkirk/codes/install/aerolab_hpc_l1/lib/x86_64-unknown-linux-gnu_opt/libmesh.so.0
MPI: #12 0x00002aaaad15420b in libMesh::NewtonSolver::solve() () from /lustre/work/benkirk/codes/install/aerolab_hpc_l1/lib/x86_64-unknown-linux-gnu_opt/libmesh.so.0
MPI: #13 0x00002aaaad19e72d in libMesh::NonlinearImplicitSystem::solve() () from /lustre/work/benkirk/codes/install/aerolab_hpc_l1/lib/x86_64-unknown-linux-gnu_opt/libmesh.so.0
MPI: #14 0x00002aaaab7818b7 in FINS::Physics::run(std::basic_string<char, std::char_traits<char>, std::allocator<char> >) () from /lustre/work/benkirk/codes/install/aerolab_hpc_l1/lib/libfins.so.0

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120725/7d47bca4/attachment.html>

More information about the petsc-users mailing list