0

0

C++怎么实现一个表达式模板(Expression Templates)_C++高性能计算与延迟求值

裘德小鎮的故事

裘德小鎮的故事

发布时间:2025-12-01 05:01:15

|

727人浏览过

|

来源于php中文网

原创

表达式模板通过CRTP构建编译期表达式树,延迟向量加法求值,避免临时对象与多次循环,实现零成本抽象与循环融合。

c++怎么实现一个表达式模板(expression templates)_c++高性能计算与延迟求值

表达式模板(Expression Templates)是一种在编译期优化数值计算的技术,常用于高性能计算库(如Eigen、Blaze)中实现延迟求值和消除临时对象。它的核心思想是:将数学表达式构建成一个模板表达式树,在运行时才真正求值,从而避免中间结果的创建。

基本动机:避免临时对象

考虑两个向量相加:

Vector a, b, c, d;
d = a + b + c;

如果没有表达式模板,每次 + 都会生成一个临时 Vector 对象,导致多次内存分配和拷贝。而表达式模板可以让编译器“看到”整个表达式结构,在最终赋值时一次性计算,跳过中间临时变量。

实现一个简单的向量加法表达式模板

我们从一个简化的例子开始,实现支持延迟加法的向量类型。

立即学习C++免费学习笔记(深入)”;

1. 定义基础向量类

template
struct Vector {
    std::vector data;

    Vector(int n) : data(n) {}

    int size() const { return data.size(); }

    T& operator[](int i) { return data[i]; }
    const T& operator[](int i) const { return data[i]; }
};

2. 定义表达式基类(CRTP)

template
struct Expr {
    const E& self() const { return static_cast(*this); }
};

这里使用了CRTP(奇异递归模板模式),让基类能访问派生类的类型信息,而无需虚函数开销。

3. 定义加法表达式类型

template
struct AddExpr : Expr> {
    const LHS& lhs;
    const RHS& rhs;

    AddExpr(const LHS& l, const RHS& r) : lhs(l), rhs(r) {}

    // 延迟求值:只在需要时计算
    auto operator[](int i) const { return lhs[i] + rhs[i]; }
    int size() const { return lhs.size(); }
};

4. 重载操作符 +

AI发型设计
AI发型设计

虚拟发型试穿工具和发型模拟器

下载
template
auto operator+(const Expr& lhs, const Expr& rhs) {
    return AddExpr(lhs.self(), rhs.self());
}

注意:所有表达式类型都继承自 Expr,所以这个 + 能接受任意组合的表达式。

5. 让 Vector 成为表达式

template
struct Vector : Expr> {
    ... // 同上
};

现在 Vector 也成了表达式类型,可以参与表达式构建。

6. 实现赋值操作以触发求值

template
Vector& operator=(Vector& vec, const Expr& expr) {
    auto& e = expr.self();
    for (int i = 0; i         vec[i] = e[i];
    }
    return vec;
}

现代C++中可写成泛型 lambda 形式,或使用具体模板参数替代 auto

实际使用与效果

Vector a(1000), b(1000), c(1000), d(1000);
// ... 初始化
d = a + b + c;

这行代码的执行过程:

  • a + b 返回一个 AddExpr
  • 再与 c 相加,返回 AddExpr, Vector>
  • 赋值时遍历一次,计算 (a[i] + b[i]) + c[i]

全程没有构造临时 Vector 对象,循环合并为一次,实现“循环融合”(loop fusion)。

扩展与注意事项

真实场景中还需考虑:

  • 支持更多操作(减法、标量乘法、点积等)
  • 处理别名问题(Aliasing),必要时强制求值
  • 表达式深度过大可能导致编译时间变长或溢出
  • 调试困难:表达式类型名非常复杂

基本上就这些。表达式模板本质是把运行时的计算调度提前到编译期建模,用模板实例化代替临时对象,是C++零成本抽象的典范之一。不复杂但容易忽略细节。

相关专题

更多
c语言const用法
c语言const用法

const是关键字,可以用于声明常量、函数参数中的const修饰符、const修饰函数返回值、const修饰指针。详细介绍:1、声明常量,const关键字可用于声明常量,常量的值在程序运行期间不可修改,常量可以是基本数据类型,如整数、浮点数、字符等,也可是自定义的数据类型;2、函数参数中的const修饰符,const关键字可用于函数的参数中,表示该参数在函数内部不可修改等等。

524

2023.09.20

string转int
string转int

在编程中,我们经常会遇到需要将字符串(str)转换为整数(int)的情况。这可能是因为我们需要对字符串进行数值计算,或者需要将用户输入的字符串转换为整数进行处理。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

318

2023.08.02

int占多少字节
int占多少字节

int占4个字节,意味着一个int变量可以存储范围在-2,147,483,648到2,147,483,647之间的整数值,在某些情况下也可能是2个字节或8个字节,int是一种常用的数据类型,用于表示整数,需要根据具体情况选择合适的数据类型,以确保程序的正确性和性能。本专题为大家提供相关的文章、下载、课程内容,供大家免费下载体验。

538

2024.08.29

c++怎么把double转成int
c++怎么把double转成int

本专题整合了 c++ double相关教程,阅读专题下面的文章了解更多详细内容。

52

2025.08.29

C++中int的含义
C++中int的含义

本专题整合了C++中int相关内容,阅读专题下面的文章了解更多详细内容。

197

2025.08.29

c++怎么把double转成int
c++怎么把double转成int

本专题整合了 c++ double相关教程,阅读专题下面的文章了解更多详细内容。

52

2025.08.29

C++中int、float和double的区别
C++中int、float和double的区别

本专题整合了c++中int和double的区别,阅读专题下面的文章了解更多详细内容。

99

2025.10.23

lambda表达式
lambda表达式

Lambda表达式是一种匿名函数的简洁表示方式,它可以在需要函数作为参数的地方使用,并提供了一种更简洁、更灵活的编码方式,其语法为“lambda 参数列表: 表达式”,参数列表是函数的参数,可以包含一个或多个参数,用逗号分隔,表达式是函数的执行体,用于定义函数的具体操作。本专题为大家提供lambda表达式相关的文章、下载、课程内容,供大家免费下载体验。

204

2023.09.15

高德地图升级方法汇总
高德地图升级方法汇总

本专题整合了高德地图升级相关教程,阅读专题下面的文章了解更多详细内容。

43

2026.01.16

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
国外Web开发全栈课程全集
国外Web开发全栈课程全集

共12课时 | 1.0万人学习

进程与SOCKET
进程与SOCKET

共6课时 | 0.3万人学习

关于我们 免责申明 举报中心 意见反馈 讲师合作 广告合作 最新更新
php中文网:公益在线php培训,帮助PHP学习者快速成长!
关注服务号 技术交流群
PHP中文网订阅号
每天精选资源文章推送

Copyright 2014-2026 https://www.php.cn/ All Rights Reserved | php.cn | 湘ICP备2023035733号