summaryrefslogtreecommitdiff
path: root/src
diff options
context:
space:
mode:
authorvrouvrea <vrouvrea@636b058d-ea47-450e-bf9e-a15bfbe3eedb>2016-04-15 13:11:36 +0000
committervrouvrea <vrouvrea@636b058d-ea47-450e-bf9e-a15bfbe3eedb>2016-04-15 13:11:36 +0000
commit6adf6b55ebd348f84b63cdc9d822df2be587d7e3 (patch)
tree90b2255f4209024033b3e16d7c62b964b9560135 /src
parentae9eb11ebccacdcef5b91cd4ca30cc6ca5984bf0 (diff)
parent91e58d98e81963ac714ec6c9508c5bc15945f20b (diff)
periodic_alpha_complex_3d_persistence.cpp bug fix
README values in consequence git-svn-id: svn+ssh://scm.gforge.inria.fr/svnroot/gudhi/trunk@1126 636b058d-ea47-450e-bf9e-a15bfbe3eedb Former-commit-id: ce435789e89e1cefe1023a6d73b8606989fca3ba
Diffstat (limited to 'src')
-rw-r--r--src/Persistent_cohomology/example/README19
-rw-r--r--src/Persistent_cohomology/example/periodic_alpha_complex_3d_persistence.cpp4
2 files changed, 12 insertions, 11 deletions
diff --git a/src/Persistent_cohomology/example/README b/src/Persistent_cohomology/example/README
index 92b80c76..f490ed9a 100644
--- a/src/Persistent_cohomology/example/README
+++ b/src/Persistent_cohomology/example/README
@@ -105,18 +105,19 @@ N.B.: alpha_complex_persistence accepts OFF files in d-Dimension.
3) 3D periodic special case
---------------------------
-./periodic_alpha_complex_3d_persistence ../../data/points/grid_10_10_10_in_0_1.off 2 0.0
+./periodic_alpha_complex_3d_persistence ../../data/points/grid_10_10_10_in_0_1.off 3 1.0
output:
Periodic Delaunay computed.
Simplex_tree dim: 3
-2 0 0.0866025 inf
-2 1 0.0866025 inf
-2 1 0.0866025 inf
-2 1 0.0866025 inf
-2 2 0.0866025 inf
-2 2 0.0866025 inf
-2 2 0.0866025 inf
+3 0 0 inf
+3 1 0.05 inf
+3 1 0.05 inf
+3 1 0.05 inf
+3 2 0.0707107 inf
+3 2 0.0707107 inf
+3 2 0.0707107 inf
+3 3 0.0866025 inf
N.B.: periodic_alpha_complex_3d_persistence accepts only OFF files in 3D dimension. In this example, the periodic cube
is hard coded to { x = [0,1]; y = [0,1]; z = [0,1] }
@@ -139,4 +140,4 @@ output:
2 0 0 inf
2 1 0 inf
-Here we retrieve the 2 entities {0,1,2,3} and {4} (Betti numbers[0] = 2) and the hole in {0,1,3} (Betti numbers[1] = 1) \ No newline at end of file
+Here we retrieve the 2 entities {0,1,2,3} and {4} (Betti numbers[0] = 2) and the hole in {0,1,3} (Betti numbers[1] = 1)
diff --git a/src/Persistent_cohomology/example/periodic_alpha_complex_3d_persistence.cpp b/src/Persistent_cohomology/example/periodic_alpha_complex_3d_persistence.cpp
index 42e6f0c6..750ebb1a 100644
--- a/src/Persistent_cohomology/example/periodic_alpha_complex_3d_persistence.cpp
+++ b/src/Persistent_cohomology/example/periodic_alpha_complex_3d_persistence.cpp
@@ -173,7 +173,7 @@ int main(int argc, char * const argv[]) {
// alpha shape construction from points. CGAL has a strange behavior in REGULARIZED mode. This is the default mode
// Maybe need to set it to GENERAL mode
- Alpha_shape_3 as(pdt, Alpha_shape_3::GENERAL);
+ Alpha_shape_3 as(pdt, 0, Alpha_shape_3::GENERAL);
// filtration with alpha values from alpha shape
std::vector<Object> the_objects;
@@ -291,7 +291,7 @@ int main(int argc, char * const argv[]) {
std::cout << "Simplex_tree dim: " << simplex_tree.dimension() << std::endl;
// Compute the persistence diagram of the complex
- Persistent_cohomology pcoh(simplex_tree);
+ Persistent_cohomology pcoh(simplex_tree, true);
// initializes the coefficient field for homology
pcoh.init_coefficients(coeff_field_characteristic);