Laplace solver interface

[bendudson]

The solve functions in class Laplacian (the interface to the user) should not be overridden by implementations. Instead, implementations should implement protected functions (not visible to the user), which are called by the interface functions.

This would help ensure a stable interface, and allow common tasks like timing and setting of result cell location to be done in one place.