【数学问题2】QR分解

问题描述

eigen

int main() {
   
	
	Matrix3f A;
	A << 1, 1, 1,
		2, -1, -1,
		2, -4, 5;

	vec3 F;
	F << 2, 2, 2;

	//Eigen::ColPivHouseholderQR<Matrix3f> qr;
	Eigen::HouseholderQR<Matrix3f> qr;
	qr.compute(A);
	MatrixXf R = qr.matrixQR().triangularView<Eigen::Upper>();
	MatrixXf Q = qr.householderQ();

	MatrixXf out;
	out = qr.solve(F);
	cout << "QR2(): HouseholderQR---------------------------------------------" << endl;
	cout << "A "