summaryrefslogtreecommitdiff
path: root/src/python/example
diff options
context:
space:
mode:
authorMarc Glisse <marc.glisse@inria.fr>2020-04-06 16:51:32 +0200
committerMarc Glisse <marc.glisse@inria.fr>2020-04-06 16:51:32 +0200
commit173506323471cf5175ea2b340abec63968c5cd5f (patch)
tree12236343583916b4478e736ac3850b6ca713250a /src/python/example
parent5eaca3ed69c564a6f44e6ff21ac33e2cc576bafa (diff)
Use compute_persistence in an example
Diffstat (limited to 'src/python/example')
-rwxr-xr-xsrc/python/example/alpha_rips_persistence_bottleneck_distance.py12
1 files changed, 4 insertions, 8 deletions
diff --git a/src/python/example/alpha_rips_persistence_bottleneck_distance.py b/src/python/example/alpha_rips_persistence_bottleneck_distance.py
index f156826d..3e12b0d5 100755
--- a/src/python/example/alpha_rips_persistence_bottleneck_distance.py
+++ b/src/python/example/alpha_rips_persistence_bottleneck_distance.py
@@ -5,6 +5,7 @@ import argparse
import math
import errno
import os
+import numpy as np
""" This file is part of the Gudhi Library - https://gudhi.inria.fr/ -
which is released under MIT.
@@ -56,7 +57,7 @@ with open(args.file, "r") as f:
message = "Number of simplices=" + repr(rips_stree.num_simplices())
print(message)
- rips_diag = rips_stree.persistence()
+ rips_stree.compute_persistence()
print("##############################################################")
print("AlphaComplex creation from points read in a OFF file")
@@ -72,18 +73,13 @@ with open(args.file, "r") as f:
message = "Number of simplices=" + repr(alpha_stree.num_simplices())
print(message)
- alpha_diag = alpha_stree.persistence()
+ alpha_stree.compute_persistence()
max_b_distance = 0.0
for dim in range(args.max_dimension):
# Alpha persistence values needs to be transform because filtration
# values are alpha square values
- funcs = [math.sqrt, math.sqrt]
- alpha_intervals = []
- for interval in alpha_stree.persistence_intervals_in_dimension(dim):
- alpha_intervals.append(
- map(lambda func, value: func(value), funcs, interval)
- )
+ alpha_intervals = np.sqrt(alpha_stree.persistence_intervals_in_dimension(dim))
rips_intervals = rips_stree.persistence_intervals_in_dimension(dim)
bottleneck_distance = gudhi.bottleneck_distance(