Date
1  5 of 5
Computing extended persistence in dionysus2
cai.507@...
Hello Dmitriy,
I am trying to compute extended persistence for graphs using dionysus2. I found in dionysus1 there is a file called 08extendedpersistence.py in the accomapnying examples . But I am not able to run it using Dionysus 2.0.6 dev0. ( from dionysus.viewer import * ImportError: No module named viewer) I assume that if extended persistence can be computed in first version of Dionysus, it should also be available in second version. I did a little search but didn't find an example. (I found relative persistence homology although). I was wondering is computing extended persistent available in Dionysus2? Best, Chen


Dmitriy Morozov
There is nothing special about extended persistence. One can compute it by coning off the domain and assigning appropriate values to all the simplices. That's how the example you are referring to works. You can either translate it to Dionysus 2, or, even better, implement the coning from scratch. If you do, I hope you contribute it back to the codebase. I don't have a readymade example. There is a C++ implementation of extended persistence in Dionysus 2, but it doesn't sound like that's what you are looking for.
On Wed, Oct 3, 2018 at 10:12 AM <cai.507@...> wrote: Hello Dmitriy,


cai.507@...
Thanks a lot for response. I tried to translate from your example to Dionysus2 and I almost get it but because API for init_diagrams in two versions are different, I fail on the last step.
I try to capture a simple cycle here. I expect the solution should be [3, 0.5] . If I understand coning strategy correctly, The filtration order seems to be right ( <1> 0 <0> 0.5 <1> 1 <3> 2 <2> 3 <0,1> 1 <0,3> 2 <1,2> 3 <2,3> 3 <1,2> 0 <1,3> 0 <1,1> 0 <1,0> 0 <1,2,3> 0 <1,1,2> 0 <1,0,1> 0 <1,0,3> 0) But the answer it gives is (3,0). I think it's due to I am lacking an analog of eval_ep function in your code. But I looked at the API for init_dionysus in Dionysus2, and didn't find where to put this eval_ep(value) in. Hopefully this makes sense to you. Thanks a lot for your time! Your code (I am having trouble installing Dionysus1 so didn't run the code yet): # Extended persistence My code: import numpy as np


Dmitriy Morozov
You can't define extended persistence for an arbitrary filtration. You need a PL function, i.e., you need values on the vertices. Here's a simple example: import numpy as np import dionysus as d w = 1 values = [1,2,3,4] simplices = [[0], [1], [2], [3], [0,1], [0,3], [1,2], [2,3]] up_simplices = [d.Simplex(s, max(values[v] for v in s)) for s in simplices] down_simplices = [d.Simplex(s + [w], min(values[v] for v in s)) for s in simplices] up_simplices.sort() down_simplices.sort(reverse = True) f = d.Filtration([d.Simplex([w], float('inf'))] + up_simplices + down_simplices) for s in f: print(s) m = d.homology_persistence(f) dgms = d.init_diagrams(m, f) for i, dgm in enumerate(dgms): print("Dimension", i) for pt in dgm: print(pt) This produces the following output. Pair (inf,inf) is there for technical reasons and should be ignored. Without it, you get the extended persistence diagrams. <1> inf <0> 1 <1> 2 <2> 3 <3> 4 <0,1> 2 <0,3> 4 <1,2> 3 <2,3> 4 <1,2,3> 3 <1,1,2> 2 <1,0,3> 1 <1,0,1> 1 <1,3> 4 <1,2> 3 <1,1> 2 <1,0> 1 Dimension 0 (inf,inf) (1,4) Dimension 1 (4,1) Dimension 2
On Wed, Oct 3, 2018 at 8:37 PM <cai.507@...> wrote: Thanks a lot for response. I tried to translate from your example to Dionysus2 and I almost get it but because API for init_diagrams in two versions are different, I fail on the last step.


cai.507@...
It's working perfectly. Many Thanks!
Best, Chen

