simplify.c (gfc_simplify_bessel_n2): Fix indention and argument type.

fortran/
2010-07-21  Steven G. Kargl  <kargl@gcc.gnu.org>

        * simplify.c (gfc_simplify_bessel_n2): Fix indention
        and argument type.

testsuite/
2010-07-21  Steven G. Kargl  <kargl@gcc.gnu.org>
            Tobias Burnus  <burnus@net-b.de>

        * gfortran.dg/bessel_7.f90: Decrease required precision.


Co-Authored-By: Tobias Burnus <burnus@net-b.de>

From-SVN: r163615
This commit is contained in:
Steven G. Kargl 2010-08-28 14:27:34 +00:00 committed by Tobias Burnus
parent e4ff7a56c6
commit 4c6e913ce8
4 changed files with 15 additions and 5 deletions

View File

@ -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> 2010-08-28 Francois-Xavier Coudert <fxcoudert@gcc.gnu.org>
PR fortran/45436 PR fortran/45436

View File

@ -1241,7 +1241,7 @@ gfc_simplify_bessel_n2 (gfc_expr *order1, gfc_expr *order2, gfc_expr *x,
if (jn && n1 == 0) if (jn && n1 == 0)
{ {
e = gfc_get_constant_expr (x->ts.type, x->ts.kind, &x->where); 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, gfc_constructor_append_expr (&result->value.constructor, e,
&x->where); &x->where);
n1++; n1++;
@ -1251,7 +1251,7 @@ gfc_simplify_bessel_n2 (gfc_expr *order1, gfc_expr *order2, gfc_expr *x,
{ {
e = gfc_get_constant_expr (x->ts.type, x->ts.kind, &x->where); e = gfc_get_constant_expr (x->ts.type, x->ts.kind, &x->where);
if (jn) if (jn)
mpfr_set_ui (e->value.real, 0.0, GFC_RND_MODE); mpfr_set_ui (e->value.real, 0, GFC_RND_MODE);
else else
mpfr_set_inf (e->value.real, -1); mpfr_set_inf (e->value.real, -1);
gfc_constructor_append_expr (&result->value.constructor, e, gfc_constructor_append_expr (&result->value.constructor, e,

View File

@ -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> 2010-08-27 Naveen H.S <naveen.S@kpitcummins.com>
Kaz Kojima <kkojima@gcc.gnu.org> Kaz Kojima <kkojima@gcc.gnu.org>

View File

@ -8,7 +8,7 @@
implicit none 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 :: 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) & 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 ! The following is sufficient for me - the values above are a bit
! more tolerant ! more tolerant
! * [0, 0, 0, 3, 3, 0, 9, 0, 2, 1, 22, 130 ] ! * [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 integer :: i
do i = 1, ubound(values,dim=1) 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 end do
contains contains