diff options
| author | tradke <tradke@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2002-06-18 13:00:08 +0000 |
|---|---|---|
| committer | tradke <tradke@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2002-06-18 13:00:08 +0000 |
| commit | 751bd44faded4ac77cb9029cdacd832b0cb8bec8 (patch) | |
| tree | e5928adc750a6d2743db33e31f4135ec6359204a /src/AHFinder_gau.F | |
| parent | c0a5e1b582bcf8b481c08c6f1f04f40fac278a3b (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.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@312 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src/AHFinder_gau.F')
| -rw-r--r-- | src/AHFinder_gau.F | 16 |
1 files changed, 8 insertions, 8 deletions
diff --git a/src/AHFinder_gau.F b/src/AHFinder_gau.F index 0381d3a..9e8ca67 100644 --- a/src/AHFinder_gau.F +++ b/src/AHFinder_gau.F @@ -809,7 +809,7 @@ do j=1,nphi circ_eq = circ_eq + dphi - . *dsqrt(dabs(half*(g22(ntheta+1,j) + g22(ntheta+1,j+1)))) + . *sqrt(abs(half*(g22(ntheta+1,j) + g22(ntheta+1,j+1)))) end do else @@ -819,12 +819,12 @@ aux = aux - int(aux) if (cartoon) then - circ_eq = 6.283185307D0*dsqrt(g22(i,0)) + circ_eq = 6.283185307D0*sqrt(g22(i,0)) else do j=1,nphi circ_eq = circ_eq + dphi - . *((one-aux)*dsqrt(dabs(half*(g22(i,j) + g22(i,j+1)))) - . + aux*dsqrt(dabs(half*(g22(i+1,j) + g22(i+1,j+1))))) + . *((one-aux)*sqrt(abs(half*(g22(i,j) + g22(i,j+1)))) + . + aux*sqrt(abs(half*(g22(i+1,j) + g22(i+1,j+1))))) end do end if end if @@ -833,7 +833,7 @@ do i=1,ntheta meri_p1 = meri_p1 + dtheta - . *dsqrt(dabs(half*(g11(i,1) + g11(i+1,1)))) + . *sqrt(abs(half*(g11(i,1) + g11(i+1,1)))) end do if (cartoon) then @@ -841,7 +841,7 @@ else if (refx.and.refy) then do i=1,ntheta meri_p2 = meri_p2 + dtheta - . *dsqrt(dabs(half*(g11(i,nphi+1) + g11(i+1,nphi+1)))) + . *sqrt(abs(half*(g11(i,nphi+1) + g11(i+1,nphi+1)))) end do else aux = half*pi/dphi @@ -849,8 +849,8 @@ aux = aux - int(aux) do i=1,ntheta meri_p2 = meri_p2 + dtheta - . *((one-aux)*dsqrt(dabs(half*(g11(i,j) + g11(i,j)))) - . + aux*dsqrt(dabs(half*(g11(i,j+1) + g11(i,j+1))))) + . *((one-aux)*sqrt(abs(half*(g11(i,j) + g11(i,j)))) + . + aux*sqrt(abs(half*(g11(i,j+1) + g11(i,j+1))))) end do end if |