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)
    # Calculated neighbors
    n0 = mesh.neighbors[0]
    n1 = mesh.neighbors[1]
    n2 = mesh.neighbors[2]
    # Expected neighbors
    en0 = [1, 2]
    en1 = [0, 2]
    en2 = [0, 1]
    assert((n0 == en0) and (n1 == en1) and (n2 == en2))