2 files changed, 39 insertions, 1 deletions
diff --git a/src/python/test/test_cubical_complex.py b/src/python/test/test_cubical_complex.py
index 5c59db8f..d0e4e9e8 100755
--- a/src/python/test/test_cubical_complex.py
+++ b/src/python/test/test_cubical_complex.py
@@ -157,3 +157,20 @@ def test_cubical_generators():
     assert np.array_equal(g[0][0], np.empty(shape=[0,2]))
     assert np.array_equal(g[0][1], np.array([[7, 4]]))
     assert np.array_equal(g[1][0], np.array([8]))
+
+def test_cubical_cofaces_of_persistence_pairs_when_pd_has_no_paired_birth_and_death():
+    cubCpx = CubicalComplex(dimensions=[1,2], top_dimensional_cells=[0.0, 1.0])
+    Diag = cubCpx.persistence(homology_coeff_field=2, min_persistence=0)
+    pairs = cubCpx.cofaces_of_persistence_pairs()
+    assert pairs[0] == []
+    assert np.array_equal(pairs[1][0], np.array([0]))
+
+def test_periodic_cofaces_of_persistence_pairs_when_pd_has_no_paired_birth_and_death():
+    perCubCpx = PeriodicCubicalComplex(dimensions=[1,2], top_dimensional_cells=[0.0, 1.0],
+                                       periodic_dimensions=[True, True])
+    Diag = perCubCpx.persistence(homology_coeff_field=2, min_persistence=0)
+    pairs = perCubCpx.cofaces_of_persistence_pairs()
+    assert pairs[0] == []
+    assert np.array_equal(pairs[1][0], np.array([0]))
+    assert np.array_equal(pairs[1][1], np.array([0, 1]))
+    assert np.array_equal(pairs[1][2], np.array([1]))
diff --git a/src/python/test/test_dtm.py b/src/python/test/test_dtm.py
index bff4c267..0a52279e 100755
--- a/src/python/test/test_dtm.py
+++ b/src/python/test/test_dtm.py
@@ -8,10 +8,11 @@
       - YYYY/MM Author: Description of the modification
 """
 
-from gudhi.point_cloud.dtm import DistanceToMeasure
+from gudhi.point_cloud.dtm import DistanceToMeasure, DTMDensity
 import numpy
 import pytest
 import torch
+import math
 
 
 def test_dtm_compare_euclidean():
@@ -66,3 +67,23 @@ def test_dtm_precomputed():
     dtm = DistanceToMeasure(2, q=2, metric="neighbors")
     r = dtm.fit_transform(dist)
     assert r == pytest.approx([2.0, 0.707, 3.5355], rel=0.01)
+
+
+def test_density_normalized():
+    sample = numpy.random.normal(0, 1, (1000000, 2))
+    queries = numpy.array([[0.0, 0.0], [-0.5, 0.7], [0.4, 1.7]])
+    expected = numpy.exp(-(queries ** 2).sum(-1) / 2) / (2 * math.pi)
+    estimated = DTMDensity(k=150, normalize=True).fit(sample).transform(queries)
+    assert estimated == pytest.approx(expected, rel=0.4)
+
+
+def test_density():
+    distances = [[0, 1, 10], [2, 0, 30], [1, 3, 5]]
+    density = DTMDensity(k=2, metric="neighbors", dim=1).fit_transform(distances)
+    expected = numpy.array([2.0, 1.0, 0.5])
+    assert density == pytest.approx(expected)
+    distances = [[0, 1], [2, 0], [1, 3]]
+    density = DTMDensity(metric="neighbors", dim=1).fit_transform(distances)
+    assert density == pytest.approx(expected)
+    density = DTMDensity(weights=[0.5, 0.5], metric="neighbors", dim=1).fit_transform(distances)
+    assert density == pytest.approx(expected)