<div>I’m sorry, I didn’t want to imply you were taking too much time.  I was just interested in trying to solve the problem myself, especially if it was a simpler fix.</div>
<div id="mb-reply"><br></div>
<div id="mb-reply">I don’t have enough of an understanding of the internals of PETSc though, if it depends on the compiler flags.</div>
<div id="mb-reply"><br></div>
<div id="mb-reply">Good luck!</div>
<div id="mb-reply"><br></div>
<div id="mb-reply">-Andrew</div>
<div class="mailbox_signature">
<br>—<br>Andrew</div>
<br><br><div class="gmail_quote"><p>On Mon, Mar 30, 2015 at 7:23 PM, Emil Constantinescu <span dir="ltr"><<a href="mailto:emconsta@mcs.anl.gov" target="_blank">emconsta@mcs.anl.gov</a>></span> wrote:<br></p><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><p>Hi Andrew,
<br><br>Sorry for the delay; things are going slower with me over the spring 
<br>break but we are working on it.
<br><br>The problem that you see may not be related to the integrator, which 
<br>runs well (we get expected convergence rates and good results) for other 
<br>problems, in particular, for problems with real variables (and no 
<br>c++11). So the high level th->endpoint implementation seems fine.
<br><br>For some reason the linear solver is misbehaving when using cn. If you 
<br>use option -snes_mf the problem disappears.
<br><br>Please do not use -ts_theta_adapt until it is fixed.
<br><br>In conclusion I don't think that looking at theta.c will be productive. 
<br>We are currently looking at what happens at the lower levels.
<br><br>Thanks,
<br>Emil
<br><br>On 3/30/15 4:25 PM, Andrew Spott wrote:
<br>> Emil:  If I want to fix this myself, do you have any idea where I should
<br>> start looking?  I assume in the th->endpoint section of TSStep_Theta in
<br>> http://www.mcs.anl.gov/petsc/petsc-current/src/ts/impls/implicit/theta/theta.c.html?
<br>>
<br>> Thanks,
<br>>
<br>> -Andrew
<br>>
<br>>
<br>>
<br></p></blockquote></div><br>