просматривая поочередно элементы диагоналей двумерного массива

#include <algorithm>
#include <iterator>
#include <iostream>
#include <iomanip>
#include <cstdlib>
#include <ctime>
 
int main()
{
    const int limit = 100;
    const int rows = 5, columns = 5;
    const int size = rows > columns ? columns : rows;
 
    int matrix [rows][columns];
    int diagonal [size];
    int npos = 0;
 
    // random fill the matrix
    srand((unsigned int)time(NULL));
    for (int i = 0; i < rows; ++i) {
        for (int j = 0; j < columns; ++j)
            matrix[i][j] = std::rand() % limit - std::rand() % limit;
    }
 
    // copy the matrix's main diagonal to the single array
    for (int i = 0; i < size; ++i) {
        if(matrix[i][i] > 0)
            diagonal[npos++] = matrix[i][i];
    }
 
    // sort 
    for (int i = 1; i < npos; ++i) {
        for (int k = i, j = i - 1; j >= 0 && diagonal[j] < diagonal[k]; --j) {
            std::swap(diagonal[k--], diagonal[j]);
        }
    }
 
    // print out the matrix
    for (int i = 0; i < rows; ++i) {
        for (int j = 0; j < columns; ++j) {
            std::cout << std::setw(3) << matrix[i][j] << " ";
        }
 
        std::cout << std::endl;
    }
 
    // print out the sorted diagonal
    std::cout << std::endl;
    std::copy(diagonal, diagonal + npos, std::ostream_iterator<int>(std::cout, " "));
    std::cout << std::endl;
 
    system("PAUSE");
    return 0;
}