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))
