% combinemattomat_fixeddt.m % A MATLAB script for summing the the files saved by running "calcsumcurr_manyareagsynmediumtau_parts_fixeddt.py 20 0.025 0.000042 10000 10000 2 myseed 200" % Expects that the variable "iseed" has been initialized. iseed=1 corresponds to myseed=1,...,10 in combinemattomat_fixeddt.m, iseed=2 corresponds to % myseed=11,...,20, etc. % Tuomo Maki-Marttunen, 2014-2016 synloctype = 2; dt = 0.1; T=10000; Nts = round(T/dt)-1; ts = dt*(0:Nts-1); Nvox=13; ica = zeros(Nvox,Nts); icap = zeros(Nvox,Nts); il = zeros(Nvox,Nts); ik = zeros(Nvox,Nts); ih = zeros(Nvox,Nts); ina = zeros(Nvox,Nts); imemb = zeros(Nvox,Nts); VtimesA = zeros(Nvox,Nts); Vsoma = cell(1,10); ts_syn = cell(1,10); part_syn = cell(1,10); for irep=1:10 disp(['Loading ' num2str(iseed)]); A=load(['currsums_parts_10000areagsynsmediumtau_fixeddt_type' num2str(synloctype) '_amp4.2e-05_tstop10000.0_nseg20_dt0.025_seed' num2str((iseed-1)*10+irep-1) '_comb200.0.mat']); ica = ica + interpolate_multidim(A.times,A.ica,ts); icap = icap + interpolate_multidim(A.times,A.icap,ts); il = il + interpolate_multidim(A.times,A.il,ts); ik = ik + interpolate_multidim(A.times,A.ik,ts); ih = ih + interpolate_multidim(A.times,A.ih,ts); ina = ina + interpolate_multidim(A.times,A.ina,ts); imemb = imemb + interpolate_multidim(A.times,A.imemb,ts); VtimesA = VtimesA + interpolate_multidim(A.times,A.VtimesA,ts); Vsoma{irep} = interpolate(A.times,A.Vsoma,ts); ts_syn{irep} = A.ts_syn; part_syn{irep} = A.part_syn; end save(['currsums_parts_10000areagsynsmediumtau_fixeddt_type' num2str(synloctype) '_amp4.2e-05_tstop10000.0_nseg20_dt0.025_seed' num2str(iseed) '_comb_summed.mat'],'ts','ica','icap','il','ik','ih','ina','imemb','VtimesA','Vsoma','ts_syn','part_syn');