diff options
author | Benoit Jacob <jacob.benoit.1@gmail.com> | 2010-06-30 10:11:55 -0400 |
---|---|---|
committer | Benoit Jacob <jacob.benoit.1@gmail.com> | 2010-06-30 10:11:55 -0400 |
commit | 4d4a23cd3e97bbc10d6307b828627c672097c5f5 (patch) | |
tree | edccff8f095690ba489a55d4fee79bf6ba78b00d /doc/examples/TutorialLinAlgExComputeSolveError.cpp | |
parent | b1741c1dc6c82737b69b20edfd817ef7e11ef28d (diff) |
nearly complete page 6 / linear algebra + examples
fix the previous/next links
Diffstat (limited to 'doc/examples/TutorialLinAlgExComputeSolveError.cpp')
-rw-r--r-- | doc/examples/TutorialLinAlgExComputeSolveError.cpp | 14 |
1 files changed, 14 insertions, 0 deletions
diff --git a/doc/examples/TutorialLinAlgExComputeSolveError.cpp b/doc/examples/TutorialLinAlgExComputeSolveError.cpp new file mode 100644 index 000000000..f362fb71a --- /dev/null +++ b/doc/examples/TutorialLinAlgExComputeSolveError.cpp @@ -0,0 +1,14 @@ +#include <iostream> +#include <Eigen/Dense> + +using namespace std; +using namespace Eigen; + +int main() +{ + MatrixXd A = MatrixXd::Random(100,100); + MatrixXd b = MatrixXd::Random(100,50); + MatrixXd x = A.fullPivLu().solve(b); + double relative_error = (A*x - b).norm() / b.norm(); // norm() is L2 norm + cout << "The relative error is:\n" << relative_error << endl; +} |