def generate(reference="SimpleNet", scale=1, format='xml'):
population_size = scale_pop_size(3, scale)
nml_doc, network = oc.generate_network(reference)
oc.add_cell_and_channels(nml_doc, 'izhikevich/RS.cell.nml', 'RS')
pop = oc.add_population_in_rectangular_region(network, 'RS_pop', 'RS',
population_size, 0, 0, 0,
100, 100, 100)
syn = oc.add_exp_two_syn(nml_doc,
id="syn0",
gbase="2nS",
erev="0mV",
tau_rise="0.5ms",
tau_decay="10ms")
pfs = oc.add_poisson_firing_synapse(nml_doc,
id="poissonFiringSyn",
average_rate="50 Hz",
synapse_id=syn.id)
oc.add_inputs_to_population(network, "Stim0", pop, pfs.id, all_cells=True)
nml_file_name = '%s.net.nml' % network.id
oc.save_network(nml_doc,
nml_file_name,
validate=(format == 'xml'),
format=format)
if format == 'xml':
oc.generate_lems_simulation(nml_doc,
network,
nml_file_name,
duration=500,
dt=0.025)
