diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-07-06 18:11:12 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-07-06 18:11:12 +0000 |
| commit | 669d5b37c19d16380ba6c63a7195cbcea8d08f35 (patch) | |
| tree | 43c35e9b941be21d80ba34711d3bbfe7825b3b6d /src/GRHydro_BvecfromAvec.F90 | |
| parent | fcbaf698714e9594ef35fbba18530a4a1ab164ba (diff) | |
GRHydro: add grid function for sqrt(detg)
* add new 1 tl grid function sdetg that stores the sqrt of the determinent
of the 3-metric.
* replace lots of re-computation of det by use of this grid function
From: Christian Ott <cott@bethe.tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@555 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_BvecfromAvec.F90')
| -rw-r--r-- | src/GRHydro_BvecfromAvec.F90 | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/src/GRHydro_BvecfromAvec.F90 b/src/GRHydro_BvecfromAvec.F90 index ea0a58d..30da906 100644 --- a/src/GRHydro_BvecfromAvec.F90 +++ b/src/GRHydro_BvecfromAvec.F90 @@ -53,7 +53,7 @@ subroutine GRHydro_BvecfromAvec(CCTK_ARGUMENTS) implicit none CCTK_REAL :: Az_y, Ay_z, Ax_z, Az_x, Ay_x, Ax_y - CCTK_REAL :: det, isdet + CCTK_REAL :: sdet, isdet CCTK_REAL :: dx, dy, dz, idx, idy, idz CCTK_INT :: i, j, k, nx, ny, nz, GRHydro_UseGeneralCoordinates @@ -102,7 +102,7 @@ subroutine GRHydro_BvecfromAvec(CCTK_ARGUMENTS) call CCTK_WARN(1,"Bvec from Avec start."); - !$OMP PARALLEL DO PRIVATE( det, isdet, local_spatial_order, & + !$OMP PARALLEL DO PRIVATE( sdet, isdet, local_spatial_order, & !$OMP Az_y, Ay_z, Ax_z, Az_x, Ay_x, Ax_y) do k = GRHydro_stencil, cctk_lsh(3)-GRHydro_stencil+1 do j = GRHydro_stencil, cctk_lsh(2)-GRHydro_stencil+1 @@ -129,8 +129,8 @@ subroutine GRHydro_BvecfromAvec(CCTK_ARGUMENTS) Ax_y = DIFF_Y_4(Avecx) end if - det = SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),gyy(i,j,k),gyz(i,j,k),gzz(i,j,k)) - isdet = 1.d0/sqrt(det) + sdet = sdetg(i,j,k) + isdet = 1.d0/sdet Bvecx(i,j,k) = isdet*( Az_y - Ay_z ) Bvecy(i,j,k) = isdet*( Ax_z - Az_x ) |