public inbox for glibc-bugs@sourceware.org
help / color / mirror / Atom feed
* [Bug math/13203] New: gfortran.dg/bessel_6.f90 in gcc testsuite failed
@ 2011-09-20 4:34 hjl.tools at gmail dot com
2011-09-20 4:34 ` [Bug math/13203] " hjl.tools at gmail dot com
` (3 more replies)
0 siblings, 4 replies; 5+ messages in thread
From: hjl.tools at gmail dot com @ 2011-09-20 4:34 UTC (permalink / raw)
To: glibc-bugs
http://sourceware.org/bugzilla/show_bug.cgi?id=13203
Bug #: 13203
Summary: gfortran.dg/bessel_6.f90 in gcc testsuite failed
Product: glibc
Version: 2.14
Status: NEW
Severity: normal
Priority: P2
Component: math
AssignedTo: aj@suse.de
ReportedBy: hjl.tools@gmail.com
Classification: Unclassified
Glibc master on 2011-09-19 failed gfortran.dg/bessel_6.f90
on Linux/x86-64:
gnu-mic-2:pts/1[83]> cat x.f90 /tmp
! { dg-do run { xfail spu-*-* } }
! { dg-add-options ieee }
!
! PR fortran/36158
! PR fortran/33197
!
! XFAILed for SPU targets since we don't have an accurate library
! implementation of the single-precision Bessel functions.
!
! Run-time tests for transformations BESSEL_JN
!
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, 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 ]
integer,parameter :: mymax(size(values)) = &
[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
do i = 11, ubound(values,dim=1)
call compare(mymax(i), values(i), myeps(i))
end do
contains
subroutine compare(mymax, X, myeps)
integer :: i, nit, mymax
real X, myeps, myeps2
rec(0:mymax) = BESSEL_JN(0, mymax, X)
lib(0:mymax) = [ (BESSEL_JN(i, X), i=0,mymax) ]
print *, 'YN for X = ', X, ' -- Epsilon = ',epsilon(x)
do i = 24, 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), abs((rec(i)-lib(i))/rec(i)) < myeps
if (rec(i) == lib(i)) CYCLE
if (abs((rec(i)-lib(i))/rec(i)) > myeps) then
call abort()
endif
end do
end
end
gnu-mic-2:pts/1[84]> bash /tmp
[hjl@gnu-mic-2 tmp]$ cat bessel_6.f90
! { dg-do run { xfail spu-*-* } }
! { dg-add-options ieee }
!
! PR fortran/36158
! PR fortran/33197
!
! XFAILed for SPU targets since we don't have an accurate library
! implementation of the single-precision Bessel functions.
!
! Run-time tests for transformations BESSEL_JN
!
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, 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 ]
integer,parameter :: mymax(size(values)) = &
[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
do i = 1, ubound(values,dim=1)
call compare(mymax(i), values(i), myeps(i))
end do
contains
subroutine compare(mymax, X, myeps)
integer :: i, nit, mymax
real X, myeps, myeps2
rec(0:mymax) = BESSEL_JN(0, mymax, X)
lib(0:mymax) = [ (BESSEL_JN(i, X), i=0,mymax) ]
!print *, 'YN for X = ', X, ' -- Epsilon = ',epsilon(x)
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), 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
end
end
[hjl@gnu-mic-2 tmp]$ gfortran bessel_6.f90
[hjl@gnu-mic-2 tmp]$ ./a.out
Aborted (core dumped)
[hjl@gnu-mic-2 tmp]$
--
Configure bugmail: http://sourceware.org/bugzilla/userprefs.cgi?tab=email
------- You are receiving this mail because: -------
You are on the CC list for the bug.
^ permalink raw reply [flat|nested] 5+ messages in thread
* [Bug math/13203] gfortran.dg/bessel_6.f90 in gcc testsuite failed
2011-09-20 4:34 [Bug math/13203] New: gfortran.dg/bessel_6.f90 in gcc testsuite failed hjl.tools at gmail dot com
@ 2011-09-20 4:34 ` hjl.tools at gmail dot com
2011-12-22 16:59 ` drepper.fsp at gmail dot com
` (2 subsequent siblings)
3 siblings, 0 replies; 5+ messages in thread
From: hjl.tools at gmail dot com @ 2011-09-20 4:34 UTC (permalink / raw)
To: glibc-bugs
http://sourceware.org/bugzilla/show_bug.cgi?id=13203
H.J. Lu <hjl.tools at gmail dot com> changed:
What |Removed |Added
----------------------------------------------------------------------------
Version|2.14 |unspecified
--
Configure bugmail: http://sourceware.org/bugzilla/userprefs.cgi?tab=email
------- You are receiving this mail because: -------
You are on the CC list for the bug.
^ permalink raw reply [flat|nested] 5+ messages in thread
* [Bug math/13203] gfortran.dg/bessel_6.f90 in gcc testsuite failed
2011-09-20 4:34 [Bug math/13203] New: gfortran.dg/bessel_6.f90 in gcc testsuite failed hjl.tools at gmail dot com
2011-09-20 4:34 ` [Bug math/13203] " hjl.tools at gmail dot com
@ 2011-12-22 16:59 ` drepper.fsp at gmail dot com
2012-02-25 13:30 ` jsm28 at gcc dot gnu.org
2014-06-27 12:02 ` fweimer at redhat dot com
3 siblings, 0 replies; 5+ messages in thread
From: drepper.fsp at gmail dot com @ 2011-12-22 16:59 UTC (permalink / raw)
To: glibc-bugs
http://sourceware.org/bugzilla/show_bug.cgi?id=13203
Ulrich Drepper <drepper.fsp at gmail dot com> changed:
What |Removed |Added
----------------------------------------------------------------------------
Status|NEW |WAITING
CC| |drepper.fsp at gmail dot
| |com
--- Comment #1 from Ulrich Drepper <drepper.fsp at gmail dot com> 2011-12-22 16:59:09 UTC ---
I'm waiting for a test case in C and then a patch.
--
Configure bugmail: http://sourceware.org/bugzilla/userprefs.cgi?tab=email
------- You are receiving this mail because: -------
You are on the CC list for the bug.
^ permalink raw reply [flat|nested] 5+ messages in thread
* [Bug math/13203] gfortran.dg/bessel_6.f90 in gcc testsuite failed
2011-09-20 4:34 [Bug math/13203] New: gfortran.dg/bessel_6.f90 in gcc testsuite failed hjl.tools at gmail dot com
2011-09-20 4:34 ` [Bug math/13203] " hjl.tools at gmail dot com
2011-12-22 16:59 ` drepper.fsp at gmail dot com
@ 2012-02-25 13:30 ` jsm28 at gcc dot gnu.org
2014-06-27 12:02 ` fweimer at redhat dot com
3 siblings, 0 replies; 5+ messages in thread
From: jsm28 at gcc dot gnu.org @ 2012-02-25 13:30 UTC (permalink / raw)
To: glibc-bugs
http://sourceware.org/bugzilla/show_bug.cgi?id=13203
Joseph Myers <jsm28 at gcc dot gnu.org> changed:
What |Removed |Added
----------------------------------------------------------------------------
Status|WAITING |RESOLVED
Resolution| |INVALID
--- Comment #2 from Joseph Myers <jsm28 at gcc dot gnu.org> 2012-02-25 13:29:02 UTC ---
The gfortran test was fixed.
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=50487
--
Configure bugmail: http://sourceware.org/bugzilla/userprefs.cgi?tab=email
------- You are receiving this mail because: -------
You are on the CC list for the bug.
^ permalink raw reply [flat|nested] 5+ messages in thread
* [Bug math/13203] gfortran.dg/bessel_6.f90 in gcc testsuite failed
2011-09-20 4:34 [Bug math/13203] New: gfortran.dg/bessel_6.f90 in gcc testsuite failed hjl.tools at gmail dot com
` (2 preceding siblings ...)
2012-02-25 13:30 ` jsm28 at gcc dot gnu.org
@ 2014-06-27 12:02 ` fweimer at redhat dot com
3 siblings, 0 replies; 5+ messages in thread
From: fweimer at redhat dot com @ 2014-06-27 12:02 UTC (permalink / raw)
To: glibc-bugs
https://sourceware.org/bugzilla/show_bug.cgi?id=13203
Florian Weimer <fweimer at redhat dot com> changed:
What |Removed |Added
----------------------------------------------------------------------------
Flags| |security-
--
You are receiving this mail because:
You are on the CC list for the bug.
^ permalink raw reply [flat|nested] 5+ messages in thread
end of thread, other threads:[~2014-06-27 12:02 UTC | newest]
Thread overview: 5+ messages (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2011-09-20 4:34 [Bug math/13203] New: gfortran.dg/bessel_6.f90 in gcc testsuite failed hjl.tools at gmail dot com
2011-09-20 4:34 ` [Bug math/13203] " hjl.tools at gmail dot com
2011-12-22 16:59 ` drepper.fsp at gmail dot com
2012-02-25 13:30 ` jsm28 at gcc dot gnu.org
2014-06-27 12:02 ` fweimer at redhat dot com
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).