[Bug fortran/36158] Transformational function BESSEL_YN(n1,n2,x) and BESSEL_JN missing

[dominiq at lps dot ens dot fr]

------- Comment #12 from dominiq at lps dot ens dot fr  2010-08-21 13:11 -------
On x86_64-apple-darwin10.4 I need the following adjustments:

--- /opt/gcc/work/gcc/testsuite/gfortran.dg/bessel_6.f90        2010-08-21
12:30:29.000000000 +0200
+++ bessel_6_db.f90     2010-08-21 15:05:20.000000000 +0200
@@ -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, 7, 5, 6, 9, 12, 12, 7, 7, 8, 75, 6 ]
+                  * [2, 7, 5, 6, 9, 12, 12, 7, 7, 8, 120, 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 ]

--- /opt/gcc/work/gcc/testsuite/gfortran.dg/bessel_7.f90        2010-08-21
12:30:30.000000000 +0200
+++ bessel_7_db.f90     2010-08-21 15:02:49.000000000 +0200
@@ -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, 2, 2, 5, 5, 2, 12, 2, 4, 3, 30, 130 ]
+                  * [2, 3, 2, 5, 7, 2, 12, 3, 7, 3, 30, 144 ]
 ! 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), myeps(i))
 end do

 contains


-- 


http://gcc.gnu.org/bugzilla/show_bug.cgi?id=36158