Fixed copy&paste error when throwing an exception.
Implemented flux calculation between segments. 1D working, some work left for 2D and 3D.

Interface is not yet final: Right now one has to pass the full PDE system to the accessor, which is too much information.
Instead, one should only pass the flux expression, which is more flexible and makes more sense from a usability point of view.