On Wed, May 30, 2012 at 7:20 PM, Jozsef Bakosi <span dir="ltr"><<a href="mailto:jbakosi@lanl.gov" target="_blank">jbakosi@lanl.gov</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hello,<br>
<br>
I have a mesh with very small cell sizes and in a Poisson equation CG solve<br>
with ML preconditioner I get the following error:<br>
<br>
[31]PETSC ERROR: --------------------- Error Message ------------------------------------<br>
[31]PETSC ERROR: Detected zero pivot in LU factorization see<br>
<a href="http://www.mcs.anl.gov/petsc/petsc-as/documentation/troubleshooting.html#ZeroPivot" target="_blank">http://www.mcs.anl.gov/petsc/petsc-as/documentation/troubleshooting.html#ZeroPivot</a>!<br>
[31]PETSC ERROR: Zero pivot row 537 value 9.41111e-13 tolerance 1e-12!<br>
<br>
I really would like PETSc to NOT treat matrix values lower than 1e-12 to<br>
be zero, so I'm trying to set the zero pivot to a lower value with<br>
PCFactorSetZeroPivot(), but I still get the same error message: as if<br>
the function call would have no effect. I also tried setting the<br>
zeropivot via PetscOptionsSetValue("-pc_factor_zeropivot", ...), but<br>
still no luck.<br>
<br>
Is there another way to set the zero pivot that I'm not aware of?<br></blockquote><div><br></div><div>1) Always look at -ksp_view</div><div><br></div><div>2) I bet this is a subsolve, so that it would be -sub_pc_factor_zeropivot <value>. You</div>
<div>    can verify by looking at -ksp_view</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Is there a way to query the zero pivot value?<br></blockquote><div><br></div><div>3) Its printed out, 1.0e-12</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">

Thanks,<br>
J<br>
<span class="HOEnZb"><font color="#888888">--<br>
Jozsef Bakosi<br>
Computational Physics Group (CCS-2)<br>
MS T006, Los Alamos National Laboratory<br>
Los Alamos, NM 87545, USA<br>
email: <a href="mailto:jbakosi@lanl.gov">jbakosi@lanl.gov</a><br>
work: <a href="tel:505-663-5607" value="+15056635607">505-663-5607</a><br>
fax : <a href="tel:505-663-5504" value="+15056635504">505-663-5504</a><br>
</font></span></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<br>