<div dir="ltr"><br><div class="gmail_extra"><div class="gmail_quote">On Tue, Oct 15, 2013 at 6:20 PM, 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:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div class="im">Christophe Ortiz <<a href="mailto:christophe.ortiz@ciemat.es">christophe.ortiz@ciemat.es</a>> writes:<br>

<br>
> Yes ! Actually, that was the problem ! Now I set ARKIMEX1BEE after<br>
> TSSetFromOptions() and it is much faster ! ARKIMEXA2 seems even faster.<br>
<br>
</div>I still recommend configuring this stuff from the command line, and I<br>
would put TSSetFromOptions last so that you can override anything at<br>
run-time.<br></blockquote><div><br></div><div>Ok.</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">

<div class="im"><br>
> Thanks Jed !<br>
><br>
> Now I see that the convergence speed of ARKIMEX depends on where I put<br>
> TSSetFromOptions().<br>
><br>
> If I do:<br>
> TSCreate();<br>
> TSSetFromOptions();<br>
> TSSetProblemType();<br>
> TSSetType();<br>
> TSARKIMEXSetType();<br>
><br>
> then ARKIMEXA2 is very fast.<br>
><br>
> Instead, if I put TSSetFromOptions() somewhere else:<br>
> TSCreate();<br>
><br>
> TSSetProblemType();<br>
> TSSetType();<br>
</div>> *TSSetFromOptions();*<br>
> TSARKIMEXSetType();<br>
<br>
I would use this order:<br>
<br>
    TSCreate();<br>
    TSSetProblemType();  // optional<br>
    TSSetType();<br>
    TSARKIMEXSetType();<br>
    TSSetFromOptions();<br>
<div class="im"><br>
> ARKIMEXA2 slows down significantly. I guess that TSSetFromOptions overrides<br>
> something...Any idea ?<br>
<br>
</div>Compare the output from -ts_view in both cases.<br></blockquote><div><br></div><div>I found where is the problem. I had set TSARKIMEXSetFullyImplicit(ts,PETSC_TRUE);</div><div>I found it produces wrong solutions, either with ARKIMEX3, 1BEE or A2. Solution becomes negative. Without this option, all ARKIMEX types work very well and very fast.</div>
<div><br></div><div>Christophe</div><div> </div></div><br></div></div>