rb-gsl 1.16.0.5 → 1.16.0.6

data/test/gsl/vector_test.rb

@@ -1,223 +0,0 @@
-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

@@ -1,130 +0,0 @@
-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

@@ -1,321 +0,0 @@
-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