Alex Jones

That would be due to `structural_simplify` simplifying away points due to the specified Neumann BCs, please use `sol[u]` to retrieve all points including those that were simplified out.

@YingboMa I'm sure that this should work, we're doing it throughout the tests. All state/observed variables look like `u[i](t)`. What is going wrong here?

Not the interpolations, @zoemcc did you get anywhere with the solution interfaces or should I spearhead this?

I'm afraid that second time derivatives are not yet supported, I will be implementing them this week however.

Complex equations are not currently supported but I will be working on this, can you reformulate as a system of equations in terms of the real and imaginary parts of...

This seems to be above board, this is a bug. I'll take a closer look at this soon.

Would the regions ever need non orthoganal/perpendicular boundaries?

@ChrisRackauckas Everything passing except invalidations, is this ok to merge?

This looks like a problem in `NonlinearSolve`, are you running with latest versions? Can you retry after running `Pkg.update()`?

Have you updated to v0.3.1? This fixed an issue with lhs advection terms like those appearing in your equation.