toggle quoted messageShow quoted text
Apologies for all the emails, but I found, what I believe to be an error in the R-implementation. It appears that when computing the wasserstein distance, the R version does not take the pth root of the value. I believe it computes the norm to the pth power.
On May 17, 2018, at 2:19 PM, Adam Spannaus <aspannau@...
Thanks for looking into it. If it helps, by looking at the distance between the same persistence diagram, even with a ‘inf’, gives the correct value of 0.
On another note, would you happen to know, of can point me in a direction, as to why the computation of Wasserstein distances with the same diagrams is different in the TDA R package as compared with your implementation?
Thanks for your assistance,
On May 15, 2018, at 7:42 PM, Dmitriy Morozov <dmitriy@...
Never mind. I can reproduce the problem:
import dionysus as d
dgm1 = d.Diagram([(1,2), (3,4), (1., float('inf'))])
dgm2 = d.Diagram([(0,2), (3,5), (2., float('inf'))])
I'll look into it. Let me meanwhile point out that you can construct a diagram out of a list of tuples of point coordinates. So you could always generate a list of points in the diagram and then generate a new diagram, tweaking the points however you like. So that could be a temporary workaround.
BTW, TDA package in R uses Dionysus 1, so it has nothing in common with the code used in Dionysus 2.