Wait, so is it for PREONLY? And what Matsolve does, solving a triangular system with both L and U available?<br>The reason that I previously guess the count be two is that I think there is a back-substituiton<br>and a forward-substitution involved in solving a linear system using factorization. If a pair of <br>
back-substitution and forward-substitution counts 1 MatSolve. Then I think we mean the same thing.<br><br>Thanks,<br><br>Yan<br><br>On Thu, Apr 7, 2011 at 10:08 AM, Jed Brown <<a href="mailto:jed@59a2.org">jed@59a2.org</a>> wrote:<br>
<div class="gmail_quote"><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div dir="ltr"><div class="im"><div class="gmail_quote">On Thu, Apr 7, 2011 at 16:06, Ryan Yan <span dir="ltr"><<a href="mailto:vyan2000@gmail.com" target="_blank">vyan2000@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div>Hi Jed,<br>How is going? :-)<br><br>PREONLY maybe also use two MatSolve? Or is there any magic I did not see. :-)</div></blockquote></div><br></div><div>Why do you say that?</div><div><br></div><div><div>$ ./ex2 -pc_type lu -ksp_type preonly -log_summary |g '^MatSolve'</div>
<div>MatSolve 1 1.0 1.1921e-05 1.0 1.22e+03 1.0 0.0e+00 0.0e+00 0.0e+00 0 22 0 0 0 0 22 0 0 0 102</div></div></div>
</blockquote></div><br>