0

0

C++怎么使用Eigen库_C++矩阵运算教程【科学】

穿越時空

穿越時空

发布时间:2026-02-23 15:51:38

|

704人浏览过

|

来源于php中文网

原创

eigen是纯头文件模板库,无需链接;只需将include目录加入编译路径,如target_include_directories(... public ${eigen3_include_dir})。

c++怎么使用eigen库_c++矩阵运算教程【科学】

怎么确认 Eigen 是头文件库,不用链接

Eigen 不是传统意义上的动态或静态库,它压根没提供 .so.lib 文件。你看到的 #include <eigen></eigen> 能直接编译通过,是因为 Eigen 全是模板头文件——所有实现都在 .h 里。编译器在实例化模板时才生成机器码。

所以:不需要 target_link_libraries,也不用配 LD_LIBRARY_PATH。常见错误是 CMake 里傻乎乎加 find_package(Eigen3 REQUIRED) 后又去 link,结果报 undefined reference ——那是因为你误以为它有符号要链接。

  • 只要把 Eigen 的 include 目录加进编译器头文件搜索路径(CMake 用 target_include_directories(... PUBLIC ${EIGEN3_INCLUDE_DIR})
  • 确保 #include 路径正确,比如 #include <eigen></eigen>,不是 #include <eigen3></eigen3>(除非你手动改了路径)
  • Mac 上用 Homebrew 装的 Eigen,默认头文件在 /opt/homebrew/include/eigen3,别漏掉末尾的 eigen3

MatrixXd 和 MatrixXf 到底该选哪个

类型选错最直接影响数值精度和运行时行为。Eigen 中矩阵类型由三部分决定:基类型(double/float)、行列数(X 表示动态)、存储顺序(默认列优先)。别光看名字缩写,得盯住底层类型。

MatrixXdMatrix<double dynamic></double>MatrixXfMatrix<float dynamic></float>。两者内存布局一样,但浮点精度差一个数量级——float 在求逆、SVD、共轭梯度等病态问题中容易溢出或丢失有效位。

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

MyMap AI
MyMap AI

使用AI将想法转化为图表

下载
  • 科研计算、需要高精度中间结果(比如优化、PDE 离散),一律用 MatrixXd
  • 嵌入式或 GPU 前端预处理、对吞吐敏感且输入本身是 float(如图像像素),可用 MatrixXf,但要自己验证误差是否可接受
  • 混用会触发隐式转换:把 MatrixXf 赋给 MatrixXd 变量,编译不报错,但运行时多一次逐元素 cast,性能白丢

为什么 resize() 后数据没了,而 conservativeResize() 有时也不行

Eigen 默认不保留旧数据,这是为了明确语义:resize() 就是“我要新尺寸”,不是“我打算扩/缩但留着原来那些”。很多人以为它像 Python 的 numpy.resize,结果发现原矩阵内容全变零。

conservativeResize() 确实会尽量保留左上角子块,但它只在“新尺寸不小于旧尺寸”且“内存连续可延展”时才安全。比如从 3x3 扩到 4x4 成功;但从 4x4 缩回 3x3,虽然能调用,但若原矩阵是表达式求值结果(如 A + B),内部没有实际分配内存,就会 UB。

  • 想保留数据,先用 resize(),再用 .block(0,0,old_rows,old_cols) = old_matrix.block(0,0,old_rows,old_cols) 显式拷贝
  • 更稳妥的做法是:一开始就用 Matrix<double dynamic>::Zero(rows, cols)</double> 预分配,避免反复 resize
  • 如果真要动态调整且频繁,考虑换 std::vector<:vector>></:vector> 或用 Eigen::Map 包装外部缓冲区

用 Map 读原始数组时对齐失败就崩溃

Eigen::Map 默认要求指针对齐(16 字节对齐用于 SSE,32 字节用于 AVX)。如果你传进去的是普通 new double[n] 或栈上数组,大概率不对齐,运行时触发 segmentation fault 或静默错误——尤其在开启 -mavx 编译时。

根本原因:Eigen 在向量化路径里直接发 movaps 指令,该指令硬性要求地址对齐。不是所有平台都报错,但行为不可靠。

  • 临时方案:构造 Map 时加 Unaligned 标签,例如 Eigen::Map<:vectorxd eigen::unaligned>(ptr, n)</:vectorxd>
  • 长期方案:用 Eigen::aligned_allocator 配合 std::vector,或直接用 Eigen::VectorXd::Map(ptr, n)(它内部自动处理对齐策略)
  • 注意:Unaligned 会降级为标量指令或非对齐向量指令,性能损失约 10%~30%,大数据量时得权衡

对齐这事看不见摸不着,但一旦出问题,core dump 里连函数名都可能被优化掉,查起来比维度错还费劲。

热门AI工具

更多
DeepSeek
DeepSeek

幻方量化公司旗下的开源大模型平台

豆包大模型
豆包大模型

字节跳动自主研发的一系列大型语言模型

通义千问
通义千问

阿里巴巴推出的全能AI助手

腾讯元宝
腾讯元宝

腾讯混元平台推出的AI助手

文心一言
文心一言

文心一言是百度开发的AI聊天机器人,通过对话可以生成各种形式的内容。

讯飞写作
讯飞写作

基于讯飞星火大模型的AI写作工具,可以快速生成新闻稿件、品宣文案、工作总结、心得体会等各种文文稿

即梦AI
即梦AI

一站式AI创作平台,免费AI图片和视频生成。

ChatGPT
ChatGPT

最最强大的AI聊天机器人程序,ChatGPT不单是聊天机器人,还能进行撰写邮件、视频脚本、文案、翻译、代码等任务。

相关专题

更多
css中float用法
css中float用法

css中float属性允许元素脱离文档流并沿其父元素边缘排列,用于创建并排列、对齐文本图像、浮动菜单边栏和重叠元素。想了解更多float的相关内容,可以阅读本专题下面的文章。

592

2024.04.28

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

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

105

2025.10.23

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

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

294

2025.08.29

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

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

105

2025.10.23

堆和栈的区别
堆和栈的区别

堆和栈的区别:1、内存分配方式不同;2、大小不同;3、数据访问方式不同;4、数据的生命周期。本专题为大家提供堆和栈的区别的相关的文章、下载、课程内容,供大家免费下载体验。

422

2023.07.18

堆和栈区别
堆和栈区别

堆(Heap)和栈(Stack)是计算机中两种常见的内存分配机制。它们在内存管理的方式、分配方式以及使用场景上有很大的区别。本文将详细介绍堆和栈的特点、区别以及各自的使用场景。php中文网给大家带来了相关的教程以及文章欢迎大家前来学习阅读。

595

2023.08.10

golang map内存释放
golang map内存释放

本专题整合了golang map内存相关教程,阅读专题下面的文章了解更多相关内容。

77

2025.09.05

golang map相关教程
golang map相关教程

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

36

2025.11.16

pixiv网页版官网登录与阅读指南_pixiv官网直达入口与在线访问方法
pixiv网页版官网登录与阅读指南_pixiv官网直达入口与在线访问方法

本专题系统整理pixiv网页版官网入口及登录访问方式,涵盖官网登录页面直达路径、在线阅读入口及快速进入方法说明,帮助用户高效找到pixiv官方网站,实现便捷、安全的网页端浏览与账号登录体验。

1127

2026.02.13

热门下载

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

精品课程

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

共94课时 | 10万人学习

C 教程
C 教程

共75课时 | 4.9万人学习

C++教程
C++教程

共115课时 | 19.1万人学习

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

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