<div dir="ltr">Hi Jed. I have a similar situation here. I am intentionally testing an iterative method. It requires me to do some handling (for example copying the current solution vector out) right after a nonlinear iteration is finished, but before anything happening of the next nonlinear iteration. Is there any way to do so?<div>
<br></div><div>Best,</div><div><br></div><div>Ling</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Nov 21, 2013 at 5:41 AM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">TAY wee-beng <<a href="mailto:zonexo@gmail.com">zonexo@gmail.com</a>> writes:<br>
<br>
> Hi Jed,<br>
><br>
> I need a solve a Poisson equation, which is part of my finite volume CFD<br>
> code.<br>
><br>
> The RHS terms are calculated from the velocity.<br>
><br>
> In the standard mtd, the RHS terms are calculated from the velocity,<br>
> which then solves the Poisson eqn to get the pressure.<br>
><br>
> Now, to improve the accuracy of the pressure obtained, a paper suggested<br>
> obtaining the pressure  from the velocity in 1 iteration. Then use the<br>
> new pressure for the momentum eqn to get a more accurate intermediate<br>
> velocity. The RHS terms are recomputed again using the new velocity,<br>
> which is then used to calculate the new pressure. This cycle goes on<br>
> until the pressure converges to a certain value.<br>
<br>
</div>This is a crude iterative method for solving a coupled problem.  I<br>
recommend that you solve the coupled problem.  See the section of the<br>
user's manual on solving block problems, and the examples in PETSc that<br>
solve Stokes-type problems.<br>
</blockquote></div><br></div>