adaptive_simulation.py

import nest
import numpy
import pylab

for J in range(0,110,10):

    nest.SetKernelStatus({"local_num_threads":8})

    neuron_population = 100
    simulation_time = 500.0

    dict_params = {"V_peak" : 0.0
    , "V_reset" : -70.0
    , "t_ref" : 2.0
    , "g_L" :  30.0
    , "C_m" : 281.0
    , "E_ex" : 0.0
    , "E_in" : -85.0
    , "E_L" : -70.0
    , "Delta_T" : 2.0
    , "tau_w" : 100.0
    , "a" : 0.0
    , "b" : 70.0
    , "V_th" : -55.0
    , "tau_syn_ex" : 0.2
    , "tau_syn_in" : 2.0
    , "I_e" : 200.0
    , "w": 0.0}

    neurons = nest.Create("aeif_cond_exp", 100)
    noise = nest.Create("poisson_generator")

    nest.SetStatus(neurons, params=dict_params)

    nest.SetStatus(noise, {"rate": 100.0})

    for neuron in neurons:
        nest.SetStatus([neuron], {"V_m": dict_params["E_L"]+(dict_params["V_th"]-dict_params["E_L"])*numpy.random.rand()})

    spike_detector = nest.Create("spike_detector")

    for n in neurons:
        nest.Connect([n], spike_detector)

    nest.Connect(noise, neurons)

    K = 10
    d = 1.0

    conn_dict = {"rule": "fixed_indegree", "indegree": K}
    syn_dict = {"delay": d, "weight": float(J)}

    nest.Connect(neurons, neurons, conn_dict, syn_dict)

    for neuron in neurons:
        nest.SetStatus([neuron], {"V_m": dict_params["E_L"]+(dict_params["V_th"]-dict_params["E_L"])*numpy.random.rand()})

    for I in range(0, 510, 10):

        nest.SetStatus(neurons, params={"I_e": float(I)})

        nest.Simulate(simulation_time)

        dSD = nest.GetStatus(spike_detector, keys='events')[0]
        evs = dSD["senders"]
        ts = dSD["times"]

        nest.ResumeSimulation()

        pylab.figure("Increasing Spike Raster J: " + str(J))
        pylab.plot(ts, evs, "b.")
        pylab.xlabel("Time: ms")
        pylab.ylabel("Neuron ID")
        pylab.title("Increasing Spike Raster J: " + str(J))
        #pylab.savefig("Decreasing Spike Raster J: " + str(J))

    nest.ResetKernel()
    nest.ResetNetwork()

pylab.show()

"""open_file = open("adaptive_decrease.txt", "w")
open_file.write(time_string)
open_file.close()"""