Actually I want to update the domain according to the global CFL number; e.g. dt = CFL * dx / Lambda; separately for each cell <br>( dx is the dimension of the cell and Lambda is the local maximum eigenvalue of the field )<br>
<br>So, here shift is the maximum time step and I should adjust other cells to this maximum value?<br><br>Thanks,<br><br><br> <br><br><div class="gmail_quote">On Thu, Oct 27, 2011 at 5:49 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov">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"><div class="gmail_quote">On Thu, Oct 27, 2011 at 06:44, behzad baghapour <span dir="ltr"><<a href="mailto:behzad.baghapour@gmail.com" target="_blank">behzad.baghapour@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>2- page 120 of PETSc manual (TS): "For location-dependent pseudo-timestepping, the interface function has not yet been created."<br>
Is there any way to control time-step locally? (This is due to variable time updating usually needed in PDE solution)<br></div></blockquote></div><br></div><div>The documentation needs to be updated. Use TSSetIJacobian() and regard the "shift" as a maximum shift (corresponds to minimum time step). If you have a cell in which you want to take a 10x longer time step, then use shift/10 in that cell. Note that you can use steady-state constraints on some fields. For example, src/ts/examples/tutorials/ex26.c (in petsc-dev) uses no shift for the velocity components (enforces incompressibility in velocity-vorticity form).</div>
</blockquote></div><br><br clear="all"><br>-- <br>==================================<br>Behzad Baghapour<br>Ph.D. Candidate, Mechecanical Engineering<br>University of Tehran, Tehran, Iran<br><a href="https://sites.google.com/site/behzadbaghapour" target="_blank">https://sites.google.com/site/behzadbaghapour</a><br>
Fax: 0098-21-88020741<br>==================================<br><br>