call_before_step_each_stage=False passes uninitialized aux to riemann solver

[GustafssonJonathan]

For sharpclaw passing to following options: solver.aux_time_de = True solver.call_before_step_each_stage = False

Seems to pass a uninitialized aux array to the Fortran Riemann solver. solver.call_before_step_each_stage = True works as expected.

[ketch]

Could you link to a full example that reproduces this behavior? Thanks.

[GustafssonJonathan]

No problem. Here is a modified version of acoustics_2d_variable The solver gives nan when solver.before_step_each_stage = False and works correctly when solver.before_step_each_stage = True

On Thu, Feb 25, 2016 at 2:46 AM, David Ketcheson [email protected] wrote:

Could you link to a full example that reproduces this behavior? Thanks.

— Reply to this email directly or view it on GitHub https://github.com/clawpack/pyclaw/issues/536#issuecomment-188657374.