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

polishing

This commit is contained in:
Martin Diehl 2016-02-02 08:45:47 +01:00
parent eaf9b41b7a
commit bf04ee60f0
1 changed files with 1 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
code/math.f90
View File

@ -2066,7 +2066,7 @@ function math_eigenvaluesSym33(m)
c(2) = - math_trace33(m) c(2) = - math_trace33(m)
c(1) = m(1,1)*m(2,2) + m(1,1)*m(3,3) + m(2,2)*m(3,3) & c(1) = m(1,1)*m(2,2) + m(1,1)*m(3,3) + m(2,2)*m(3,3) &
-(m(1,2)**2 + m(2,3)**2 + m(1,3)**2) -(m(1,2)**2 + m(1,3)**2 + m(2,3)**2)
c(0) = m(1,1)*m(2,3)**2 + m(2,2)*m(1,3)**2 + m(3,3)*m(1,2)**2 & c(0) = m(1,1)*m(2,3)**2 + m(2,2)*m(1,3)**2 + m(3,3)*m(1,2)**2 &
-(m(1,1)*m(2,2)*m(3,3) + 2.0_pReal * m(1,3)*m(1,2)*m(2,3)) -(m(1,1)*m(2,2)*m(3,3) + 2.0_pReal * m(1,3)*m(1,2)*m(2,3))
@ -2081,8 +2081,6 @@ function math_eigenvaluesSym33(m)
-cos(phi)+sqrt(1.0_pReal/3.0_pReal)*sin(phi) & -cos(phi)+sqrt(1.0_pReal/3.0_pReal)*sin(phi) &
] * sqrt(p) -c(2))/3.0_pReal ] * sqrt(p) -c(2))/3.0_pReal
end function math_eigenvaluesSym33 end function math_eigenvaluesSym33
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------