diff options
| author | pollney <pollney@f75ba9e5-694f-0410-ac2c-87ea7ce7132b> | 2002-06-18 13:17:46 +0000 |
|---|---|---|
| committer | pollney <pollney@f75ba9e5-694f-0410-ac2c-87ea7ce7132b> | 2002-06-18 13:17:46 +0000 |
| commit | 4d5305354d28689d8a55f4ccc14c28d06a08d4d0 (patch) | |
| tree | 92f9742d3154cb613a0bedbd5ba37c39ba0dfa39 | |
| parent | 37053cfb2b313548a3952475b6989858ef7db38a (diff) | |
Use the generic routines sqrt(), abs(), log(), etc. instead of double-precision
ones dsqrt(), dabs(), dlog(), detc.
This makes the code compile in single precision.
(Fix courtesy of ThomasR)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@33 f75ba9e5-694f-0410-ac2c-87ea7ce7132b
| -rw-r--r-- | src/findnormals.F90 | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/findnormals.F90 b/src/findnormals.F90 index 4c963cf..518679c 100644 --- a/src/findnormals.F90 +++ b/src/findnormals.F90 @@ -82,7 +82,7 @@ subroutine excision_findnormals (ierr, mask, dirx, diry, dirz, ni, nj, nk) ! Find magnitude of (sx,sy,sz). - smag = dsqrt(sx**2 + sy**2 + sz**2) + smag = sqrt(sx**2 + sy**2 + sz**2) ! Normalize (sx,sy,sz). @@ -125,7 +125,7 @@ subroutine excision_findnormals (ierr, mask, dirx, diry, dirz, ni, nj, nk) vz = dble(kk) dot = (sx*vx + sy*vy + sz*vz) & - / dsqrt(vx**2 + vy**2 + vz**2) + / sqrt(vx**2 + vy**2 + vz**2) ! If the new product is larger than the ! largest so far, redifine the normal. |