<div dir="ltr">H<div style="font-family:arial,helvetica,sans-serif;display:inline" class="gmail_default">i all,</div><div><div style="font-family:arial,helvetica,sans-serif;display:inline" class="gmail_default"></div></div>

<div><div style="font-family:arial,helvetica,sans-serif;display:inline" class="gmail_default">Here's a bit of a doozy of a problem.  I noticed the other day a sharp increase in the number of iterations aspect required for the Stokes solver to converge.  As in, the cheap solver never converges, and the expensive one takes ~50-200 iterations.  After a bunch of experimentation, I think it may be due to some recent changes in deal.ii.  If I run simple_2d_shells.prm on 8 cores with deal.ii 7.3.0, I get the following representative output: </div>
</div><div><div style="font-family:arial,helvetica,sans-serif;display:inline" class="gmail_default"><br></div></div><div><div style="display:inline" class="gmail_default"><div><font face="arial, helvetica, sans-serif">*** Timestep 148:  t=4.6211e+06 years</font></div>
<div><font face="arial, helvetica, sans-serif">   Solving temperature system... 13 iterations.</font></div><div><font face="arial, helvetica, sans-serif">   Solving Stokes system... 13 iterations.</font></div><div style="font-family:arial,helvetica,sans-serif">
<br></div><div style="font-family:arial,helvetica,sans-serif">If, however, I run it using deal.ii r29892 I get the following</div><div style="font-family:arial,helvetica,sans-serif"><br></div><div><div><font face="arial, helvetica, sans-serif">*** Timestep 27:  t=2.68009e+06 years</font></div>
<div><font face="arial, helvetica, sans-serif">   Solving temperature system... 13 iterations.</font></div><div><font face="arial, helvetica, sans-serif">   Solving Stokes system... 30+219 iterations.</font></div><div style="font-family:arial,helvetica,sans-serif">
<br></div></div><div style="font-family:arial,helvetica,sans-serif">I have no idea what may be causing this, or how reproducible it is across different systems (though I have noticed the same increase in iterations on TACC Stampede).  In both cases, I am using trilinos 11.2.3 and p4est 0.3.4.  Everything is built with icc-12.1, mkl math libraries, and using myrinet MPI.  Any insight would be greatly appreciated.</div>
<div style="font-family:arial,helvetica,sans-serif">Best,</div><div style="font-family:arial,helvetica,sans-serif">Ian</div><div style="font-family:arial,helvetica,sans-serif"><br></div></div></div></div>