summaryrefslogtreecommitdiff
path: root/src/Kernels/example/kernel_basic_example.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'src/Kernels/example/kernel_basic_example.cpp')
-rw-r--r--src/Kernels/example/kernel_basic_example.cpp49
1 files changed, 49 insertions, 0 deletions
diff --git a/src/Kernels/example/kernel_basic_example.cpp b/src/Kernels/example/kernel_basic_example.cpp
new file mode 100644
index 00000000..0a8d83b3
--- /dev/null
+++ b/src/Kernels/example/kernel_basic_example.cpp
@@ -0,0 +1,49 @@
+/* This file is part of the Gudhi Library. The Gudhi library
+ * (Geometric Understanding in Higher Dimensions) is a generic C++
+ * library for computational topology.
+ *
+ * Authors: Francois Godi, small modifications by Pawel Dlotko
+ *
+ * Copyright (C) 2015 INRIA
+ *
+ * This program is free software: you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation, either version 3 of the License, or
+ * (at your option) any later version.
+ *
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program. If not, see <http://www.gnu.org/licenses/>.
+ */
+
+#include <gudhi/SW.h>
+//#include <gudhi/PSS.h>
+//#include <gudhi/PWG.h>
+
+#include <iostream>
+#include <vector>
+#include <utility> // for pair
+#include <limits> // for numeric_limits
+
+int main() {
+ std::vector< std::pair<double, double> > v1, v2;
+
+ v1.emplace_back(2.7, 3.7);
+ v1.emplace_back(9.6, 14.);
+ v1.emplace_back(34.2, 34.974);
+
+ v2.emplace_back(2.8, 4.45);
+ v2.emplace_back(9.5, 14.1);
+
+
+ double b1 = Gudhi::sliced_wasserstein::compute_approximate_SW (v1, v2);
+ double b2 = Gudhi::sliced_wasserstein::compute_exact_SW (v1, v2);
+
+ std::cout << "Approximate Sliced Wasserstein distance = " << b1 << std::endl;
+ std::cout << "Exact Sliced Wasserstein distance = " << b2 << std::endl;
+
+}