eigen/doc/snippets/FullPivLU_image.cpp

Matrix3d m;
m << 1,1,0,
     1,3,2,
     0,1,1;
cout << "Here is the matrix m:" << endl << m << endl;
cout << "Notice that the middle column is the sum of the two others, so the "
     << "columns are linearly dependent." << endl;
cout << "Here is a matrix whose columns have the same span but are linearly independent:"
     << endl << m.fullPivLu().image(m) << endl;
* fix in IO.h, a useless copy was made because of assignment from Derived to MatrixBase. * the optimization of eval() for Matrix now consists in a partial specialization of ei_eval, which returns a reference type for Matrix. No overriding of eval() in Matrix anymore. Consequence: careful, ei_eval is no longer guaranteed to give a plain matrix type! For that, use ei_plain_matrix_type, or the PlainMatrixType typedef. * so lots of changes to adapt to that everywhere. Hope this doesn't break (too much) MSVC compilation. * add code examples for the new image() stuff. * lower a bit the precision for floats in the unit tests as we were already doing some workarounds in inverse.cpp and we got some failed tests. 2008-12-19 04:36:25 +08:00			`Matrix3d m;`
			`m << 1,1,0,`
			`1,3,2,`
			`0,1,1;`
			`cout << "Here is the matrix m:" << endl << m << endl;`
			`cout << "Notice that the middle column is the sum of the two others, so the "`
			`<< "columns are linearly dependent." << endl;`
			`cout << "Here is a matrix whose columns have the same span but are linearly independent:"`
simplifications in the ei_solve_impl system, factor out some boilerplate code 2009-11-09 05:51:41 +08:00			`<< endl << m.fullPivLu().image(m) << endl;`