| Message ID | 4C791D2D.30006@net-b.de |
|---|---|
| State | New |
| Headers | show |
On Sat, Aug 28, 2010 at 04:29:01PM +0200, Tobias Burnus wrote: > Committed as obvious. (Rev. 163615). > Thanks.
On Sat, Aug 28, 2010 at 7:29 AM, Tobias Burnus <burnus@net-b.de> wrote: > Committed as obvious. (Rev. 163615). > I still got http://gcc.gnu.org/ml/gcc-testresults/2010-08/msg03035.html: FAIL: gfortran.dg/bessel_7.f90 -O0 execution test FAIL: gfortran.dg/bessel_7.f90 -O1 execution test FAIL: gfortran.dg/bessel_7.f90 -O2 execution test FAIL: gfortran.dg/bessel_7.f90 -O3 -fomit-frame-pointer execution test FAIL: gfortran.dg/bessel_7.f90 -O3 -fomit-frame-pointer -funroll-loops execution test FAIL: gfortran.dg/bessel_7.f90 -O3 -fomit-frame-pointer -funroll-all-loops -finline-functions execution test FAIL: gfortran.dg/bessel_7.f90 -O3 -g execution test FAIL: gfortran.dg/bessel_7.f90 -Os execution test
H.J. Lu wrote: > I still got > http://gcc.gnu.org/ml/gcc-testresults/2010-08/msg03035.html: > > FAIL: gfortran.dg/bessel_7.f90 -O0 execution test As written before: Can you uncomment the print lines (all!) in bessel_7 - and comment the if/abort lines? That gives a long output, but one can use it to find the values for which it fails. Thanks! Tobias
On Thu, Sep 2, 2010 at 12:53 PM, Tobias Burnus <burnus@net-b.de> wrote: > H.J. Lu wrote: >> >> I still got >> http://gcc.gnu.org/ml/gcc-testresults/2010-08/msg03035.html: >> >> FAIL: gfortran.dg/bessel_7.f90 -O0 execution test > > As written before: Can you uncomment the print lines (all!) in bessel_7 - > and comment the if/abort lines? That gives a long output, but one can use it > to find the values for which it fails. > > Thanks! > > Tobias > I got YN for X = 0.0000000 -- Epsilon = 1.19209290E-07 YN for X = 0.50000000 -- Epsilon = 1.19209290E-07 YN for X = 1.0000000 -- Epsilon = 1.19209290E-07 YN for X = 0.89999998 -- Epsilon = 1.19209290E-07 YN for X = 1.8000000 -- Epsilon = 1.19209290E-07 33 -0.277982435E+37 -0.277982181E+37 -0.25E+31 7.6507163048 T F 34 -0.101848647E+39 -0.101848556E+39 -0.91E+32 7.5173835754 T F YN for X = 2.0000000 -- Epsilon = 1.19209290E-07 YN for X = 3.0000000 -- Epsilon = 1.19209290E-07 YN for X = 4.0000000 -- Epsilon = 1.19209290E-07 41 -0.130542104E+36 -0.130542174E+36 0.69E+29 -4.4547867775 T F 42 -0.266956917E+37 -0.266957075E+37 0.16E+31 -4.9791855812 T F 43 -0.559304076E+38 -0.559304431E+38 0.35E+32 -5.3235278130 T F YN for X = 4.2500000 -- Epsilon = 1.19209290E-07 YN for X = 8.0000000 -- Epsilon = 1.19209290E-07 YN for X = 34.529999 -- Epsilon = 1.19209290E-07 YN for X = 475.78000 -- Epsilon = 1.19209290E-07
Index: testsuite/ChangeLog =================================================================== --- testsuite/ChangeLog (Revision 163613) +++ testsuite/ChangeLog (Arbeitskopie) @@ -1,3 +1,8 @@ +2010-07-21 Steven G. Kargl <kargl@gcc.gnu.org> + Tobias Burnus <burnus@net-b.de> + + * gfortran.dg/bessel_7.f90: Decrease required precision. + 2010-08-27 Naveen H.S <naveen.S@kpitcummins.com> Kaz Kojima <kkojima@gcc.gnu.org> Index: testsuite/gfortran.dg/bessel_7.f90 =================================================================== --- testsuite/gfortran.dg/bessel_7.f90 (Revision 163613) +++ testsuite/gfortran.dg/bessel_7.f90 (Arbeitskopie) @@ -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, 3, 3, 5, 7, 2, 12, 4, 7, 3, 30, 168 ] + * [2, 3, 4, 5, 7, 2, 12, 4, 7, 6, 30, 168 ] ! 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 ] @@ -19,7 +19,7 @@ real :: rec(0:Nmax), lib(0:Nmax) integer :: i do i = 1, ubound(values,dim=1) - call compare(values(i), myeps(i), nit(i), 3*epsilon(0.0)) + call compare(values(i), myeps(i), nit(i), 6*epsilon(0.0)) end do contains Index: fortran/ChangeLog =================================================================== --- fortran/ChangeLog (Revision 163613) +++ fortran/ChangeLog (Arbeitskopie) @@ -1,3 +1,8 @@ +2010-07-21 Steven G. Kargl <kargl@gcc.gnu.org> + + * simplify.c (gfc_simplify_bessel_n2): Fix indention + and argument type. + 2010-08-28 Francois-Xavier Coudert <fxcoudert@gcc.gnu.org> PR fortran/45436 Index: fortran/simplify.c =================================================================== --- fortran/simplify.c (Revision 163613) +++ fortran/simplify.c (Arbeitskopie) @@ -1241,7 +1241,7 @@ gfc_simplify_bessel_n2 (gfc_expr *order1 if (jn && n1 == 0) { e = gfc_get_constant_expr (x->ts.type, x->ts.kind, &x->where); - mpfr_set_ui (e->value.real, 1.0, GFC_RND_MODE); + mpfr_set_ui (e->value.real, 1, GFC_RND_MODE); gfc_constructor_append_expr (&result->value.constructor, e, &x->where); n1++; @@ -1251,7 +1251,7 @@ gfc_simplify_bessel_n2 (gfc_expr *order1 { e = gfc_get_constant_expr (x->ts.type, x->ts.kind, &x->where); if (jn) - mpfr_set_ui (e->value.real, 0.0, GFC_RND_MODE); + mpfr_set_ui (e->value.real, 0, GFC_RND_MODE); else mpfr_set_inf (e->value.real, -1); gfc_constructor_append_expr (&result->value.constructor, e, @@ -1315,7 +1315,7 @@ gfc_simplify_bessel_n2 (gfc_expr *order1 return &gfc_bad_expr; } if (jn) - gfc_constructor_insert_expr (&result->value.constructor, e, &x->where, -2); + gfc_constructor_insert_expr (&result->value.constructor, e, &x->where, -2); else gfc_constructor_append_expr (&result->value.constructor, e, &x->where);