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

equivalent, but more systematic

This commit is contained in:
Martin Diehl 2020-03-15 15:03:53 +01:00
parent 65decfc48a
commit 253ee0013d
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

@ -1006,8 +1006,8 @@ pure function math_eigenvectorBasisSym33(m)
EB(1:3,1:3,3) = 0.0_pReal EB(1:3,1:3,3) = 0.0_pReal
elseif (abs(v(3)-v(1)) < TOL) then twoSimilarEigVals elseif (abs(v(3)-v(1)) < TOL) then twoSimilarEigVals
EB(1:3,1:3,2) = matmul(N(1:3,1:3,3),N(1:3,1:3,1))/((v(2)-v(3))*(v(2)-v(1))) EB(1:3,1:3,2) = matmul(N(1:3,1:3,3),N(1:3,1:3,1))/((v(2)-v(3))*(v(2)-v(1)))
EB(1:3,1:3,1) = math_I3-EB(1:3,1:3,2) EB(1:3,1:3,3) = math_I3-EB(1:3,1:3,2)
EB(1:3,1:3,3) = 0.0_pReal EB(1:3,1:3,1) = 0.0_pReal
else twoSimilarEigVals else twoSimilarEigVals
EB(1:3,1:3,1) = matmul(N(1:3,1:3,2),N(1:3,1:3,3))/((v(1)-v(2))*(v(1)-v(3))) EB(1:3,1:3,1) = matmul(N(1:3,1:3,2),N(1:3,1:3,3))/((v(1)-v(2))*(v(1)-v(3)))
EB(1:3,1:3,2) = matmul(N(1:3,1:3,1),N(1:3,1:3,3))/((v(2)-v(3))*(v(2)-v(1))) EB(1:3,1:3,2) = matmul(N(1:3,1:3,1),N(1:3,1:3,3))/((v(2)-v(3))*(v(2)-v(1)))