<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sat, Apr 11, 2015 at 7:03 PM, TAY wee-beng <span dir="ltr"><<a href="mailto:zonexo@gmail.com" target="_blank">zonexo@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hi,<br>
<br>
I'm using DMDA. The LHS of the poisson eqn is fixed, only the RHS changes.<br>
<br>
I thought I can't find ex29.c in the current ver of PETSc, so I thought it's no longer valid. Do I always have to get the jacobian as in ex29 to solve the Poisson eqn?<br></blockquote><div><br></div><div>It is there:</div><div><br></div><div>  <a href="https://bitbucket.org/petsc/petsc/src/0dbc604a2409e3ecaa129deb24062cf2be4d41a7/src/ksp/ksp/examples/tutorials/ex29.c?at=master">https://bitbucket.org/petsc/petsc/src/0dbc604a2409e3ecaa129deb24062cf2be4d41a7/src/ksp/ksp/examples/tutorials/ex29.c?at=master</a></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
What about a good ex. using KSPBCGS to solve my momentum eqn with DMDA?<br></blockquote><div><br></div><div>All examples are independent of the precise Krylov method. You can use that anywhere:</div><div><br></div><div>  cd src/snes/examples/tutorials</div><div>  make ex5</div><div>  ./ex5 -ksp_type bcgs</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
Thank you.<br>
<br>
Yours sincerely,<br>
<br>
TAY wee-beng<br>
<br>
On 4/11/2015 6:43 AM, Barry Smith wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
   Are you using DMDA or DMPLEX or some other DM? Does the Poisson equation change (variable coefficient) or is it a constant coefficient with the same values for all solves? ex29.c could still possibly serve as a model for the Poisson subproblem, why do you say it is no longer valid?<br>
<br>
   Barry<br>
<br>
<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
On Apr 8, 2015, at 10:18 PM, TAY wee-beng <<a href="mailto:zonexo@gmail.com" target="_blank">zonexo@gmail.com</a>> wrote:<br>
<br>
Hi,<br>
<br>
I was following example ex29 in the past to write the code to solve my Poisson eqn using DM. Now it seems that ex29 is no longer valid. Which example should I follow now?<br>
<br>
Also, I am using KSPBCGS to solve my momentum eqn in the past. Do I need to change much if I'm now using DM? If so, what is a good example?<span class=""><font color="#888888"><br>
<br>
-- <br>
Thank you<br>
<br>
Yours sincerely,<br>
<br>
TAY wee-beng<br>
<br>
</font></span></blockquote></blockquote>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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</div>
</div></div>