def METHOD_NAME(cls, geom, face_centers: np.ndarray[tuple[int, Literal[2]], np.uint32]) -> THybridGeometry:
    """
    Create from existing facets and just select the ones that lie inside this polygon.
    """
    if isinstance(geom, (LineString, MultiLineString)):
        return HybridGeometry(geom, ())
    faces = tuple(
        set(np.argwhere(shapely_to_mpl(subgeom).contains_points(face_centers)).flatten())
        for subgeom in assert_multipolygon(geom)
    )
    faces = tuple(reduce(operator.or_, faces, set()))
    return HybridGeometry(geom, faces)  # old code had wrong typing
    # return HybridGeometry(geom, tuple(f for f in faces if f))  # old code had wrong typing