| // SPDX-FileCopyrightText: The Eigen Authors | |
| // SPDX-License-Identifier: MPL-2.0 | |
| MatrixXd X = MatrixXd::Random(4, 4); | |
| MatrixXd A = X * X.transpose(); | |
| cout << "Here is a random positive-definite matrix, A:" << endl << A << endl << endl; | |
| SelfAdjointEigenSolver<MatrixXd> es(A); | |
| cout << "The inverse square root of A is: " << endl; | |
| cout << es.operatorInverseSqrt() << endl; | |
| cout << "We can also compute it with operatorSqrt() and inverse(). That yields: " << endl; | |
| cout << es.operatorSqrt().inverse() << endl; |