public inbox for glibc-bugs@sourceware.org help / color / mirror / Atom feed
From: "hjl.tools at gmail dot com" <sourceware-bugzilla@sourceware.org> To: glibc-bugs@sources.redhat.com Subject: [Bug math/13203] New: gfortran.dg/bessel_6.f90 in gcc testsuite failed Date: Tue, 20 Sep 2011 04:34:00 -0000 [thread overview] Message-ID: <bug-13203-131@http.sourceware.org/bugzilla/> (raw) 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.
next reply other threads:[~2011-09-20 4:34 UTC|newest] Thread overview: 5+ messages / expand[flat|nested] mbox.gz Atom feed top 2011-09-20 4:34 hjl.tools at gmail dot com [this message] 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
Reply instructions: You may reply publicly to this message via plain-text email using any one of the following methods: * Save the following mbox file, import it into your mail client, and reply-to-all from there: mbox Avoid top-posting and favor interleaved quoting: https://en.wikipedia.org/wiki/Posting_style#Interleaved_style * Reply using the --to, --cc, and --in-reply-to switches of git-send-email(1): git send-email \ --in-reply-to=bug-13203-131@http.sourceware.org/bugzilla/ \ --to=sourceware-bugzilla@sourceware.org \ --cc=glibc-bugs@sources.redhat.com \ /path/to/YOUR_REPLY https://kernel.org/pub/software/scm/git/docs/git-send-email.html * If your mail client supports setting the In-Reply-To header via mailto: links, try the mailto: linkBe sure your reply has a Subject: header at the top and a blank line before the message body.
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).