降低基于本本的模板

我正在尝试基于特征矩阵进行减少。

#include <iostream>
#include <Eigen/Dense>
#include <type_traits>
template<typename T1, typename T2, int n1, int n2>
auto reduction(Eigen::Matrix<T1, n1, n2> &a1,
               Eigen::Matrix<T2, n1, n2> &a2)
     -> decltype(T1{}*T2{})
{
  using BaseT3 = 
    typename std::remove_cv<typename std::remove_reference<decltype(T1{}*T2{})>::type>::type;
  BaseT3 res = a1(0, 0)*a2(0, 0);
  for (int i=0; i<n1; ++i)
    for (int j=0; j<n2; ++j)
      if (i+j)
        res = res + a1(i, j)*a2(i, j);
  return res;
}
int main()
{
  Eigen::Matrix<double, 3, 3> m;
  Eigen::Matrix<Eigen::Vector3d, 3, 3> n;
  std::cout << reduction(m, n) << std::endl;
}

基本上，我试图获得sum_{i, j} a1[i, j] * a2[i, j]，其中a1和a2是一些特征Mathix，但我会遇到编译错误。我得到的错误是

error: no match for ‘operator=’ (operand types are ‘BaseT3 {aka 
Eigen::CwiseUnaryOp<Eigen::internal::scalar_multiple_op<double>, 
const Eigen::Matrix<double, 3, 1> >}’ 
and 
‘const Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double>, 
const Eigen::CwiseUnaryOp<Eigen::internal::scalar_multiple_op<double>, 
const Eigen::Matrix<double, 3, 1> >, 
const Eigen::CwiseUnaryOp<Eigen::internal::scalar_multiple_op<double>, 
const Eigen::Matrix<double, 3, 1> > >’)
         res = res + a1(i, j)*a2(i, j);
             ^

如果我没记错的话，对于给定的main，类型BaseT3应该是Eigen::Vector3d。我也尝试静态铸造，以便operator=不应该失败，但我会遇到其他错误。

这是C 11，我使用eigen3，编译器为G 5.4.1。