gsl 1.12.109 → 1.14.5

data/tests/gsl_test2.rb

@@ -0,0 +1,107 @@
+#!/usr/bin/env ruby
+require("gsl")
+include GSL
+
+module GSL
+  module Test
+    Verbose = true
+    $tests = 0
+    $passed = 0
+    $failed = 0
+
+# PASS if status == 0, FAIL otherwise
+    def test(status, desc)
+      $tests += 1
+      if status == 0
+        $passed += 1
+        printf("PASS: #{desc}\n")
+      else
+        $failed += 1
+        printf("FAIL: #{desc}\n")
+      end
+    end
+
+# PASS if status == true, FAIL otherwise
+    def test2(status, desc)
+      if status == true
+        printf("PASS: #{desc}\n")
+      else
+        printf("FAIL: #{desc}\n")
+      end
+    end
+
+    def test_factor(result, expected, factor, desc)
+      status = nil
+      if result == expected
+        status = 0
+      elsif expected == 0.0
+        status = (result > expected or result < expected) ? 1 : 0
+      else
+        u = result/expected
+        status = (u > factor or u < 1.0/factor) ? 1 : 0
+      end
+      $tests += 1
+      if status == 0
+        $passed += 1
+        printf("PASS: #{desc} (%g observed vs %g expected)\n", result, expected) 
+      else
+        $failed += 1
+        printf("FAIL: #{desc} (%.18g observed vs %.18g expected)\n", result, expected) 
+      end
+    end
+
+    def test_rel(result, expected, relerr, desc)
+      status = nil
+      if isnan?(result) or isnan?(expected)
+        status = isnan?(result) != isnan?(expected) ? 1 : 0
+      elsif isinf?(result) or isinf?(expected)
+        status = isinf?(result) != isinf?(expected) ? 1 : 0
+      elsif expected.to_f != 0.0
+        status = (result - expected).abs/expected.abs > relerr ? 1 : 0
+      else
+        status = result.abs > relerr ? 1 : 0
+      end
+      $tests += 1
+      if status == 0
+        $passed += 1
+        printf("PASS: #{desc} (%g observed vs %g expected)\n", result, expected) 
+      else
+        $failed += 1
+        printf("FAIL: #{desc} (%.18g observed vs %.18g expected)\n", result, expected) 
+      end
+    end
+
+    def test_abs(result, expected, abserr, desc)
+      status = nil
+      if isnan?(result) or isnan?(expected)
+        status = isnan?(result) != isnan?(expected) ? 1 : 0
+      elsif isinf?(result) or isinf?(expected)
+        status = isinf?(result) != isinf?(expected) ? 1 : 0
+      else
+        status = (result - expected).abs > abserr ? 1 : 0
+      end
+      $tests += 1
+      if status == 0
+        $passed += 1
+        printf("PASS: #{desc} (%g observed vs %g expected)\n", result, expected) 
+      else
+        $failed += 1
+        printf("FAIL: #{desc} (%.18g observed vs %.18g expected)\n", result, expected) 
+      end      
+    end
+
+    def test_int(result, expected, desc)
+      status = (result != expected) ? 1 : 0
+      $tests += 1
+      if status == 0
+        $passed += 1
+        printf("PASS: #{desc} (%d observed vs %d expected)\n", result, expected) 
+      else
+        $failed += 1
+        printf("FAIL: #{desc} (%d observed vs %d expected)\n", result, expected) 
+      end
+    end
+
+  end
+end
+

data/tests/histo.rb

@@ -0,0 +1,12 @@
+#!/usr/bin/env ruby
+require("gsl")
+
+h = GSL::Histogram.alloc(10, [0, 10])
+p h
+
+p h.get_range(2)
+
+p h.range
+p h.bin
+
+

data/tests/integration/integration1.rb

@@ -0,0 +1,72 @@
+#!/usr/bin/env ruby
+require("gsl")
+include Math
+
+f = GSL::Function::alloc{ |x|
+  exp(x)*cos(x)
+}
+
+xmin = 0.0
+xmax = 1.0
+limit = 1000
+
+puts("QNG")
+p f.integration_qng(xmin, xmax, 0.0, 1.0e-7)
+p f.integration_qng(xmin, xmax)
+p f.integration_qng([xmin, xmax])
+p f.integration_qng([xmin, xmax], [0.0, 1.0e-7])
+p f.integration_qng([xmin, xmax], 0.0, 1.0e-7)
+p f.integration_qng(xmin, xmax, [0.0, 1.0e-7])
+
+puts("QAG")
+p f.integration_qag(xmin, xmax, 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
+p f.integration_qag([xmin, xmax], 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
+p f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
+p f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
+p f.integration_qag(xmin, xmax, GSL::Integration::GAUSS15)
+p f.integration_qag([xmin, xmax], GSL::Integration::GAUSS15)
+w = GSL::Integration::Workspace.alloc
+p f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
+p f.integration_qag(xmin, xmax, 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
+p f.integration_qag([xmin, xmax], 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
+p f.integration_qag([xmin, xmax], [0.0, 1.0e-7], GSL::Integration::GAUSS15, w)
+p f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
+p f.integration_qag([xmin, xmax], w)
+
+puts("QAGS")
+p f.integration_qags(xmin, xmax)
+p f.integration_qags([xmin, xmax])
+p f.integration_qags(xmin, xmax, 0.0, 1e-7)
+p f.integration_qags(xmin, xmax, [0.0, 1e-7])
+p f.integration_qags([xmin, xmax], [0.0, 1e-7])
+p f.integration_qags([xmin, xmax], 0.0, 1e-7)
+p f.integration_qags([xmin, xmax], 0.0, 1e-7, limit)
+p f.integration_qags([xmin, xmax], [0.0, 1e-7], limit)
+p f.integration_qags(xmin, xmax, [0.0, 1e-7], limit)
+p f.integration_qags([xmin, xmax], 0.0, 1e-7, limit, w)
+p f.integration_qags(xmin, xmax, 0.0, 1e-7, limit, w)
+p f.integration_qags(xmin, xmax, [0.0, 1e-7], w)
+p f.integration_qags([xmin, xmax], [0.0, 1e-7], w)
+
+p f.integration_qags([xmin, xmax], limit)
+p f.integration_qags(xmin, xmax, limit)
+p f.integration_qags([xmin, xmax], w)
+p f.integration_qags(xmin, xmax, w)
+p f.integration_qags(xmin, xmax, limit, w)
+p f.integration_qags([xmin, xmax], limit, w)
+
+puts("QAGP")
+p f.integration_qagp([xmin, xmax])
+p f.integration_qagp([xmin, xmax], [0.0, 1e-7])
+p f.integration_qagp([xmin, xmax], 0.0, 1e-7)
+p f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit)
+p f.integration_qagp([xmin, xmax], [0.0, 1e-7], limit)
+p f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit, w)
+p f.integration_qagp([xmin, xmax], [0.0, 1e-7], w)
+
+p f.integration_qagp([xmin, xmax], limit)
+p f.integration_qagp([xmin, xmax], w)
+p f.integration_qagp([xmin, xmax], limit, w)
+
+
+__END__

data/tests/integration/integration2.rb

@@ -0,0 +1,71 @@
+#!/usr/bin/env ruby
+require("gsl")
+include Math
+
+f = GSL::Function::alloc{ |x|
+  sin(x)/x
+}
+
+xmin = 0.0
+xmax = 2.0*PI
+limit = 1000
+
+puts("QNG")
+p f.integration_qng(xmin, xmax, 0.0, 1.0e-7)
+p f.integration_qng(xmin, xmax)
+p f.integration_qng([xmin, xmax])
+p f.integration_qng([xmin, xmax], [0.0, 1.0e-7])
+p f.integration_qng([xmin, xmax], 0.0, 1.0e-7)
+p f.integration_qng(xmin, xmax, [0.0, 1.0e-7])
+
+puts("QAG")
+p f.integration_qag(xmin, xmax, 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
+p f.integration_qag([xmin, xmax], 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
+p f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
+p f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
+p f.integration_qag(xmin, xmax, GSL::Integration::GAUSS15)
+p f.integration_qag([xmin, xmax], GSL::Integration::GAUSS15)
+w = GSL::Integration::Workspace.alloc(2000)
+p f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
+p f.integration_qag(xmin, xmax, 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
+p f.integration_qag([xmin, xmax], 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
+p f.integration_qag([xmin, xmax], [0.0, 1.0e-7], GSL::Integration::GAUSS15, w)
+p f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
+
+puts("QAGS")
+p f.integration_qags(xmin, xmax)
+p f.integration_qags([xmin, xmax])
+p f.integration_qags(xmin, xmax, 0.0, 1e-7)
+p f.integration_qags(xmin, xmax, [0.0, 1e-7])
+p f.integration_qags([xmin, xmax], [0.0, 1e-7])
+p f.integration_qags([xmin, xmax], 0.0, 1e-7)
+p f.integration_qags([xmin, xmax], 0.0, 1e-7, limit)
+p f.integration_qags([xmin, xmax], [0.0, 1e-7], limit)
+p f.integration_qags(xmin, xmax, [0.0, 1e-7], limit)
+p f.integration_qags([xmin, xmax], 0.0, 1e-7, limit, w)
+p f.integration_qags(xmin, xmax, 0.0, 1e-7, limit, w)
+p f.integration_qags(xmin, xmax, [0.0, 1e-7], w)
+p f.integration_qags([xmin, xmax], [0.0, 1e-7], w)
+
+p f.integration_qags([xmin, xmax], limit)
+p f.integration_qags(xmin, xmax, limit)
+p f.integration_qags([xmin, xmax], w)
+p f.integration_qags(xmin, xmax, w)
+p f.integration_qags(xmin, xmax, limit, w)
+p f.integration_qags([xmin, xmax], limit, w)
+
+puts("QAGP")
+p f.integration_qagp([xmin, xmax])
+p f.integration_qagp([xmin, xmax], [0.0, 1e-7])
+p f.integration_qagp([xmin, xmax], 0.0, 1e-7)
+p f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit)
+p f.integration_qagp([xmin, xmax], [0.0, 1e-7], limit)
+p f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit, w)
+p f.integration_qagp([xmin, xmax], [0.0, 1e-7], w)
+
+p f.integration_qagp([xmin, xmax], limit)
+p f.integration_qagp([xmin, xmax], w)
+p f.integration_qagp([xmin, xmax], limit, w)
+
+
+__END__

data/tests/integration/integration3.rb

@@ -0,0 +1,71 @@
+#!/usr/bin/env ruby
+require("gsl")
+include Math
+
+f = GSL::Function::alloc{ |x|
+  exp(-x)/sqrt(x)
+}
+
+xmin = 0.0
+xmax = 1.0
+limit = 1000
+
+=begin
+puts("QNG")
+p f.integration_qng(xmin, xmax, 0.0, 1.0e-7)
+p f.integration_qng(xmin, xmax)
+p f.integration_qng([xmin, xmax])
+p f.integration_qng([xmin, xmax], [0.0, 1.0e-7])
+p f.integration_qng([xmin, xmax], 0.0, 1.0e-7)
+p f.integration_qng(xmin, xmax, [0.0, 1.0e-7])
+=end
+
+puts("QAG")
+p f.integration_qag(xmin, xmax, 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
+p f.integration_qag([xmin, xmax], 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
+p f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
+p f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
+p f.integration_qag(xmin, xmax, GSL::Integration::GAUSS15)
+p f.integration_qag([xmin, xmax], GSL::Integration::GAUSS15)
+w = GSL::Integration::Workspace.alloc(2000)
+p f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
+p f.integration_qag(xmin, xmax, 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
+p f.integration_qag([xmin, xmax], 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
+p f.integration_qag([xmin, xmax], [0.0, 1.0e-7], GSL::Integration::GAUSS15, w)
+p f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
+
+puts("QAGS")
+p f.integration_qags(xmin, xmax)
+p f.integration_qags([xmin, xmax])
+p f.integration_qags(xmin, xmax, 0.0, 1e-7)
+p f.integration_qags(xmin, xmax, [0.0, 1e-7])
+p f.integration_qags([xmin, xmax], [0.0, 1e-7])
+p f.integration_qags([xmin, xmax], 0.0, 1e-7)
+p f.integration_qags([xmin, xmax], 0.0, 1e-7, limit)
+p f.integration_qags([xmin, xmax], [0.0, 1e-7], limit)
+p f.integration_qags(xmin, xmax, [0.0, 1e-7], limit)
+p f.integration_qags([xmin, xmax], 0.0, 1e-7, limit, w)
+p f.integration_qags(xmin, xmax, 0.0, 1e-7, limit, w)
+p f.integration_qags(xmin, xmax, [0.0, 1e-7], w)
+p f.integration_qags([xmin, xmax], [0.0, 1e-7], w)
+
+p f.integration_qags([xmin, xmax], limit)
+p f.integration_qags(xmin, xmax, limit)
+p f.integration_qags([xmin, xmax], w)
+p f.integration_qags(xmin, xmax, w)
+p f.integration_qags(xmin, xmax, limit, w)
+p f.integration_qags([xmin, xmax], limit, w)
+
+puts("QAGP")
+p f.integration_qagp([xmin, xmax])
+p f.integration_qagp([xmin, xmax], [0.0, 1e-7])
+p f.integration_qagp([xmin, xmax], 0.0, 1e-7)
+p f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit)
+p f.integration_qagp([xmin, xmax], [0.0, 1e-7], limit)
+p f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit, w)
+p f.integration_qagp([xmin, xmax], [0.0, 1e-7], w)
+
+p f.integration_qagp([xmin, xmax], limit)
+p f.integration_qagp([xmin, xmax], w)
+p f.integration_qagp([xmin, xmax], limit, w)
+

data/tests/integration/integration4.rb

@@ -0,0 +1,71 @@
+#!/usr/bin/env ruby
+require("gsl")
+include Math
+
+f = GSL::Function::alloc{ |x|
+  (1.0 - exp(-x))/sqrt(x)
+}
+
+xmin = 0.0
+xmax = 0.2
+limit = 1000
+
+=begin
+puts("QNG")
+p f.integration_qng(xmin, xmax, 0.0, 1.0e-7)
+p f.integration_qng(xmin, xmax)
+p f.integration_qng([xmin, xmax])
+p f.integration_qng([xmin, xmax], [0.0, 1.0e-7])
+p f.integration_qng([xmin, xmax], 0.0, 1.0e-7)
+p f.integration_qng(xmin, xmax, [0.0, 1.0e-7])
+=end
+
+puts("QAG")
+p f.integration_qag(xmin, xmax, 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
+p f.integration_qag([xmin, xmax], 0.0, 1.0e-7, limit, GSL::Integration::GAUSS15)
+p f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
+p f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15)
+p f.integration_qag(xmin, xmax, GSL::Integration::GAUSS15)
+p f.integration_qag([xmin, xmax], GSL::Integration::GAUSS15)
+w = GSL::Integration::Workspace.alloc(2000)
+p f.integration_qag(xmin, xmax, [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
+p f.integration_qag(xmin, xmax, 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
+p f.integration_qag([xmin, xmax], 0.0, 1.0e-7, GSL::Integration::GAUSS15, w)
+p f.integration_qag([xmin, xmax], [0.0, 1.0e-7], GSL::Integration::GAUSS15, w)
+p f.integration_qag([xmin, xmax], [0.0, 1.0e-7], limit, GSL::Integration::GAUSS15, w)
+
+puts("QAGS")
+p f.integration_qags(xmin, xmax)
+p f.integration_qags([xmin, xmax])
+p f.integration_qags(xmin, xmax, 0.0, 1e-7)
+p f.integration_qags(xmin, xmax, [0.0, 1e-7])
+p f.integration_qags([xmin, xmax], [0.0, 1e-7])
+p f.integration_qags([xmin, xmax], 0.0, 1e-7)
+p f.integration_qags([xmin, xmax], 0.0, 1e-7, limit)
+p f.integration_qags([xmin, xmax], [0.0, 1e-7], limit)
+p f.integration_qags(xmin, xmax, [0.0, 1e-7], limit)
+p f.integration_qags([xmin, xmax], 0.0, 1e-7, limit, w)
+p f.integration_qags(xmin, xmax, 0.0, 1e-7, limit, w)
+p f.integration_qags(xmin, xmax, [0.0, 1e-7], w)
+p f.integration_qags([xmin, xmax], [0.0, 1e-7], w)
+
+p f.integration_qags([xmin, xmax], limit)
+p f.integration_qags(xmin, xmax, limit)
+p f.integration_qags([xmin, xmax], w)
+p f.integration_qags(xmin, xmax, w)
+p f.integration_qags(xmin, xmax, limit, w)
+p f.integration_qags([xmin, xmax], limit, w)
+
+puts("QAGP")
+p f.integration_qagp([xmin, xmax])
+p f.integration_qagp([xmin, xmax], [0.0, 1e-7])
+p f.integration_qagp([xmin, xmax], 0.0, 1e-7)
+p f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit)
+p f.integration_qagp([xmin, xmax], [0.0, 1e-7], limit)
+p f.integration_qagp([xmin, xmax], 0.0, 1e-7, limit, w)
+p f.integration_qagp([xmin, xmax], [0.0, 1e-7], w)
+
+p f.integration_qagp([xmin, xmax], limit)
+p f.integration_qagp([xmin, xmax], w)
+p f.integration_qagp([xmin, xmax], limit, w)
+

data/tests/interp.rb

@@ -0,0 +1,45 @@
+#!/usr/bin/env ruby
+# Ruby/GSL implementation of GSL "interpolation/test.c"
+require("gsl")
+require("./gsl_test2.rb")
+include GSL::Test
+include Math
+
+def test_bsearch()
+  status = 0
+  x_array = GSL::Vector.alloc(0.0, 1.0, 2.0, 3.0, 4.0)
+  index_result = Interp.bsearch(x_array, 1.5, 0, 4)
+  s = (index_result != 1) ? 1 : 0
+  status += 1
+  test(s, "simple bsearch")
+
+  index_result = x_array.bsearch(4.0, 0, 4)
+  s = (index_result != 3) ? 1 : 0
+  status += s;
+  test(s, "upper endpoint bsearch");
+
+  index_result = Interp.bsearch(x_array, 0.0, 0, 4);
+  s = (index_result != 0) ? 1 : 0
+  status += s;
+  test(s, "lower endpoint bsearch");
+
+  index_result = Interp.bsearch(x_array, 2.0, 0, 4);
+  s = (index_result != 2) ? 1 : 0
+  status += s;
+  test(s, "degenerate bsearch");
+
+  index_result = Interp.bsearch(x_array, 10.0, 0, 4);
+  s = (index_result != 3) ? 1 : 0
+  status += s;
+  test(s, "out of bounds bsearch +");
+
+  index_result = Interp.bsearch(x_array, -10.0, 0, 4);
+  s = (index_result != 0) ? 1 : 0
+  status += s;
+  test(s, "out of bounds bsearch -");
+
+end
+
+
+
+test_bsearch()