rb-gsl 1.16.0.5 → 1.16.0.6

data/test/gsl/err_test.rb

@@ -1,23 +0,0 @@
-require 'test_helper'
-
-class ErrTest < GSL::TestCase
-
-  MAX_ERRS = 64
-
-  ERRORS = %w[
-    SUCCESS FAILURE CONTINUE EDOM ERANGE EFAULT EINVAL EFAILED
-    EFACTOR ESANITY ENOMEM EBADFUNC ERUNAWAY EMAXITER EZERODIV
-    EBADTOL ETOL EUNDRFLW EOVRFLW ELOSS EROUND EBADLEN ENOTSQR
-    ESING EDIVERGE EUNSUP EUNIMPL ECACHE ETABLE ENOPROG ENOPROGJ
-    ETOLF ETOLX ETOLG EOF
-  ].map { |name| GSL.const_get(name) }
-
-  def test_number
-    assert ERRORS.uniq == ERRORS
-  end
-
-  def test_message
-    assert ERRORS.map { |e| GSL.strerror(e) }.uniq.size == ERRORS.size
-  end
-
-end

data/test/gsl/fit_test.rb

@@ -1,101 +0,0 @@
-require 'test_helper'
-
-class FitTest < GSL::TestCase
-
-  def setup
-    @x = GSL::Vector.alloc(1000)
-    @y = GSL::Vector.alloc(1000)
-    @w = GSL::Vector.alloc(1000)
-
-    @x.stride = 2
-    @w.stride = 3
-    @y.stride = 5
-  end
-
-  def test_fit_1
-    n = 36
-
-    x = GSL::Vector.alloc(
-      0.2, 337.4, 118.2, 884.6, 10.1, 226.5, 666.3, 996.3,
-      448.6, 777.0, 558.2, 0.4, 0.6, 775.5, 666.9, 338.0,
-      447.5, 11.6, 556.0, 228.1, 995.8, 887.6, 120.2, 0.3,
-      0.3, 556.8, 339.1, 887.2, 999.0, 779.0, 11.1, 118.3,
-      229.2, 669.1, 448.9, 0.5
-    )
-
-    y = GSL::Vector.alloc(
-      0.1, 338.8, 118.1, 888.0, 9.2, 228.1, 668.5, 998.5,
-      449.1, 778.9, 559.2, 0.3, 0.1, 778.1, 668.8, 339.3,
-      448.9, 10.8, 557.7, 228.3, 998.0, 888.8, 119.6, 0.3,
-      0.6, 557.6, 339.3, 888.0, 998.5, 778.9, 10.2, 117.6,
-      228.9, 668.4, 449.2, 0.2
-    )
-
-    n.times { |i|
-      @x.set(i, x[i])
-      @w.set(i, 1.0)
-      @y.set(i, y[i])
-    }
-
-    c0, c1, cov00, cov01, cov11, sumsq = GSL::Fit.linear(@x, @y, n)
-    assert_rel c0,           -0.262323073774029,           1e-10, 'norris gsl_fit_linear c0'
-    assert_rel c1,            1.00211681802045,            1e-10, 'norris gsl_fit_linear c1'
-    assert_rel cov00, GSL.pow(0.232818234301152, 2.0),     1e-10, 'norris gsl_fit_linear cov00'
-    assert_rel cov01,        -7.74327536339570e-05,        1e-10, 'norris gsl_fit_linear cov01'
-    assert_rel cov11, GSL.pow(0.429796848199937E-03, 2.0), 1e-10, 'norris gsl_fit_linear cov11'
-    assert_rel sumsq,        26.6173985294224,             1e-10, 'norris gsl_fit_linear sumsq'
-
-    c0, c1, cov00, cov01, cov11, sumsq = GSL::Fit.wlinear(@x, @w, @y, n)
-    assert_rel c0,    -0.262323073774029,     1e-10, 'norris gsl_fit_wlinear c0'
-    assert_rel c1,     1.00211681802045,      1e-10, 'norris gsl_fit_wlinear c1'
-    assert_rel cov00,  6.92384428759429e-02,  1e-10, 'norris gsl_fit_wlinear cov00'
-    assert_rel cov01, -9.89095016390515e-05,  1e-10, 'norris gsl_fit_wlinear cov01'
-    assert_rel cov11,  2.35960747164148e-07,  1e-10, 'norris gsl_fit_wlinear cov11'
-    assert_rel sumsq, 26.6173985294224,       1e-10, 'norris gsl_fit_wlinear sumsq'
-  end
-
-  def test_fit_2
-    n = 11
-    x = GSL::Vector.alloc(60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70)
-    y = GSL::Vector.alloc(130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140)
-
-    n.times { |i|
-      @x.set(i, x[i])
-      @w.set(i, 1.0)
-      @y.set(i, y[i])
-    }
-
-    c1, cov11, sumsq = GSL::Fit.mul(@x, @y, n)
-    assert_rel c1,            2.07438016528926,            1e-10, 'noint1 gsl_fit_mul c1'
-    assert_rel cov11, GSL.pow(0.165289256198347E-01, 2.0), 1e-10, 'noint1 gsl_fit_mul cov11'
-    assert_rel sumsq,       127.272727272727,              1e-10, 'noint1 gsl_fit_mul sumsq'
-
-    c1, cov11, sumsq = GSL::Fit.wmul(@x, @w, @y, n)
-    assert_rel c1,      2.07438016528926,     1e-10, 'noint1 gsl_fit_wmul c1'
-    assert_rel cov11,   2.14661371686165e-05, 1e-10, 'noint1 gsl_fit_wmul cov11'
-    assert_rel sumsq, 127.272727272727,       1e-10, 'noint1 gsl_fit_wmul sumsq'
-  end
-
-  def test_fit_3
-    n = 3
-    x = GSL::Vector.alloc(4, 5, 6)
-    y = GSL::Vector.alloc(3, 4, 4)
-
-    n.times { |i|
-      @x.set(i, x[i])
-      @w.set(i, 1.0)
-      @y.set(i, y[i])
-    }
-
-    c1, cov11, sumsq = GSL::Fit.mul(@x, @y, n)
-    assert_rel c1,            0.727272727272727,           1e-10, 'noint2 gsl_fit_mul c1'
-    assert_rel cov11, GSL.pow(0.420827318078432E-01, 2.0), 1e-10, 'noint2 gsl_fit_mul cov11'
-    assert_rel sumsq,         0.272727272727273,           1e-10, 'noint2 gsl_fit_mul sumsq'
-
-    c1, cov11, sumsq = GSL::Fit.wmul(@x, @w, @y, n)
-    assert_rel c1,    0.727272727272727,    1e-10, 'noint2 gsl_fit_wmul c1'
-    assert_rel cov11, 1.29870129870130e-02, 1e-10, 'noint2 gsl_fit_wmul cov11'
-    assert_rel sumsq, 0.272727272727273,    1e-10, 'noint2 gsl_fit_wmul sumsq'
-  end
-
-end

data/test/gsl/histo_test.rb

@@ -1,14 +0,0 @@
-require 'test_helper'
-
-class HistoTest < GSL::TestCase
-
-  def test_histo
-    h = GSL::Histogram.alloc(10, [0, 10])
-
-    assert h
-    assert h.get_range(2)
-    assert h.range
-    assert h.bin
-  end
-
-end

data/test/gsl/index_test.rb

@@ -1,61 +0,0 @@
-require 'test_helper'
-
-class IndexTest < GSL::TestCase
-
-  # helper(s)
-
-  def _create_index array
-    i = GSL::Index.alloc(array.size)
-    array.each_with_index { |e, j| i[j] = e }
-    i
-  end
-
-  # tests
-
-  def test_get_int
-    i = GSL::Index.alloc(5)
-    assert_equal 2, i.get(2)
-    assert_equal 4, i.get(4)
-    assert_equal 3, i.get(-2)
-    assert_equal 0, i.get(-5)
-  end
-
-  def test_get_array
-    i = GSL::Index.alloc(5)
-    assert_equal _create_index([2, 3]), i.get([2, 3])
-    assert_equal _create_index([4, 2]), i.get([-1, 2])
-    assert_equal _create_index([4, 3, 1]), i.get([4, -2, 1])
-    assert_equal _create_index([0, 4, 4, 0]), i.get([0, 4, -1, -5])
-  end
-
-  def test_get_range
-    i = GSL::Index.alloc(5)
-    assert_equal _create_index([2, 3]), i.get(2..3)
-    assert_equal _create_index([0, 1, 2, 3, 4]), i.get(0..4)
-    assert_equal _create_index([0, 1, 2, 3, 4]), i.get(0..5)
-    assert_equal _create_index([0, 1, 2, 3, 4]), i.get(-5..-1)
-    assert_equal _create_index([0]), i.get(0..-5)
-    assert_equal _create_index([0]), i.get(-5..0)
-    assert_equal _create_index([3]), i.get(-2..3)
-    assert_equal _create_index([3, 4]), i.get(-2..4)
-    assert_equal _create_index([1, 2, 3]), i.get(-4...-1) # note the exclusive range operator!
-  end
-
-  def test_get_failure
-    i = GSL::Index.alloc(5)
-    assert_raises(RangeError) { i.get(5) }
-    assert_raises(RangeError) { i.get(-6) }
-    assert_raises(RangeError) { i.get(1_000_000) }
-    assert_raises(ArgumentError) { i.get(10**100) }
-
-    assert_raises(RangeError) { i.get([5]) }
-    assert_raises(RangeError) { i.get([-6]) }
-    assert_raises(RangeError) { i.get([-6, 0, 5]) }
-
-    assert_raises(RangeError) { i.get(-6..-1) }
-    assert_raises(RangeError) { i.get(3..-5) }
-    assert_raises(RangeError) { i.get(-2..-3) }
-    assert_raises(RangeError) { i.get(-2..1) }
-  end
-
-end

data/test/gsl/integration_test.rb

@@ -1,274 +0,0 @@
-require 'test_helper'
-
-class IntegrationTest < GSL::TestCase
-
-  def test_integration1
-    f = GSL::Function.alloc { |x| Math.exp(x) * Math.cos(x) }
-
-    xmin = 0.0
-    xmax = 1.0
-    limit = 1000
-
-    # QNG
-    assert f.integration_qng(xmin, xmax, 0.0, 1.0e-7)
-    assert f.integration_qng(xmin, xmax)
-    assert f.integration_qng([xmin, xmax])
-    assert f.integration_qng([xmin, xmax], [0.0, 1.0e-7])
-    assert f.integration_qng([xmin, xmax], 0.0, 1.0e-7)
-    assert f.integration_qng(xmin, xmax, [0.0, 1.0e-7])
-
-    # QAG
-    assert f.integration_qag(xmin, xmax, 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag([xmin, xmax], 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag(xmin, xmax, GSL::Integration::GAUSS15)
-    assert f.integration_qag([xmin, xmax], GSL::Integration::GAUSS15)
-
-    w = GSL::Integration::Workspace.alloc
-    assert f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
-    assert f.integration_qag(xmin, xmax, 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
-    assert f.integration_qag([xmin, xmax], 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
-    assert f.integration_qag([xmin, xmax], [0.0, 1.0e-7], GSL::Integration::GAUSS15, w)
-    assert f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
-    assert f.integration_qag([xmin, xmax], w)
-
-    # QAGS
-    assert f.integration_qags(xmin, xmax)
-    assert f.integration_qags([xmin, xmax])
-    assert f.integration_qags(xmin, xmax, 0.0, 1e-7)
-    assert f.integration_qags(xmin, xmax, [0.0, 1e-7])
-    assert f.integration_qags([xmin, xmax], [0.0, 1e-7])
-    assert f.integration_qags([xmin, xmax], 0.0, 1e-7)
-    assert f.integration_qags([xmin, xmax], 0.0, 1e-7, limit)
-    assert f.integration_qags([xmin, xmax], [0.0, 1e-7], limit)
-    assert f.integration_qags(xmin, xmax, [0.0, 1e-7], limit)
-    assert f.integration_qags([xmin, xmax], 0.0, 1e-7, limit, w)
-    assert f.integration_qags(xmin, xmax, 0.0, 1e-7, limit, w)
-    assert f.integration_qags(xmin, xmax, [0.0, 1e-7], w)
-    assert f.integration_qags([xmin, xmax], [0.0, 1e-7], w)
-
-    assert f.integration_qags([xmin, xmax], limit)
-    assert f.integration_qags(xmin, xmax, limit)
-    assert f.integration_qags([xmin, xmax], w)
-    assert f.integration_qags(xmin, xmax, w)
-    assert f.integration_qags(xmin, xmax, limit, w)
-    assert f.integration_qags([xmin, xmax], limit, w)
-
-    # QAGP
-    assert f.integration_qagp([xmin, xmax])
-    assert f.integration_qagp([xmin, xmax], [0.0, 1e-7])
-    assert f.integration_qagp([xmin, xmax], 0.0, 1e-7)
-    assert f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit)
-    assert f.integration_qagp([xmin, xmax], [0.0, 1e-7], limit)
-    assert f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit, w)
-    assert f.integration_qagp([xmin, xmax], [0.0, 1e-7], w)
-
-    assert f.integration_qagp([xmin, xmax], limit)
-    assert f.integration_qagp([xmin, xmax], w)
-    assert f.integration_qagp([xmin, xmax], limit, w)
-  end
-
-  def test_integration2
-    f = GSL::Function.alloc { |x| Math.sin(x) / x }
-
-    xmin = 0.0
-    xmax = 2.0 * Math::PI
-    limit = 1000
-
-    # QNG
-    assert f.integration_qng(xmin, xmax, 0.0, 1.0e-7)
-    assert f.integration_qng(xmin, xmax)
-    assert f.integration_qng([xmin, xmax])
-    assert f.integration_qng([xmin, xmax], [0.0, 1.0e-7])
-    assert f.integration_qng([xmin, xmax], 0.0, 1.0e-7)
-    assert f.integration_qng(xmin, xmax, [0.0, 1.0e-7])
-
-    # QAG
-    assert f.integration_qag(xmin, xmax, 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag([xmin, xmax], 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag(xmin, xmax, GSL::Integration::GAUSS15)
-    assert f.integration_qag([xmin, xmax], GSL::Integration::GAUSS15)
-
-    w = GSL::Integration::Workspace.alloc(2000)
-    assert f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
-    assert f.integration_qag(xmin, xmax, 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
-    assert f.integration_qag([xmin, xmax], 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
-    assert f.integration_qag([xmin, xmax], [0.0, 1.0e-7], GSL::Integration::GAUSS15, w)
-    assert f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
-
-    # QAGS
-    assert f.integration_qags(xmin, xmax)
-    assert f.integration_qags([xmin, xmax])
-    assert f.integration_qags(xmin, xmax, 0.0, 1e-7)
-    assert f.integration_qags(xmin, xmax, [0.0, 1e-7])
-    assert f.integration_qags([xmin, xmax], [0.0, 1e-7])
-    assert f.integration_qags([xmin, xmax], 0.0, 1e-7)
-    assert f.integration_qags([xmin, xmax], 0.0, 1e-7, limit)
-    assert f.integration_qags([xmin, xmax], [0.0, 1e-7], limit)
-    assert f.integration_qags(xmin, xmax, [0.0, 1e-7], limit)
-    assert f.integration_qags([xmin, xmax], 0.0, 1e-7, limit, w)
-    assert f.integration_qags(xmin, xmax, 0.0, 1e-7, limit, w)
-    assert f.integration_qags(xmin, xmax, [0.0, 1e-7], w)
-    assert f.integration_qags([xmin, xmax], [0.0, 1e-7], w)
-
-    assert f.integration_qags([xmin, xmax], limit)
-    assert f.integration_qags(xmin, xmax, limit)
-    assert f.integration_qags([xmin, xmax], w)
-    assert f.integration_qags(xmin, xmax, w)
-    assert f.integration_qags(xmin, xmax, limit, w)
-    assert f.integration_qags([xmin, xmax], limit, w)
-
-    # QAGP
-    assert f.integration_qagp([xmin, xmax])
-    assert f.integration_qagp([xmin, xmax], [0.0, 1e-7])
-    assert f.integration_qagp([xmin, xmax], 0.0, 1e-7)
-    assert f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit)
-    assert f.integration_qagp([xmin, xmax], [0.0, 1e-7], limit)
-    assert f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit, w)
-    assert f.integration_qagp([xmin, xmax], [0.0, 1e-7], w)
-
-    assert f.integration_qagp([xmin, xmax], limit)
-    assert f.integration_qagp([xmin, xmax], w)
-    assert f.integration_qagp([xmin, xmax], limit, w)
-  end
-
-  def test_integration3
-    f = GSL::Function.alloc { |x| Math.exp(-x) / Math.sqrt(x) }
-
-    xmin = 0.0
-    xmax = 1.0
-    limit = 1000
-
-    # QNG
-    # XXX GSL::ERROR::ETOL: Ruby/GSL error code 14, failed to reach tolerance with
-    # highest-order rule (file qng.c, line 189), failed to reach the specified tolerance
-    #assert f.integration_qng(xmin, xmax, 0.0, 1.0e-7)
-    #assert f.integration_qng(xmin, xmax)
-    #assert f.integration_qng([xmin, xmax])
-    #assert f.integration_qng([xmin, xmax], [0.0, 1.0e-7])
-    #assert f.integration_qng([xmin, xmax], 0.0, 1.0e-7)
-    #assert f.integration_qng(xmin, xmax, [0.0, 1.0e-7])
-
-    # QAG
-    assert f.integration_qag(xmin, xmax, 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag([xmin, xmax], 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag(xmin, xmax, GSL::Integration::GAUSS15)
-    assert f.integration_qag([xmin, xmax], GSL::Integration::GAUSS15)
-
-    w = GSL::Integration::Workspace.alloc(2000)
-    assert f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
-    assert f.integration_qag(xmin, xmax, 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
-    assert f.integration_qag([xmin, xmax], 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
-    assert f.integration_qag([xmin, xmax], [0.0, 1.0e-7], GSL::Integration::GAUSS15, w)
-    assert f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
-
-    # QAGS
-    assert f.integration_qags(xmin, xmax)
-    assert f.integration_qags([xmin, xmax])
-    assert f.integration_qags(xmin, xmax, 0.0, 1e-7)
-    assert f.integration_qags(xmin, xmax, [0.0, 1e-7])
-    assert f.integration_qags([xmin, xmax], [0.0, 1e-7])
-    assert f.integration_qags([xmin, xmax], 0.0, 1e-7)
-    assert f.integration_qags([xmin, xmax], 0.0, 1e-7, limit)
-    assert f.integration_qags([xmin, xmax], [0.0, 1e-7], limit)
-    assert f.integration_qags(xmin, xmax, [0.0, 1e-7], limit)
-    assert f.integration_qags([xmin, xmax], 0.0, 1e-7, limit, w)
-    assert f.integration_qags(xmin, xmax, 0.0, 1e-7, limit, w)
-    assert f.integration_qags(xmin, xmax, [0.0, 1e-7], w)
-    assert f.integration_qags([xmin, xmax], [0.0, 1e-7], w)
-
-    assert f.integration_qags([xmin, xmax], limit)
-    assert f.integration_qags(xmin, xmax, limit)
-    assert f.integration_qags([xmin, xmax], w)
-    assert f.integration_qags(xmin, xmax, w)
-    assert f.integration_qags(xmin, xmax, limit, w)
-    assert f.integration_qags([xmin, xmax], limit, w)
-
-    # QAGP
-    assert f.integration_qagp([xmin, xmax])
-    assert f.integration_qagp([xmin, xmax], [0.0, 1e-7])
-    assert f.integration_qagp([xmin, xmax], 0.0, 1e-7)
-    assert f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit)
-    assert f.integration_qagp([xmin, xmax], [0.0, 1e-7], limit)
-    assert f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit, w)
-    assert f.integration_qagp([xmin, xmax], [0.0, 1e-7], w)
-
-    assert f.integration_qagp([xmin, xmax], limit)
-    assert f.integration_qagp([xmin, xmax], w)
-    assert f.integration_qagp([xmin, xmax], limit, w)
-  end
-
-  def test_integration4
-    f = GSL::Function.alloc { |x| (1.0 - Math.exp(-x)) / Math.sqrt(x) }
-
-    xmin = 0.0
-    xmax = 0.2
-    limit = 1000
-
-    # QNG
-    # XXX GSL::ERROR::ETOL: Ruby/GSL error code 14, failed to reach tolerance with
-    # highest-order rule (file qng.c, line 189), failed to reach the specified tolerance
-    #assert f.integration_qng(xmin, xmax, 0.0, 1.0e-7)
-    #assert f.integration_qng(xmin, xmax)
-    #assert f.integration_qng([xmin, xmax])
-    #assert f.integration_qng([xmin, xmax], [0.0, 1.0e-7])
-    #assert f.integration_qng([xmin, xmax], 0.0, 1.0e-7)
-    #assert f.integration_qng(xmin, xmax, [0.0, 1.0e-7])
-
-    # QAG
-    assert f.integration_qag(xmin, xmax, 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag([xmin, xmax], 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
-    assert f.integration_qag(xmin, xmax, GSL::Integration::GAUSS15)
-    assert f.integration_qag([xmin, xmax], GSL::Integration::GAUSS15)
-
-    w = GSL::Integration::Workspace.alloc(2000)
-    assert f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
-    assert f.integration_qag(xmin, xmax, 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
-    assert f.integration_qag([xmin, xmax], 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
-    assert f.integration_qag([xmin, xmax], [0.0, 1.0e-7], GSL::Integration::GAUSS15, w)
-    assert f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
-
-    # QAGS
-    assert f.integration_qags(xmin, xmax)
-    assert f.integration_qags([xmin, xmax])
-    assert f.integration_qags(xmin, xmax, 0.0, 1e-7)
-    assert f.integration_qags(xmin, xmax, [0.0, 1e-7])
-    assert f.integration_qags([xmin, xmax], [0.0, 1e-7])
-    assert f.integration_qags([xmin, xmax], 0.0, 1e-7)
-    assert f.integration_qags([xmin, xmax], 0.0, 1e-7, limit)
-    assert f.integration_qags([xmin, xmax], [0.0, 1e-7], limit)
-    assert f.integration_qags(xmin, xmax, [0.0, 1e-7], limit)
-    assert f.integration_qags([xmin, xmax], 0.0, 1e-7, limit, w)
-    assert f.integration_qags(xmin, xmax, 0.0, 1e-7, limit, w)
-    assert f.integration_qags(xmin, xmax, [0.0, 1e-7], w)
-    assert f.integration_qags([xmin, xmax], [0.0, 1e-7], w)
-
-    assert f.integration_qags([xmin, xmax], limit)
-    assert f.integration_qags(xmin, xmax, limit)
-    assert f.integration_qags([xmin, xmax], w)
-    assert f.integration_qags(xmin, xmax, w)
-    assert f.integration_qags(xmin, xmax, limit, w)
-    assert f.integration_qags([xmin, xmax], limit, w)
-
-    # QAGP
-    assert f.integration_qagp([xmin, xmax])
-    assert f.integration_qagp([xmin, xmax], [0.0, 1e-7])
-    assert f.integration_qagp([xmin, xmax], 0.0, 1e-7)
-    assert f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit)
-    assert f.integration_qagp([xmin, xmax], [0.0, 1e-7], limit)
-    assert f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit, w)
-    assert f.integration_qagp([xmin, xmax], [0.0, 1e-7], w)
-
-    assert f.integration_qagp([xmin, xmax], limit)
-    assert f.integration_qagp([xmin, xmax], w)
-    assert f.integration_qagp([xmin, xmax], limit, w)
-  end
-
-end