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.