------- Comment #1 from dominiq at lps dot ens dot fr 2010-08-21 18:08 ------- Could you try the following patch
--- /opt/gcc/_clean/gcc/testsuite/gfortran.dg/bessel_6.f90 2010-08-21 12:18:37.000000000 +0200 +++ bessel_6_db.f90 2010-08-21 15:05:20.000000000 +0200 @@ -8,7 +8,7 @@ implicit none real,parameter :: values(*) = [0.0, 0.5, 1.0, 0.9, 1.8,2.0,3.0,4.0,4.25,8.0,34.53, 475.78] real,parameter :: myeps(size(values)) = epsilon(0.0) & - * [2, 7, 5, 6, 9, 12, 12, 7, 7, 8, 75, 6 ] + * [2, 7, 5, 6, 9, 12, 12, 7, 7, 8, 120, 15 ] ! The following is sufficient for me - the values above are a bit ! more tolerant ! * [0, 5, 3, 4, 6, 7, 7, 5, 5, 6, 66, 4 ] --- /opt/gcc/_clean/gcc/testsuite/gfortran.dg/bessel_7.f90 2010-08-21 12:18:37.000000000 +0200 +++ bessel_7_db.f90 2010-08-21 20:04:16.000000000 +0200 @@ -8,12 +8,12 @@ implicit none real,parameter :: values(*) = [0.0, 0.5, 1.0, 0.9, 1.8,2.0,3.0,4.0,4.25,8.0,34.53, 475.78] real,parameter :: myeps(size(values)) = epsilon(0.0) & - * [2, 2, 2, 5, 5, 2, 12, 2, 4, 3, 30, 130 ] + * [2, 3, 2, 5, 7, 2, 12, 3, 7, 3, 30, 144 ] ! The following is sufficient for me - the values above are a bit ! more tolerant ! * [0, 0, 0, 3, 3, 0, 9, 0, 2, 1, 22, 130 ] integer,parameter :: nit(size(values)) = & - [100, 100, 100, 100, 100, 100, 10, 100, 100, 100, 10, 25 ] + [100, 100, 100, 25, 15, 100, 10, 100, 15, 100, 7, 25 ] integer, parameter :: Nmax = 100 real :: rec(0:Nmax), lib(0:Nmax) integer :: i -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=45367
