<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, May 5, 2014 at 7:25 AM, linjing bo <span dir="ltr"><<a href="mailto:francium87@hotmail.com" target="_blank">francium87@hotmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">


<div><div dir="ltr">


<div dir="ltr">Hi, I'm trying to use PETSc's ksp method to solve a linear system. When running, Error is reported by VecValidValues() that NaN or Inf is found with error message listed below <br><br><br>[3]PETSC ERROR: --------------------- Error Message ------------------------------------           <br>
[3]PETSC ERROR: Floating point exception!                                                          <br>[3]PETSC ERROR: Vec entry at local location 0 is not-a-number or infinite at beginning of function: Parameter number 2!                                                                               <br>
[3]PETSC ERROR: ------------------------------------------------------------------------           <br>[3]PETSC ERROR: Petsc Release Version 3.4.4, Mar, 13, 2014                                         <br>[3]PETSC ERROR: See docs/changes/index.html for recent updates.                                    <br>
[3]PETSC ERROR: See docs/faq.html for hints about trouble shooting.                                <br>[3]PETSC ERROR: See docs/index.html for manual pages.                                              <br>[3]PETSC ERROR: ------------------------------------------------------------------------           <br>
[3]PETSC ERROR: ./gtc on a arch-linux2-c-debug named <a href="http://node2.indac.info" target="_blank">node2.indac.info</a> by jlin Mon May  5 20:03:20 2014                                                                                                <br>
[3]PETSC ERROR: Libraries linked from /opt/sfw/petsc/3.4.4/intel/openmpi/lib                       <br>[3]PETSC ERROR: Configure run at Sat Apr 26 20:19:41 2014                                          <br>[3]PETSC ERROR: Configure options --prefix=/opt/sfw/petsc/3.4.4/intel/openmpi --with-mpi-dir=/opt/sfw/openmpi/1.6.3/intel --with-blas-lapack-dir=/opt/sfw/intel/composer_xe_2011_sp1.7.256/mkl/lib/intel64 --with-mpiexec=mpiexec                                                                        <br>
[3]PETSC ERROR: ------------------------------------------------------------------------           <br>[3]PETSC ERROR: VecValidValues() line 28 in /tmp/petsc-3.4.4/src/vec/vec/interface/rvector.c  </div></div></div></blockquote>
<div><br></div><div>It looks like the vector after preconditioner application is bad. What is the preconditioner?</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><div dir="ltr"><div dir="ltr">     <br>[3]PETSC ERROR: PCApply() line 436 in /tmp/petsc-3.4.4/src/ksp/pc/interface/precon.c               <br>[3]PETSC ERROR: KSP_PCApply() line 227 in /tmp/petsc-3.4.4/include/petsc-private/kspimpl.h         <br>
[3]PETSC ERROR: KSPInitialResidual() line 64 in /tmp/petsc-3.4.4/src/ksp/ksp/interface/itres.c     <br>[3]PETSC ERROR: KSPSolve_GMRES() line 239 in /tmp/petsc-3.4.4/src/ksp/ksp/impls/gmres/gmres.c      <br>[3]PETSC ERROR: KSPSolve() line 441 in /tmp/petsc-3.4.4/src/ksp/ksp/interface/itfunc.c             <br>
<br><br>After read the source code shown by backtrack informations, I realize the problem is in the right hand side vector. So I make a trial of set right hand side vector to ONE by VecSet, But the program still shows error message above, and using VecView or VecGetValue to investigate the first value of rhs vec shows the value is 1.0 as I set it to. Hope I clearly describe the problem. The code related is listed below<br>
<br>---------------------------Solver section--------------------------<br><br>  call VecSet( pet_bp_b, one, ierr)<br><br>  vecidx=[0,1]<br>  call VecGetValues( pet_bp_b, 2, vecidx, first, ierr)<br>  write(*,*) ' first two values ', first(1), first(2)<br>
<br>  call KSPSetInitialGuessNonzero(solver_bp,Petsc_True,ierr)<br>  call KSPSolve(solver_bp,pet_bp_b,pet_bp_x,ierr)<br>  call KSPView(solver_bp, PETSC_VIEWER_STDOUT_SELF,ierr)<br>  CHKERRQ(ierr)<br><br><br></div>
                                          </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>