This MR attempts to simplify the boundary condition interface by providing:

• an abstract interface
• more documentation
• clearer function names
• deduction for some template arguments

The introduction of virtual functions requires additional template arguments Mat, Sol, Rhs. Those may be removed later with the change to local application of constraints.