Calculate solution errors

In addition to setting up initial conditions and setting the solution all the time, offer to calculate the solution error as analysis quantity.

The implementation is straightforward, only two implementation choices are non-trivial:
- Calculating the solution error is implemented by transforming (via pattern matching) the calculation that sets up the solution
- By default, the solution error is not calculated; a keyword parameter "error_method" needs to be set to activate this.