From 7677bc876b4ce3ae04fc77cf86355bb2e7fa47a0 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Sun, 15 Mar 2020 15:58:25 +0100 Subject: [PATCH] more precise documentation symmetric is not enough, could result in taking the square root of negative eigenvalues --- src/math.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/math.f90 b/src/math.f90 index f1dda4f79..20ca9aaef 100644 --- a/src/math.f90 +++ b/src/math.f90 @@ -962,12 +962,12 @@ end subroutine math_eigh33 !-------------------------------------------------------------------------------------------------- -!> @brief eigenvector basis of symmetric 3x3 matrix +!> @brief eigenvector basis of positive-definite 3x3 matrix !-------------------------------------------------------------------------------------------------- pure function math_eigenvectorBasisSym33(m) real(pReal), dimension(3,3) :: math_eigenvectorBasisSym33 - real(pReal), dimension(3,3), intent(in) :: m !< quadratic matrix of which the eigenvector basis is computed + real(pReal), dimension(3,3), intent(in) :: m !< positive-definite matrix of which the basis is computed real(pReal), dimension(3) :: I, v real(pReal) :: P, Q, rho, phi