def METHOD_NAME():
    # Construct a single triangle
    vertices = np.array([[0.5, 0, 0.5], [0, 0.5, 0.5], [0.5, 0.5, 0.5]])
    faces = np.array([[0, 1, 2]])
    mesh = triangular_mesh.TriangularMesh(vertices, faces)
    mesh.add_vertex(0)
    expected_faces = np.array([[0,3,2],[0,1,3],[1,2,3]])
    np.testing.assert_array_almost_equal(mesh.faces, expected_faces, 4)