def METHOD_NAME(surfs_1, surfs_2, multiplicity1, multiplicity2):
    """
    Given surface points and multiplicities of the surfaces this returns the maximal overlap fraction of the sites
    Attr:
    surfs : lists of the surface points on each side of the interface.
    multiplicity : lists of the multiplicity of the lattice vectors of u and v for each side of the interface.
    Returns:
    max_csl : float, the maximum fraction overlap found.
    """
    csl_values = []
    for surface_1 in surfs_1:
        if len(surface_1) > 0:
            surf_1_super = super_surface(
                np.asarray(surface_1), np.asarray(multiplicity1)
            )
            for surface_2 in surfs_2:
                if len(surface_2) > 0:
                    surf_2_super = super_surface(
                        np.asarray(surface_2), np.asarray(multiplicity2)
                    )
                    for i in np.arange(0, 1, 0.1):
                        for j in np.arange(0, 1, 0.1):
                            t_surf = translate(surf_2_super, [i, j])
                            csl_values.append(
                                csl(surf_1_super, t_surf, multiplicity1)
                            )
    return max(csl_values)