Torque calculation from body forces is complicated by discontinuities in field gradients across material boundaries

Fix for Mesh::insert_internal_boundaries with small airgaps
    This method is probably not 100% fixed, but it is more stable than before (SalientPoleSynchrel regression tests run with CylindricalAirgap<2>)