# -*- coding:utf-8 -*-
from pylab import *
import tables

"""
TODO
- Faire tous les MPS
- Faire tous les STS

- Faire le déphasage

"""

DB = tables.openFile('db.h5')

# Get data
DATA = []
for g in DB.walkGroups():
    try:
        pset = g.paramset._v_attrs
        res = g.results._v_attrs
    except tables.NoSuchNodeError:
        pass
    else:
        new_data_item = {}

        # Interconnection rate and strength
        common = pset['Common']
        interco_strength = common['inter_conn_strength'][0][1]
        interco_rate = common['inter_conn_rate'][0][1]
        new_data_item['interco'] = {}
        new_data_item['interco']['strength'] = interco_strength
        new_data_item['interco']['rate'] = interco_rate

        # MPS
        new_data_item['MPS'] = res['MPS']

        # MPS
        new_data_item['STS'] = res['STS']

        # Add data item
        DATA.append(new_data_item)



# Plot data

# MPS
MPS = []
for simu in DATA:
    rate = simu['interco']['rate']
    strength = simu['interco']['strength']
    MPS.append((rate, strength, simu['MPS']))
MPS.sort()

# plot MPS whole
X = linspace(0, 1, 10)
Y = linspace(0, 1, 10)
MESH = zeros((len(X), len(Y)))
for indx in xrange(len(X)):
    for indy in xrange(len(Y)):
        MESH[indx][indy] = MPS[indy*len(X) + indx][2]['whole']

CS = contourf(X, Y, MESH)
colorbar(CS)  # TODO mauvaise échelle...
show()

DB.close()

Add an analysis file for 10x10 simulations
Add an analysis file for 10x10 simulations

Not finished, MPS Whole figure is quite ok.
