summaryrefslogtreecommitdiffhomepage
path: root/eigen/doc/snippets/RealQZ_compute.cpp
diff options
context:
space:
mode:
authorStanislaw Halik <sthalik@misaki.pl>2019-03-03 21:09:10 +0100
committerStanislaw Halik <sthalik@misaki.pl>2019-03-03 21:10:13 +0100
commitf0238cfb6997c4acfc2bd200de7295f3fa36968f (patch)
treeb215183760e4f615b9c1dabc1f116383b72a1b55 /eigen/doc/snippets/RealQZ_compute.cpp
parent543edd372a5193d04b3de9f23c176ab439e51b31 (diff)
don't index Eigen
Diffstat (limited to 'eigen/doc/snippets/RealQZ_compute.cpp')
-rw-r--r--eigen/doc/snippets/RealQZ_compute.cpp17
1 files changed, 0 insertions, 17 deletions
diff --git a/eigen/doc/snippets/RealQZ_compute.cpp b/eigen/doc/snippets/RealQZ_compute.cpp
deleted file mode 100644
index a18da42..0000000
--- a/eigen/doc/snippets/RealQZ_compute.cpp
+++ /dev/null
@@ -1,17 +0,0 @@
-MatrixXf A = MatrixXf::Random(4,4);
-MatrixXf B = MatrixXf::Random(4,4);
-RealQZ<MatrixXf> qz(4); // preallocate space for 4x4 matrices
-qz.compute(A,B); // A = Q S Z, B = Q T Z
-
-// print original matrices and result of decomposition
-cout << "A:\n" << A << "\n" << "B:\n" << B << "\n";
-cout << "S:\n" << qz.matrixS() << "\n" << "T:\n" << qz.matrixT() << "\n";
-cout << "Q:\n" << qz.matrixQ() << "\n" << "Z:\n" << qz.matrixZ() << "\n";
-
-// verify precision
-cout << "\nErrors:"
- << "\n|A-QSZ|: " << (A-qz.matrixQ()*qz.matrixS()*qz.matrixZ()).norm()
- << ", |B-QTZ|: " << (B-qz.matrixQ()*qz.matrixT()*qz.matrixZ()).norm()
- << "\n|QQ* - I|: " << (qz.matrixQ()*qz.matrixQ().adjoint() - MatrixXf::Identity(4,4)).norm()
- << ", |ZZ* - I|: " << (qz.matrixZ()*qz.matrixZ().adjoint() - MatrixXf::Identity(4,4)).norm()
- << "\n";