如何在Python中使用np.einsum爱因斯坦求和约定优化计算

大萱姑娘_5653

大萱姑娘_5653

2026-09-15

421人浏览

原创

该用np.einsum而非np.dot或广播时,是需执行多维“带下标规则的乘加”操作,如批量矩阵乘、迹计算、跨轴归约或避免中间数组分配;它控制计算图结构,非单纯语法糖。

如何在python中使用np.einsum爱因斯坦求和约定优化计算

什么时候该用 np.einsum 而不是常规的 np.dot 或广播?

当你要做多维数组之间的“带下标规则的乘加”时,np.einsum 才真正不可替代。比如批量矩阵乘(batch@matrix)、对角线提取+求和、跨轴归约合并,或者需要避免中间大数组分配的场景。它不是语法糖,而是控制计算图结构的工具——如果你只是两个二维矩阵相乘,np.dot@ 更快也更可读。

  • 典型适用场景:"ij,ijk->ik"(每个二维切片与对应三维张量相乘)、"ii->"(迹)、"ij,jk->ik"(等价于 dot,但显式可控)
  • 常见误用:拿 np.einsum("ij,jk->ik", a, b) 替代 a @ b——性能反而略差,且可读性下降
  • 关键优势在于:能跳过冗余维度展开,比如 "i,i->"np.sum(a * b) 少一次临时数组分配

np.einsum 的字符串格式怎么写才不出错?

核心是“输入下标 + 箭头 + 输出下标”,每个字母代表一个轴,重复字母表示该轴要被求和(收缩),未出现在输出中的字母即为求和轴。大小写敏感,且最多支持 52 个不同字母(a–z, A–Z),但别真用满——可读性会崩。

  • 错误示例:"ij,jk->i" —— 左边 j 出现在两个输入中,右边没出现,按规则应求和;但 i 在左边只出现一次,右边却只写了 i,意味着你丢掉了 k 维,实际会报 ValueError: output shape not matched
  • 正确对照:"ij,jk->ik"(标准矩阵乘)、"ij,ij->i"(每行点积)、"...i,...i->..."(支持广播的批量内积)
  • 省略号 ... 表示“前面任意数量的批处理维度”,但必须两边一致;"...ij,...jk->...ik" 是安全的批量矩阵乘写法

为什么开了 optimize=True 反而变慢或出错?

optimize 参数控制是否启用路径优化,默认 False。设为 True 时,np.einsum 会调用 opt_einsum 的逻辑预计算最优收缩顺序(类似动态规划),对三路及以上张量运算收益明显;但对两路运算,开销常大于收益,甚至因浮点误差累积导致结果微小偏差。

Shadows Python Sensei
Shadows Python Sensei

Python 最佳实践助手——代码规范、设计模式、性能优化、测试与类型注解。适用于编写或审查 Python 代码。

下载
  • 建议策略:两输入时保持 optimize=False;三输入及以上(如 "ij,jk,kl->il")务必设 optimize=True,并用 np.einsum_path 预览路径
  • 陷阱:optimize="greedy""optimal" 在大数组上可能卡住——内存爆掉或耗时分钟级,生产环境慎用 "optimal"
  • 验证方式:np.allclose(np.einsum(..., optimize=False), np.einsum(..., optimize=True)) 必须为 True,否则说明数值不稳定,应回退

如何确认 np.einsum 真的比手写循环/广播更快?

不能只看单次执行时间。np.einsum 的优势在避免中间数组、融合多个操作;但如果输入很小(

  • 实测建议:用 timeit 对比,输入尺寸至少到 (100, 100) 以上;关注内存占用(用 memory_profiler),有时速度差不多,但 einsum 内存少 3 倍
  • 典型胜出案例:"ij,j->i"(向量加权和)比 np.sum(a * b[None, :], axis=1) 稳定快 1.2–1.5×,且不生成 (i,j) 临时数组
  • 容易翻车的点:含大量 ... 的高维表达式,解释器解析成本陡增;此时拆成 reshape + 简单 einsum 更稳

最常被忽略的是:np.einsum 不自动处理类型提升,float32 输入默认仍按 float32 计算,但某些硬件上 float64 路径反而更快——得自己用 dtype 参数指定。

Python免费学习笔记(深入):立即使用
在学习笔记中,你将探索 Python 的核心概念和高级技巧!

相关文章

PHP速学视频免费教程(入门到精通)
PHP速学视频免费教程(入门到精通)

PHP怎么学习?PHP怎么入门?PHP在哪学?PHP怎么学才快?不用担心,这里为大家提供了PHP速学教程(入门到精通),有需要的小伙伴保存下载就能学习啦!

下载

本站声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn

相关专题

更多
python打包成可执行文件
python打包成可执行文件

本专题为大家带来python打包成可执行文件相关的文章,大家可以免费的下载体验。

2023.07.20

1571

4

python能做什么
python能做什么

python能做的有:可用于开发基于控制台的应用程序、多媒体部分开发、用于开发基于Web的应用程序、使用python处理数据、系统编程等等。本专题为大家提供python相关的各种文章、以及下载和课程。

2023.07.25

3704

7

format在python中的用法
format在python中的用法

Python中的format是一种字符串格式化方法,用于将变量或值插入到字符串中的占位符位置。通过format方法,我们可以动态地构建字符串,使其包含不同值。php中文网给大家带来了相关的教程以及文章,欢迎大家前来阅读学习。

2023.07.31

1569

3

python教程
python教程

Python已成为一门网红语言,即使是在非编程开发者当中,也掀起了一股学习的热潮。本专题为大家带来python教程的相关文章,大家可以免费体验学习。

2023.08.03

21237

23

python环境变量的配置
python环境变量的配置

Python是一种流行的编程语言,被广泛用于软件开发、数据分析和科学计算等领域。在安装Python之后,我们需要配置环境变量,以便在任何位置都能够访问Python的可执行文件。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

2627

5

python eval
python eval

eval函数是Python中一个非常强大的函数,它可以将字符串作为Python代码进行执行,实现动态编程的效果。然而,由于其潜在的安全风险和性能问题,需要谨慎使用。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

2687

5

scratch和python区别
scratch和python区别

scratch和python的区别:1、scratch是一种专为初学者设计的图形化编程语言,python是一种文本编程语言;2、scratch使用的是基于积木的编程语法,python采用更加传统的文本编程语法等等。本专题为大家提供scratch和python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

1083

5

python合并两个列表
python合并两个列表

Python是一种强大的编程语言,具有许多方便的功能和工具。在Python中,有多种方法可以合并两个列表。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.10

576

4

python是前端还是后端
python是前端还是后端

Python属于前端也属于后端,其灵活性和丰富的生态系统使得开发人员能够在不同的领域中灵活运用。本专题为大家提供python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

2063

5

热门下载

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

精品课程

更多
相关推荐
/
热门推荐
/
最新课程