
<html><head></head><body><div class="yahoo-style-wrap" style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:13px;"><div dir="ltr" data-setdir="false">Hello,</div><div dir="ltr" data-setdir="false"><br></div><div dir="ltr" data-setdir="false">I am using the Petsc native direct solver <font face="Helvetica, Arial, sans-serif">(with KSPPREONLY and PCLU) to solver a small Linear equation system. The matrix has a couple of diagonal terms with value 0, which makes the solver fail. Then I set a very small value (e.g., 1e-10) at these terms and the solver works. Can this cause noticeable errors in the solution? Are there any options in Petsc to make the direct solver do row/column pivoting so that all the diagonal terms are non-zero?</font></div><div dir="ltr" data-setdir="false"><font face="Helvetica, Arial, sans-serif"><br></font></div><div dir="ltr" data-setdir="false"><font face="Helvetica, Arial, sans-serif">Thanks for any suggestions.</font></div><div dir="ltr" data-setdir="false"><font face="Helvetica, Arial, sans-serif"><br></font></div><div dir="ltr" data-setdir="false"><font face="Helvetica, Arial, sans-serif">Regards,</font></div><div dir="ltr" data-setdir="false"><font face="Helvetica, Arial, sans-serif">Qin</font></div></div></body></html>