Еще раз матрицы

#include <iostream>
#include "Eigen/Core"
#include "Eigen/Array" // Matrix.setRandom()

USING_PART_OF_NAMESPACE_EIGEN
using namespace std;

static const int MATR_SIZE = 5;

int main(int argc, const char *argv[])
{
    // Создаем матрицу и заполняем случайными числами
  MatrixXf m(MATR_SIZE, MATR_SIZE);
  m.setRandom(MATR_SIZE, MATR_SIZE);
  cout << "input matrix: " << endl << m << endl;

  double product = 1; // произведение элементов на главной диагонали
  int i;
  for (i = 0; i < MATR_SIZE; i++) {
      product *= abs(m(i, i));
  }
  product = pow( product, 1.0/i ); // среднее геометрическое

  cout << "Geometric mean on main diag: " << product << endl;

  // заменяем эл-ты с четной суммой коэфф-в на среднее геометрическое
  for (int i = 0; i < MATR_SIZE; i++) {
      for (int j = 0; j < MATR_SIZE; j++) {
          if ( !((i + j) & 1) ) {
              m(i,j) = product;
          }
      }
  }

  cout << "output matrix: " << endl << m << endl;
  return 0;
}