diff options
| -rw-r--r-- | src/rotation.F90 | 5 |
1 files changed, 3 insertions, 2 deletions
diff --git a/src/rotation.F90 b/src/rotation.F90 index 2d8d80b..de3c511 100644 --- a/src/rotation.F90 +++ b/src/rotation.F90 @@ -47,11 +47,12 @@ contains subroutine make_euler (phi, theta, psi, tt) CCTK_REAL, intent(in) :: phi, theta, psi CCTK_REAL, intent(out) :: tt(0:3,0:3) - CCTK_REAL :: tt1(0:3,0:3), tt2(0:3,0:3), tt3(0:3,0:3) + CCTK_REAL :: tt1(0:3,0:3), tt2(0:3,0:3), tt3(0:3,0:3), tt4(0:3,0:3) call make_euler1 (phi, tt1) call make_euler2 (theta, tt2) call make_euler3 (psi, tt3) - tt = matmul(tt3, matmul(tt2, tt1)) + tt4 = matmul(tt2, tt1) + tt = matmul(tt3, tt4) end subroutine make_euler end module rotation |