Introduced problem_description class holding the mesh pointer and associated quantities.

This way the PDE solver is decoupled from the quantities, so one can plug&play different PDE solvers on top of a the same quantities.