diff options
| author | schnetter <schnetter@e296648e-0e4f-0410-bd07-d597d9acff87> | 2010-04-03 04:41:35 +0000 |
|---|---|---|
| committer | schnetter <schnetter@e296648e-0e4f-0410-bd07-d597d9acff87> | 2010-04-03 04:41:35 +0000 |
| commit | e12d4fce536f8bf53758c99075011f2760919199 (patch) | |
| tree | 8a074403608636f3f2bcc4d9b88c68a52b827b67 | |
| parent | 25601b8fa42d5918201fb9f32abec5f109affcd2 (diff) | |
Correct OpenMP parallelisation of initialisation of boost and rotation
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/Exact/trunk@253 e296648e-0e4f-0410-bd07-d597d9acff87
| -rw-r--r-- | src/boost.F77 | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/src/boost.F77 b/src/boost.F77 index dc26340..babe398 100644 --- a/src/boost.F77 +++ b/src/boost.F77 @@ -53,9 +53,6 @@ c output arguments $ guxx, guyy, guzz, guxy, guyz, guxz, $ psi -c intrinsic functions called - CCTK_REAL sqrt - c static local variables describing Lorentz transformation logical firstcall data firstcall /.true./ @@ -68,9 +65,8 @@ c static local variables describing Lorentz transformation CCTK_REAL R(0:3,0:3) save firstcall save gamma - save vv, nn + save vv save parallel, perp - save Cx_par, Cx_perp save partial_Mx_wrt_Cx save partial_Cx_wrt_Mx save R @@ -145,7 +141,8 @@ c c compute Lorentz transformation information on first call c if (firstcall) then - firstcall = .false. +c$omp critical + if (firstcall) then c boost velocity vv(1) = boost_vx @@ -295,6 +292,9 @@ c Mbetax(3) = Nbetax(1) c Mbetax(i) + R(i,j) Nbetax(j) [wrong] c Mbetax(i) + R(j,i) Nbetax(j) [correct] + firstcall = .false. + end if +c$omp end critical end if cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |