GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-07-26 06:44:27 +08:00
Code Issues Packages Projects Releases Wiki Activity

fix compilation of LU class example

Browse Source
This commit is contained in:
Benoit Jacob 2011-01-24 07:41:47 -05:00
parent 1dabd133cc
commit 5331fa3033
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
doc/snippets/class_FullPivLU.cpp
View File

@ -7,10 +7,10 @@ cout << "Here is, up to permutations, its LU decomposition matrix:"
<< endl << lu.matrixLU() << endl; << endl << lu.matrixLU() << endl;
cout << "Here is the L part:" << endl; cout << "Here is the L part:" << endl;
Matrix5x5 l = Matrix5x5::Identity(); Matrix5x5 l = Matrix5x5::Identity();
l.block<5,3>(0,0).part<StrictlyLower>() = lu.matrixLU(); l.block<5,3>(0,0).triangularView<StrictlyLower>() = lu.matrixLU();
cout << l << endl; cout << l << endl;
cout << "Here is the U part:" << endl; cout << "Here is the U part:" << endl;
Matrix5x3 u = lu.matrixLU().part<Upper>(); Matrix5x3 u = lu.matrixLU().triangularView<Upper>();
cout << u << endl; cout << u << endl;
cout << "Let us now reconstruct the original matrix m:" << endl; cout << "Let us now reconstruct the original matrix m:" << endl;
cout << lu.permutationP().inverse() * l * u * lu.permutationQ().inverse() << endl; cout << lu.permutationP().inverse() * l * u * lu.permutationQ().inverse() << endl;