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_fun.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_fun.F')
| -rw-r--r-- | src/AHFinder_fun.F | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/src/AHFinder_fun.F b/src/AHFinder_fun.F index e565782..0af446c 100644 --- a/src/AHFinder_fun.F +++ b/src/AHFinder_fun.F @@ -99,7 +99,7 @@ do l=1,lmax - dsqrt_factor2(l) = c0(l)*dsqrt(dble(2*l+1)) + dsqrt_factor2(l) = c0(l)*sqrt(dble(2*l+1)) lm_factor1(l,0) = dble(2*l-1)/dble(l) lm_factor2(l,0) = dble(l-1)/dble(l) @@ -113,11 +113,11 @@ do l=1,lmax - dsqrt_factor(l,1) = dsqrt(2.0D0*dble(2*l+1)/dble(l*(l+1))) + dsqrt_factor(l,1) = sqrt(2.0D0*dble(2*l+1)/dble(l*(l+1))) do m=2,l dsqrt_factor(l,m) = dsqrt_factor(l,m-1) - . / dsqrt(dble((l-m+1)*(l+m))) + . / sqrt(dble((l-m+1)*(l+m))) end do end do @@ -237,7 +237,7 @@ ! Find rp. - rp = dsqrt(xp**2 + yp**2 + zp**2) + rp = sqrt(xp**2 + yp**2 + zp**2) ! Monopole term. @@ -281,7 +281,7 @@ ! if (nonaxi) then - aux2 = dsqrt((one-cost)*(one+cost)) + aux2 = sqrt((one-cost)*(one+cost)) ! This general loop can only be used until l=lmax-2 since ! the recursion relations requires 2 starting values. |