Нахождение максимального элемента среди элементов, лежащих выше побочной диагонали матрицы

#include <cstdlib>
#include <iostream>
#include <iomanip>
#include <numeric>
 
template < typename T, size_t size >
T max_matrix_element(T matrix[size][size])
{
    T max = std::numeric_limits<T>::min();
    for (size_t i = 0; i < size; i++)
        for (size_t j = i + 1; j < size; j++)
            if (matrix[i][j] > max)
                max = matrix[i][j];
    return max;
}
 
 
int main()
{
    int matrix[7][7];
    std::srand(std::time(NULL));
    for (size_t i = 0; i < 7; i++)
    {
        for (size_t j = 0; j < 7; j++)
        {
            std::cout << std::setw(4) 
                      << (matrix[i][j] = std::rand() % 20);
        }
        std::cout << "\n";
    }
    std::cout << "Max: " << max_matrix_element<int, 7>(matrix) << std::endl;
#ifdef _MSC_VER || _WIN32
    system("pause");
#endif
    return 0;
}