Network = gpuArray.randi(3, 1, M);

NetworkReplicated = repmat(Network, [N 1]);

[signalPowerSorted, I] = sort(signalPowerDecibels, 2, 'descend');

NetworkSorted = NetworkReplicated(sub2ind([N M], RowIndex(:), I(:)));

[~, J] = sort(reshape(NetworkSorted, [N M]), 2);

signalPowerGrouped = signalPowerSorted(J);

NetworkGrouped = sort(Network);

groupStartColumn = find(diff([0, NetworkGrouped]));

signalPowerCum = cumsum([zeros(N,1), signalPowerGrouped], 2);

signalMap = (signalPowerCum(:,groupStartColumn+3) - ...

signalPowerCum(:,groupStartColumn)) / 3;

drawMultipleSignalMaps(map, masts, AntennaDirection, H, signalMap)