#np.savetxt(save_dir+'/seed_'+str(trial)+'.txt', np.array([seed]))
prng = Random()
prng.seed(seed)
# generate random models
for i in range(n_models):
# create a model randomly from all channels
weights_model = problem.generator(prng, None)
for j, dt in enumerate(dts):
print 'run '+'model: ' + str(i) + ' dt: ' +str(dt)
# change dt
data_newdt = change_dt(dt, problem.data)
problem.simulation_params = extract_simulation_params(data_newdt)
# run simulation
problem.update_cell(weights_model)
currents = [problem.cell.soma.record_from(channel_list[k], 'i'+ion_list[k], pos=.5) for k in range(len(channel_list))]
v_newdt, t_newdt = run_simulation(problem.cell, **problem.simulation_params)
#pl.figure()
#pl.plot(t_newdt, v_newdt)
#pl.show()
# compute parameter
dvdt_newdt = np.concatenate((np.array([(v_newdt[1]-v_newdt[0])/dt]), np.diff(v_newdt)/dt))
i_newdt = data_newdt.i.values
celsius = problem.simulation_params['celsius']
