bessel_6.f90: Reduce mymax for Tru64.

2011-04-03  Tobias Burnus  <burnus@net-b.de>

       * gfortran.dg/bessel_6.f90: Reduce mymax for Tru64.

From-SVN: r171912
This commit is contained in:
Tobias Burnus 2011-04-03 18:04:25 +02:00 committed by Tobias Burnus
parent 99ebea0382
commit 3aa8b384e3
2 changed files with 8 additions and 3 deletions

View File

@ -1,3 +1,7 @@
2011-04-03 Tobias Burnus <burnus@net-b.de>
* gfortran.dg/bessel_6.f90: Reduce mymax for Tru64.
2011-04-03 Richard Guenther <rguenther@suse.de>
Ira Rosen <ira.rosen@linaro.org>

View 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