背景

---

最近接触到了矩阵库Eigen，在此记录一下使用中遇到的一些问题。

1. 块操作

---

MatrixXd M = MatrixXd::Zero(50, 50);
cout << M.block(2, 3, 5, 6) << endl;  //指的是M矩阵从(2,3)位置开始，下数5个元素，右数6个元素

M.block(a, b, m ,n)指的是从下标(a,b)的元素开始，向下m个元素，向右n个元素。

2. 矩阵和C++数组转换

---

MatrxXd M = MatrixXd::Zero(50,50);
//矩阵转数组
double* test = M.data();
//数组转矩阵
Map<MatrixXd>tM(test, 50, 50);

注意，这两种转换方式都是浅拷贝，改变源数据的值拷贝后对象的值也会改变。

3. 矩阵输出到文件

---

Eigen重载了<<，直接输出到文件即可。

MatrixXd M = MatrixXd::Zero(5, 5);
ofstream fout("test.txt");
fout << M;
fout.close();

测试结果：

测试结果