Aggregate matrix cells (lower the resolution by a factor)
aggregate_matrix(M, fact)
void
internal function to be documented