def test_automatic_dimension():
st = SimplexTree()
assert st.__is_defined() == True
assert st.__is_persistence_defined() == False
# insert test
assert st.insert([0, 1, 3], filtration=0.5) == True
assert st.insert([0, 1, 2], filtration=1.0) == True
assert st.num_vertices() == 4
assert st.num_simplices() == 11
assert st.dimension() == 2
assert st.upper_bound_dimension() == 2
assert st.prune_above_filtration(0.6) == True
assert st.dimension() == 2
assert st.upper_bound_dimension() == 2
st.assign_filtration([0, 1, 3], 0.7)
assert st.filtration([0, 1, 3]) == 0.7
st.remove_maximal_simplex([0, 1, 3])
assert st.upper_bound_dimension() == 2
assert st.dimension() == 1
assert st.upper_bound_dimension() == 1
def test_extend_filtration():
# Inserted simplex:
# 5 4
# o o
# / \ /
# o o
# /2\ /3
# o o
# 1 0
st = SimplexTree()
st.insert([0, 2])
st.insert([1, 2])
st.insert([0, 3])
st.insert([2, 5])
st.insert([3, 4])
st.insert([3, 5])
st.assign_filtration([0], 1.)
st.assign_filtration([1], 2.)
st.assign_filtration([2], 3.)
st.assign_filtration([3], 4.)
st.assign_filtration([4], 5.)
st.assign_filtration([5], 6.)
assert list(st.get_filtration()) == [([0, 2], 0.0), ([1, 2], 0.0),
([0, 3], 0.0), ([3, 4], 0.0),
([2, 5], 0.0), ([3, 5], 0.0),
([0], 1.0), ([1], 2.0), ([2], 3.0),
([3], 4.0), ([4], 5.0), ([5], 6.0)]
st.extend_filtration()
assert list(st.get_filtration()) == [([6], -3.0), ([0], -2.0), ([1], -1.8),
([2], -1.6), ([0, 2], -1.6),
([1, 2], -1.6), ([3], -1.4),
([0, 3], -1.4), ([4], -1.2),
([3, 4], -1.2), ([5], -1.0),
([2, 5], -1.0), ([3, 5], -1.0),
([5, 6], 1.0), ([4, 6], 1.2),
([3, 6], 1.4), ([3, 4, 6], 1.4),
([3, 5, 6], 1.4), ([2, 6], 1.6),
([2, 5, 6], 1.6), ([1, 6], 1.8),
([1, 2, 6], 1.8), ([0, 6], 2.0),
([0, 2, 6], 2.0), ([0, 3, 6], 2.0)]
dgms = st.extended_persistence(min_persistence=-1.)
assert dgms[0][0][1][0] == pytest.approx(2.)
assert dgms[0][0][1][1] == pytest.approx(3.)
assert dgms[1][0][1][0] == pytest.approx(5.)
assert dgms[1][0][1][1] == pytest.approx(4.)
assert dgms[2][0][1][0] == pytest.approx(1.)
assert dgms[2][0][1][1] == pytest.approx(6.)
assert dgms[3][0][1][0] == pytest.approx(6.)
assert dgms[3][0][1][1] == pytest.approx(1.)
def test_make_filtration_non_decreasing():
st = SimplexTree()
assert st.__is_defined() == True
assert st.__is_persistence_defined() == False
# Inserted simplex:
# 1
# o
# /X\
# o---o---o---o
# 2 0 3\X/4
# o
# 5
assert st.insert([2, 1, 0], filtration=2.0) == True
assert st.insert([3, 0], filtration=2.0) == True
assert st.insert([3, 4, 5], filtration=2.0) == True
assert st.make_filtration_non_decreasing() == False
# Because of non decreasing property of simplex tree, { 0 } , { 1 } and
# { 0, 1 } are going to be set from value 2.0 to 1.0
st.insert([0, 1, 6, 7], filtration=1.0)
assert st.make_filtration_non_decreasing() == False
# Modify specific values to test make_filtration_non_decreasing
st.assign_filtration([0, 1, 6, 7], 0.8)
st.assign_filtration([0, 1, 6], 0.9)
st.assign_filtration([0, 6], 0.6)
st.assign_filtration([3, 4, 5], 1.2)
st.assign_filtration([3, 4], 1.1)
st.assign_filtration([4, 5], 1.99)
assert st.make_filtration_non_decreasing() == True
assert st.filtration([0, 1, 6, 7]) == 1.0
assert st.filtration([0, 1, 6]) == 1.0
assert st.filtration([0, 1]) == 1.0
assert st.filtration([0]) == 1.0
assert st.filtration([1]) == 1.0
assert st.filtration([3, 4, 5]) == 2.0
assert st.filtration([3, 4]) == 2.0
assert st.filtration([4, 5]) == 2.0