diff options
author | Jitse Niesen <jitse@maths.leeds.ac.uk> | 2011-05-10 09:59:58 +0100 |
---|---|---|
committer | Jitse Niesen <jitse@maths.leeds.ac.uk> | 2011-05-10 09:59:58 +0100 |
commit | 0c463a21c4a87e9ff3daf7df365082f7df55e3dd (patch) | |
tree | ecd959faa0a2201d03726b09f8fd4ab44f34460f | |
parent | d7e3c949be71dd26e4554547ebb0aa3673a1ec47 (diff) |
Forgot to 'hg add' example file in last commit.
-rw-r--r-- | unsupported/Eigen/MatrixFunctions | 2 | ||||
-rw-r--r-- | unsupported/doc/examples/MatrixSquareRoot.cpp | 16 |
2 files changed, 17 insertions, 1 deletions
diff --git a/unsupported/Eigen/MatrixFunctions b/unsupported/Eigen/MatrixFunctions index d49350f67..ecdf23cef 100644 --- a/unsupported/Eigen/MatrixFunctions +++ b/unsupported/Eigen/MatrixFunctions @@ -310,7 +310,7 @@ cut. The computation is the same as in the real case, except that the complex Schur decomposition is used to reduce the matrix to a triangular matrix. The theoretical cost is the same. Details are in: -Åke Björck and Scen Hammarling, "A Schur method for the +Åke Björck and Sven Hammarling, "A Schur method for the square root of a matrix", <em>Linear Algebra Appl.</em>, 52/53:127–140, 1983. diff --git a/unsupported/doc/examples/MatrixSquareRoot.cpp b/unsupported/doc/examples/MatrixSquareRoot.cpp new file mode 100644 index 000000000..88e7557d7 --- /dev/null +++ b/unsupported/doc/examples/MatrixSquareRoot.cpp @@ -0,0 +1,16 @@ +#include <unsupported/Eigen/MatrixFunctions> +#include <iostream> + +using namespace Eigen; + +int main() +{ + const double pi = std::acos(-1.0); + + MatrixXd A(2,2); + A << cos(pi/3), -sin(pi/3), + sin(pi/3), cos(pi/3); + std::cout << "The matrix A is:\n" << A << "\n\n"; + std::cout << "The matrix square root of A is:\n" << A.sqrt() << "\n\n"; + std::cout << "The square of the last matrix is:\n" << A.sqrt() * A.sqrt() << "\n"; +} |