diff options
author | schnetter <schnetter@b716e942-a2de-43ad-8f52-f3dfe468e4e7> | 2003-11-21 13:29:20 +0000 |
---|---|---|
committer | schnetter <schnetter@b716e942-a2de-43ad-8f52-f3dfe468e4e7> | 2003-11-21 13:29:20 +0000 |
commit | 0ae33d296c4c7f2db69729f8410cf80d81e2b7a5 (patch) | |
tree | 6556d1507996797730721a203ddda730b2a5b253 | |
parent | a63f47a0af8c4f78cf3ae4d2b6d4e0940caaea59 (diff) |
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
-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 |