fm_components function

Compute connected mesh subsets