<div class="gmail_quote">On Mon, Apr 23, 2012 at 20:40, tribur <span dir="ltr"><<a href="mailto:burckhardt@itis.ethz.ch">burckhardt@itis.ethz.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Jed,<br>
<br>
Thanks for the ultra fast response.<br>
<br>
Apart from the nonzero initial guess, is -ksp_type preonly -pc_type sor -pc_sor_its 10 right?<br></blockquote><div><br></div><div>Yes.</div><div><br></div><div>You can use -ksp_initial_guess_nonzero -ksp_type richardson</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
Am <a href="tel:24.04.2012%2003" value="+12404201203" target="_blank">24.04.2012 03</a>:29, schrieb Jed Brown:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
On Mon, Apr 23, 2012 at 20:27, tribur <<a href="mailto:burckhardt@itis.ethz.ch" target="_blank">burckhardt@itis.ethz.ch</a> [1]><div class="im"><br>
wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I would like to do 10 iterations of a point-by-point Gauss-Seidel<br>
procedure for the purpose of testing.<br>
<br>
Is it correct if I type "mpiexec -n 1 ... -ksp_type preonly<br>
-pc_type sor -pc_sor_its 10"?<br>
<br>
How can I define the initial guess nonzero?<br>
</blockquote>
<br>
With a nonzero initial guess, you are doing Richardson (not preonly).<br>
You can set the nonzero initial guess as usual<br>
(KSPSetInitialGuessNonzero()).<br>
<br></div>
Links:<br>
------<br>
[1] mailto:<a href="mailto:burckhardt@itis.ethz.ch" target="_blank">burckhardt@itis.ethz.ch</a><br>
</blockquote>
<br>
</blockquote></div><br>