<div class="gmail_quote">On Fri, Jul 8, 2011 at 08:26, TAY wee-beng <span dir="ltr">&lt;<a href="mailto:zonexo@gmail.com">zonexo@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div bgcolor="#ffffff" text="#000000">[0]PETSC ERROR: KSPSolve_BCGS() line 75 in
    src/ksp/ksp/impls/bcgs/C:\Codes\PETSC<br>
    -~1.1-P\src\ksp\ksp\impls\bcgs\bcgs.c</div></blockquote></div><div><br></div><div><div>    ierr = KSP_PCApplyBAorAB(ksp,P,V,T);CHKERRQ(ierr);  /*   v &lt;- K p           */</div><div>    ierr = VecDot(V,RP,&amp;d1);CHKERRQ(ierr);</div>
<div>    if (d1 == 0.0) SETERRQ(PETSC_ERR_PLIB,&quot;Divide by zero&quot;);</div></div><div><br></div><div>I suspect your preconditioner is singular. What options are you running with? What happens if you use -ksp_type gmres -ksp_gmres_restart 1000 -ksp_gmres_modifiedgramschmidt -ksp_monitor_true_residual? Also do the same run with -ksp_type fgmres.</div>