[aspect-devel] Problems with recent deal.ii

[Ian Rose]

Hi Timo,
I just tried undoing that change, and am still getting large iteration
counts.
Ian



On Mon, Jul 1, 2013 at 9:06 AM, Timo Heister <heister at math.tamu.edu> wrote:

> Dear Ian,
>
> Could this be related to r29539 ?
> http://dealii.org/websvn/revision.php?repname=deal.II+Repository&path=%2F&rev=29539&peg=29539
> On Jul 1, 2013 5:50 PM, "Ian Rose" <ian.rose at berkeley.edu> wrote:
>
>> H
>> i all,
>> 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:
>>
>> *** Timestep 148:  t=4.6211e+06 years
>>    Solving temperature system... 13 iterations.
>>    Solving Stokes system... 13 iterations.
>>
>> If, however, I run it using deal.ii r29892 I get the following
>>
>> *** Timestep 27:  t=2.68009e+06 years
>>    Solving temperature system... 13 iterations.
>>    Solving Stokes system... 30+219 iterations.
>>
>> 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.
>> Best,
>> Ian
>>
>>
>> _______________________________________________
>> Aspect-devel mailing list
>> Aspect-devel at geodynamics.org
>> http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>>
>
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geodynamics.org/pipermail/aspect-devel/attachments/20130701/1c38a203/attachment.html>