FullPivLU_kernel.cpp revision c981c48f5bc9aefeffc0bcb0cc3934c2fae179dd
1MatrixXf m = MatrixXf::Random(3,5); 2cout << "Here is the matrix m:" << endl << m << endl; 3MatrixXf ker = m.fullPivLu().kernel(); 4cout << "Here is a matrix whose columns form a basis of the kernel of m:" 5 << endl << ker << endl; 6cout << "By definition of the kernel, m*ker is zero:" 7 << endl << m*ker << endl; 8