Commit 60f22d5d by Tobias Burnus Committed by Tobias Burnus

re PR fortran/36158 (Transformational function BESSEL_YN(n1,n2,x) and BESSEL_JN missing)

2010-08-22  Tobias Burnus  <burnus@net-b.de>

        PR fortran/36158
        * gfortran.dg/bessel_7.f90: Disable accidently enabled debug
        output.

From-SVN: r163455
parent 8af7a073
2010-08-22 Tobias Burnus <burnus@net-b.de>
PR fortran/36158
* gfortran.dg/bessel_7.f90: Disable accidently enabled debug output.
2010-08-22 Tobias Burnus <burnus@net-b.de>
Dominique d'Humieres <dominiq@lps.ens.fr>
PR fortran/45367
......
......@@ -32,13 +32,13 @@ real X, myeps, myeps2
rec = BESSEL_YN(0, Nmax, X)
lib = [ (BESSEL_YN(i, X), i=0,Nmax) ]
print *, 'YN for X = ', X, ' -- Epsilon = ',epsilon(x)
!print *, 'YN for X = ', X, ' -- Epsilon = ',epsilon(x)
do i = 0, Nmax
print '(i2,2e17.9,e12.2,f14.10,2l3)', i, rec(i), lib(i), &
rec(i)-lib(i), ((rec(i)-lib(i))/rec(i))/epsilon(x), &
i > nit .or. rec(i) == lib(i) &
.or. abs((rec(i)-lib(i))/rec(i)) < myeps2, &
rec(i) == lib(i) .or. abs((rec(i)-lib(i))/rec(i)) < myeps
! print '(i2,2e17.9,e12.2,f14.10,2l3)', i, rec(i), lib(i), &
! rec(i)-lib(i), ((rec(i)-lib(i))/rec(i))/epsilon(x), &
! i > nit .or. rec(i) == lib(i) &
! .or. abs((rec(i)-lib(i))/rec(i)) < myeps2, &
! rec(i) == lib(i) .or. abs((rec(i)-lib(i))/rec(i)) < myeps
if (.not. (i > nit .or. rec(i) == lib(i) &
.or. abs((rec(i)-lib(i))/rec(i)) < myeps2)) &
call abort ()
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment