def distance_residues(r1, r2):
    """Return the distance between two residues."""

    center_1 = r1.get_mass_center()
    center_2 = r2.get_mass_center()
    return (sum(map(lambda t: (t[0] - t[1])**2, zip(center_1, center_2))))**0.5