add a method to return the gradient of the solution, also let this and the get_solution()
method optionally store the result onto a different grid
use the new interfaces for get_solution() and get_solution_gradient()
use the new interfaces for get_solution() and get_solution_gradient()