<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Oct 14, 2014 at 2:07 PM, Sharp Stone <span dir="ltr"><<a href="mailto:thronesf@gmail.com" target="_blank">thronesf@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear all,<div><br></div><div>I'm trying to solve a linear system with petsc ksp solver. But after I called the kspsolve routine, I always got the "Vec entry at local location 0 is not-a-number or infinite at end of function" errors below. I checked the matrix and the vectors, and they are all numbers (and zeros). Is this a bug or I miss something? Thank you very much in advance!</div></div></blockquote><div><br></div><div>Start simpler. Just use -pc_type lu -ksp_type gmres and see if you are getting a believable solution. Then back off on the</div><div>PC and KSP.</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>
<p>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</p><p>[0]PETSC ERROR: Floating point exception</p><p>[0]PETSC ERROR: Vec entry at local location 0 is not-a-number or infinite at end of function: Parameter number 3</p><p>[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.</p><p>[0]PETSC ERROR: Petsc Release Version 3.5.2, Sep, 08, 2014 </p><p>[0]PETSC ERROR: ./pStreamer on a arch-darwin-cxx-debug named <a href="http://host206-47.physci.fit.edu" target="_blank">host206-47.physci.fit.edu</a> by GPL Tue Oct 14 14:58:05 2014</p><p>[0]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ --with-fc=gfortran --with-clanguage=cxx --download-fblaslapack --download-mpich</p><p>[0]PETSC ERROR: #12 VecValidValues() line 34 in /Users/GPL/local/PETSc/petsc-3.5.2/src/vec/vec/interface/rvector.c</p><p>[0]PETSC ERROR: #13 PCApply() line 442 in /Users/GPL/local/PETSc/petsc-3.5.2/src/ksp/pc/interface/precon.c</p><p>[0]PETSC ERROR: #14 KSP_PCApply() line 230 in /Users/GPL/local/PETSc/petsc-3.5.2/include/petsc-private/kspimpl.h</p><p>[0]PETSC ERROR: #15 KSPSolve_PREONLY() line 26 in /Users/GPL/local/PETSc/petsc-3.5.2/src/ksp/ksp/impls/preonly/preonly.c</p><p>[0]PETSC ERROR: #16 KSPSolve() line 459 in /Users/GPL/local/PETSc/petsc-3.5.2/src/ksp/ksp/interface/itfunc.c</p><p>[0]PETSC ERROR: #17 PCApply_BJacobi_Singleblock() line 675 in /Users/GPL/local/PETSc/petsc-3.5.2/src/ksp/pc/impls/bjacobi/bjacobi.c</p><p>[0]PETSC ERROR: #18 PCApply() line 440 in /Users/GPL/local/PETSc/petsc-3.5.2/src/ksp/pc/interface/precon.c</p><p>[0]PETSC ERROR: #19 KSP_PCApply() line 230 in /Users/GPL/local/PETSc/petsc-3.5.2/include/petsc-private/kspimpl.h</p><p>[0]PETSC ERROR: #20 KSPInitialResidual() line 56 in /Users/GPL/local/PETSc/petsc-3.5.2/src/ksp/ksp/interface/itres.c</p><p>[0]PETSC ERROR: #21 KSPSolve_BCGS() line 50 in /Users/GPL/local/PETSc/petsc-3.5.2/src/ksp/ksp/impls/bcgs/bcgs.c</p><p>
</p><p>[0]PETSC ERROR: #22 KSPSolve() line 459 in /Users/GPL/local/PETSc/petsc-3.5.2/src/ksp/ksp/interface/itfunc.c</p><span class="HOEnZb"><font color="#888888"><div><div><br></div>-- <br><div dir="ltr"><div>Best regards,</div><div><br></div></div>
</div></font></span></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener
</div></div>