function [accN,accP] = grt_acc(M);

[r,c,m] = size(M);

corr = 0;
for mi = 1:m
    for i = 1:c
        for j = i
            corr = corr + M(i,j,mi);
        end
    end
end
accN = corr;
accP = corr/sum(sum(sum(M)));