def METHOD_NAME():
    cell = triangle
    x = SpatialCoordinate(cell)
    I = Identity(cell.geometric_dimension())
    i, j = indices(2)
    s = (x[i] * x[j]) * I[i, j]
    e = s((5, 7))
    # v = sum_i sum_j x_i x_j delta_ij = x_0 x_0 + x_1 x_1
    v = 5 ** 2 + 7 ** 2
    assert e == v