From 0ae33d296c4c7f2db69729f8410cf80d81e2b7a5 Mon Sep 17 00:00:00 2001 From: schnetter Date: Fri, 21 Nov 2003 13:29:20 +0000 Subject: Simplify routine to make it work with gfortran. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinUtils/TGRtensor/trunk@10 b716e942-a2de-43ad-8f52-f3dfe468e4e7 --- src/rotation.F90 | 5 +++-- 1 file 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 -- cgit v1.2.3