diff --git a/src/math.f90 b/src/math.f90
index de27519e8..7e35ca390 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(m,n,o,p) &
-                                      + rot_tensor(m,i) * rot_tensor(n,j) &
-                                      * rot_tensor(o,k) * rot_tensor(p,l) * tensor(i,j,k,l)
+     math_rotate_forward3333(i,j,k,l) = math_rotate_forward3333(i,j,k,l) &
+                                      + rot_tensor(i,m) * rot_tensor(j,n) &
+                                      * rot_tensor(k,o) * rot_tensor(l,p) * tensor(m,n,o,p)
  enddo; enddo; enddo; enddo; enddo; enddo; enddo; enddo
 
 end function math_rotate_forward3333