diff options
author | Marc Glisse <marc.glisse@inria.fr> | 2020-04-24 18:13:20 +0200 |
---|---|---|
committer | GitHub <noreply@github.com> | 2020-04-24 18:13:20 +0200 |
commit | d895af4897cb4169486197dec925a8a16d4de9cf (patch) | |
tree | 500c55f4b47a33b933f03ffad16b8c26df121830 /src/python/example | |
parent | e3f276ab5b7503ba7ce278fffbf73ebe66d6351c (diff) | |
parent | 031e6879f94503e5250c005f8cb71e581799d2f3 (diff) |
Merge pull request #277 from mglisse/compute_persistence
compute_persistence in python
Diffstat (limited to 'src/python/example')
-rwxr-xr-x | src/python/example/alpha_rips_persistence_bottleneck_distance.py | 12 |
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( |