Eigen 取三角矩阵（对应matlab tril triu函数）

#include <Eigen/Dense>
#include <iostream>
using namespace Eigen;
using namespace std;
int main(){
	MatrixXf zzz(4, 4);
	Eigen::MatrixXf l;
	zzz << 1, 2, 3, 4, 5, 6, 7, 8, 9, 7, 8, 9, 8, 8, 4, 5;
	cout << "MatrixXf" << endl << zzz << endl;
	l = zzz.triangularView<Eigen::Lower>();
	cout << "Lower" << endl << l << endl;
	l = zzz.triangularView<Eigen::Upper>();
	cout << "Upper" << endl << l << endl;
	l = zzz.triangularView<Eigen::ZeroDiag>();  //后来发现在某些情况下不对，找不到原因，只能改用StrictlyUpper
	cout << "ZeroDiag" << endl << l << endl;
	l = zzz.triangularView<Eigen::UnitLower>();
	cout << "UnitLower" << endl << l << endl;
	l = zzz.triangularView<Eigen::UnitUpper>();
	cout << "UnitUpper" << endl << l << endl;
	l = zzz.triangularView<Eigen::StrictlyLower>();
	cout << "StrictlyLower" << endl << l << endl;
	l = zzz.triangularView<Eigen::StrictlyUpper>();
	cout << "StrictlyUpper" << endl << l << endl;
	l = zzz.triangularView<Eigen::SelfAdjoint>();
	cout << "SelfAdjoint" << endl << l << endl;
	l = zzz.triangularView<Eigen::Symmetric>();
	cout << "Symmetric" << endl << l << endl;
}

结果：

链接：（官网）

triangularView

各个功能详细介绍