WIP: Extend `ProximalNewtonSolver` to more general `ProximalSolver`
With this change we can use the structure of the Proximal Newton solver for problem with less differential order of the smooth part.
Things that need some further input:
- I want to provide shortcuts to
ProximalNewtonSolverandProximalGradientSolverwith aliases usingusing. - However, this is a
C++20feature if I want to omit the explicit template arguments - The Methods
setHessianetc. are only meaningful fordiffOrder > 1. Is there a smooth SFINAE method to provide these methods only in this case?