使用模板类作为参数专门化模板结构

Question

我正在尝试通过创建一个包含矩阵和对这些矩阵的操作的库来提高我的模板技能（我知之甚少）。基本上，我希望我的矩阵具有非常强的类型（在编译时已知的数据类型和大小），并且我还希望能够自动推断转置矩阵的类型。

模板<类型名TDataType，size_t行，size_t列> class MyMatrix

矩阵可以嵌套，因此 TDataType 可以是整型，但它本身也是一个 MyMatrix<...>，导致数据类型为转置矩阵不一定与原始矩阵相同，例如： Transpose( MyMatrix, 4, 6 > ) ==>我的矩阵我的矩阵char, 3, 2 >, 6, 4 > （外部矩阵的数据类型已更改）

我对转置类型推导的第一次尝试是：

template< typename TDataType >
struct Transpose
  {
  typedef TDataType type;
  };

template<>
struct Transpose< MyMatrix<TDataType, rows, cols> >
  {
  typedef MyMatrix<typename Transpose<TDataType>::type, cols, rows> type;
  };

我找不到这样做的方法，因为我不能' t 似乎专门使用 MyMatrix 转置模板（TDataType 未知和类似错误）。

我想出的唯一可编译的解决方案（我什至不知道它是否有效）是这样的：

template< typename TMatrixType, typename TDataType, size_t rows, size_t cols >
struct Transpose
  {
  typedef TMatrixType type;
  };

template< typename TDataType, size_t rows, size_t cols >
struct Transpose< MyMatrix<TDataType, rows, cols>, TDataType, rows, cols >
  {
  typedef MyMatrix< typename Transpose<TDataType,TDataType,rows,cols>::type, cols, rows > type;
  };

我相信我让事情变得太复杂了；有没有更简单的解决方案来实现我想要的？

---

我的问题的答案（我在没有帐户的情况下发布了问题，所以我没有足够的代表以正常方式做事）。已经非常感谢了！

@Bo Persson @Will A：我不打算将其用作通用矩阵库，我想对特定（事先已知）大小的矩阵执行操作，并想看看使用这种方法可以得到什么。它可以让我优化矩阵的内存布局（例如，在 32 字节边界上对齐行向量）并执行其他类型的时髦操作。我希望这样做会多次搬起石头砸自己的脚，但我在这里尝试的主要事情是积累经验，找出什么有效，什么无效（以及什么很难做，什么不可以） 't）。

@Bo Perrson：我知道为什么第一个版本无法编译，但我想知道我的第二次尝试是否有一个更简单的版本可以工作。主要问题是 MyMatrix 本身就是一个类模板，我需要以某种方式将它的模板参数获取到 Transpose 结构。

@VJo：我认为这行不通。如果T是MyMatrix<...>那么转置矩阵应该将 Transpose 作为数据类型，而不是 T 本身。对于所有基本类型（char、int、double...），这当然是正确的并且更简单。

Answer 1

是的，你让事情变得复杂了。

如果你有这样的声明：

template< typename TDataType, size_t rows, size_t cols > class MyMatrix

那么转置函数应该是这样的：

template< typename T, size_t rows, size_t cols >
MyMatrix< T, cols, rows > Transpose( const MyMatrix< T, rows, cols > & m )
{
  MyMatrix< T, cols, rows > res;
  // implementation
  return res;
}

Answer 2

将行/列作为模板定义的一部分可以得到什么？

我个人的感觉是，这让事情变得太复杂了。无论如何，在模板中包含数据类型，但包含维度似乎是错误的。

Answer 3

第一次尝试失败，因为特化是一种单独的类型，并且基本模板的模板参数在那里未知。

根据语言，第二个版本是正确的，但就像 Will A 所说的那样 - 你真的希望 row 和 col 的每个组合都创建一个新类型吗？

Answer 4

我会这样写，以允许任意深度的递归（矩阵的矩阵的矩阵......）

template<typename T, unsigned rows, unsigned cols>
struct MyMatrix
{
  typedef T value_type;
  T stuff[rows][cols];  // or whatever                                      
};

// For basic types, transpose is identity.                                  
template<typename T>
struct Transpose {
  typedef T result_type;
  result_type operator()(const T & in) {
    return in;
  }
};

// For more complex types, specialize and invoke recursively.
template<typename T, unsigned rows, unsigned cols>
struct Transpose<MyMatrix<T, rows, cols> > {
  typedef MyMatrix<Transpose<T>, cols, rows> result_type;
  result_type operator()(const MyMatrix<T, rows, cols> & in) {
    Transpose<T> transposer;
    // (invoke transposer on each element of in and build result)           
  }
};

这里，Transpose 是一个函子；您创建它的一个实例，但将其作为函数调用。为了获得额外的积分，您可以让它继承 unary_function 并免费获得 result_type typedef...