1 files changed, 98 insertions, 0 deletions
diff --git a/src/python/test/test_cubical_complex.py b/src/python/test/test_cubical_complex.py
new file mode 100755
index 00000000..68f54fbe
--- /dev/null
+++ b/src/python/test/test_cubical_complex.py
@@ -0,0 +1,98 @@
+from gudhi import CubicalComplex
+
+""" This file is part of the Gudhi Library - https://gudhi.inria.fr/ - which is released under MIT.
+    See file LICENSE or go to https://gudhi.inria.fr/licensing/ for full license details.
+    Author(s):       Vincent Rouvreau
+
+    Copyright (C) 2016 Inria
+
+    Modification(s):
+      - YYYY/MM Author: Description of the modification
+"""
+
+__author__ = "Vincent Rouvreau"
+__copyright__ = "Copyright (C) 2016 Inria"
+__license__ = "MIT"
+
+
+def test_empty_constructor():
+    # Try to create an empty CubicalComplex
+    cub = CubicalComplex()
+    assert cub.__is_defined() == False
+    assert cub.__is_persistence_defined() == False
+
+
+def test_non_existing_perseus_file_constructor():
+    # Try to open a non existing file
+    cub = CubicalComplex(perseus_file="pouetpouettralala.toubiloubabdou")
+    assert cub.__is_defined() == False
+    assert cub.__is_persistence_defined() == False
+
+
+def test_dimension_or_perseus_file_constructor():
+    # Create test file
+    test_file = open("CubicalOneSphere.txt", "w")
+    test_file.write("2\n3\n3\n0\n0\n0\n0\n100\n0\n0\n0\n0\n")
+    test_file.close()
+    # CubicalComplex can be constructed from dimensions and
+    # top_dimensional_cells OR from a Perseus-style file name.
+    cub = CubicalComplex(
+        dimensions=[3, 3],
+        top_dimensional_cells=[1, 2, 3, 4, 5, 6, 7, 8, 9],
+        perseus_file="CubicalOneSphere.txt",
+    )
+    assert cub.__is_defined() == False
+    assert cub.__is_persistence_defined() == False
+
+    cub = CubicalComplex(
+        top_dimensional_cells=[1, 2, 3, 4, 5, 6, 7, 8, 9],
+        perseus_file="CubicalOneSphere.txt",
+    )
+    assert cub.__is_defined() == False
+    assert cub.__is_persistence_defined() == False
+
+    cub = CubicalComplex(dimensions=[3, 3], perseus_file="CubicalOneSphere.txt")
+    assert cub.__is_defined() == False
+    assert cub.__is_persistence_defined() == False
+
+
+def test_dimension_simple_constructor():
+    cub = CubicalComplex(
+        dimensions=[3, 3], top_dimensional_cells=[1, 2, 3, 4, 5, 6, 7, 8, 9]
+    )
+    assert cub.__is_defined() == True
+    assert cub.__is_persistence_defined() == False
+    assert cub.persistence() == [(0, (1.0, float("inf")))]
+    assert cub.__is_persistence_defined() == True
+    assert cub.betti_numbers() == [1, 0, 0]
+    assert cub.persistent_betti_numbers(0, 1000) == [0, 0, 0]
+
+
+def test_user_case_simple_constructor():
+    cub = CubicalComplex(
+        dimensions=[3, 3],
+        top_dimensional_cells=[float("inf"), 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0],
+    )
+    assert cub.__is_defined() == True
+    assert cub.__is_persistence_defined() == False
+    assert cub.persistence() == [(1, (0.0, 1.0)), (0, (0.0, float("inf")))]
+    assert cub.__is_persistence_defined() == True
+    other_cub = CubicalComplex(
+        dimensions=[3, 3],
+        top_dimensional_cells=[1000.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0],
+    )
+    assert other_cub.persistence() == [(1, (0.0, 1.0)), (0, (0.0, float("inf")))]
+
+
+def test_dimension_file_constructor():
+    # Create test file
+    test_file = open("CubicalOneSphere.txt", "w")
+    test_file.write("2\n3\n3\n0\n0\n0\n0\n100\n0\n0\n0\n0\n")
+    test_file.close()
+    cub = CubicalComplex(perseus_file="CubicalOneSphere.txt")
+    assert cub.__is_defined() == True
+    assert cub.__is_persistence_defined() == False
+    assert cub.persistence() == [(1, (0.0, 100.0)), (0, (0.0, float("inf")))]
+    assert cub.__is_persistence_defined() == True
+    assert cub.betti_numbers() == [1, 0, 0]
+    assert cub.persistent_betti_numbers(0, 1000) == [1, 0, 0]