[aspect-devel] Aspect hangs after several time steps

[Wolfgang Bangerth]

On 02/23/2016 04:58 PM, Rene Gassmoeller wrote:
>
> The problem is that the 'Simpler' material model is not using compositional
> fields, and does not consider that the model might still be solved with
> passive compositional fields. Therefore the 'reaction_terms' variable, which
> describes the change of the compositional fields over one timestep remains
> unfilled and contains the uninitialized value 'NaN'.

Ah, wow, good detective work. I think that the material models should actually 
set these reaction terms to zero, since that's the material behavior they 
describe. This is likely an issue for not just the Simpler material model.