Refactor and modularize the penalty resolution classes
The two almost identical implementations of the linear and quadratic penalty method (resolution_penalty.cc
and resolution_penalty_quadratic.cc
) are merged into a single class templated with PenaltyFunction
objects of which new ones can easily be added now.