From 4d5305354d28689d8a55f4ccc14c28d06a08d4d0 Mon Sep 17 00:00:00 2001 From: pollney Date: Tue, 18 Jun 2002 13:17:46 +0000 Subject: 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 --- src/findnormals.F90 | 4 ++-- 1 file 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. -- cgit v1.2.3