rb-gsl 1.15.3.1 → 1.15.3.2

data/test/gsl/vector_test.rb

@@ -0,0 +1,223 @@
+require 'test_helper'
+
+class VectorTest < GSL::TestCase
+
+  def test_get
+    v = GSL::Vector::Int.indgen(5)
+    assert_equal GSL::Vector::Int[3, 1, 2], v.get([3, 1, 2])
+  end
+
+  def test_addsub
+    a = GSL::Vector::Int[2, 5, 4]
+    b = GSL::Vector::Int[10, 30, 20]
+    c = GSL::Vector::Int[12, 35, 24]
+    d = GSL::Vector::Int[8, 25, 16]
+
+    assert_equal c, a + b
+    assert_equal d, b - a
+  end
+
+  def test_collect
+    v = GSL::Vector::Int.indgen(5)
+    u = GSL::Vector::Int[0, 1, 4, 9, 16]
+
+    assert_equal u, v.collect { |val| val * val }
+  end
+
+  def test_ispos_neg
+    v = GSL::Vector::Int.indgen(5)
+    assert_equal 0,     v.ispos
+    assert_equal false, v.ispos?
+    assert_equal 0,     v.isneg
+    assert_equal false, v.isneg?
+
+    v += 1
+    assert_equal 1,     v.ispos
+    assert_equal true,  v.ispos?
+    assert_equal 0,     v.isneg
+    assert_equal false, v.isneg?
+
+    v -= 100
+    assert_equal 0,     v.ispos
+    assert_equal false, v.ispos?
+    assert_equal 1,     v.isneg
+    assert_equal true,  v.isneg?
+  end
+
+  def test_isnonneg
+    v = GSL::Vector::Int.indgen(5)
+    assert_equal 1,     v.isnonneg
+    assert_equal true,  v.isnonneg?
+    assert_equal 0,     v.isneg
+    assert_equal false, v.isneg?
+
+    v -= 100
+    assert_equal 0,     v.isnonneg
+    assert_equal false, v.isnonneg?
+    assert_equal 1,     v.isneg
+    assert_equal true,  v.isneg?
+
+    v += 200
+    assert_equal 1,     v.isnonneg
+    assert_equal true,  v.isnonneg?
+    assert_equal 1,     v.ispos
+    assert_equal true,  v.ispos?
+  end
+
+  def test_subvector
+    v = GSL::Vector::Int.indgen(12)
+
+    vv = v.subvector
+    assert_not_equal v.object_id, vv.object_id
+    assert_equal     v.subvector, v
+
+    vv = v.subvector(3)
+    assert_equal [0, 1, 2], vv.to_a
+    assert_nothing_raised('subvector(-1)') { v.subvector(-1) }
+
+    vv = v.subvector(-1)
+    assert_equal [11], vv.to_a
+
+    vv = v.subvector(-2)
+    assert_equal [10, 11], vv.to_a
+    assert_raises(RangeError) { v.subvector(-13) }
+
+    vv = v.subvector(2, 3)
+    assert_equal [2, 3, 4], vv.to_a
+
+    vv = v.subvector(-4, 3)
+    assert_equal [8, 9, 10], vv.to_a
+    assert_nothing_raised('subvector(-4, -3)') { v.subvector(-4, -3) }
+
+    vv = v.subvector(-4, -3)
+    assert_equal [8, 7, 6], vv.to_a
+    assert_raises(GSL::ERROR::EINVAL) { v.subvector(-11, -3) }
+
+    vv = v.subvector(1, 3, 4)
+    assert_equal [1, 4, 7, 10], vv.to_a
+
+    {
+    # ( range ) => [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]
+      ( 1..  3) => [   1, 2, 3                          ],
+      ( 1... 3) => [   1, 2                             ],
+      ( 3..  1) => [   3, 2, 1                          ],
+      ( 3... 1) => [      3, 2                          ],
+      (-7..  9) => [               5, 6, 7, 8, 9        ],
+      (-7... 9) => [               5, 6, 7, 8           ],
+      ( 4.. -3) => [            4, 5, 6, 7, 8, 9        ],
+      ( 4...-3) => [            4, 5, 6, 7, 8           ],
+      ( 2.. -2) => [      2, 3, 4, 5, 6, 7, 8, 9, 10    ],
+      ( 2...-2) => [      2, 3, 4, 5, 6, 7, 8, 9        ],
+      (-2..  2) => [     10, 9, 8, 7, 6, 5, 4, 3,  2    ],
+      (-2... 2) => [     10, 9, 8, 7, 6, 5, 4, 3        ],
+      (-3.. -1) => [                           9, 10, 11],
+      (-3...-1) => [                           9, 10    ],
+      (-1.. -3) => [                          11, 10,  9],
+      (-1...-3) => [                          11, 10    ]
+    }.each { |r, x|
+      assert_nothing_raised("subvector(#{r})") { v.subvector(r) }
+      assert_equal x, v.subvector(r).to_a, "subvector(#{r})"
+    }
+
+    {
+    # [( range ), s] => [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]
+      [( 1..  6), 2] => [   1,    3,    5                    ],
+      [( 1... 6), 2] => [   1,    3,    5                    ],
+      [( 0..  6), 3] => [0,       3,      6                  ],
+      [( 0... 6), 3] => [0,    3                             ]
+    }.each { |(r, s), x|
+      assert_nothing_raised("subvector(#{r},#{s})") { v.subvector(r) }
+      assert_equal x, v.subvector(r,s).to_a, "subvector(#{r},#{s})"
+    }
+  end
+
+  def test_complex_get
+    v = GSL::Vector::Complex.indgen(5)
+    assert_equal GSL::Vector::Complex[[3, 0], [1, 0], [2, 0]], v.get([3, 1, 2])
+  end
+
+  def test_complex_addsub
+    a = GSL::Vector::Complex[[-2,  5], [ 4, -1]]
+    b = GSL::Vector::Complex[[10, 30], [20, 40]]
+    c = GSL::Vector::Complex[[ 8, 35], [24, 39]]
+    d = GSL::Vector::Complex[[12, 25], [16, 41]]
+
+    assert_equal c, a + b
+    assert_equal d, b - a
+  end
+
+  def test_complex_collect
+    v = GSL::Vector::Complex.indgen(5)
+    u = GSL::Vector::Complex[[0, 0], [1, 0], [4, 0], [9, 0], [16, 0]]
+
+    assert_equal u, v.collect { |val| val * val }
+  end
+
+  def test_complex_subvector
+    v = GSL::Vector::Complex.indgen(12)
+
+    vv = v.subvector
+    assert_not_equal v.object_id, vv.object_id
+    assert_equal     v.subvector, v
+
+    vv = v.subvector(3)
+    assert_equal [0, 0, 1, 0, 2, 0], vv.to_a
+    assert_nothing_raised('subvector(-1)') { v.subvector(-1) }
+
+    vv = v.subvector(-1)
+    assert_equal [11, 0], vv.to_a
+
+    vv = v.subvector(-2)
+    assert_equal [10, 0, 11, 0], vv.to_a
+    assert_raises(RangeError) { v.subvector(-13) }
+
+    vv = v.subvector(2, 3)
+    assert_equal [2, 0, 3, 0, 4, 0], vv.to_a
+
+    vv = v.subvector(-4, 3)
+    assert_equal [8, 0, 9, 0, 10, 0], vv.to_a
+    assert_nothing_raised('subvector(-4, -3)') { v.subvector(-4, -3) }
+
+    vv = v.subvector(-4, -3)
+    assert_equal [8, 0, 7, 0, 6, 0], vv.to_a
+    assert_raises(GSL::ERROR::EINVAL) { v.subvector(-11, -3) }
+
+    vv = v.subvector(1, 3, 4)
+    assert_equal [1, 0, 4, 0, 7, 0, 10, 0], vv.to_a
+
+    {
+    # ( range ) => [0, 0, 1, 0, 2, 0, 3, 0, 4, 0, 5, 0, 6, 0, 7, 0, 8, 0, 9, 0, 10, 0, 11, 0]
+      ( 1..  3) => [      1, 0, 2, 0, 3, 0                                                  ],
+      ( 1... 3) => [      1, 0, 2, 0,                                                       ],
+      ( 3..  1) => [      3, 0, 2, 0, 1, 0                                                  ],
+      ( 3... 1) => [            3, 0, 2, 0                                                  ],
+      (-7..  9) => [                              5, 0, 6, 0, 7, 0, 8, 0, 9, 0              ],
+      (-7... 9) => [                              5, 0, 6, 0, 7, 0, 8, 0                    ],
+      ( 4.. -3) => [                        4, 0, 5, 0, 6, 0, 7, 0, 8, 0, 9, 0              ],
+      ( 4...-3) => [                        4, 0, 5, 0, 6, 0, 7, 0, 8, 0                    ],
+      ( 2.. -2) => [            2, 0, 3, 0, 4, 0, 5, 0, 6, 0, 7, 0, 8, 0, 9, 0, 10, 0       ],
+      ( 2...-2) => [            2, 0, 3, 0, 4, 0, 5, 0, 6, 0, 7, 0, 8, 0, 9, 0              ],
+      (-2..  2) => [           10, 0, 9, 0, 8, 0, 7, 0, 6, 0, 5, 0, 4, 0, 3, 0,  2, 0       ],
+      (-2... 2) => [           10, 0, 9, 0, 8, 0, 7, 0, 6, 0, 5, 0, 4, 0, 3, 0              ],
+      (-3.. -1) => [                                                      9, 0, 10, 0, 11, 0],
+      (-3...-1) => [                                                      9, 0, 10, 0       ],
+      (-1.. -3) => [                                                     11, 0, 10, 0,  9, 0],
+      (-1...-3) => [                                                     11, 0, 10, 0       ]
+    }.each { |r, x|
+      assert_nothing_raised("subvector(#{r})") { v.subvector(r) }
+      assert_equal x, v.subvector(r).to_a, "subvector(#{r})"
+    }
+
+    {
+    # [( range ), s] => [0, 0, 1, 0, 2, 0, 3, 0, 4, 0, 5, 0, 6, 0, 7, 0, 8, 0, 9, 0, 10, 0, 11, 0]
+      [( 1..  6), 2] => [      1, 0,       3, 0,       5, 0                                      ],
+      [( 1... 6), 2] => [      1, 0,       3, 0,       5, 0                                      ],
+      [( 0..  6), 3] => [0, 0,             3, 0,             6, 0                                ],
+      [( 0... 6), 3] => [0, 0,             3, 0                                                  ]
+    }.each { |(r, s), x|
+      assert_nothing_raised("subvector(#{r},#{s})") { v.subvector(r) }
+      assert_equal x, v.subvector(r,s).to_a, "subvector(#{r},#{s})"
+    }
+  end
+
+end

data/test/gsl/wavelet_test.rb

@@ -0,0 +1,130 @@
+require 'test_helper'
+
+class WaveletTest < GSL::TestCase
+
+  MEMBERS = [309, 307, 305, 303, 301, 208, 206, 204, 202, 105, 103]
+
+  def _urand
+    x = 1
+    x = (1103515245 * x + 12345) & 0x7fffffff
+    x / 2147483648.0
+  end
+
+  def _test_1d(n, stride, type, member)
+    nn = n * stride
+    data = GSL::Vector.alloc(nn)
+    nn.times { |i| data[i] = 12345.0 + i }
+
+    v1 = data.view_with_stride(0, stride, n)
+    n.times { |i| v1[i] = _urand }
+
+    v2 = GSL::Vector.alloc(n)
+    GSL::Vector.memcpy(v2, v1)
+
+    vdelta = GSL::Vector.alloc(n)
+
+    work = GSL::Wavelet::Workspace.alloc(n)
+    w = GSL::Wavelet.alloc(type, member)
+    w.transform_forward(v2, work)
+    w.transform_inverse(v2, work)
+
+    n.times { |i| vdelta[i] = (v1[i] - v2[i]).abs }
+
+    i = vdelta.max_index
+    x1, x2 = v1[i], v2[i]
+
+    refute((x2 - x1).abs > n * 1e-15,
+      "#{w.name}(#{member}), n = #{n}, stride = #{stride}, maxerr = #{(x2 - x1).abs}")
+
+    assert((0...nn).all? { |j| j % stride == 0 || data[j] == 12345.0 + j },
+      "#{w.name}(#{member}) other data untouched, n = #{n}, stride = #{stride}") if stride > 1
+  end
+
+  def _test_2d(n, tda, t, member, type)
+    nn = n * tda
+    data = GSL::Vector.alloc(nn)
+    nn.times { |i| data[i] = 12345.0 + i }
+
+    m1 = data.matrix_view_with_tda(n, n, tda)
+    n.times { |i| n.times { |j| m1.set(i, j, _urand) } }
+
+    m2 = GSL::Matrix.alloc(n, n)
+    GSL::Matrix.memcpy(m2, m1)
+
+    mdelta = GSL::Matrix.alloc(n, n)
+
+    work = GSL::Wavelet::Workspace.alloc(n)
+    w = GSL::Wavelet.alloc(t, member)
+
+    typename = case type
+      when 1
+        GSL::Wavelet2d.transform_matrix_forward(w, m2, work)
+        GSL::Wavelet2d.transform_matrix_inverse(w, m2, work)
+        'standard'
+      when 2
+        GSL::Wavelet2d.nstransform_matrix_forward(w, m2, work)
+        GSL::Wavelet2d.nstransform_matrix_inverse(w, m2, work)
+        'nonstd'
+    end
+
+    n.times { |i| n.times { |j| mdelta.set(i, j, (m1[i, j] - m2[i, j]).abs) } }
+
+    i, j = mdelta.max_index
+    x1, x2 = m1[i, j], m1[i, j]
+
+    refute((x2 - x1).abs > n * 1e-15,
+      "#{w.name}(#{member})-2d #{typename}, n = #{n}, tda = #{tda}, maxerr = #{(x2 - x1).abs}")
+
+    assert((0...n).to_a.product((n...tda).to_a).all? { |k, l| data[k * tda + l] == 12345.0 + k * tda + l },
+      "#{w.name}(#{member})-2d #{typename} other data untouched, n = #{n}, tda = #{tda}") if tda > n
+  end
+
+  def _each_pow(n = 14)
+    n.times { |i| yield 2 ** i }
+  end
+
+  def _each_n(n = 9, m = 4)
+    n.times { |i| yield m + 2 * i }
+  end
+
+  def test_1d_bspline
+    _each_pow { |n| MEMBERS.each { |m| _test_1d(n, 1, 'bspline', m) } }
+  end
+
+  def test_1d_bspline_centered
+    _each_pow { |n| MEMBERS.each { |m| _test_1d(n, 1, 'bspline_centered', m) } }
+  end
+
+  def test_1d_daubechies
+    _each_pow { |n| _each_n { |i| _test_1d(n, 1, 'daubechies', i) } }
+  end
+
+  def test_1d_daubechies_centered
+    _each_pow { |n| _each_n { |i|_test_1d(n, 1, 'daubechies_centered', i) } }
+  end
+
+  def test_1d_haar
+    _each_pow { |n| _test_1d(n, 1, 'haar', 2) }
+  end
+
+  def test_1d_haar_centered
+    _each_pow { |n| _test_1d(n, 1, 'haar_centered', 2) }
+  end
+
+  def test_2d_bspline_standard
+    _each_pow(6) { |n| MEMBERS.each { |m| _test_2d(n, n, 'bspline', m, 1) } }
+  end
+
+  def test_2d_bspline_centered_standard
+    _each_pow(6) { |n| MEMBERS.each { |m| _test_2d(n, n, 'bspline_centered', m, 1) } }
+  end
+
+  def test_2d_bspline_nonstd
+    _each_pow(6) { |n| MEMBERS.each { |m| _test_2d(n, n, 'bspline', m, 2) } }
+  end
+
+  def test_2d_bspline_centered_nonstd
+    _each_pow(6) { |n| MEMBERS.each { |m| _test_2d(n, n, 'bspline_centered', m, 2) } }
+  end
+
+end

data/test/gsl_test.rb

@@ -0,0 +1,321 @@
+require 'test_helper'
+
+class GSLTest < GSL::TestCase
+
+  def test_expm1
+    y = GSL.expm1(0.0)
+    y_expected = 0.0
+    assert_rel y, y_expected, 1e-15, 'GSL.expm1(0.0)'
+
+    y = GSL.expm1(1e-10)
+    y_expected = 1.000000000050000000002e-10
+    assert_rel y, y_expected, 1e-15, 'GSL.expm1(1e-10)'
+
+    y = GSL.expm1(-1e-10)
+    y_expected = -9.999999999500000000017e-11
+    assert_rel y, y_expected, 1e-15, 'GSL.expm1(-1e-10)'
+
+    y = GSL.expm1(0.1)
+    y_expected = 0.1051709180756476248117078264902
+    assert_rel y, y_expected, 1e-15, 'GSL.expm1(0.1)'
+
+    y = GSL.expm1(-0.1)
+    y_expected = -0.09516258196404042683575094055356
+    assert_rel y, y_expected, 1e-15, 'GSL.expm1(-0.1)'
+
+    y = GSL.expm1(10.0)
+    y_expected = 22025.465794806716516957900645284
+    assert_rel y, y_expected, 1e-15, 'GSL.expm1(10.0)'
+
+    y = GSL.expm1(-10.0)
+    y_expected = -0.99995460007023751514846440848444
+    assert_rel y, y_expected, 1e-15, 'GSL.expm1(-10.0)'
+  end
+
+  def test_log1p
+    y = GSL.log1p(0.0)
+    y_expected = 0.0
+    assert_rel y, y_expected, 1e-15, 'GSL.log1p(0.0)'
+
+    y = GSL.log1p(1e-10)
+    y_expected = 9.9999999995000000000333333333308e-11
+    assert_rel y, y_expected, 1e-15, 'GSL.log1p(1e-10)'
+
+    y = GSL.log1p(0.1)
+    y_expected = 0.095310179804324860043952123280765
+    assert_rel y, y_expected, 1e-15, 'GSL.log1p(0.1)'
+
+    y = GSL.log1p(10.0)
+    y_expected = 2.3978952727983705440619435779651
+    assert_rel y, y_expected, 1e-15, 'GSL.log1p(10.0)'
+  end
+
+  def test_hypot
+    y = GSL.hypot(0.0, 0.0)
+    y_expected = 0.0
+    assert_rel y, y_expected, 1e-15, 'GSL.hypot(0.0, 0.0)'
+
+    y = GSL.hypot(1e-10, 1e-10)
+    y_expected = 1.414213562373095048801688e-10
+    assert_rel y, y_expected, 1e-15, 'GSL.hypot(1e-10, 1e-10)'
+
+    y = GSL.hypot(1e-38, 1e-38)
+    y_expected = 1.414213562373095048801688e-38
+    assert_rel y, y_expected, 1e-15, 'GSL.hypot(1e-38, 1e-38)'
+
+    y = GSL.hypot(1e-10, -1.0)
+    y_expected = 1.000000000000000000005
+    assert_rel y, y_expected, 1e-15, 'GSL.hypot(1e-10, -1)'
+
+    y = GSL.hypot(-1.0, 1e-10)
+    y_expected = 1.000000000000000000005
+    assert_rel y, y_expected, 1e-15, 'GSL.hypot(-1, 1e-10)'
+
+    #y = GSL.hypot(1e307, 1e301)
+    #y_expected = 1.000000000000499999999999e307
+    #assert_rel y, y_expected, 1e-15, 'GSL.hypot(1e307, 1e301)'
+
+    #y = GSL.hypot(1e301, 1e307)
+    #y_expected = 1.000000000000499999999999e307
+    #assert_rel y, y_expected, 1e-15, 'GSL.hypot(1e301, 1e307)'
+
+    #y = GSL.hypot(1e307, 1e307)
+    #y_expected = 1.414213562373095048801688e307
+    #assert_rel y, y_expected, 1e-15, 'GSL.hypot(1e307, 1e307)'
+  end
+
+  def test_acosh
+    y = GSL.acosh(1.0)
+    y_expected = 0.0
+    assert_rel y, y_expected, 1e-15, 'GSL.acosh(1.0)'
+
+    y = GSL.acosh(1.1)
+    y_expected = 4.435682543851151891329110663525e-1
+    assert_rel y, y_expected, 1e-15, 'GSL.acosh(1.1)'
+
+    y = GSL.acosh(10.0)
+    y_expected = 2.9932228461263808979126677137742e0
+    assert_rel y, y_expected, 1e-15, 'GSL.acosh(10.0)'
+
+    y = GSL.acosh(1e10)
+    y_expected = 2.3718998110500402149594646668302e1
+    assert_rel y, y_expected, 1e-15, 'GSL.acosh(1e10)'
+  end
+
+  def test_asinh
+    y = GSL.asinh(0.0)
+    y_expected = 0.0
+    assert_rel y, y_expected, 1e-15, 'GSL.asinh(0.0)'
+
+    y = GSL.asinh(1e-10)
+    y_expected = 9.9999999999999999999833333333346e-11
+    assert_rel y, y_expected, 1e-15, 'GSL.asinh(1e-10)'
+
+    y = GSL.asinh(-1e-10)
+    y_expected = -9.9999999999999999999833333333346e-11
+    assert_rel y, y_expected, 1e-15, 'GSL.asinh(1e-10)'
+
+    y = GSL.asinh(0.1)
+    y_expected = 9.983407889920756332730312470477e-2
+    assert_rel y, y_expected, 1e-15, 'GSL.asinh(0.1)'
+
+    y = GSL.asinh(-0.1)
+    y_expected = -9.983407889920756332730312470477e-2
+    assert_rel y, y_expected, 1e-15, 'GSL.asinh(-0.1)'
+
+    y = GSL.asinh(1.0)
+    y_expected = 8.8137358701954302523260932497979e-1
+    assert_rel y, y_expected, 1e-15, 'GSL.asinh(1.0)'
+
+    y = GSL.asinh(-1.0)
+    y_expected = -8.8137358701954302523260932497979e-1
+    assert_rel y, y_expected, 1e-15, 'GSL.asinh(-1.0)'
+
+    y = GSL.asinh(10.0)
+    y_expected = 2.9982229502979697388465955375965e0
+    assert_rel y, y_expected, 1e-15, 'GSL.asinh(10)'
+
+    y = GSL.asinh(-10.0)
+    y_expected = -2.9982229502979697388465955375965e0
+    assert_rel y, y_expected, 1e-15, 'GSL.asinh(-10)'
+
+    y = GSL.asinh(1e10)
+    y_expected = 2.3718998110500402149599646668302e1
+    assert_rel y, y_expected, 1e-15, 'GSL.asinh(1e10)'
+
+    y = GSL.asinh(-1e10)
+    y_expected = -2.3718998110500402149599646668302e1
+    assert_rel y, y_expected, 1e-15, 'GSL.asinh(-1e10)'
+  end
+
+  def test_atanh
+    y = GSL.atanh(0.0)
+    y_expected = 0.0
+    assert_rel y, y_expected, 1e-15, 'GSL.atanh(0.0)'
+
+    y = GSL.atanh(1e-20)
+    y_expected = 1e-20
+    assert_rel y, y_expected, 1e-15, 'GSL.atanh(1e-20)'
+
+    y = GSL.atanh(-1e-20)
+    y_expected = -1e-20
+    assert_rel y, y_expected, 1e-15, 'GSL.atanh(-1e-20)'
+
+    y = GSL.atanh(0.1)
+    y_expected = 1.0033534773107558063572655206004e-1
+    assert_rel y, y_expected, 1e-15, 'GSL.atanh(0.1)'
+
+    y = GSL.atanh(-0.1)
+    y_expected = -1.0033534773107558063572655206004e-1
+    assert_rel y, y_expected, 1e-15, 'GSL.atanh(-0.1)'
+
+    y = GSL.atanh(0.9)
+    y_expected = 1.4722194895832202300045137159439e0
+    assert_rel y, y_expected, 1e-15, 'GSL.atanh(0.9)'
+
+    y = GSL.atanh(-0.9)
+    y_expected = -1.4722194895832202300045137159439e0
+    assert_rel y, y_expected, 1e-15, 'GSL.atanh(0.9)'
+  end
+
+  def test_pow_int
+    y = GSL.pow_2(-3.14)
+    y_expected = GSL.pow(-3.14, 2.0)
+    assert_rel y, y_expected, 1e-15, 'GSL.pow_2(-3.14)'
+
+    y = GSL.pow_3(-3.14)
+    y_expected = GSL.pow(-3.14, 3.0)
+    assert_rel y, y_expected, 1e-15, 'GSL.pow_3(-3.14)'
+
+    y = GSL.pow_4(-3.14)
+    y_expected = GSL.pow(-3.14, 4.0)
+    assert_rel y, y_expected, 1e-15, 'GSL.pow_4(-3.14)'
+
+    y = GSL.pow_5(-3.14)
+    y_expected = GSL.pow(-3.14, 5.0)
+    assert_rel y, y_expected, 1e-15, 'GSL.pow_5(-3.14)'
+
+    y = GSL.pow_6(-3.14)
+    y_expected = GSL.pow(-3.14, 6.0)
+    assert_rel y, y_expected, 1e-15, 'GSL.pow_6(-3.14)'
+
+    y = GSL.pow_7(-3.14)
+    y_expected = GSL.pow(-3.14, 7.0)
+    assert_rel y, y_expected, 1e-15, 'GSL.pow_7(-3.14)'
+
+    y = GSL.pow_8(-3.14)
+    y_expected = GSL.pow(-3.14, 8.0)
+    assert_rel y, y_expected, 1e-15, 'GSL.pow_8(-3.14)'
+
+    y = GSL.pow_9(-3.14)
+    y_expected = GSL.pow(-3.14, 9.0)
+    assert_rel y, y_expected, 1e-15, 'GSL.pow_9(-3.14)'
+
+    -9.upto(9) { |n|
+      y = GSL.pow_int(-3.14, n)
+      y_expected = GSL.pow(-3.14, n)
+      assert_rel y, y_expected, 1e-15, "GSL.pow_n(-3.14,#{n})"
+    }
+  end
+
+  def test_ldexp
+    y = GSL.ldexp(GSL::M_PI, -2)
+    y_expected = GSL::M_PI_4
+    assert_rel y, y_expected, 1e-15, 'GSL.ldexp(pi,-2)'
+
+    y = GSL.ldexp(1.0, 2)
+    y_expected = 4.000000
+    assert_rel y, y_expected, 1e-15, 'GSL.ldexp(1.0,2)'
+
+    y = GSL.ldexp(0.0, 2)
+    y_expected = 0.0
+    assert_rel y, y_expected, 1e-15, 'GSL.ldexp(0.0,2)'
+  end
+
+  def test_frexp
+    y, e = GSL.frexp(GSL::M_PI)
+    y_expected = GSL::M_PI_4
+    e_expected = 2
+    assert_rel y, y_expected, 1e-15, 'GSL.frexp(pi) fraction'
+    assert_int e, e_expected, 'GSL.frexp(pi) exponent'
+
+    y, e = GSL.frexp(2.0)
+    y_expected = 0.5
+    e_expected = 2
+    assert_rel y, y_expected, 1e-15, 'GSL.frexp(2.0) fraction'
+    assert_int e, e_expected, 'GSL.frexp(2.0) exponent'
+
+    y, e = GSL.frexp(1.0 / 4.0)
+    y_expected = 0.5
+    e_expected = -1
+    assert_rel y, y_expected, 1e-15, 'GSL.frexp(0.25) fraction'
+    assert_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
+    assert_rel y, y_expected, 1e-15, 'GSL.frexp(0.25-eps) fraction'
+    assert_int e, e_expected, 'GSL.frexp(0.25-eps) exponent'
+  end
+
+  def test_gsl
+    x = GSL::M_PI
+    y = 22.0 / 7.0
+
+    10.times { |i|
+      tol = GSL.pow(10, -i)
+      res = GSL.fcmp(x, y, tol)
+      assert_int res, i >= 4 ? -1 : 0, "GSL.fcmp(#{x},#{y},#{tol})"
+
+      res = GSL.fcmp(y, x, tol)
+      assert_int res, i >= 4 ? 1 : 0, "GSL.fcmp(#{y},#{x},#{tol})"
+    }
+
+    zero = 0.0
+    one = 1.0
+    inf = Math.exp(1.0e10)
+    nan = inf / inf
+
+    s = GSL.isinf(zero)
+    assert_int s, 0, 'GSL.isinf(0)'
+
+    s = GSL.isinf(one)
+    assert_int s, 0, 'GSL.isinf(1)'
+
+    s = GSL.isinf(inf)
+    assert_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)
+    #assert_int s, -1, 'GSL.isinf(-inf)'
+
+    s = GSL.isinf(nan)
+    assert_int s, 0, 'GSL.isinf(nan)'
+
+    s = GSL.isnan(zero)
+    assert_int s, 0, 'GSL.isnan(0)'
+
+    s = GSL.isnan(one)
+    assert_int s, 0, 'GSL.isnan(1)'
+    s = GSL.isnan(inf)
+    assert_int s, 0, 'GSL.isnan(inf)'
+
+    s = GSL.isnan(nan)
+    assert_int s, 1, 'GSL.isnan(nan)'
+
+    s = GSL.finite(zero)
+    assert_int s, 1, 'GSL.finite(0)'
+
+    s = GSL.finite(one)
+    assert_int s, 1, 'GSL.finite(1)'
+
+    s = GSL.finite(inf)
+    assert_int s, 0, 'GSL.finite(inf)'
+
+    s = GSL.finite(nan)
+    assert_int s, 0, 'GSL.finite(nan)'
+  end
+
+end