Faster way of counting spikes

[dimokaramanlis]

No description provided.

[marius10p]

Thanks, give me a few days until I can test it.

[alexmorley]

LGTM. The following runs fine:

Nk = 10;
st3 = randsample(1:Nk, 1000, true);

nspk_1 = accumarray(st3(:), 1, [Nk, 1], @sum);


nspk_2 = zeros(Nk, 1);
for j = 1:Nk
    nspk_2(j) = sum(st3==j); % determine total number of spikes in each neuron
end

assert(all(nspk_1 == nspk_2))