From cd4115bada88783012b6cd6014457768cb95465f Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Mon, 26 Feb 2018 11:14:02 +0100 Subject: [PATCH] now using the correct rotation formula --- src/math.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/math.f90 b/src/math.f90 index 055ec6c62..de27519e8 100644 --- a/src/math.f90 +++ b/src/math.f90 @@ -2789,9 +2789,9 @@ pure function math_rotate_forward3333(tensor,rot_tensor) do i = 1_pInt,3_pInt; do j = 1_pInt,3_pInt; do k = 1_pInt,3_pInt; do l = 1_pInt,3_pInt do m = 1_pInt,3_pInt; do n = 1_pInt,3_pInt; do o = 1_pInt,3_pInt; do p = 1_pInt,3_pInt - math_rotate_forward3333(m,n,o,p) = math_rotate_forward3333(i,j,k,l) & + math_rotate_forward3333(m,n,o,p) = math_rotate_forward3333(m,n,o,p) & + rot_tensor(m,i) * rot_tensor(n,j) & - * rot_tensor(o,k) * rot_tensor(p,l) * tensor(m,n,o,p) + * rot_tensor(o,k) * rot_tensor(p,l) * tensor(i,j,k,l) enddo; enddo; enddo; enddo; enddo; enddo; enddo; enddo end function math_rotate_forward3333