TutorialLinAlgSVDSolve.cpp revision c981c48f5bc9aefeffc0bcb0cc3934c2fae179dd
191bca4244286fb519c93fe92329da96b0e6f32eejvr#include <iostream> 291bca4244286fb519c93fe92329da96b0e6f32eejvr#include <Eigen/Dense> 305b4b4a27160e90307372f85dd99be69a9d972ffjvr 405b4b4a27160e90307372f85dd99be69a9d972ffjvrusing namespace std; 510de259aec75d3ac0c34b444b2f0423fa86a7709jvrusing namespace Eigen; 691bca4244286fb519c93fe92329da96b0e6f32eejvr 791bca4244286fb519c93fe92329da96b0e6f32eejvrint main() 891bca4244286fb519c93fe92329da96b0e6f32eejvr{ 991bca4244286fb519c93fe92329da96b0e6f32eejvr MatrixXf A = MatrixXf::Random(3, 2); 1091bca4244286fb519c93fe92329da96b0e6f32eejvr cout << "Here is the matrix A:\n" << A << endl; 1191bca4244286fb519c93fe92329da96b0e6f32eejvr VectorXf b = VectorXf::Random(3); 1291bca4244286fb519c93fe92329da96b0e6f32eejvr cout << "Here is the right hand side b:\n" << b << endl; 1391bca4244286fb519c93fe92329da96b0e6f32eejvr cout << "The least-squares solution is:\n" 1491bca4244286fb519c93fe92329da96b0e6f32eejvr << A.jacobiSvd(ComputeThinU | ComputeThinV).solve(b) << endl; 1510de259aec75d3ac0c34b444b2f0423fa86a7709jvr} 1605b4b4a27160e90307372f85dd99be69a9d972ffjvr