Add an error-estimation kernel.

This can be used as part of the implementation of an adaptive
time integration scheme.