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_mask.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_mask.F')
| -rw-r--r-- | src/AHFinder_mask.F | 16 |
1 files changed, 8 insertions, 8 deletions
diff --git a/src/AHFinder_mask.F b/src/AHFinder_mask.F index 194f13c..cef8bbb 100644 --- a/src/AHFinder_mask.F +++ b/src/AHFinder_mask.F @@ -84,7 +84,7 @@ dd = max(dx,dy,dz) aux = ahf_maskbuffer*dd - if (dabs(buffer).lt.aux) buffer = - aux + if (abs(buffer).lt.aux) buffer = - aux ! ***************************************** @@ -101,7 +101,7 @@ do i=1,nx if (ahfgrid(i,j,k).ge.buffer) then - aux = dsqrt((x(i,j,k)-xc)**2 + (y(i,j,k)-yc)**2 + aux = sqrt((x(i,j,k)-xc)**2 + (y(i,j,k)-yc)**2 . + (z(i,j,k)-zc)**2) if (rhor.gt.aux) rhor = aux end if @@ -167,9 +167,9 @@ aux = 0.57D0*rhor - xa = dabs(x(i,j,k)-xc) - ya = dabs(y(i,j,k)-yc) - za = dabs(z(i,j,k)-zc) + xa = abs(x(i,j,k)-xc) + ya = abs(y(i,j,k)-yc) + za = abs(z(i,j,k)-zc) if ((xa.lt.aux).and.(ya.lt.aux).and.(za.lt.aux)) then ahmask(i,j,k) = zero @@ -202,9 +202,9 @@ aux = 0.7D0*rhor - xa = dabs(x(i,j,k)-xc) - ya = dabs(y(i,j,k)-yc) - za = dabs(z(i,j,k)-zc) + xa = abs(x(i,j,k)-xc) + ya = abs(y(i,j,k)-yc) + za = abs(z(i,j,k)-zc) if ((xa.lt.aux).and.(ya.lt.aux).and.(za.lt.aux).and. . ((xa+ya+za).lt.2.0D0*aux)) then |