Написать функцию для вычисления произведения значений матрицы вещественных или целых чисел

#include <iostream>
#include <algorithm>
#include <functional>
#include <random>
#include <iterator>
#include <type_traits>
 
template <typename T>
void f()
{
    if(!std::is_same<int, T>::value && !std::is_same<double, T>::value)
    {
        std::cerr << "Wrong template type" << std::endl;
        return ;
    }
 
    const std::size_t N = 2;
    const std::size_t M = 4;
    const std::size_t LEN = M * N;
    T arr[M][N];
 
    std::mt19937 gen { std::random_device()() };
    std::uniform_real_distribution<> dist(static_cast<T>(1),static_cast<T>(9));
    //or if(std::is_same<int, T>::value) int_distribution else real_distribution
 
    std::generate(&arr[0][0], &arr[0][0] + LEN, [&]{ return dist(gen); });
 
    for(std::size_t i = 0; i < M; ++i)
    {
        for(std::size_t j = 0; j < N; ++j)
            std::cout << arr[i][j] << ' ';
        std::cout << std::endl;
    }
 
    std::cout << "Mul: "
              << std::accumulate(&arr[0][0],&arr[0][0] + LEN,static_cast<T>(1),std::multiplies<T>())
              << std::endl;
}
 
 
int main(int argc, char *argv[])
{
    f<int>(); //or f<double>();
}