tools/plot-values/scatter.py
author Dmitriy Morozov <dmitriy@mrzv.org>
Thu Jul 09 02:42:47 2009 -0700 (2 years ago)
branchdev
changeset 138 96030f8d2f2c
child 141cda0b85ffc50
permissions -rwxr-xr-x
Rips pairwise cohomology produces output necessary for 1-cocycle parametrization
     1 #!/usr/bin/env python
     2 
     3 from    pylab           import scatter, show, cm, colorbar, axes, savefig
     4 from    itertools       import izip
     5 from    sys             import argv, exit
     6 import  os.path         as     osp
     7 
     8 
     9 def plot(val1_fn, val2_fn):
    10     values1 = []
    11     with open(val1_fn) as fp:
    12         for line in fp.xreadlines():
    13             values1.append(float(line.split()[1]))
    14     
    15     values2 = []
    16     with open(val2_fn) as fp:
    17         for line in fp.xreadlines():
    18             values2.append(float(line.split()[1]))
    19     
    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)
    23 
    24     scatter(values1, values2)
    25     axes().set_aspect('equal')
    26     show()
    27 
    28 if __name__ == '__main__':
    29     if len(argv) < 3:
    30         print "Usage: %s VALUES1 VALUES2" % argv[0]
    31         exit()
    32 
    33     val1_fn = argv[1]
    34     val2_fn  = argv[2]
    35     plot(val1_fn, val2_fn)