obtained gcc-4.6.4.tar.bz2 from upstream website;upstream

verified gcc-4.6.4.tar.bz2.sig;
imported gcc-4.6.4 source tree from verified upstream tarball.

downloading a git-generated archive based on the 'upstream' tag
should provide you with a source tree that is binary identical
to the one extracted from the above tarball.

if you have obtained the source via the command 'git clone',
however, do note that line-endings of files in your working
directory might differ from line-endings of the respective
files in the upstream repository.

1 files changed, 49 insertions, 0 deletions

diff --git a/gcc/testsuite/gfortran.dg/bessel_6.f90 b/gcc/testsuite/gfortran.dg/bessel_6.f90
new file mode 100644
index 000000000..1671d1118
--- /dev/null
+++ b/gcc/testsuite/gfortran.dg/bessel_6.f90
@@ -0,0 +1,49 @@
+! { 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, 92, 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, 36, 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) .or. abs((rec(i)-lib(i))/rec(i)) < myeps
+if (.not. (rec(i) == lib(i) .or. abs((rec(i)-lib(i))/rec(i)) < myeps)) &
+  call abort()
+end do
+
+end
+end