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

Merge remote-tracking branch 'origin/development' into polishing

This commit is contained in:
Martin Diehl 2022-06-15 20:13:39 +02:00
parent 703798b080 0e4a8f284d
commit 0b8473d7ec
3 changed files with 120 additions and 758 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
VERSION
View File

@ -1 +1 @@
3.0.0-alpha6-458-g132823519 3.0.0-alpha6-464-g540eeb4a6

6
src/math.f90
View File

@ -1129,7 +1129,7 @@ pure function math_rotationalPart(F) result(R)
- I_U(1)*I_F(1) * transpose(F) & - I_U(1)*I_F(1) * transpose(F) &
+ I_U(1) * transpose(matmul(F,F)) & + I_U(1) * transpose(matmul(F,F)) &
- matmul(F,C) - matmul(F,C)
R = R /(I_U(1)*I_U(2)-I_U(3)) R = R*math_det33(R)**(-1.0_pReal/3.0_pReal)
end function math_rotationalPart end function math_rotationalPart
@ -1422,7 +1422,9 @@ subroutine selfTest()
t33_2 = math_rotationalPart(transpose(t33)) t33_2 = math_rotationalPart(transpose(t33))
t33 = math_rotationalPart(t33) t33 = math_rotationalPart(t33)
if (any(dNeq0(matmul(t33_2,t33) - math_I3,tol=1.0e-10_pReal))) & if (any(dNeq0(matmul(t33_2,t33) - math_I3,tol=1.0e-10_pReal))) &
error stop 'math_rotationalPart' error stop 'math_rotationalPart (forward-backward)'
if (dNeq(1.0_pReal,math_det33(math_rotationalPart(t33)),tol=1.0e-10_pReal)) &
error stop 'math_rotationalPart (determinant)'
call random_number(r) call random_number(r)
d = int(r*5.0_pReal) + 1 d = int(r*5.0_pReal) + 1

870
src/rotations.f90
View File

File diff suppressed because it is too large Load Diff