gsl 1.12.109 → 1.14.5

data/tests/stats.rb

@@ -0,0 +1,174 @@
+#!/usr/bin/env ruby
+require("gsl")
+require("gsl_test.rb")
+include GSL::Test
+
+lew = GSL::Vector.alloc(    -213, -564,  -35,  -15,  141,  115, -420, -360,  203, -338, -431,  194,
+    -220, -513,  154, -125, -559,   92,  -21, -579,  -52,   99, -543, -175,
+     162, -457, -346,  204, -300, -474,  164, -107, -572,   -8,   83, -541,
+    -224,  180, -420, -374,  201, -236, -531,   83,   27, -564, -112,  131,
+    -507, -254,  199, -311, -495,  143,  -46, -579,  -90,  136, -472, -338,
+     202, -287, -477,  169, -124, -568,   17,   48, -568, -135,  162, -430,
+    -422,  172,  -74, -577,  -13,   92, -534, -243,  194, -355, -465,  156,
+     -81, -578,  -64,  139, -449, -384,  193, -198, -538,  110,  -44, -577,
+      -6,   66, -552, -164,  161, -460, -344,  205, -281, -504,  134,  -28,
+    -576, -118,  156, -437, -381,  200, -220, -540,   83,   11, -568, -160,
+     172, -414, -408,  188, -125, -572,  -32,  139, -492, -321,  205, -262,
+    -504,  142,  -83, -574,    0,   48, -571, -106,  137, -501, -266,  190,
+    -391, -406,  194, -186, -553,   83,  -13, -577,  -49,  103, -515, -280,
+     201,  300, -506,  131,  -45, -578,  -80,  138, -462, -361,  201, -211,
+    -554,   32,   74, -533, -235,  187, -372, -442,  182, -147, -566,   25,
+      68, -535, -244,  194, -351, -463,  174, -125, -570,   15,   72, -550,
+    -190,  172, -424, -385,  198, -218, -536,   96 )
+
+expected_mean = -177.435000000000;
+expected_sd = 277.332168044316;
+expected_lag1 = -0.307304800605679;
+
+mean = lew.mean()
+sd = lew.sd()
+lag1 = lew.lag1_autocorrelation()
+
+GSL::Test::test_rel(mean, expected_mean, 1e-15, "lew gsl_stats_mean") ;
+GSL::Test::test_rel(sd, expected_sd, 1e-15, "lew gsl_stats_sd") ;
+GSL::Test::test_rel(lag1, expected_lag1, 1e-14, "lew autocorrelation") ;
+
+rel = 1e-10
+
+rawa = GSL::Vector.alloc(0.0421, 0.0941, 0.1064, 0.0242, 0.1331,0.0773, 0.0243, 0.0815, 0.1186, 0.0356,0.0728, 0.0999, 0.0614, 0.0479)
+
+rawb = GSL::Vector.alloc(0.1081, 0.0986, 0.1566, 0.1961, 0.1125,
+ 0.1942, 0.1079, 0.1021, 0.1583, 0.1673,
+ 0.1675, 0.1856, 0.1688, 0.1512)
+ 
+raww = GSL::Vector.alloc(0.0000, 0.0000, 0.0000, 3.000, 0.0000,
+ 1.000, 1.000, 1.000, 0.000, 0.5000,
+ 7.000, 5.000, 4.000, 0.123)
+ 
+mean = rawa.mean
+expected = 0.0728
+GSL::Test::test_rel(mean, expected, rel, "gsl_stats_mean")
+
+var = rawa.variance_with_fixed_mean(mean)
+expected = 0.00113837428571429
+GSL::Test::test_rel(var, expected, rel, "gsl_stats_variance_with_fixed_mean")
+
+sd = rawa.sd_with_fixed_mean(mean)
+expected = 0.0337398026922845
+GSL::Test::test_rel(sd, expected, rel, "gsl_stats_sd_with_fixed_mean")
+
+
+ 
+var = rawb.variance()
+expected = 0.00124956615384615
+GSL::Test::test_rel(var, expected, rel, "gsl_stats_variance")
+
+sd = rawa.sd()
+expected = 0.0350134479659107
+GSL::Test::test_rel(sd, expected, rel, "gsl_stats_sd")
+
+absdev = rawa.absdev()
+expected = 0.0287571428571429
+GSL::Test::test_rel(absdev, expected, rel, "gsl_stats_absdev")
+
+skew = rawa.skew()
+expected = 0.0954642051479004
+GSL::Test::test_rel(skew, expected, rel, "gsl_stats_skew")
+
+kurtosis = rawa.kurtosis()
+expected = -1.38583851548909
+GSL::Test::test_rel(kurtosis, expected, rel, "gsl_stats_kurtosis")
+
+wmean = rawa.wmean(raww)
+expected = 0.0678111523670601
+GSL::Test::test_rel(wmean, expected, rel, "gsl_stats_wmean")
+
+wvariance = rawa.wvariance(raww)
+expected = 0.000769562962860317
+GSL::Test::test_rel(wvariance, expected, rel, "gsl_stats_wvariance")
+
+wsd = rawa.wsd(raww)
+expected = 0.0277409978706664
+GSL::Test::test_rel(wsd, expected, rel, "gsl_stats_wsd")
+
+wabsdev = rawa.wabsdev(raww)
+expected = 0.0193205027504008
+GSL::Test::test_rel(wabsdev, expected, rel, "gsl_stats_wabsdev")
+
+wskew = rawa.wskew(raww)
+expected = -0.373631000307076
+GSL::Test::test_rel(wskew, expected, rel, "gsl_stats_wskew")
+
+wkurtosis = rawa.wkurtosis(raww)
+expected =  -1.48114233353963
+GSL::Test::test_rel(wkurtosis, expected, rel, "gsl_stats_wkurtosis")
+
+c = GSL::Stats::covariance(rawa, rawb)
+expected = -0.000139021538461539
+GSL::Test::test_rel(c, expected, rel, "gsl_stats_covariance")
+
+if GSL_VERSION >= "1.9.90"
+	r = GSL::Stats::correlation(rawa, rawb)
+	expected = -0.112322712666074171
+	GSL::Test::test_rel(r, expected, rel, "gsl_stats_correlation")
+	
+	pv = GSL::Stats::pvariance(rawa, rawb)
+	expected = 0.00123775384615385
+	GSL::Test::test_rel(pv, expected, rel, "gsl_stats_pvariance")	
+	
+	t = GSL::Stats::ttest(rawa, rawb)
+	expected =  -5.67026326985851
+	GSL::Test::test_rel(t, expected, rel, "gsl_stats_ttest")		
+end
+
+expected = 0.1331
+GSL::Test::test_rel(rawa.max, expected, rel, "gsl_stats_max")
+
+expected = 0.0242
+GSL::Test::test_rel(rawa.min, expected, rel, "gsl_stats_min")
+min, max = rawa.minmax
+expected = 0.1331
+GSL::Test::test_rel(max, expected, rel, "gsl_stats_minmax: max")
+expected = 0.0242
+GSL::Test::test_rel(min, expected, rel, "gsl_stats_minmax: min")
+
+expected = 4
+max_index = rawa.max_index
+GSL::Test::test(max_index != expected, "gsl_stats_max_index")
+
+expected = 3
+min_index = rawa.min_index
+GSL::Test::test(min_index != expected, "gsl_stats_min_index")
+
+min_index, max_index = rawa.minmax_index
+expected = 4
+GSL::Test::test(max_index != expected, "gsl_stats_minmax_index: max")
+
+expected = 3
+GSL::Test::test(min_index != expected, "gsl_stats_minmax_index: min")
+
+rawa.sort!
+
+median = rawa.median_from_sorted_data()
+expected = 0.07505
+GSL::Test::test_rel(median, expected, rel, "gsl_stats_median_from_sorted_data")
+
+median = rawa.subvector(0, rawa.size-1).median_from_sorted_data()
+expected = 0.0728
+GSL::Test::test_rel(median, expected, rel, "gsl_stats_median_from_sorted_data")
+
+quantile = rawa.quantile_from_sorted_data(0.0)
+expected = 0.0242
+GSL::Test::test_rel(quantile, expected, rel, "gsl_stats_quantile_from_sorted_data")
+
+quantile = rawa.quantile_from_sorted_data(1.0)
+expected = 0.1331
+GSL::Test::test_rel(quantile, expected, rel, "gsl_stats_quantile_from_sorted_data (100)")
+
+quantile = rawa.quantile_from_sorted_data(0.5)
+expected = 0.07505
+GSL::Test::test_rel(quantile, expected, rel, "gsl_stats_quantile_from_sorted_data (50even)")
+
+quantile = rawa.subvector(0, rawa.size-1).quantile_from_sorted_data(0.5)
+expected = 0.0728
+GSL::Test::test_rel(quantile, expected, rel, "gsl_stats_quantile_from_sorted_data (50odd)")

data/tests/sum.rb

@@ -0,0 +1,98 @@
+#!/usr/bin/env ruby
+# Ruby/GSL implementation of GSL "sum/test.c"
+require("gsl")
+require("./gsl_test2.rb")
+include GSL::Test
+include Math
+
+N = 50
+
+def check_trunc(t, expected, desc)
+  w = GSL::Sum::Levin_utrunc.alloc(N)
+  sum_accel, err, = w.accel(t)
+  desc2 = sprintf("trunc result, %s", desc)
+  GSL::Test::test_rel(sum_accel, expected, 1e-8, desc2)
+end
+
+def check_full(t, expected, desc)
+  w = GSL::Sum::Levin_u.alloc(N)
+  sum_accel, err_est, = w.accel(t)
+  desc2 = sprintf("full result, %s", desc)
+  GSL::Test::test_rel(sum_accel, expected, 1e-8, desc2)
+
+  sd_est = -log10(err_est/sum_accel.abs);
+  sd_actual = -log10(GSL::DBL_EPSILON + ((sum_accel - expected)/expected).abs);
+  desc2 = sprintf("full significant digits, %s (%g vs %g)", 
+                  desc, sd_est, sd_actual)
+  GSL::Test::test((sd_est > sd_actual + 1.0) ? 1 : 0, desc2)
+end
+
+GSL::IEEE::env_setup()
+
+Zeta_2 = M_PI*M_PI/6.0
+t = GSL::Vector.alloc(N)
+for n in 0...N
+  np1 = n + 1.0
+  t[n] = 1.0/(np1*np1)
+end
+check_trunc(t, Zeta_2, "zeta(2)")
+check_full(t, Zeta_2, "zeta(2)")
+
+x = 10.0
+y = exp(x)
+t[0] = 1.0
+for n in 1...N
+  t[n] = t[n-1]*(x/n)
+end
+check_trunc(t, y, "exp(10)");
+check_full(t, y, "exp(10)");
+
+x = -10.0
+y = exp(x)
+t[0] = 1.0
+for n in 1...N
+  t[n] = t[n-1]*(x/n)
+end
+check_trunc(t, y, "exp(-10)");
+check_full(t, y, "exp(-10)");
+
+x = 0.5
+y = -log(1-x)
+t[0] = x
+for n in 1...N
+  t[n] = t[n-1]*(x*n)/(n + 1.0)
+end
+check_trunc(t, y, "-log(1/2)")
+check_full(t, y, "-log(1/2)")
+
+x = -1.0
+y = -log(1-x)
+t[0] = x
+for n in 1...N
+  t[n] = t[n-1]*(x*n)/(n + 1.0)
+end
+check_trunc(t, y, "-log(2)")
+check_full(t, y, "-log(2)")
+
+result = 0.192594048773
+t[0] = 3.0 / (M_PI * M_PI)
+for n in 1...N
+  t[n] = -t[n - 1] * (4.0 * (n + 1.0) - 1.0) / (M_PI * M_PI)
+end
+check_trunc(t, result, "asymptotic series")
+check_full(t, result, "asymptotic series")
+
+result = 0.5772156649015328606065120900824;
+t[0] = 1.0;
+for n in 1...N
+  t[n] = 1/(n+1.0) + log(n/(n+1.0))
+end
+check_trunc(t, result, "Euler's constant")
+check_full(t, result, "Euler's constant")
+
+result = 0.6048986434216305 
+for n in 0...N
+  t[n] = (n%2 == 1 ? -1 : 1) * 1.0 /sqrt(n + 1.0)
+end
+check_trunc(t, result, "eta(1/2)")
+check_full(t, result, "eta(1/2)")

data/tests/sys.rb

@@ -0,0 +1,323 @@
+#!/usr/bin/env ruby
+# Ruby/GSL implementation of GSL "sum/test.c"
+require("gsl")
+require("./gsl_test2.rb")
+include GSL::Test
+
+GSL::IEEE::env_setup()
+
+y = GSL::expm1(0.0);
+y_expected = 0.0;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::expm1(0.0)");
+
+y = GSL::expm1(1e-10);
+y_expected = 1.000000000050000000002e-10;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::expm1(1e-10)");
+
+y = GSL::expm1(-1e-10);
+y_expected = -9.999999999500000000017e-11;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::expm1(-1e-10)");
+
+y = GSL::expm1(0.1);
+y_expected = 0.1051709180756476248117078264902;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::expm1(0.1)");
+
+y = GSL::expm1(-0.1);
+y_expected = -0.09516258196404042683575094055356;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::expm1(-0.1)");
+
+y = GSL::expm1(10.0);
+y_expected = 22025.465794806716516957900645284;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::expm1(10.0)");
+
+y = GSL::expm1(-10.0);
+y_expected = -0.99995460007023751514846440848444;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::expm1(-10.0)");
+
+# Test for log1p 
+
+y = GSL::log1p(0.0);
+y_expected = 0.0;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::log1p(0.0)");
+
+y = GSL::log1p(1e-10);
+y_expected = 9.9999999995000000000333333333308e-11;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::log1p(1e-10)");
+
+y = GSL::log1p(0.1);
+y_expected = 0.095310179804324860043952123280765;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::log1p(0.1)");
+
+y = GSL::log1p(10.0);
+y_expected = 2.3978952727983705440619435779651;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::log1p(10.0)");
+
+# Test for GSL::hypot 
+
+y = GSL::hypot(0.0, 0.0);
+y_expected = 0.0;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::hypot(0.0, 0.0)");
+
+y = GSL::hypot(1e-10, 1e-10);
+y_expected = 1.414213562373095048801688e-10;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::hypot(1e-10, 1e-10)");
+
+y = GSL::hypot(1e-38, 1e-38);
+y_expected = 1.414213562373095048801688e-38;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::hypot(1e-38, 1e-38)");
+
+y = GSL::hypot(1e-10, -1.0);
+y_expected = 1.000000000000000000005;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::hypot(1e-10, -1)");
+
+y = GSL::hypot(-1.0, 1e-10);
+y_expected = 1.000000000000000000005;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::hypot(-1, 1e-10)");
+
+#y = GSL::hypot(1e307, 1e301);
+#y_expected = 1.000000000000499999999999e307;
+#GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::hypot(1e307, 1e301)");
+
+#y = GSL::hypot(1e301, 1e307);
+#y_expected = 1.000000000000499999999999e307;
+#GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::hypot(1e301, 1e307)");
+
+#y = GSL::hypot(1e307, 1e307);
+#y_expected = 1.414213562373095048801688e307;
+#GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::hypot(1e307, 1e307)");
+
+y = GSL::acosh(1.0);
+y_expected = 0.0;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::acosh(1.0)");
+
+y = GSL::acosh(1.1);
+y_expected = 4.435682543851151891329110663525e-1;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::acosh(1.1)");
+
+y = GSL::acosh(10.0);
+y_expected = 2.9932228461263808979126677137742e0;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::acosh(10.0)");
+
+y = GSL::acosh(1e10);
+y_expected = 2.3718998110500402149594646668302e1;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::acosh(1e10)");
+
+# Test for asinh 
+
+y = GSL::asinh(0.0);
+y_expected = 0.0;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::asinh(0.0)");
+
+y = GSL::asinh(1e-10);
+y_expected = 9.9999999999999999999833333333346e-11;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::asinh(1e-10)");
+
+y = GSL::asinh(-1e-10);
+y_expected = -9.9999999999999999999833333333346e-11;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::asinh(1e-10)");
+
+y = GSL::asinh(0.1);
+y_expected = 9.983407889920756332730312470477e-2;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::asinh(0.1)");
+
+y = GSL::asinh(-0.1);
+y_expected = -9.983407889920756332730312470477e-2;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::asinh(-0.1)");
+
+y = GSL::asinh(1.0);
+y_expected = 8.8137358701954302523260932497979e-1;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::asinh(1.0)");
+
+y = GSL::asinh(-1.0);
+y_expected = -8.8137358701954302523260932497979e-1;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::asinh(-1.0)");
+
+y = GSL::asinh(10.0);
+y_expected = 2.9982229502979697388465955375965e0;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::asinh(10)");
+
+y = GSL::asinh(-10.0);
+y_expected = -2.9982229502979697388465955375965e0;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::asinh(-10)");
+
+y = GSL::asinh(1e10);
+y_expected = 2.3718998110500402149599646668302e1;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::asinh(1e10)");
+
+y = GSL::asinh(-1e10);
+y_expected = -2.3718998110500402149599646668302e1;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::asinh(-1e10)");
+
+y = GSL::atanh(0.0);
+y_expected = 0.0;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::atanh(0.0)");
+
+y = GSL::atanh(1e-20);
+y_expected = 1e-20;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::atanh(1e-20)");
+
+y = GSL::atanh(-1e-20);
+y_expected = -1e-20;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::atanh(-1e-20)");
+
+y = GSL::atanh(0.1);
+y_expected = 1.0033534773107558063572655206004e-1;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::atanh(0.1)");
+
+y = GSL::atanh(-0.1);
+y_expected = -1.0033534773107558063572655206004e-1;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::atanh(-0.1)");
+
+y = GSL::atanh(0.9);
+y_expected = 1.4722194895832202300045137159439e0;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::atanh(0.9)");
+
+y = GSL::atanh(-0.9);
+y_expected = -1.4722194895832202300045137159439e0;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::atanh(0.9)");
+
+# Test for pow_int 
+
+y = GSL::pow_2(-3.14);
+y_expected = pow(-3.14, 2.0);
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::pow_2(-3.14)");
+
+y = GSL::pow_3(-3.14);
+y_expected = pow(-3.14, 3.0);
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::pow_3(-3.14)");
+
+y = GSL::pow_4(-3.14);
+y_expected = pow(-3.14, 4.0);
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::pow_4(-3.14)");
+
+y = GSL::pow_5(-3.14);
+y_expected = pow(-3.14, 5.0);
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::pow_5(-3.14)");
+
+y = GSL::pow_6(-3.14);
+y_expected = pow(-3.14, 6.0);
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::pow_6(-3.14)");
+
+y = GSL::pow_7(-3.14);
+y_expected = pow(-3.14, 7.0);
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::pow_7(-3.14)");
+
+y = GSL::pow_8(-3.14);
+y_expected = pow(-3.14, 8.0);
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::pow_8(-3.14)");
+
+y = GSL::pow_9(-3.14);
+y_expected = pow(-3.14, 9.0);
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::pow_9(-3.14)");
+
+for n in -9...10
+  y = GSL::pow_int(-3.14, n);
+  y_expected = pow(-3.14, n);
+  GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::pow_n(-3.14,#{n})")
+end
+  
+# Test for ldexp 
+
+y = GSL::ldexp(M_PI, -2);
+y_expected = M_PI_4;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::ldexp(pi,-2)");
+
+y = GSL::ldexp(1.0, 2);
+y_expected = 4.000000;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::ldexp(1.0,2)");
+
+y = GSL::ldexp(0.0, 2);
+y_expected = 0.0;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::ldexp(0.0,2)");
+
+# Test for frexp 
+
+y, e = GSL::frexp(M_PI);
+y_expected = M_PI_4;
+e_expected = 2;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::frexp(pi) fraction");
+GSL::Test::test_int(e, e_expected, "GSL::frexp(pi) exponent");
+
+y, e = GSL::frexp(2.0);
+y_expected = 0.5;
+e_expected = 2;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::frexp(2.0) fraction");
+GSL::Test::test_int(e, e_expected, "GSL::frexp(2.0) exponent");
+
+y, e = GSL::frexp(1.0 / 4.0);
+y_expected = 0.5;
+e_expected = -1;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::frexp(0.25) fraction");
+GSL::Test::test_int(e, e_expected, "GSL::frexp(0.25) exponent");
+
+y, e = GSL::frexp(1.0 / 4.0 - 4.0 * GSL::DBL_EPSILON);
+y_expected = 0.999999999999996447;
+e_expected = -2;
+GSL::Test::test_rel(y, y_expected, 1e-15, "GSL::frexp(0.25-eps) fraction");
+GSL::Test::test_int(e, e_expected, "GSL::frexp(0.25-eps) exponent");
+
+
+x = M_PI;
+y = 22.0 / 7.0;
+
+# test the basic function 
+
+for i in 0...10
+  tol = pow(10, -i);
+  res = GSL::fcmp(x, y, tol);
+  GSL::Test::test_int(res, -((i >= 4) ? 1 : 0), "GSL::fcmp(#{x},#{y},#{tol})")
+end
+
+for i in 0...10
+  tol = pow(10, -i);
+  res = GSL::fcmp(y, x, tol);
+  GSL::Test::test_int(res,(i >= 4) ? 1 : 0, "GSL::fcmp(#{y},#{x},#{tol})")
+end
+
+zero = 0.0;
+one = 1.0;
+inf = Math::exp(1.0e10)
+nan = inf / inf;
+
+s = GSL::isinf(zero);
+GSL::Test::test_int(s, 0, "GSL::isinf(0)");
+
+s = GSL::isinf(one);
+GSL::Test::test_int(s, 0, "GSL::isinf(1)");
+
+s = GSL::isinf(inf);
+GSL::Test::test_int(s, 1, "GSL::isinf(inf)");
+
+# Commented out 2008/Oct/17 by YT
+# This test fails in (Darwin 9.5.0, gcc4.0.1):
+#  gsl_isinf() returns 1 for -inf 
+#s = GSL::isinf(-inf);
+#GSL::Test::test_int(s, -1, "GSL::isinf(-inf)");
+
+s = GSL::isinf(nan);
+GSL::Test::test_int(s, 0, "GSL::isinf(nan)");
+
+
+s = GSL::isnan(zero);
+GSL::Test::test_int(s, 0, "GSL::isnan(0)");
+
+s = GSL::isnan(one);
+GSL::Test::test_int(s, 0, "GSL::isnan(1)");
+s = GSL::isnan(inf);
+GSL::Test::test_int(s, 0, "GSL::isnan(inf)");
+
+s = GSL::isnan(nan);
+GSL::Test::test_int(s, 1, "GSL::isnan(nan)");
+
+
+s = GSL::finite(zero);
+GSL::Test::test_int(s, 1, "GSL::finite(0)");
+
+s = GSL::finite(one);
+GSL::Test::test_int(s, 1, "GSL::finite(1)");
+
+s = GSL::finite(inf);
+GSL::Test::test_int(s, 0, "GSL::finite(inf)");
+
+s = GSL::finite(nan);
+GSL::Test::test_int(s, 0, "GSL::finite(nan)");