Add projection equation system to ADRSolver

This MR adds a projection operator to the ADRSolver, which can be used to impose C^0 connectivity on an input forcing field and is useful in e.g. the creation of boundary conditions on complex geometries.