Problem with setting initial conditions and history variable

General
2

Hi everyone,

I managed to resolve the issues in my code. The NaN error (“Does your Newton solver fail on the first step producing NaNs?”) was due to the need to interpolate the solution ‘up‘ between the previous and current timesteps, as described here.

Additionally, I discovered a simple mistake in how I was storing some values, which led to unexpected (weird) results. The correct order should be something like:

    tm_sl = up.sub(0).collapse()
    hyddeg_sl = up.sub(1).collapse()

    tm_sl.x.scatter_forward()
    hyddeg_sl.x.scatter_forward()

    tm_sol.interpolate(tm_sl)
    hyddeg_sol.interpolate(hyddeg_sl)

    # and then:
    tm_n.x.array[:] = tm_sol.x.array
    hyddeg_n.x.array[:] = hyddeg_sol.x.array

I hope this helps anyone encountering similar problems early on!

Best regards,
Silvio