I was trying to be general and allow for the case where you know the flux vector, but not the normal. But maybe this was a bad idea, particularly since I had to add the API for the case where you only know the normal component (see CatalyticWall).

And you're right - I've not at all thought about 4th order problems. I'm definitely open to updating API to handle that.