<span id="mailbox-conversation"><div>Emil:  If I want to fix this myself, do you have any idea where I should start looking?  I assume in the th->endpoint section of TSStep_Theta in http://www.mcs.anl.gov/petsc/petsc-current/src/ts/impls/implicit/theta/theta.c.html?</div>
<div><br></div>
<div>Thanks,</div>
<div><br></div>
<div>-Andrew</div></span><div class="mailbox_signature"><br></div>
<br><br><div class="gmail_quote"><p>On Tue, Mar 24, 2015 at 10:52 AM, Lisandro Dalcin <span dir="ltr"><<a href="mailto:dalcinl@gmail.com" target="_blank">dalcinl@gmail.com</a>></span> wrote:<br></p><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><p>On 24 March 2015 at 18:40, Emil Constantinescu <emconsta@mcs.anl.gov> wrote:
<br>>> Of course, though my B-D is order 2, and as you said, it is the just
<br>>> for the estimator.
<br>>>
<br>>>> if the time step varies wildly (can be easily constrained),
<br>>>> it may create some problems with the estimator.
<br>>>>
<br>>>
<br>>> Indeed, the default clipping 0.1,10 is usually too wide, I usually do
<br>>> -ts_adapt_basic_clip 0.5,2.0
<br>>
<br>>
<br>> The limits are a bit more strict for BDF2, but that's for the worst case and
<br>> no filtering. So that looks good too.
<br><br>I'm not sure I'm following you 100%. Please note I'm not using BDF2 in
<br>the usual sense, what I do is to compute an A POSTERIORI error
<br>estimator using a BD of order 2 to estimate LTE^{n+1} with the
<br>solution vectors U^{n+1}, U^{n} and U^{n-1}, and this is AFTER
<br>computing U^{n-1}.
<br><br><br>-- 
<br>Lisandro Dalcin
<br>============
<br>Research Scientist
<br>Computer, Electrical and Mathematical Sciences & Engineering (CEMSE)
<br>Numerical Porous Media Center (NumPor)
<br>King Abdullah University of Science and Technology (KAUST)
<br>http://numpor.kaust.edu.sa/
<br><br>4700 King Abdullah University of Science and Technology
<br>al-Khawarizmi Bldg (Bldg 1), Office # 4332
<br>Thuwal 23955-6900, Kingdom of Saudi Arabia
<br>http://www.kaust.edu.sa
<br><br>Office Phone: +966 12 808-0459
<br></p></blockquote></div><br>