<HTML>
<HEAD>
<TITLE>Re: [petsc-users] KSP no longer returning NaNs?</TITLE>
</HEAD>
<BODY>
<FONT FACE="Lucida Grande"><SPAN STYLE='font-size:10pt'><BR>
On 7/25/12 5:04 PM, "Jed Brown" <<a href="jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>> wrote:<BR>
<BR>
<FONT COLOR="#0000FF">> <BR>
> This is probably the main complaint:<BR>
> <BR>
> <a href="http://petsc.cs.iit.edu/petsc/petsc-dev/rev/7d6f5cbe67bc">http://petsc.cs.iit.edu/petsc/petsc-dev/rev/7d6f5cbe67bc</a><BR>
</FONT><BR>
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.<BR>
<BR>
Is there an alternative, recommended way to ask a PETSc vector if it has NaNs without causing an abort?<BR>
<BR>
<BR>
<BR>
MPI: (gdb) #0 0x000000310240effe in waitpid () from /lib64/libpthread.so.0<BR>
MPI: #1 0x00002aaaaafaa3fc in mpi_sgi_system (header=<value optimized out>) at sig.c:89<BR>
MPI: #2 MPI_SGI_stacktraceback (header=<value optimized out>) at sig.c:272<BR>
MPI: #3 0x00002aaaaafaae8e in first_arriver_handler (signo=6, stack_trace_sem=0x2aaab2fc0500) at sig.c:415<BR>
MPI: #4 0x00002aaaaafab000 in slave_sig_handler (signo=6, siginfo=<value optimized out>, extra=<value optimized out>) at sig.c:494<BR>
MPI: #5 <signal handler called><BR>
MPI: #6 0x0000003101c32885 in raise () from /lib64/libc.so.6<BR>
MPI: #7 0x0000003101c34065 in abort () from /lib64/libc.so.6<BR>
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<BR>
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<BR>
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<BR>
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<BR>
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<BR>
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<BR>
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<BR>
<BR>
</SPAN></FONT>
</BODY>
</HTML>