diff options
author | ROUVREAU Vincent <vincent.rouvreau@inria.fr> | 2021-04-21 10:10:32 +0200 |
---|---|---|
committer | ROUVREAU Vincent <vincent.rouvreau@inria.fr> | 2021-04-21 10:10:32 +0200 |
commit | fb2aca88ba373868fc3e67683258b1b98467d10d (patch) | |
tree | 842a59d76bb3aba5d6af1602100f6218e4bb6036 /src/python/example/alpha_complex_diagram_persistence_from_off_file_example.py | |
parent | 29fa56b6768d7314fcddf39d7681aa6a82417b4e (diff) | |
parent | 73fa5b763a53179444304ccbe0583b616403bb0a (diff) |
Merge branch 'master' into coxeter_integration
Diffstat (limited to 'src/python/example/alpha_complex_diagram_persistence_from_off_file_example.py')
-rwxr-xr-x | src/python/example/alpha_complex_diagram_persistence_from_off_file_example.py | 23 |
1 files changed, 11 insertions, 12 deletions
diff --git a/src/python/example/alpha_complex_diagram_persistence_from_off_file_example.py b/src/python/example/alpha_complex_diagram_persistence_from_off_file_example.py index 1e0273b3..fe03be31 100755 --- a/src/python/example/alpha_complex_diagram_persistence_from_off_file_example.py +++ b/src/python/example/alpha_complex_diagram_persistence_from_off_file_example.py @@ -25,12 +25,12 @@ parser = argparse.ArgumentParser( description="AlphaComplex creation from " "points read in a OFF file.", epilog="Example: " "example/alpha_complex_diagram_persistence_from_off_file_example.py " - "-f ../data/points/tore3D_300.off -a 0.6" + "-f ../data/points/tore3D_300.off" "- Constructs a alpha complex with the " "points from the given OFF file.", ) parser.add_argument("-f", "--file", type=str, required=True) -parser.add_argument("-a", "--max_alpha_square", type=float, default=0.5) +parser.add_argument("-a", "--max_alpha_square", type=float, required=False) parser.add_argument("-b", "--band", type=float, default=0.0) parser.add_argument( "--no-diagram", @@ -47,21 +47,20 @@ with open(args.file, "r") as f: print("##############################################################") print("AlphaComplex creation from points read in a OFF file") - message = "AlphaComplex with max_edge_length=" + repr(args.max_alpha_square) - print(message) - alpha_complex = gudhi.AlphaComplex(off_file=args.file) - simplex_tree = alpha_complex.create_simplex_tree( - max_alpha_square=args.max_alpha_square - ) + if args.max_alpha_square is not None: + print("with max_edge_length=", args.max_alpha_square) + simplex_tree = alpha_complex.create_simplex_tree( + max_alpha_square=args.max_alpha_square + ) + else: + simplex_tree = alpha_complex.create_simplex_tree() - message = "Number of simplices=" + repr(simplex_tree.num_simplices()) - print(message) + print("Number of simplices=", simplex_tree.num_simplices()) diag = simplex_tree.persistence() - print("betti_numbers()=") - print(simplex_tree.betti_numbers()) + print("betti_numbers()=", simplex_tree.betti_numbers()) if args.no_diagram == False: import matplotlib.pyplot as plot |