代码之家 › 专栏 › 技术社区 › Tom de Geus

部分模板专门化以展开特定大小的循环

c++

Tom de Geus · 技术社区 · 7 年前

我写了一个矩阵类,它可以有不同的大小。现在我要展开特定大小的循环。我该怎么做?

唯一能让我工作的方法是一个2-D的儿童类,但是我想避免这种情况,因为它会导致很多重复的代码。

例如:

#include <iostream>

template<class T, size_t M, size_t N>
class matrix
{
  matrix<T,M,N>& operator*= (const matrix<T,M,N> &B);
};

template<class T, size_t M, size_t N>
matrix<T,M,N>& matrix<T,M,N>::operator*= (const matrix<T,M,N> &B)
{
  // ...
  return *this;
}

int main()
{
  return 0;
}

现在我想为这个案例添加一个实现 M = 2 和 N = 2 我展开所有循环以获得效率。

(我在前面的一个示例中对展开进行了计时,这似乎确实有意义,特别是对于更复杂的操作,然后在本例中进行了介绍。)

1 回复 | 直到 7 年前

Maxim Egorushkin 7 年前

你可以委派 operator*= 到重载的函数模板。例如。:

template<class T, size_t M, size_t N>
class matrix
{
public:
    matrix<T,M,N>& operator*=(const matrix<T,M,N>&);
};

// Generic version.
template<class T, size_t M, size_t N>
void inplace_dot(matrix<T,M,N>& a, matrix<T,M,N> const& b); 

// Overload for 2x2 matrix.
template<class T>
void inplace_dot(matrix<T,2,2>& a, matrix<T,2,2> const& b);

template<class T, size_t M, size_t N>
matrix<T,M,N>& matrix<T,M,N>::operator*=(const matrix<T,M,N>& b)
{
    inplace_dot(*this, b);
    return *this;
}

推荐文章

AstralHex · 矩阵乘法代码工作不正常

6 月前

Baba Dan Constantin · SSE4.1在矩阵4x4乘法上比SSE3慢?

6 月前

Giogre · 为包含许多数值字段的简单“struct”重载比较运算符

6 月前

einpoklum · 定义一个并不真正提供now()函数的std::chrono Clock是“合法的”吗?

6 月前

Fishie · 作为类成员的智能指针是否仍然自动释放?[关闭]

6 月前

Die4Toast · 递归调用成员箭头运算符->

6 月前

Angle.Bracket · 如何用C++将UTF-8文件名写入MS Windows控制台?

6 月前

Anka HanÄ±m · 关于结构和动态数组地址的问题

6 月前

Adam Barnes · 我如何定义一个基于constexpr函数返回值进行限制的概念?

6 月前

user2138149 · 为什么我不能获取包含多个元素的结构体中某些元素的地址?[副本]

6 月前