Rips pairwise cohomology produces output necessary for 1-cocycle parametrization
3 from pylab import scatter, show, cm, colorbar, axes, savefig
4 from itertools import izip
5 from sys import argv, exit
9 def plot(val1_fn, val2_fn):
11 with open(val1_fn) as fp:
12 for line in fp.xreadlines():
13 values1.append(float(line.split()[1]))
16 with open(val2_fn) as fp:
17 for line in fp.xreadlines():
18 values2.append(float(line.split()[1]))
20 values1 = [v % 1. for v in values1]
21 values2 = [v % 1. for v in values2]
22 print min(values1), max(values2), min(values1), min(values2)
24 scatter(values1, values2)
25 axes().set_aspect('equal') 28 if __name__ == '__main__':
30 print "Usage: %s VALUES1 VALUES2" % argv[0]
35 plot(val1_fn, val2_fn)