Mixed-type boundary conditions for box
The current implementation in
assembly/boxlocalresidual.hh forbids boundary conditions of mixed type, i.e. where one subset of equations gets Dirichlet conditions while the complementary subset gets Neumann conditions. This is too restrictive. The box method can handle such mixed-type boundary conditions easily. They are frequently encountered in, e.g., solid mechanics problems.