Note that this will need to be revisited if we add support for stress boundary conditions, as it is essentially the same problem