Alex Jones

I'd recommend to do this with an interpolation, There is an idea in the works to to do this automatically in some cases but in needs changes in Interpolations.jl

Ah ok I have an idea, yes I will implement a helper that will discretize a new u0 on to the grid, or you can use it directly if you...

If you wanted to solve steady state, you need to write your problem independent of time and solve with a NonlinearProblem, see the docs: https://docs.sciml.ai/MethodOfLines/stable/tutorials/heatss/

I'll add that feature then, thanks for the heads up!

I see the problem, you also need to supply boundary conditions for c, this needs a better error message though

There is no hook to wrap the solution object in `step!`, this is by design as it isn't an optimal store for ODE solving. If you can construct an ODESolution...

Unfortunately this straight up will not work, until the generated ODEFunction is size invariant. Generated functions are not currently invariant under domain size so there is no way to get...

I'm afraid that'd be what I recommend, best of luck! Could you consider posting your model when you get it running so I can add a tutorial for this?

This is an unknown problem, thank you for bringing this to my attention. I'll add a note to the docs. Can you try without aux variables? System transformation should handle...

Yes, to my knowledge this should be equivalent -- this is definitely unexpected. From experience the cause is likely upstream, but this is tough to verify. How does it do...