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_pow.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_pow.F')
| -rw-r--r-- | src/AHFinder_pow.F | 18 |
1 files changed, 9 insertions, 9 deletions
diff --git a/src/AHFinder_pow.F b/src/AHFinder_pow.F index 149a474..fe53079 100644 --- a/src/AHFinder_pow.F +++ b/src/AHFinder_pow.F @@ -137,8 +137,8 @@ ! And record it if it is the largest decrease so far. - if (dabs(FPTT-FRET).gt.DEL) then - DEL = dabs(FPTT-FRET) + if (abs(FPTT-FRET).gt.DEL) then + DEL = abs(FPTT-FRET) IBIG = I end if @@ -148,12 +148,12 @@ ! Termination criterion. - TOL1 = 0.5D0*FTOL*(dabs(FP)+dabs(FRET)) + ZEPS + TOL1 = 0.5D0*FTOL*(abs(FP)+abs(FRET)) + ZEPS - if (dabs(FP-FRET).le.TOL1) then + if (abs(FP-FRET).le.TOL1) then found = .true. return - else if (dabs(FRET).lt.ZEPS) then + else if (abs(FRET).lt.ZEPS) then found = .true. return end if @@ -423,7 +423,7 @@ S = (XB-XA)*(FB-FC) Q = (XB-XC)*(FB-FA) U = XB - ((XB-XC)*Q - (XB-XA)*S) - . / (2.0D0*SIGN(MAX(DABS(Q-S),TINY),Q-S)) + . / (2.0D0*SIGN(MAX(ABS(Q-S),TINY),Q-S)) ! We will not go further than this. @@ -584,8 +584,8 @@ ! Test for done here. - TOL1 = 0.5D0*TOL*(dabs(FP)+dabs(FX)) + ZEPS - if (dabs(FP-FX).le.TOL1) goto 3 + TOL1 = 0.5D0*TOL*(abs(FP)+abs(FX)) + ZEPS + if (abs(FP-FX).le.TOL1) goto 3 ! Construct a trial parabolic fit. @@ -612,7 +612,7 @@ ! Check if X was already the minimum. - if (dabs(U-XX).lt.ZEPS) goto 3 + if (abs(U-XX).lt.ZEPS) goto 3 ! This is the one function evaluation per iteration. |