numerically better normalization

slightly more expensive, but guarantees a det(R) = 1.0
2022-06-11 23:21:34 +02:00 · 2022-06-11 23:21:34 +02:00 · f60c950a7c
parent 9f4e354b12
commit f60c950a7c
1 changed files with 1 additions and 1 deletions
--- 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