4 files changed, 22 insertions, 3 deletions
diff --git a/bottleneck/tests/data/test_880_A b/bottleneck/tests/data/test_880_A
new file mode 100644
index 0000000..5ad7454
--- /dev/null
+++ b/bottleneck/tests/data/test_880_A
@@ -0,0 +1,4 @@
+0.4217142857142855 inf
+0.8868393268494345 1.8
+0.45290401176484574 1.8
+0.9154285714285715 1.8
diff --git a/bottleneck/tests/data/test_880_B b/bottleneck/tests/data/test_880_B
new file mode 100644
index 0000000..22129a2
--- /dev/null
+++ b/bottleneck/tests/data/test_880_B
@@ -0,0 +1,4 @@
+0.8742857142857143 1.8
+0.46285714285714286 inf
+0.4722695736819835 1.8
+0.8474648251910301 1.8
diff --git a/bottleneck/tests/data/test_list.txt b/bottleneck/tests/data/test_list.txt
index 24b462a..91c03a4 100644
--- a/bottleneck/tests/data/test_list.txt
+++ b/bottleneck/tests/data/test_list.txt
@@ -790,3 +790,4 @@ test_876_A test_876_B 0.01 165
 test_877_A test_877_B 0.01 190
 test_878_A test_878_B 0.01 242
 test_879_A test_879_B 0.01 187
+test_880_A test_880_B 0.0 0.0411428571428574
diff --git a/bottleneck/tests/test_hera_bottleneck.cpp b/bottleneck/tests/test_hera_bottleneck.cpp
index f22e415..c74e800 100644
--- a/bottleneck/tests/test_hera_bottleneck.cpp
+++ b/bottleneck/tests/test_hera_bottleneck.cpp
@@ -159,6 +159,9 @@ TEST_CASE("infinity points", "bottleneckDistApprox")
         double d = hera::bottleneckDistApprox<>(diagram_A, diagram_B, delta);
         double corr_answer = 1.0;
         REQUIRE(  fabs(d - corr_answer) <= delta * corr_answer);
+
+        double exact_d = hera::bottleneckDistExact<>(diagram_A, diagram_B);
+        REQUIRE(exact_d == corr_answer);
     }
 
     SECTION("two points at infinity") {
@@ -559,15 +562,23 @@ TEST_CASE("file cases", "bottleneck_dist")
 
     SECTION("from file:") {
 
-        for (const auto& ts : test_params) {
+        for (auto& ts : test_params) {
             bool read_file_A = hera::readDiagramPointSet(dir_prefix + ts.file_1, diagram_A);
             bool read_file_B = hera::readDiagramPointSet(dir_prefix + ts.file_2, diagram_B);
             REQUIRE(read_file_A);
             REQUIRE(read_file_B);
 
-            double hera_answer = hera::bottleneckDistApprox(diagram_A, diagram_B, ts.delta, longest_edge, true);
+            double hera_answer;
+            if (ts.delta > 0)
+                hera_answer = hera::bottleneckDistApprox(diagram_A, diagram_B, ts.delta, longest_edge, true);
+            else
+                hera_answer = hera::bottleneckDistExact(diagram_A, diagram_B);
             std::pair<int, int> hera_le { longest_edge.first.get_user_id(), longest_edge.second.get_user_id() };
 
+            // cannot store exact answer in test_list.txt, but need to make sure that Exact is called
+            if (ts.delta == 0.0)
+                ts.delta = 0.00000001;
+
             REQUIRE((hera_answer == ts.answer or fabs(hera_answer - ts.answer) <= ts.delta * hera_answer));
             REQUIRE((ts.longest_edges.empty() or
                      std::find(ts.longest_edges.begin(), ts.longest_edges.end(), hera_le) != ts.longest_edges.end()));
@@ -604,7 +615,6 @@ TEST_CASE("file cases", "bottleneck_dist")
 //                          << hera_answer_exact << std::endl;
 //            }
             REQUIRE( (not check_longest_edge_cost or fabs(hera_le_cost - hera_answer_exact) < 0.0001 * hera_answer_exact) );
-            std::cout << ts << " PASSED " << std::endl;
         }
     }