223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

more precise documentation 

symmetric is not enough, could result in taking the square root of negative eigenvalues
This commit is contained in:
Martin Diehl 2020-03-15 15:58:25 +01:00
parent 253ee0013d
commit 7677bc876b
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/math.f90
View File

@ -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) pure function math_eigenvectorBasisSym33(m)
real(pReal), dimension(3,3) :: math_eigenvectorBasisSym33 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), dimension(3) :: I, v
real(pReal) :: P, Q, rho, phi real(pReal) :: P, Q, rho, phi