main.py

import network

import argparse
import datetime
import scipy.stats as stat
import matplotlib.pyplot as plot

def main():
    k = 100
    time_max = 1000

    spikes = []
    times = []

    fig = plot.figure(figsize=(6,6))

    ax1 = fig.add_subplot(111)

    t0 = datetime.datetime.now()
    net = network.Network(7 * k, 3 * k)
    t1 = datetime.datetime.now()

    print('configured network in:', t1 - t0)

    t0 = datetime.datetime.now()

    for t in range(time_max):
        net.input[0 : net.numEx] = 5.0 * stat.norm.rvs(size=net.numEx)
        net.input[net.numEx : ] = 2.0 * stat.norm.rvs(size=net.numIn)

        tmp = net.update()
        for n in tmp:
            spikes.append(n)
            times.append(t)

    t1 = datetime.datetime.now()

    print('simulated', time_max, 'steps in: ', t1 - t0)

    ax1.plot(times, spikes, ',k')

    xl, xr = ax1.get_xlim()
    yb, yt = ax1.get_ylim()

    ax1.set_aspect(abs((xr - xl) / (yb - yt))  * 1.0)
    ax1.axhline(color='r', y=net.numEx - 0.5, xmax=time_max)

    plot.show()

if __name__ == "__main__":
    main()