1 files changed, 55 insertions, 0 deletions
diff --git a/geom_bottleneck/example/bottleneck_dist.cpp b/geom_bottleneck/example/bottleneck_dist.cpp
new file mode 100644
index 0000000..9a50ce2
--- /dev/null
+++ b/geom_bottleneck/example/bottleneck_dist.cpp
@@ -0,0 +1,55 @@
+#include <iomanip>
+#include "bottleneck.h"
+
+// any container of pairs of doubles can be used,
+// we use vector in this example.
+
+typedef std::vector<std::pair<double, double>>  PairVector;
+
+int main(int argc, char* argv[])
+{
+    if (argc < 3 ) {
+        std::cerr << "Usage: " << argv[0] << " file1 file2 [relative_error]. Without relative_error calculate the exact distance." << std::endl;
+        return 1;
+    }
+
+    PairVector diagramA, diagramB;
+    if (!geom_bt::readDiagramPointSet(argv[1], diagramA)) {
+        std::exit(1);
+    }
+
+    if (!geom_bt::readDiagramPointSet(argv[2], diagramB)) {
+        std::exit(1);
+    }
+
+    double res;
+    if (argc >= 4) {
+        // the third parameter is epsilon,
+        // return approximate distance (faster)
+        double approxEpsilon =  atof(argv[3]);
+        if (approxEpsilon > 0.0) {
+            res = geom_bt::bottleneckDistApprox(diagramA, diagramB, approxEpsilon);
+        } else if (approxEpsilon == 0.0) {
+            res = geom_bt::bottleneckDistExact(diagramA, diagramB);
+        } else {
+            std::cerr << "The third parameter (relative error) must be positive!" << std::endl;
+            std::exit(1);
+        }
+    } else {
+        // only filenames have been supplied, return exact distance
+        res = geom_bt::bottleneckDistExact(diagramA, diagramB);
+    }
+    std::cout << std::setprecision(15) << res << std::endl;
+
+    // Alternative could be to construct DiagramPointSet
+    // using the constructor with iterators.
+    // May be useful if the same diagram is used multiple times
+    // to avoid copying data from user's container each time.
+    
+    //geom_bt::DiagramPointSet dA(diagramA.begin(), diagramA.end());
+    //geom_bt::DiagramPointSet dB(diagramB.begin(), diagramB.end());
+    //double result1 = geom_bt::bottleneckDistExact(dA, dB);
+    //std::cout << std::setprecision(15) << result1 << std::endl;
+
+    return 0;
+}