parameters['ISYN_MAXCONDUCTANCE']=float(rates[y])/100.0 parameters['BRAINNAME']='vb-exp'+str(x)+'-'+str(y) brainsim8.simBrain(parameters) # For speed: # - move associated Report files to separate folder # - remove associated files y=y+1 y=0 #reset x=x+1 x=0 y=0 while x < len(rates): while y < len(rates): print rates[x],rates[y] lJ=brainsim8.loadBrain('vb-exp'+str(x)+'-'+str(y)+'-JReport.txt') lM=brainsim8.loadBrain('vb-exp'+str(x)+'-'+str(y)+'-MReport.txt') lX=brainsim8.loadBrain('vb-exp'+str(x)+'-'+str(y)+'-XReport.txt') lY=brainsim8.loadBrain('vb-exp'+str(x)+'-'+str(y)+'-YReport.txt') lAP=brainsim8.loadBrain('vb-exp'+str(x)+'-'+str(y)+'-APReport.txt') lOP=brainsim8.loadBrain('vb-exp'+str(x)+'-'+str(y)+'-OPReport.txt') lAQ=brainsim8.loadBrain('vb-exp'+str(x)+'-'+str(y)+'-AQReport.txt') lOQ=brainsim8.loadBrain('vb-exp'+str(x)+'-'+str(y)+'-OQReport.txt') lI=brainsim8.loadBrain('vb-exp'+str(x)+'-'+str(y)+'-IReport.txt') _,spikesJ=brainalyze.countSpikes(lJ,10000*(parameters['ENDSIM']/2)) _,spikesM=brainalyze.countSpikes(lM,10000*(parameters['ENDSIM']/2)) _,spikesX=brainalyze.countSpikes(lX,10000*(parameters['ENDSIM']/2)) _,spikesY=brainalyze.countSpikes(lY,10000*(parameters['ENDSIM']/2)) _,spikesAP=brainalyze.countSpikes(lAP,10000*(parameters['ENDSIM']/2)) _,spikesOP=brainalyze.countSpikes(lOP,10000*(parameters['ENDSIM']/2)) _,spikesAQ=brainalyze.countSpikes(lAQ,10000*(parameters['ENDSIM']/2))
plt.subplot(914) plt.plot(sp,lY,'b') plt.subplot(915) plt.plot(sp,lAP,'b') plt.subplot(916) plt.plot(sp,lOP,'b') plt.subplot(917) plt.plot(sp,lAQ,'b') plt.subplot(918) plt.plot(sp,lOQ,'b') plt.subplot(919) plt.plot(sp,lI,'b') plt.show() brainsim8.simBrain(parameters) lJ=brainsim8.loadBrain('vb-exp-JReport.txt') lM=brainsim8.loadBrain('vb-exp-MReport.txt') lX=brainsim8.loadBrain('vb-exp-XReport.txt') lY=brainsim8.loadBrain('vb-exp-YReport.txt') lAP=brainsim8.loadBrain('vb-exp-APReport.txt') lOP=brainsim8.loadBrain('vb-exp-OPReport.txt') lAQ=brainsim8.loadBrain('vb-exp-AQReport.txt') lOQ=brainsim8.loadBrain('vb-exp-OQReport.txt') lI=brainsim8.loadBrain('vb-exp-IReport.txt') _,spikesJ=brainalyze.countSpikes(lJ,10000*(parameters['ENDSIM']/2)) _,spikesM=brainalyze.countSpikes(lM,10000*(parameters['ENDSIM']/2)) _,spikesX=brainalyze.countSpikes(lX,10000*(parameters['ENDSIM']/2)) _,spikesY=brainalyze.countSpikes(lY,10000*(parameters['ENDSIM']/2)) _,spikesAP=brainalyze.countSpikes(lAP,10000*(parameters['ENDSIM']/2)) _,spikesOP=brainalyze.countSpikes(lOP,10000*(parameters['ENDSIM']/2)) _,spikesAQ=brainalyze.countSpikes(lAQ,10000*(parameters['ENDSIM']/2))