Changes for page Code description

Last modified by galluzziandrea on 2022/06/20 12:33
From version 8.1
edited by galluzziandrea
on 2021/12/09 14:58
Change comment: There is no comment for this version
To version 29.1
edited by galluzziandrea
on 2022/06/20 12:33
Change comment: There is no comment for this version

Summary

Details

Page properties
Content
... ... @@ -235,6 +235,8 @@
235 235  # [.....],[],...]
236 236  {{/code}}
237 237  
238 +=== [[image:image-20220127173048-1.png||height="508" width="948"]] ===
239 +
238 238  === Defining general and nest.kernel parameters ===
239 239  
240 240  {{code language="python"}}
... ... @@ -323,6 +323,207 @@
323 323  endbuild = time.time()
324 324  {{/code}}
325 325  
326 -=== Results ===
328 +=== [[image:image-20220127165908-2.png||height="659" width="1149"]] ===
327 327  
330 +=== Connecting the network nodes: neuronal populations, Poisson processes and spike detectors ===
331 +
332 +{{code language="python"}}
333 +#############################------------------------------------------------------------------------
334 +print("Connecting ")
335 +#############################------------------------------------------------------------------------
336 +
337 +startconnect = time.time()
338 +Connessioni=[]
339 +Medie=[]
340 +
341 +#create and define the connections between the populations of neurons and the poisson generators
342 +#and between the populations of neurons and the spike detectors with the parameters extracted from the.ini files
343 +
344 +for i in range(0,int(InfoBuild[0])):
345 + nest.Connect(NoisePop[i], NeuronPop[i], syn_spec={'synapse_model': 'static_synapse_hpc',
346 + 'delay': dt,
347 + 'weight': nest.math.redraw(nest.random.normal(mean=float(InfoConnectNoise[i+1][0]),
348 + std=(float(InfoConnectNoise[i+1][1])*float(InfoConnectNoise[i+1][0]))),
349 + min=0., max=float('Inf'))
350 + })
351 + nest.Connect(NeuronPop[i][:int(InfoBuild[i+1][0])], DetectorPop[i], syn_spec={"weight": 1.0, "delay": dt})
352 +
353 +#create and define the connections between the populations of neurons with the parameters extracted from the.ini files
354 +
355 +for i in range(0,len(InfoConnectPop[1:])):
356 +
357 + conn=nest.Connect(NeuronPop[int(InfoConnectPop[i+1][1])], NeuronPop[int(InfoConnectPop[i+1][0])],
358 + {'rule': 'pairwise_bernoulli',
359 + 'p':float(InfoConnectPop[i+1][2]) },
360 + syn_spec={'synapse_model': 'static_synapse_hpc',
361 + 'delay':nest.math.redraw(nest.random.exponential(beta=float(1./(2.99573227355/(float(InfoConnectPop[i+1][4])-float(InfoConnectPop[i+1][3]))))),
362 + min= numpy.max([dt,float(1./float(InfoConnectPop[i+1][4]))]),
363 + max= float(1./(float(InfoConnectPop[i+1][3])-dt/2))),
364 +
365 + 'weight':nest.random.normal(mean=float(InfoConnectPop[i+1][6]),
366 + std=math.fabs(float(InfoConnectPop[i+1][6])*float(InfoConnectPop[i+1][7])))})
367 +
368 +
369 +endconnect = time.time()
370 +{{/code}}
371 +
372 +=== ===
373 +
374 +=== ===
375 +
376 +=== [[image:image-20220127170722-1.png]] ===
377 +
378 +=== Simulating: neuronal time evolution. ===
379 +
380 +=== ===
381 +
382 +{{code language="python"}}
383 + #############################------------------------------------------------------------------------
384 + print("Simulating")
385 + #############################------------------------------------------------------------------------
386 + ###################################################################################################################################################################
387 + if Salva:
388 + print("I m going to save the data")
389 + #x=str(iterazioni)
390 + f = open(FileName,"w")
391 + if len(InfoProtocol):
392 + print("I m going to split the simulation")
393 + tempo=0
394 + for contatore in range(0,len(InfoProtocol)):
395 + appoggio1=int((tempo+InfoProtocol[contatore][0])/1000.)
396 + appoggio2=int(tempo/1000.)
397 + appoggio3=tempo+InfoProtocol[contatore][0]
398 + if (appoggio1-appoggio2)>=1:
399 + T1=(1+appoggio2)*1000-tempo
400 + nest.Simulate(T1)
401 + #Save the Data!!!!
402 + ###########################################################
403 + Equilibri=[]
404 + for i in range(0,int(InfoBuild[0])):
405 + Equilibri.append([])
406 + a=nest.GetStatus(DetectorPop[i])[0]["events"]["times"]
407 + if len(a)>0:
408 + Trange=(1000*int(numpy.min(a)/1000.),1000*int(numpy.min(a)/1000.)+1000)
409 + hist,Tbin=numpy.histogram(a,200,(Trange[0],Trange[1]))
410 + Equilibri[i]=hist*1000./(5.*int(InfoBuild[i+1][0]))
411 + else:
412 + Trange=(1000*int(tempo/1000.),1000*int(tempo/1000.)+1000)
413 + hist=numpy.zeros(200)
414 + Tbin=numpy.linspace(Trange[0],Trange[1],num=201)
415 + Equilibri[i]=hist
416 + nest.SetStatus(DetectorPop[i],{'n_events':0})
417 + for j in range(0,len(hist)):
418 + f.write(str(Tbin[j])+" ")
419 + for i in range(0,int(InfoBuild[0])):
420 + f.write(str(Equilibri[i][j])+" ")
421 + f.write("\n ")
422 + ###########################################################
423 + tempo=tempo+T1
424 + for contatore2 in range(1,(appoggio1-appoggio2)):
425 + nest.Simulate(1000.)
426 + #Save the Data!!!!
427 + ###########################################################
428 + Equilibri=[]
429 + for i in range(0,int(InfoBuild[0])):
430 + Equilibri.append([])
431 + a=nest.GetStatus(DetectorPop[i])[0]["events"]["times"]
432 + if len(a)>0:
433 + Trange=(1000*int(numpy.min(a)/1000.),1000*int(numpy.min(a)/1000.)+1000)
434 + hist,Tbin=numpy.histogram(a,200,(Trange[0],Trange[1]))
435 + Equilibri[i]=hist*1000./(5.*int(InfoBuild[i+1][0]))
436 + else:
437 + Trange=(1000*int(tempo/1000.),1000*int(tempo/1000.)+1000)
438 + hist=numpy.zeros(200)
439 + Tbin=numpy.linspace(Trange[0],Trange[1],num=201)
440 + Equilibri[i]=hist
441 + nest.SetStatus(DetectorPop[i],{'n_events':0})
442 + for j in range(0,len(hist)):
443 + f.write(str(Tbin[j])+" ")
444 + for i in range(0,int(InfoBuild[0])):
445 + f.write(str(Equilibri[i][j])+" ")
446 + f.write("\n ")
447 + tempo=tempo+1000.
448 + T2=appoggio3-tempo
449 + nest.Simulate(T2);
450 + tempo=tempo+T2;
451 + else:
452 + nest.Simulate(InfoProtocol[contatore][0])
453 + temp=InfoProtocol[contatore][0]
454 + tempo=tempo+temp
455 + if InfoProtocol[contatore][2]==4:
456 + nest.SetStatus(NoisePop[InfoProtocol[contatore][1]],params={"rate": float(InfoBuild[1+InfoProtocol[contatore][1]][2]*InfoProtocol[contatore][3])})
457 + if InfoProtocol[contatore][2]==12:
458 + nest.SetStatus(NeuronPop[InfoProtocol[contatore][1]], params={"b": float(InfoProtocol[contatore][3])})
459 + else:
460 + nest.Simulate(simtime)
461 + tempo=simtime
462 + if (simtime-tempo)>0.:
463 + nest.Simulate(simtime-tempo)
464 +
465 +
466 + endsimulate = time.time()
467 + f.close()
468 + else:
469 + if len(InfoProtocol):
470 + tempo=0
471 + for contatore in range(0,len(InfoProtocol)):
472 + nest.Simulate(InfoProtocol[contatore][0])
473 + temp=InfoProtocol[contatore][0]
474 + tempo=tempo+temp
475 + if InfoProtocol[contatore][2]==4:
476 + nest.SetStatus(NoisePop[InfoProtocol[contatore][1]],params={"rate": float(InfoBuild[1+InfoProtocol[contatore][1]][2]*InfoProtocol[contatore][3])})
477 + #print "Population:", InfoProtocol[contatore][1] ,";Parameter:", InfoProtocol[contatore][2] ,"; Value: ",InfoProtocol[contatore][3]
478 + if InfoProtocol[contatore][2]==12:
479 + nest.SetStatus(NeuronPop[InfoProtocol[contatore][1]], params={"b": float(InfoProtocol[contatore][3])})
480 + #print "Population:", InfoProtocol[contatore][1] ,";Parameter:", InfoProtocol[contatore][2] ,"; Value: ",InfoProtocol[contatore][3]
481 +
482 + else:
483 + nest.Simulate(simtime)
484 + tempo=simtime
485 + if (simtime-tempo)>0.:
486 + nest.Simulate(simtime-tempo)
487 + endsimulate = time.time()
488 +
489 +
490 + ###################################################################################################################################################################
491 +
492 + #############################------------------------------------------------------------------------
493 + #print some information from the simulation
494 + #############################------------------------------------------------------------------------
495 +
496 + num_synapses = nest.GetDefaults('static_synapse_hpc')["num_connections"]
497 + build_time = endbuild - startbuild
498 + connect_time = endconnect - startconnect
499 + sim_time = endsimulate - endconnect
500 +
501 + N_neurons=0
502 + for i in range(0,int(InfoBuild[0])):
503 + N_neurons=N_neurons+int(InfoBuild[i+1][0])
504 +
505 + print(" Network simulation (Python) neuron type:",InfoPerseo[0])
506 + print("Number of neurons : {0}".format(N_neurons))
507 + print("Number of synapses: {0}".format(num_synapses))
508 + print("Building time : %.2f s" % build_time)
509 + print("Connecting time : %.2f s" % connect_time)
510 + print("Simulation time : %.2f s" % sim_time)
511 +
512 +Fine=time.time()
513 +print ("Total Simulation time : %.2f s" % (Fine-Inizio))
514 +{{/code}}
515 +
516 +=== ===
517 +
518 +=== ===
519 +
520 +=== [[image:image-20220127171242-1.png]] ===
521 +
522 +=== Results: ===
523 +
524 +[[the output of this simulation is...>>https://drive.ebrains.eu/smart-link/215f8213-17e3-468b-b573-e6eaf49d315e/]]
525 +
526 +
527 +
528 +
529 +
530 +
328 328  ==== ====
CorticalField_t=80_DeepSpontPlanar.mp4
Author
... ... @@ -1,0 +1,1 @@
1 +XWiki.galluzziandrea
Size
... ... @@ -1,0 +1,1 @@
1 +3.4 MB
Content
image-20220127165822-1.png
Author
... ... @@ -1,0 +1,1 @@
1 +XWiki.galluzziandrea
Size
... ... @@ -1,0 +1,1 @@
1 +237.6 KB
Content
image-20220127165908-2.png
Author
... ... @@ -1,0 +1,1 @@
1 +XWiki.galluzziandrea
Size
... ... @@ -1,0 +1,1 @@
1 +106.6 KB
Content
image-20220127170104-1.png
Author
... ... @@ -1,0 +1,1 @@
1 +XWiki.galluzziandrea
Size
... ... @@ -1,0 +1,1 @@
1 +149.6 KB
Content
image-20220127170155-2.png
Author
... ... @@ -1,0 +1,1 @@
1 +XWiki.galluzziandrea
Size
... ... @@ -1,0 +1,1 @@
1 +280.0 KB
Content
image-20220127170722-1.png
Author
... ... @@ -1,0 +1,1 @@
1 +XWiki.galluzziandrea
Size
... ... @@ -1,0 +1,1 @@
1 +220.5 KB
Content
image-20220127171242-1.png
Author
... ... @@ -1,0 +1,1 @@
1 +XWiki.galluzziandrea
Size
... ... @@ -1,0 +1,1 @@
1 +301.3 KB
Content
image-20220127173048-1.png
Author
... ... @@ -1,0 +1,1 @@
1 +XWiki.galluzziandrea
Size
... ... @@ -1,0 +1,1 @@
1 +71.5 KB
Content