def METHOD_NAME():
    dtype_float = np.float32
    dtype_int = np.int64
    vertices = np.array([[0, 0, 0], [1, 0, 0], [0, 1, 0], [0, 0, 1]],
                        dtype=dtype_float)
    faces = np.array([[0, 2, 1], [0, 3, 2], [0, 1, 3]], dtype=dtype_int)
    mesh = MeshData(vertices=vertices, faces=faces)
    vertex_normals_unnormalized = np.array(
        [[-1, -1, -1], [0, -1, -1], [-1, 0, -1], [-1, -1, 0]],
        dtype=dtype_float)
    norms = np.sqrt((vertex_normals_unnormalized**2).sum(axis=1,
                                                         keepdims=True))
    vertex_normals = vertex_normals_unnormalized / norms
    expected_vertex_normals = vertex_normals[faces]
    computed_vertex_normals = mesh.get_vertex_normals(indexed="faces")
    assert_array_equal(expected_vertex_normals, computed_vertex_normals)