| Benoit Jacob | f04c1cb | 2008-08-11 21:26:37 +0000 | [diff] [blame] | 1 | MatrixXf m = MatrixXf::Random(3,5); |
| 2 | cout << "Here is the matrix m:" << endl << m << endl; | ||||
| 3 | MatrixXf ker = m.lu().kernel(); | ||||
| 4 | cout << "Here is a matrix whose columns form a basis of the kernel of m:" | ||||
| 5 | << endl << ker << endl; | ||||
| 6 | cout << "By definition of the kernel, m*ker is zero:" | ||||
| 7 | << endl << m*ker << endl; | ||||