diff --git a/src/math.f90 b/src/math.f90 index 25c90ccf4..b92e1870b 100644 --- a/src/math.f90 +++ b/src/math.f90 @@ -1129,7 +1129,7 @@ pure function math_rotationalPart(F) result(R) - I_U(1)*I_F(1) * transpose(F) & + I_U(1) * transpose(matmul(F,F)) & - matmul(F,C) - R = R /(I_U(1)*I_U(2)-I_U(3)) + R = R*math_det33(R)**(-1.0/3.0) end function math_rotationalPart