aboutsummaryrefslogtreecommitdiff
path: root/src/AHFinder_fun.F
diff options
context:
space:
mode:
authormiguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c>2000-10-13 08:44:14 +0000
committermiguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c>2000-10-13 08:44:14 +0000
commit09b19f9e686db7098d93ae961066b14273e0be5f (patch)
tree639a5d9b75a52e47a299f394701cede2d4d27c91 /src/AHFinder_fun.F
parent04bde353cfefd2239092c24eb8234d7f11066483 (diff)
Optimizing.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@126 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src/AHFinder_fun.F')
-rw-r--r--src/AHFinder_fun.F14
1 files changed, 1 insertions, 13 deletions
diff --git a/src/AHFinder_fun.F b/src/AHFinder_fun.F
index d09455f..a6a426e 100644
--- a/src/AHFinder_fun.F
+++ b/src/AHFinder_fun.F
@@ -100,13 +100,12 @@
! Non-axisymmetric terms. Notice how the sum over M is first.
! This will allow me to use the recursion relations to avoid
! having to start from scratch every time. Also, I sum over
-! all l's even if I don't want some terms. This is because
+! all l even if I do not want some terms. This is because
! in order to use the recursion relations I need all polynomials.
if (nonaxi) then
do m=1,lmax
aux2 = dble(m)*phi
-
sina = sin(aux2)
cosa = cos(aux2)
do l=m,lmax
@@ -121,17 +120,6 @@
ahfgrid(i,j,k) = rp - aux1
-! Check for NaNs.
-
- if (ahfgrid(i,j,k).ge.zero) then
- else if (ahfgrid(i,j,k).lt.zero) then
- else
- if (myproc.eq.0) then
- write(*,*) 'NaN in ahfgrid at point',i,j,k
- end if
- STOP
- end if
-
end do
end do
end do