diff options
Diffstat (limited to 'geom_bottleneck/example/bottleneck_dist.cpp')
| -rw-r--r-- | geom_bottleneck/example/bottleneck_dist.cpp | 11 |
1 files changed, 10 insertions, 1 deletions
diff --git a/geom_bottleneck/example/bottleneck_dist.cpp b/geom_bottleneck/example/bottleneck_dist.cpp index 655de44..b27bdbb 100644 --- a/geom_bottleneck/example/bottleneck_dist.cpp +++ b/geom_bottleneck/example/bottleneck_dist.cpp @@ -6,6 +6,11 @@ typedef std::vector<std::pair<double, double>> PairVector; +// estimate initial guess on sampled diagram? +constexpr bool useSamplingHeur = false; +// if diagrams contain fewer points, don't use heuristic +constexpr int heurThreshold = 30000; + int main(int argc, char* argv[]) { if (argc < 3 ) { @@ -29,7 +34,11 @@ int main(int argc, char* argv[]) // return approximate distance (faster) double approxEpsilon = atof(argv[3]); if (approxEpsilon > 0.0) { - res = geom_bt::bottleneckDistApprox(diagramA, diagramB, approxEpsilon); + if (useSamplingHeur && diagramA.size() > heurThreshold && diagramB.size() > heurThreshold) { + res = geom_bt::bottleneckDistApproxHeur(diagramA, diagramB, approxEpsilon); + } else { + res = geom_bt::bottleneckDistApprox(diagramA, diagramB, approxEpsilon); + } } else if (approxEpsilon == 0.0) { res = geom_bt::bottleneckDistExact(diagramA, diagramB, decPrecision); } else { |