rb-gsl 1.15.3.1

data/tests/roots.rb

@@ -0,0 +1,76 @@
+#!/usr/bin/env ruby
+# Ruby/GSL implementation of GSL "roots/test.c"
+require("gsl")
+require("./gsl_test2.rb")
+include GSL::Test
+include Math
+
+EPSREL = 10.0*GSL::DBL_EPSILON
+EPSABS = 10.0*GSL::DBL_EPSILON
+MAX_ITERATIONS = 150
+
+def WITHIN_TOL(a, b, epsrel, epsabs)
+  (a - b).abs < epsrel*GSL::MIN(a.abs, b.abs) + epsabs
+end
+
+def test_f(type, desc, f, lower, upper, correct)
+  s = GSL::Root::FSolver.alloc(type)
+  s.set(f, lower, upper)
+  iter = 0
+  begin
+    iter += 1
+    s.iterate
+    r = s.root
+    a = s.x_lower
+    b = s.x_upper
+    if a > b
+      GSL::Test::test(GSL::FAILURE, "interval is invalid (#{a},#{b})")
+    end
+    if r < a or r > b
+      GSL::Test::test(GSL::FAILURE, "r lies outside interval #{r} (#{a},#{b})")
+    end
+    status = GSL::Root::test_interval(a, b, EPSABS, EPSREL)
+  end while status == GSL::CONTINUE and iter < MAX_ITERATIONS
+  GSL::Test::test(status, "#{s.name}, #{desc} (#{s.root} obs vs #{correct} expected)")
+  if iter == MAX_ITERATIONS
+    GSL::Test::test(GSL::FAILURE, "exceeded maximum number of iterations")
+  end
+  if !WITHIN_TOL(r, correct, EPSREL, EPSABS)
+    GSL::Test::test(GSL::FAILURE, "incorrect precision (#{r} obs vs #{correct} expected)")
+  end
+end
+
+def test_fdf(type, desc, fdf, root, correct)
+  s = GSL::Root::FdfSolver.alloc(type)
+  s.set(fdf, root)
+  iter = 0
+  begin
+    iter += 1
+    prev = s.root
+    s.iterate
+    status = GSL::Root::test_delta(s.root, prev, EPSABS, EPSREL)
+  end while status == GSL::CONTINUE and iter < MAX_ITERATIONS
+  GSL::Test::test(status, "#{s.name} #{desc} (#{s.root} obs vs #{correct} expected)")
+  if iter == MAX_ITERATIONS
+    GSL::Test::test(GSL::FAILURE, "exceeded maximum number of iterations")
+  end
+  if !WITHIN_TOL(s.root, correct, EPSREL, EPSABS)
+    GSL::Test::test(GSL::FAILURE, "incorrect precision (#{r} obs vs #{correct} expected)")
+  end    
+end
+
+func1 = GSL::Function.alloc { |x| pow(x, 20.0) - 1 }
+func1_fdf = Proc.new { |x| 20.0*pow(x, 19.0) }
+
+fdf1 = GSL::Function_fdf.alloc(func1.f, func1_fdf)
+
+fsolvers = ["bisection", "brent", "falsepos"]
+fsolvers.each do |type|
+  test_f(type, "x^20 - 1 [0.1, 2]", func1, 0.1, 2.0, 1.0)
+end
+
+fdfsolvers = ["newton", "secant", "steffenson"]
+fdfsolvers = ["newton", "steffenson"]
+fdfsolvers.each do |type|
+  test_fdf(type, "x^{20} - 1 {0.9}", fdf1, 0.9, 1.0)
+end

data/tests/run-test.sh

@@ -0,0 +1,17 @@
+#!/bin/sh
+
+topdir=".."
+for f in `find . -name '*.rb'`; do
+	if test "$f" = "./sf/make_test.rb"; then continue; fi
+	if test "$f" = "./speed/bench.rb"; then continue; fi
+	if test "$f" = "./tensor.rb"; then continue; fi
+	pushd `dirname $f` > /dev/null
+#	ruby -w -I $topdir/lib -I $topdir/ext `basename $f`
+	ruby -w `basename $f`
+	result=$?
+	popd > /dev/null
+	if test "$result" != "0"; then
+		echo $f failed
+		exit $result
+	fi
+done

data/tests/sf/gsl_test_sf.rb

@@ -0,0 +1,249 @@
+#!/usr/bin/env ruby
+require("gsl")
+require("../gsl_test2.rb")
+include GSL::Test
+include Math
+
+module GSL
+  module Test
+    module Sf
+      TEST_TOL0 = 2.0*GSL::DBL_EPSILON
+      TEST_TOL1 = 16.0*GSL::DBL_EPSILON
+      TEST_TOL2 = 256.0*GSL::DBL_EPSILON
+      TEST_TOL3 = 2048.0*GSL::DBL_EPSILON
+      TEST_TOL4 = 16384.0*GSL::DBL_EPSILON
+      TEST_TOL5 = 131072.0*GSL::DBL_EPSILON
+      TEST_TOL6 = 1048576.0*GSL::DBL_EPSILON
+      TEST_SQRT_TOL0 = 2.0*GSL::SQRT_DBL_EPSILON
+      TEST_SNGL = 1.0e-06
+      
+      TEST_SF_INCONS = 1
+      TEST_SF_ERRNEG = 2
+      TEST_SF_TOLBAD = 4
+      TEST_SF_RETBAD = 8
+      
+      TEST_FACTOR = 100.0
+      TEST_SIGMA = 1.5
+      
+      def test_sf_frac_diff(x1, x2)
+        if x1.zero? and x2.zero?
+          return 0.0
+        elsif x1 <= DBL_MAX and x2 < DBL_MAX and (x1 + x1 != 0.0)
+          return ((x1 - x2)/(x1 + x2)).abs
+        else
+          return 1.0
+        end
+      end
+      
+      def test_sf_check_result(mbuf, r, val, tol)
+        s = 0
+        f = 0.0
+        if isnan?(r.val) or isnan?(val)
+          s = isnan(r.val) != isnan(val) ? TEST_SF_INCONS : s
+        elsif isinf?(r.val) or isinf?(val)
+          s = isinf(r.val) != isinf(val) ? TEST_SF_INCONS : s
+        else
+          f = test_sf_frac_diff(val, r.val)
+          if (val - r.val).abs > 2.0*TEST_SIGMA*r.err
+            s |= TEST_SF_INCONS
+          end
+          if r.err < 0.0
+            s |= TEST_SF_ERRNEG
+          end
+          if f > TEST_FACTOR * tol
+            s |= TEST_SF_TOLBAD
+          end
+        end
+        if s != 0
+          message = sprintf("  expected: %20.16g\n", val)
+          mbuf += message
+          message = sprintf("  obtained: %20.16g   %20.16g  %g\n", 
+                            r.val, r.err, r.err/(r.val.abs + r.err))
+          mbuf += message
+          message = sprintf("  fracdiff: %20.16g\n", f)
+          mbuf += message
+        end
+        
+        if s & TEST_SF_INCONS
+          mbuf += "  value/expected not consistent within reported error\n"
+        end
+        if s &  TEST_SF_ERRNEG
+          mbuf += "  reported error negative\n"
+        end
+        if s & TEST_SF_TOLBAD
+          mbuf += "  value not within tolerance of expected value\n"
+        end
+        return s, mbuf
+      end
+      
+      def test_sf_check_val(mbuf, rval, val, tol)
+        s = 0;
+        f = test_sf_frac_diff(val, rval)
+        if f > TEST_FACTOR * tol
+          s |= TEST_SF_TOLBAD
+        end
+        if s != 0
+          buf = sprintf("  expected: %20.16g\n", val)
+          mbuf += buf
+          buf = sprintf("  obtained: %20.16g\n", rval)
+          mbuf += buf
+          buf = sprintf("  fracdiff: %20.16g\n", f)
+          mbuf += buf
+        end
+        if s & TEST_SF_TOLBAD
+          mbuf +=  "  value not within tolerance of expected value\n"
+        end
+        return s, mbuf
+      end
+      
+      def test_sf_check_result_relax(mbuf, r, val, tol)
+        s = 0;
+        f = test_sf_frac_diff(val, r.val)
+        
+        if f > GSL_MAX_DBL(TEST_SNGL, TEST_FACTOR * tol)
+          s |= TEST_SF_INCONS
+        end
+        if r.err < 0.0
+          s |= TEST_SF_ERRNEG
+        end
+        if f > TEST_FACTOR * tol    
+          s |= TEST_SF_TOLBAD
+        end
+        if s != 0
+          buf = sprintf("  expected: %20.16g\n", val)
+          mbuf += buf
+          buf = sprintf("  obtained: %20.16g   %20.16g  %g\n", r.val, r.err, r.err/(r.val.abs + r.err))
+          mbuf += buf
+          buf = sprintf("  fracdiff: %20.16g\n", f)
+          mbuf += buf
+        end
+        if s & TEST_SF_INCONS
+          mbuf += "  value/expected not consistent MAX(tol,SINGLE_PREC)\n"
+        end
+        if s & TEST_SF_ERRNEG
+          mbuf += "  reported error negative\n"
+        end
+        if s & TEST_SF_TOLBAD
+          mbuf += "  value not within tolerance of expected value\n"
+        end
+        return s, mbuf
+      end
+      
+      def test_sf_check_return(mbuf, val, expected)
+        if val != expected
+          buf = sprintf("  unexpected return code: %d\n", val)
+          mbuf += buf
+          return TEST_SF_RETBAD, mbuf
+        else
+          return 0, mbuf
+        end
+      end
+      
+      def test_sf(r, val, tol, status, expect, desc)
+        local_s = 0
+        mbuf = ""
+        s, mbuf = test_sf_check_result(mbuf, r, val, tol)
+        local_s |= s
+        s, mbuf = test_sf_check_return(mbuf, status, expect)
+        local_s |= s
+        GSL::Test::test(local_s, desc)
+        if local_s != 0
+          print(mbuf)
+          printf("  %22.18g  %22.18g\n", r.val, r.err)
+        end
+        return local_s
+      end
+      
+      def test_sf_val(val, val_in, tol, desc)
+        local_s = 0
+        mbuf = ""
+        s, mbuf = test_sf_check_val(mbuf, val, val_in, tol)
+        local_s |= s
+        GSL::Test::test(local_s, desc)
+        if local_s != 0
+          printf("%s", mbuf)
+          printf("  %22.18g\n", val)
+        end
+        return local_s
+      end
+      
+      def test_sf_rlx(r, val_in, tol, status, expect_return, desc)
+        local_s = 0
+        message_buff = ""
+        s, message_buff = test_sf_check_result_relax(message_buff, r, val_in, tol)
+        local_s |= s
+        s, message_buff = test_sf_check_return(message_buff, status, expect_return)
+        local_s |= s
+        
+        GSL::Test::test(local_s, desc)
+        if local_s != 0
+          printf("%s", message_buff)
+          printf("  %22.18g  %22.18g\n", r.val, r.err)
+        end
+        return local_s
+      end
+      
+      def test_sf_2(r1, val1, tol1, r2, val2, tol2, status, expect_return, desc)
+        int local_s = 0
+        message_buff = ""
+        s, message_buff = test_sf_check_result(message_buff, r1, val1, tol1)
+        local_s |= s
+        s, message_buff = test_sf_check_result(message_buff, r2, val2, tol2)
+        local_s |= s
+        s, message_buff = test_sf_check_return(message_buff, status, expect_return)
+        local_s |= s
+        
+        GSL::Test::test(local_s, desc)
+        if local_s != 0
+          printf("%s", message_buff)
+          printf("  %22.18g  %22.18g\n", r1.val, r1.err)
+          printf("  %22.18g  %22.18g\n", r2.val, r2.err)
+        end
+        return local_s
+      end
+      
+      def test_sf_sgn(r, sgn, val_in, tol, expect_sgn, status, expect_return, desc)
+        local_r
+        local_s = 0
+        message_buff = ""
+        local_r.val = sgn
+        local_r.err = 0.0
+        s, message_buff = test_sf_check_result(message_buff, r, val_in, tol)
+        local_s |= s
+        s, message_buff = test_sf_check_result(message_buff, local_r, expect_sgn, 0.0)
+        local_s |= s
+        s, message_buff = test_sf_check_return(message_buff, status, expect_return)
+        local_s |= s
+        GSL::Test::test(local_s, desc)
+        if local_s != 0
+          printf("%s", message_buff)
+          printf("  %22.18g  %22.18g\n", r.val, r.err)
+        end
+        return local_s
+      end
+      
+      def TEST_SF(stat, func, args, val_in,  tol, expect_return)
+        r, = eval("#{func}#{args}")
+        #      p r
+        #      p val_in
+        status = 0
+        stat += test_sf(r, val_in, tol, status, expect_return, "#{func}#{args}")
+        return stat
+      end
+      def TEST_SF_2(stat, func, args, val1, tol1, val2, tol2, expect_return) 
+        r, = eval("#{func}#{args}") 
+        status = 0
+        stat += test_sf_2(r1, val1, tol1, r2, val2, tol2, status, 
+                          expect_return, "#{func}#{args}")
+        return stat
+      end
+      def TEST_SF_SGN(stat, func, args, val_in, tol, expect_sgn, expect_return)
+        r, = eval("#{func}#{args}")
+        status = 0
+        stat += test_sf_sgn(r, sgn, val_in, tol, expect_sgn, status, 
+                            expect_return, "#{func}#{args}")
+        return stat
+      end
+    end
+  end
+end

data/tests/sf/test_airy.rb

@@ -0,0 +1,83 @@
+#!/usr/bin/env ruby
+require("gsl")
+require("../gsl_test2.rb")
+require("./gsl_test_sf.rb")
+include GSL::Test
+include Math
+include GSL::Test::Sf
+s = 0
+m = GSL::MODE_DEFAULT
+  TEST_SF(s, "GSL::Sf::airy_Ai_e", "(-500.0, GSL::MODE_DEFAULT)",              0.0725901201040411396, TEST_TOL4, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_e", "(-5.0, GSL::MODE_DEFAULT)",                0.3507610090241142,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_e", "(-0.3000000000000094, GSL::MODE_DEFAULT)", 0.4309030952855831,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_e", "(0.6999999999999907, GSL::MODE_DEFAULT)",  0.1891624003981519,    TEST_TOL0, GSL::SUCCESS)
+    TEST_SF(s, "GSL::Sf::airy_Ai_e", "(1.649999999999991, GSL::MODE_DEFAULT)",   0.05831058618720882,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_e", "(1.649999999999991, GSL::MODE_DEFAULT)",   0.0583105861872088521,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_e", "(2.54999999999999, GSL::MODE_DEFAULT)",    0.01446149513295428,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_e", "(3.499999999999987, GSL::MODE_DEFAULT)",   0.002584098786989702,  TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_e", "(5.39999999999998, GSL::MODE_DEFAULT)",    4.272986169411866e-05, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_scaled_e", "(-5.0, GSL::MODE_DEFAULT)",                  0.3507610090241142, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_scaled_e", "(0.6999999999999907, GSL::MODE_DEFAULT)", 0.2795125667681217, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_scaled_e", "(1.649999999999991, GSL::MODE_DEFAULT)",  0.2395493001442741, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_scaled_e", "(2.54999999999999, GSL::MODE_DEFAULT)",   0.2183658595899388, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_scaled_e", "(3.499999999999987, GSL::MODE_DEFAULT)",  0.2032920808163519, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_scaled_e", "(5.39999999999998, GSL::MODE_DEFAULT)",   0.1836050093282229, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_e", "(-500.0, GSL::MODE_DEFAULT)",             -0.094688570132991028, TEST_TOL4, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_e", "(-5.0, GSL::MODE_DEFAULT)",               -0.1383691349016005,   TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_e", "(0.6999999999999907, GSL::MODE_DEFAULT)",  0.9733286558781599,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_e", "(1.649999999999991, GSL::MODE_DEFAULT)",   2.196407956850028,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_e", "(2.54999999999999, GSL::MODE_DEFAULT)",    6.973628612493443,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_e", "(3.499999999999987, GSL::MODE_DEFAULT)",   33.05550675461069,    TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_e", "(5.39999999999998, GSL::MODE_DEFAULT)",    1604.476078241272,    TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_scaled_e", "(-5.0, GSL::MODE_DEFAULT)",                  -0.1383691349016005, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_scaled_e", "(0.6999999999999907, GSL::MODE_DEFAULT)",  0.6587080754582302, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_scaled_e", "(1.649999999999991, GSL::MODE_DEFAULT)",   0.5346449995597539, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_scaled_e", "(2.54999999999999, GSL::MODE_DEFAULT)",    0.461835455542297,  TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_scaled_e", "(3.499999999999987, GSL::MODE_DEFAULT)",   0.4201771882353061, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_scaled_e", "(5.39999999999998, GSL::MODE_DEFAULT)",    0.3734050675720473, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_deriv_e", "(-5.0, GSL::MODE_DEFAULT)",                 0.3271928185544435,       TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_deriv_e", "(-0.5500000000000094, GSL::MODE_DEFAULT)", -0.1914604987143629,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_deriv_e", "(0.4999999999999906, GSL::MODE_DEFAULT)",  -0.2249105326646850,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_deriv_e", "(1.899999999999992, GSL::MODE_DEFAULT)",   -0.06043678178575718,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_deriv_e", "(3.249999999999988, GSL::MODE_DEFAULT)",   -0.007792687926790889,  TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_deriv_e", "(5.199999999999981, GSL::MODE_DEFAULT)",   -0.0001589434526459543, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_deriv_scaled_e", "(-5.0, GSL::MODE_DEFAULT)",                0.3271928185544435, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_deriv_scaled_e", "(0.5499999999999906, GSL::MODE_DEFAULT)", -0.2874057279170166, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_deriv_scaled_e", "(1.499999999999991, GSL::MODE_DEFAULT)",  -0.3314199796863637, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_deriv_scaled_e", "(2.49999999999999, GSL::MODE_DEFAULT)",   -0.3661089384751620, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_deriv_scaled_e", "(3.649999999999986, GSL::MODE_DEFAULT)",  -0.3974033831453963, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Ai_deriv_scaled_e", "(6.299999999999977, GSL::MODE_DEFAULT)",  -0.4508799189585947, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_deriv_e", "(-5.0, GSL::MODE_DEFAULT)",                0.778411773001899,  TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_deriv_e", "(-0.5500000000000094, GSL::MODE_DEFAULT)", 0.5155785358765014, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_deriv_e", "(0.4999999999999906, GSL::MODE_DEFAULT)",  0.5445725641405883, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_deriv_e", "(1.899999999999992, GSL::MODE_DEFAULT)",   3.495165862891568,  TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_deriv_e", "(3.249999999999988, GSL::MODE_DEFAULT)",   36.55485149250338,  TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_deriv_e", "(5.199999999999981, GSL::MODE_DEFAULT)",   2279.748293583233,  TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_deriv_scaled_e", "(-5.0, GSL::MODE_DEFAULT)",               0.778411773001899,  TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_deriv_scaled_e", "(0.5499999999999906, GSL::MODE_DEFAULT)", 0.4322811281817566, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_deriv_scaled_e", "(1.499999999999991, GSL::MODE_DEFAULT)",  0.5542307563918037, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_deriv_scaled_e", "(2.49999999999999, GSL::MODE_DEFAULT)",   0.6755384441644985, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_deriv_scaled_e", "(3.649999999999986, GSL::MODE_DEFAULT)",  0.7613959373000228, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_Bi_deriv_scaled_e", "(6.299999999999977, GSL::MODE_DEFAULT)",  0.8852064139737571, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Ai_e", "(2)",  -4.087949444130970617, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Ai_e", "(50)", -38.02100867725525443, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Ai_e", "(100)",  -60.45555727411669871, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Ai_e", "(110)",  -64.43135670991324811, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Bi_e", "(2)", -3.271093302836352716, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Bi_e", "(50)", -37.76583438165180116, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Bi_e", "(100)", -60.25336482580837088, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Bi_e", "(110)", -64.2355167606561537,  TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Bi_e", "(111)", -64.6268994819519378,  TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Bi_e", "(200)", -95.88699147356682665, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Ai_deriv_e", "(2)", -3.248197582179836561, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Ai_deriv_e", "(50)", -37.76565910053887108, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Ai_deriv_e", "(100)", -60.25329596442479317, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Ai_deriv_e", "(110)", -64.23545617243546956, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Ai_deriv_e", "(1000)", -280.9378080358935071, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Bi_deriv_e", "(2)", -4.073155089071828216, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Bi_deriv_e", "(50)", -38.02083574095788210, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Bi_deriv_e", "(100)", -60.45548887257140819, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Bi_deriv_e", "(110)", -64.43129648944845060, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Bi_deriv_e", "(111)", -64.82208737584206093, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Bi_deriv_e", "(200)", -96.04731050310324450, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::airy_zero_Bi_deriv_e", "(1000)", -281.0315164471118527, TEST_TOL0, GSL::SUCCESS)

data/tests/sf/test_bessel.rb

@@ -0,0 +1,306 @@
+#!/usr/bin/env ruby
+require("gsl")
+require("../gsl_test2.rb")
+require("./gsl_test_sf.rb")
+include GSL::Test
+include Math
+include GSL::Test::Sf
+s = 0
+m = GSL::MODE_DEFAULT
+  TEST_SF(s, "GSL::Sf::bessel_J0_e", "(0.1)",     0.99750156206604003230,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_J0_e", "(2.0)",     0.22389077914123566805,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_J0_e", "(100.0)",   0.019985850304223122424,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_J0_e", "(1.0e+10)", 2.1755917502468917269e-06, TEST_SQRT_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_J1_e", "(0.1)",      0.04993752603624199756,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_J1_e", "(2.0)",      0.57672480775687338720,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_J1_e", "(100.0)",   -0.07714535201411215803,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_J1_e", "(1.0e+10)", -7.676508175684157103e-06, TEST_TOL4, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jn_e", "(4, 0.1)",     2.6028648545684032338e-07, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jn_e", "(5, 2.0)",     0.007039629755871685484,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jn_e", "(10, 20.0)",   0.18648255802394508321,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jn_e", "(100, 100.0)", 0.09636667329586155967,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jn_e", "(2, 900.0)", -0.019974345269680646400, TEST_TOL4, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jn_e", "(2, 15000.0)", -0.0020455820181216382666, TEST_TOL4, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jn_e", "(0, 1.0e+10)", 2.1755917502468917269e-06, TEST_SQRT_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jn_e", "(1, 1.0e+10)", -7.676508175684157103e-06, TEST_TOL4, GSL::SUCCESS)
+#  TEST_SF(s, "GSL::Sf::bessel_Jn_e", "(0, 20000)", -7.676508175684157103e-06, TEST_TOL4, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Y0_e", "(0.1)",         -1.5342386513503668441,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Y0_e", "(2)",            0.5103756726497451196,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Y0_e", "(256.0)",       -0.03381290171792454909 ,  TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Y0_e", "(4294967296.0)", 3.657903190017678681e-06, TEST_SQRT_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Y1_e", "(0.1)",         -6.45895109470202698800,     TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Y1_e", "(2)",           -0.10703243154093754689,     TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Y1_e", "(100.0)",       -0.020372312002759793305,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Y1_e", "(4294967296.0)", 0.000011612249378370766284, TEST_TOL4, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Yn_e", "(4, 0.1)",            -305832.29793353160319,    TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Yn_e", "(5, 2)",              -9.935989128481974981,     TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Yn_e", "(100, 100.0)",        -0.16692141141757650654,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Yn_e", "(100, 4294967296.0)",  3.657889671577715808e-06, TEST_SQRT_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Yn_e", "(1000, 4294967296.0)", 3.656551321485397501e-06, 2.0e-05, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Yn_e", "(2, 15000.0)", -0.006185217273358617849, TEST_TOL4, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_I0_scaled_e", "(1e-10)",   0.99999999990000000001,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_I0_scaled_e", "(0.1)",     0.90710092578230109640,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_I0_scaled_e", "(2)",       0.30850832255367103953,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_I0_scaled_e", "(100.0)",   0.03994437929909668265,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_I0_scaled_e", "(65536.0)", 0.0015583712551952223537, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_I1_scaled_e", "(0.1)",     0.04529844680880932501,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_I1_scaled_e", "(2)",       0.21526928924893765916,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_I1_scaled_e", "(100.0)",   0.03974415302513025267,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_I1_scaled_e", "(65536.0)", 0.0015583593657207350452, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_In_scaled_e", "(  -4,    0.1)", 2.3575258620054605307e-07, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_In_scaled_e", "(   4,    0.1)", 2.3575258620054605307e-07, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_In_scaled_e", "(   5,    2.0)", 0.0013297610941881578142, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_In_scaled_e", "( 100,  100.0)", 1.7266862628167695785e-22, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_I0_e", "(0.1)", 1.0025015629340956014, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_I0_e", "(2.0)", 2.2795853023360672674, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_I0_e", "(100.0)", 1.0737517071310738235e+42, TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_I1_e", "(0.1)", 0.05006252604709269211,      TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_I1_e", "(2.0)", 1.59063685463732906340,      TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_I1_e", "(100.0)", 1.0683693903381624812e+42, TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_In_e", "(   4,    0.1)", 2.6054690212996573677e-07, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_In_e", "(   5,    2.0)", 0.009825679323131702321,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_In_e", "( 100,  100.0)", 4.641534941616199114e+21,  TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_K0_scaled_e", "(0.1)", 2.6823261022628943831, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_K0_scaled_e", "(2.0)", 0.8415682150707714179, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_K0_scaled_e", "(100.0)", 0.1251756216591265789, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_K1_scaled_e", "(0.1)", 10.890182683049696574, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_K1_scaled_e", "(2.0)", 1.0334768470686885732, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_K1_scaled_e", "(100.0)", 0.1257999504795785293, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Kn_scaled_e", "(   4,    0.1)", 530040.2483725626207, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Kn_scaled_e", "(   5,    2.0)", 69.68655087607675118, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Kn_scaled_e", "( 100,  100.0)", 2.0475736731166756813e+19, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_K0_e", "(0.1)", 2.4270690247020166125, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_K0_e", "(2.0)", 0.11389387274953343565, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_K0_e", "(100.0)", 4.656628229175902019e-45, TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_K1_e", "(0.1)", 9.853844780870606135,       TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_K1_e", "(2.0)", 0.13986588181652242728,     TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_K1_e", "(100.0)", 4.679853735636909287e-45, TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Kn_e", "(   4,    0.1)", 479600.2497925682849,     TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Kn_e", "(   5,    2.0)", 9.431049100596467443,     TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Kn_e", "( 100,  100.0)", 7.617129630494085416e-25, TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j0_e", "(-10.0)", -0.05440211108893698134, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j0_e", "(0.001)", 0.9999998333333416667, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j0_e", "(  1.0)", 0.84147098480789650670, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j0_e", "( 10.0)", -0.05440211108893698134, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j0_e", "(100.0)", -0.005063656411097587937, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j0_e", "(1048576.0)", 3.1518281938718287624e-07, TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j1_e", "(-10.0)", -0.07846694179875154709, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j1_e", "(0.01)", 0.003333300000119047399, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j1_e", "(  1.0)", 0.30116867893975678925, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j1_e", "( 10.0)", 0.07846694179875154709, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j1_e", "(100.0)", -0.008673825286987815220, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j1_e", "(1048576.0)", -9.000855242905546158e-07, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j2_e", "(-10.0)", 0.07794219362856244547, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j2_e", "(0.01)", 6.666619047751322551e-06, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j2_e", "(  1.0)", 0.06203505201137386110, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j2_e", "( 10.0)", 0.07794219362856244547, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j2_e", "(100.0)", 0.004803441652487953480, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_j2_e", "(1048576.0)", -3.1518539455252413111e-07, TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_jl_e", "(0, 0.0)", 1.0, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_jl_e", "(1,       10.0)",   0.07846694179875154709000, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_jl_e", "(5,        1.0)",   0.00009256115861125816357, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_jl_e", "(10,      10.0)",   0.06460515449256426427,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_jl_e", "(100,    100.0)",   0.010880477011438336539,   TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_jl_e", "(2000, 1048576.0)", 7.449384239168568534e-07,  TEST_SQRT_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_jl_e", "(2, 900.0)",   -0.0011089115568832940086,  TEST_TOL4, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_jl_e", "(2, 15000.0)", -0.00005955592033075750554, TEST_TOL4, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y0_e", "(0.001)", -999.99950000004166670, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y0_e", "(  1.0)", -0.5403023058681397174, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y0_e", "( 10.0)", 0.08390715290764524523, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y0_e", "(100.0)", -0.008623188722876839341, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y0_e", "(65536.0)", 0.000011014324202158573930, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y0_e", "(4294967296.0)", 2.0649445131370357007e-10, TEST_SQRT_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y1_e", "( 0.01)", -10000.499987500069444, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y1_e", "(  1.0)", -1.3817732906760362241, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y1_e", "( 10.0)", 0.06279282637970150586, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y1_e", "(100.0)", 0.004977424523868819543, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y1_e", "(4294967296.0)", 1.0756463271573404688e-10, TEST_SQRT_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y2_e", "( 0.01)", -3.0000500012499791668e+06, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y2_e", "(  1.0)", -3.605017566159968955, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y2_e", "( 10.0)", -0.06506930499373479347, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y2_e", "(100.0)", 0.008772511458592903927, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_y2_e", "(4294967296.0)", -2.0649445123857054207e-10, TEST_SQRT_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_yl_e", "(0,        0.01)", -99.995000041666528,    TEST_TOL0, GSL::SUCCESS) 
+  TEST_SF(s,  "GSL::Sf::bessel_yl_e", "(0,        1.0)",  -0.54030230586813972,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_yl_e", "(1,       10.0)",   0.062792826379701506,   TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_yl_e", "(5,        1.0)",  -999.44034339223641,     TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_yl_e", "(10,       0.01)", -6.5473079797378378e+30, TEST_TOL0, GSL::SUCCESS) 
+  TEST_SF(s,  "GSL::Sf::bessel_yl_e", "(10,      10.0)",  -0.172453672088057849,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_yl_e", "(100,      1.0)",  -6.6830794632586775e+186, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_yl_e", "(100,    100.0)",  -0.0229838504915622811,   TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s,  "GSL::Sf::bessel_yl_e", "(2000, 1048576.0)", 5.9545201447146155e-07,  TEST_SQRT_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_i0_scaled_e", "(0.1)", 0.9063462346100907067, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_i0_scaled_e", "(2.0)", 0.24542109027781645493, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_i0_scaled_e", "(100.0)", 0.005000000000000000000, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_i1_scaled_e", "(0.1)", 0.030191419289002226846, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_i1_scaled_e", "(2.0)", 0.131868364583275317610, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_i1_scaled_e", "(100.0)", 0.004950000000000000000, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_i2_scaled_e", "(0.1)", 0.0006036559400239012567, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_i2_scaled_e", "(2.0)", 0.0476185434029034785100, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_i2_scaled_e", "(100.0)", 0.0048515000000000000000, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_il_scaled_e", "(   4, 0.001)", 1.0571434341190365013e-15, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_il_scaled_e", "(   4,   0.1)", 9.579352242057134927e-08,  TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_il_scaled_e", "(   5,   2.0)", 0.0004851564602127540059,  TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_il_scaled_e", "(   5, 100.0)", 0.004300446777500000000,   TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_il_scaled_e", "( 100, 100.0)", 1.3898161964299132789e-23, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_k0_scaled_e", "(0.1)", 15.707963267948966192, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_k0_scaled_e", "(2.0)", 0.7853981633974483096, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_k0_scaled_e", "(100.0)", 0.015707963267948966192, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_k1_scaled_e", "(0.1)", 172.78759594743862812, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_k1_scaled_e", "(2.0)", 1.1780972450961724644, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_k1_scaled_e", "(100.0)", 0.015865042900628455854, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_k2_scaled_e", "(0.1)", 5199.335841691107810, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_k2_scaled_e", "(2.0)", 2.5525440310417070063, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_k2_scaled_e", "(100.0)", 0.016183914554967819868, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_kl_scaled_e", "(   4, 1.0/256.0)", 1.8205599816961954439e+14, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_kl_scaled_e", "(   4, 1.0/8.0)",   6.1173217814406597530e+06, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_kl_scaled_e", "(   5,   2.0)",     138.10735829492005119,     TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_kl_scaled_e", "( 100, 100.0)",     3.985930768060258219e+18,  TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jnu_e", "(0.0001, 1.0)",         0.7652115411876708497,  TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jnu_e", "(0.0001, 10.0)",       -0.2459270166445205,     TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jnu_e", "(0.0009765625, 10.0)", -0.2458500798634692,     TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jnu_e", "(0.75, 1.0)",           0.5586524932048917478,  TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jnu_e", "(0.75, 10.0)",         -0.04968928974751508135, TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jnu_e", "( 1.0, 0.001)", 0.0004999999375000026,     TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jnu_e", "( 1.0,   1.0)", 0.4400505857449335160,     TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jnu_e", "( 1.75,  1.0)", 0.168593922545763170103,     TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jnu_e", "(30.0,   1.0)", 3.482869794251482902e-42,  TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jnu_e", "(30.0, 100.0)", 0.08146012958117222297,    TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jnu_e", "(10.0,   1.0)", 2.6306151236874532070e-10, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jnu_e", "(10.0, 100.0)", -0.05473217693547201474,   TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jnu_e", "(10.2, 100.0)", -0.03548919161046526864,   TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jnu_e", "(2.0, 900.0)",   -0.019974345269680646400,  TEST_TOL4, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Jnu_e", "(2.0, 15000.0)", -0.0020455820181216382666, TEST_TOL4, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Ynu_e", "(0.0001, 1.0)",  0.08813676933044478439,    TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Ynu_e", "(0.0001,10.0)",  0.05570979797521875261,    TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Ynu_e", "( 0.75,  1.0)", -0.6218694174429746383,     TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Ynu_e", "( 0.75, 10.0)",  0.24757063446760384953,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Ynu_e", "( 1.0, 0.001)", -636.6221672311394281,      TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Ynu_e", "( 1.0,   1.0)", -0.7812128213002887165,     TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Ynu_e", "(30.0,   1.0)", -3.0481287832256432162e+39, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Ynu_e", "(30.0, 100.0)",  0.006138839212010033452,   TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Ynu_e", "(10.0,   1.0)", -1.2161801427868918929e+08, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Ynu_e", "(10.0, 100.0)",  0.05833157423641492875,    TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Ynu_e", "(10.2, 100.0)",  0.07169383985546287091,    TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_scaled_e", "(0.0001,10.0)", 0.12783333709581669672,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_scaled_e", "( 1.0, 0.001)", 0.0004995003123542213370,  TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_scaled_e", "( 1.0,   1.0)", 0.20791041534970844887,    TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_scaled_e", "(30.0,   1.0)", 1.3021094983785914437e-42, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_scaled_e", "(30.0, 100.0)", 0.0004486987756920986146,  TEST_TOL3, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_scaled_e", "(10.0,   1.0)", 1.0127529864692066036e-10, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_scaled_e", "(10.0, 100.0)", 0.024176682718258828365,   TEST_TOL3, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_scaled_e", "(10.2, 100.0)", 0.023691628843913810043,   TEST_TOL3, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_e", "(0.0001,10.0)", 2815.7166269770030352,     TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_e", "( 1.0, 0.001)", 0.0005000000625000026042,  TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_e", "( 1.0,   1.0)", 0.5651591039924850272,     TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_e", "(30.0,   1.0)", 3.539500588106447747e-42,  TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_e", "(30.0, 100.0)", 1.2061548704498434006e+40, TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_e", "(10.0,   1.0)", 2.7529480398368736252e-10, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_e", "(10.0, 100.0)", 6.498975524720147799e+41,  TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Inu_e", "(10.2, 100.0)", 6.368587361287030443e+41,  TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_scaled_e", "(0.0001,10.0)", 0.3916319346235421817, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_scaled_e", "( 1.0, 0.001)", 1000.9967345590684524, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_scaled_e", "( 1.0,   1.0)", 1.6361534862632582465, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_scaled_e", "(30.0,   1.0)", 1.2792629867539753925e+40, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_scaled_e", "(30.0, 100.0)", 10.673443449954850040, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_scaled_e", "(10.0,   1.0)", 4.912296520990198599e+08, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_scaled_e", "(10.0, 100.0)", 0.20578687173955779807, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_scaled_e", "(10.0, 1000.0)", 0.04165905142800565788, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_scaled_e", "(10.0, 1.0e+8)", 0.00012533147624060789938, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_scaled_e", "(10.2, 100.0)", 0.20995808355244385075, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_e", "(0.0001,0.001)", 7.023689431812884141,      TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_e", "(0.0001,10.0)", 0.000017780062324654874306, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_e", "( 1.0, 0.001)", 999.9962381560855743,      TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_e", "( 1.0,   1.0)", 0.6019072301972345747,     TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_e", "(10.0, 0.001)", 1.8579455483904008064e+38, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_e", "(10.0,   1.0)", 1.8071328990102945469e+08, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_e", "(10.0, 100.0)", 7.655427977388100611e-45,  TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_e", "(10.2, 100.0)", 7.810600225948217841e-45,  TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_e", "(30.0,   1.0)", 4.706145526783626883e+39,  TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_Knu_e", "(30.0, 100.0)", 3.970602055959398739e-43,  TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_lnKnu_e", "(0.0001,1.0e-100)", 5.439794449319847, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_lnKnu_e", "(0.0001,0.0001)", 2.232835507214331, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_lnKnu_e", "(0.0001,10.0)", -10.93743282256098, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_lnKnu_e", "( 1.0, 1.0e-100)", 230.2585092994045, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_lnKnu_e", "( 1.0, 1.0e-10)", 23.025850929940456840, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_lnKnu_e", "( 1.0, 0.001)", 6.907751517131146, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_lnKnu_e", "( 1.0,   1.0)", -0.5076519482107523309, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_lnKnu_e", "(30.0, 1.0e-100)", 6999.113586185543475, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_lnKnu_e", "(30.0,   1.0)", 91.34968784026325464, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_lnKnu_e", "(30.0, 100.0)", -97.63224126416760932, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_lnKnu_e", "(100.0, 1.0e-100)", 23453.606706185466825, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_lnKnu_e", "(100.0, 1.0)", 427.7532510250188083, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_lnKnu_e", "(100.0, 100.0)", -55.53422771502921431, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_lnKnu_e", "(1000.0, 1.0e-100)", 236856.183755993135, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_lnKnu_e", "(10000.0, 1.0e-100)", 2.39161558914890695e+06, TEST_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_J0_e", "( 1)",  2.404825557695771, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_J0_e", "( 2)",  5.520078110286304, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_J0_e", "(20)", 62.048469190227081, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_J0_e", "(25)", 77.756025630388058, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_J0_e", "(100)", 313.37426607752784, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_J1_e", "( 1)", 3.831705970207512, TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_J1_e", "( 2)", 7.015586669815619, TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_J1_e", "(20)", 63.61135669848124, TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_J1_e", "(25)", 79.32048717547630, TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_J1_e", "(100)", 314.9434728377672, TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(0.0,  1)",  2.404825557695771, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(0.0,  2)",  5.520078110286304, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(0.0, 20)", 62.048469190227081, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(0.0, 25)", 77.756025630388058, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(0.0, 100)", 313.37426607752784, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 1.5, 1)",  4.4934094579090641, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 1)",  8.7714838159599540, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 1.5, 2)",  7.7252518369377072, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 2)",  12.338604197466944, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 1.5, 3)",  10.904121659428900, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 3)",  15.700174079711671, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 1.5, 4)",  14.066193912831473, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 4)",  18.980133875179921, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 1.5, 5)",  17.220755271930768, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 5)",  22.217799896561268, TEST_SQRT_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 8.0, 5)",  26.266814641176644, TEST_SQRT_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(20.0, 5)",  41.413065513892636, TEST_SQRT_TOL0, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 1.5, 6)",  20.371302959287563, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 6)",  25.430341154222704, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 8.0, 6)",  29.545659670998550, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 1.5, 7)",  23.519452498689007, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 7)",  28.626618307291138, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 8.0, 7)",  32.795800037341462, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 1.5, 8)",  26.666054258812674, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 8)",  31.811716724047763, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(10.0, 8)",  38.761807017881651, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 1.5, 9)",  29.811598790892959, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 9)",  34.988781294559295, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(10.0, 9)",  42.004190236671805, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 1.5, 10)",  32.956389039822477, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 10)",  38.159868561967132, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(15.0, 10)",  52.017241278881633, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 11)", 41.326383254047406, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(15.0, 11)", 55.289204146560061, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 12)", 44.4893191232197314, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(15.0, 12)", 58.5458289043850856, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 13)", 47.6493998066970948, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(15.0, 13)", 61.7897598959450550, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 14)", 50.8071652030063595, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(15.0, 14)", 65.0230502510422545, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 15)", 53.9630265583781707, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(15.0, 15)", 68.2473219964207837, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 16)", 57.1173027815042647, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(15.0, 16)", 71.4638758850226630, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 17)", 60.2702450729428077, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(15.0, 17)", 74.6737687121404241, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 18)", 63.4220540458757799, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(15.0, 18)", 77.8778689734863729, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 19)", 66.5728918871182703, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(15.0, 19)", 81.0768977206328326, TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 5.0, 20)", 69.722891161716742,  TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(15.0, 20)", 84.271459069716442,  TEST_TOL1, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 23.0, 11)", 65.843393469524653, TEST_TOL6, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 30.0, 11)", 74.797306585175426, TEST_TOL6, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 32.0, 15)", 90.913637691861741, TEST_TOL6, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 50.0, 15)", 113.69747988073942, TEST_TOL6, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(  5.0, 22)", 76.020793430591605, TEST_TOL2, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 10.0, 22)", 83.439189796105756, TEST_TOL3, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "( 12.0, 22)", 86.345496520534055, TEST_TOL6, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(100.0, 22)", 199.82150220122519, TEST_TOL4, GSL::SUCCESS)
+  TEST_SF(s, "GSL::Sf::bessel_zero_Jnu_e", "(500.0, 22)", 649.34132440891735, TEST_TOL2, GSL::SUCCESS)