OSDN Git Service

2011-04-03 Tobias Burnus <burnus@net-b.de>
[pf3gnuchains/gcc-fork.git] / gcc / testsuite / gfortran.dg / bessel_6.f90
index e760d4c..bfbf1e1 100644 (file)
@@ -16,7 +16,7 @@ real,parameter :: myeps(size(values)) = epsilon(0.0) &
 ! more tolerant
 !                  * [0, 5, 3, 4, 6, 7, 7, 5, 5, 6, 66, 4 ]
 integer,parameter :: mymax(size(values)) =  &
-                 [100, 17, 23, 21, 27, 28, 32, 35, 36, 41, 47, 37 ]
+                 [100, 17, 23, 21, 27, 28, 32, 35, 31, 41, 47, 37 ]
 integer, parameter :: Nmax = 100
 real :: rec(0:Nmax), lib(0:Nmax)
 integer :: i
@@ -39,8 +39,9 @@ lib(0:mymax) = [ (BESSEL_JN(i, X), i=0,mymax) ]
 do i = 0, mymax
 !  print '(i2,2e17.9,e12.2,f18.10,2l3)', i, rec(i), lib(i), &
 !        rec(i)-lib(i),           ((rec(i)-lib(i))/rec(i))/epsilon(x), &
-!        rec(i) == lib(i) .or. abs((rec(i)-lib(i))/rec(i)) < myeps
-if (.not. (rec(i) == lib(i) .or. abs((rec(i)-lib(i))/rec(i)) < myeps)) &
+!        rec(i) == lib(i), abs((rec(i)-lib(i))/rec(i)) < myeps
+if (rec(i) == lib(i)) CYCLE
+if (abs((rec(i)-lib(i))/rec(i)) > myeps) &
   call abort()
 end do