[aspect-devel] Composition gradient refinement strategy

[Lev Karatun]

Hi Wolfgang,

given the initial global refinement level, one of the layers was too thin
and didn't show in the output (I guess not a single quadrature point was
inside it?).  When I removed the layer from the setup (and made the
neighboring layer thicker), the error disappeared; changing the refinement
from 3+3 to 6+0 also helped. Also, switching to "composition" refinement
strategy without removing the layer helps, so I do believe there is a
problem with the plug-in.

Best regards,
Lev Karatun.

2017-06-07 13:48 GMT-04:00 Wolfgang Bangerth <bangerth at tamu.edu>:

> On 05/22/2017 11:35 PM, Lev Karatun wrote:
>
>>
>> I tried to use the composition gradient refinement strategy (using the
>> latest development version), but on the 1st refinement iteration the region
>> with the lowest gradient was refined. At the next iteration the simulation
>> just hangs even for the simple setup (prm file attached). Am I doing
>> something wrong or is the plugin not working properly?
>>
>
> Lev -- are you running this in debug mode? Because when I do, I run into
> assertions that say that some of the mesh refinement criteria are infinite.
> I haven't tried to figure out why that may be, but it's worth looking at
> the first error first.
>
> Best
>  W.
>
>
> --
> ------------------------------------------------------------------------
> Wolfgang Bangerth          email:                 bangerth at colostate.edu
>                            www: http://www.math.colostate.edu/~bangerth/
>
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20170611/044d8acf/attachment.html>