def face_with_verts(bm, verts, default=None):
    """ Find a face in the bmesh with the given verts"""
    for face in bm.faces:
        if len(set(list(face.verts) + verts)) == len(verts):
            return face
    return default