Динамический массив: вычислить среднее геометрическое чисел, которые входят в заданный интервал

#include <iostream>
#include <iomanip>
#include <cmath>
 
void input_matrix(int **matrix, const size_t matrix_size){
  for(int i = 0; i < matrix_size; ++i)
    for(int j = 0; j < matrix_size; ++j){
      std::cout << "\nMtarix[" << i << "][" << j << "] = ";
      std::cin >> matrix[i][j];
    }
}
 
void output_matrix(int **matrix, const size_t matrix_size){
  for(int i = 0; i < matrix_size; ++i){
    for(int j = 0; j < matrix_size; ++j)
      std::cout << std::setw(3) << matrix[i][j] << ' ';
  }
  std::cout << '\n';
}
 
double get_mid_geo(int **matrix, const size_t matrix_size){
  int multy = 1, cnt = 0;
  for(int i = 0; i < matrix_size; ++i)
    for(int j = 0; j < matrix_size; ++j)
      if( matrix[i][j] >= -2 && matrix[i][j] < 5 ){
        multy *= matrix[i][j];
        ++cnt;
      }
  return pow( multy, static_cast<double>(cnt) );
}
 
int main(){
  size_t matrix_size;
  std::cout << "\n\nInput matrix size -> ";
  std::cin >> matrix_size;
  int **matrix = new int*[matrix_size];
  for(int i = 0; i < matrix_size; ++i)
    matrix[i] = new int[matrix_size];
  input_matrix(matrix, matrix_size);
  std::cout << '\n';
  output_matrix(matrix, matrix_size);
  std::cout << '\n' << "Result -> " << get_mid_geo(matrix, matrix_size);
  for(int i = 0; i < matrix_size; ++i)
    delete[] matrix[i];
  delete[] matrix;
  return 0;
}