gsl 1.15.3 → 1.16.0.6

data/Gemfile

@@ -0,0 +1,4 @@
+source 'https://rubygems.org'
+gemspec
+
+gem 'narray' if ENV['NARRAY']

data/README.md

@@ -0,0 +1,77 @@
+# Ruby/GSL, a Ruby interface to GSL (GNU Scientific library)
+
+Permission is granted to copy, distribute and/or modify this document under
+the terms of the GNU Free Documentation License.
+
+## Description
+
+Ruby/GSL is a Ruby interface to the [GNU Scientific Library](https://gnu.org/software/gsl/)
+(GSL), for numerical computing with [Ruby](http://www.ruby-lang.org/).
+
+
+## Installation
+
+Ruby/GSL may be installed as a Ruby Gem by simply running
+
+  gem install rb-gsl
+
+Note that the GSL libraries must already be installed before Ruby/GSL
+can be installed:
+
+Debian/Ubuntu: +libgsl0-dev+
+Fedora/SuSE:   +gsl-devel+
+Gentoo:        +sci-libs/gsl+
+OS X:          <tt>brew install gsl</tt>
+
+It is recommended to install the [GNU plotutils](https://gnu.org/software/plotutils/plotutils.html)
+package. Some of the example scripts in the +examples/+ directory use the
++graph+ utility included in the package to plot the results. Windows cygwin
+binaries of <tt>GNU plotutils</tt> and related packages are available
+[here](http://gnuwin32.sourceforge.net/packages/plotutils.htm).
+
+
+## Reference
+
+The [Ruby/GSL reference manual](link:rdoc/ref_rdoc.html) follows and borrows
+large parts of the GSL reference manual.
+
+
+## Examples
+
+See scripts in +examples/+ and +test/+ directories.
+
+
+## Related Projects
+
+* [ruby-gsl](http://ruby-gsl.sourceforge.net/): Another Ruby binding,
+  developed by Arno Erpenbeck.
+
+
+## Licence
+
+Ruby/GSL is free software: you can redistribute it and/or modify it
+under the terms of the GNU General Public License.
+This library is distributed in the hope that it will be useful, but
+WITHOUT ANY WARRANTY.
+
+
+## Bug Reports
+
+Any bug reports are welcome. If you encounter bugs in Ruby/GSL, please
+report them on GitHub(https://github.com/SciRuby/rb-gsl/issues).
+
+
+## Links
+
+Documentation: https://sciruby.github.com/rb-gsl
+Source code:   https://github.com/SciRuby/rb-gsl
+RubyGem:       https://rubygems.org/gems/rb-gsl
+Bug tracker:   https://github.com/SciRuby/rb-gsl/issues
+Travis CI:     https://travis-ci.org/SciRuby/rb-gsl
+
+
+## Authors
+
+* Yoshiki Tsunesada <y-tsunesada at mm dot em-net dot ne dot jp> (July, 2004)
+* David MacMahon <davidm@astro.berkeley.edu> (November, 2010)
+* Jens Wille <mailto:jens.wille@gmail.com> (November, 2013)

data/Rakefile

@@ -1,99 +1,20 @@
-require 'rubygems'
-require 'rubygems/package_task'
-require 'rake/rdoctask'
-
-RB_GSL_VERSION = File.readlines('VERSION')[0].chomp
-
-spec = Gem::Specification.new do |s|
-  # Basics
-  s.name = 'gsl'
-  s.version = RB_GSL_VERSION
-  s.summary = 'Ruby interface to GNU Scientific Library'
-  s.description = 'Ruby/GSL is a Ruby interface to the GNU Scientific Library, for numerical computing with Ruby'
-  #s.platform = Gem::Platform::Ruby
-  s.required_ruby_version = '>= 1.8.1'
-  s.requirements << 'GSL (http://www.gnu.org/software/gsl/)'
-  # plotlib?
-  s.add_dependency('narray', '>= 0.5.9')
-
-  # About
-  s.authors = ['Yoshiki Tsunesada', 'David MacMahon']
-  s.email = 'y-tsunesada@mm.em-net.ne.jp'
-  s.homepage = 'http://rb-gsl.rubyforge.org/'
-  s.rubyforge_project = 'rb-gsl' 
-
-  # Files, Libraries, and Extensions
-  s.files = FileList[
-    'AUTHORS',
-    'COPYING',
-    'ChangeLog',
-    'README',
-    'Rakefile',
-    'setup.rb',
-    'THANKS',
-    'VERSION',
-    'examples/**/*',
-    'ext/extconf.rb',
-    'ext/*.c',
-    'lib/**/*',
-    'include/*',
-    'rdoc/*',
-    'tests/**/*'
-  ].to_a
-  s.require_paths = ['lib', 'lib/gsl', 'lib/ool', 'ext']
-  #s.autorequire = nil
-  #s.bindir = 'bin'
-  #s.executables = []
-  #s.default_executable = nil
-
-  # C compilation
-  s.extensions = %w[ ext/extconf.rb ]
-
-  # Documentation
-  s.has_rdoc = true
-  s.rdoc_options = [
-    '--title', 'Ruby/GSL',
-    '--main', 'rdoc/index.rdoc',
-    '--exclude', 'ext/',
-    '--exclude', 'include/',
-    '--exclude', 'lib/',
-  ]
-  s.extra_rdoc_files = FileList['rdoc/*'].to_a
+require File.expand_path(%q{../lib/gsl/version}, __FILE__)
 
-  # Testing TODO
-  #s.test_files = []
-end
+require 'bundler/setup'
+require 'rubygems/package_task'
+require 'rake/extensiontask'
+require 'rake/testtask'
 
-Rake::PackageTask.new('rb-gsl', RB_GSL_VERSION) do |pkg|
-  pkg.need_zip = true
-  pkg.need_tar = true
-  pkg.package_files = spec.files
-end
+Bundler::GemHelper.install_tasks
 
-Gem::PackageTask.new(spec) do |pkg|
-  pkg.need_zip = false
-  pkg.need_tar = false
+Rake::TestTask.new do |t|
+  t.libs << 'test'
+  t.libs << 'test/gsl'
+  t.test_files = FileList['test/*.rb', 'test/gsl/*.rb']
 end
 
-task :default => [:package, :gem]
+spec = eval(IO.read('gsl.gemspec'))
+Gem::PackageTask.new(spec).define
+Rake::ExtensionTask.new(:gsl_native, spec)
 
-# --------------------------------------------------------------------
-# Create a task to build the RDOC documentation tree.
-
-desc "Create the RDoc html files"
-Rake::RDocTask.new("rdoc") { |rdoc|
-  rdoc.rdoc_dir = 'html'
-  rdoc.title    = 'Ruby/GSL'
-  rdoc.main     = 'rdoc/index.rdoc'
-  rdoc.options << '--exclude' << 'ext/'
-  rdoc.options << '--exclude' << 'include/'
-  rdoc.options << '--exclude' << 'lib/'
-  rdoc.rdoc_files.include('rdoc/*.rdoc')
-}
-
-desc "Publish the RDoc files on RubyForge"
-task :pub_rdoc => ["html/index.html"] do
-  mkdir_p "emptydir"
-  sh "scp -rq html/* www.rubyforge.org:/var/www/gforge-projects/rb-gsl/."
-  rm_r "emptydir"
-end
+task :default => [:compile, :test]

data/THANKS

@@ -1,5 +1,5 @@
 Takeshi Horinouchi
-Masahide Kikkawa 
+Masahide Kikkawa
 Seiya Nishizawa (Patches, NArray supports)
 Arno Erpenbeck (ruby-gsl project)
 Austin Newton

data/examples/blas/blas.rb

@@ -10,4 +10,4 @@ B = b.matrix_view(3, 2)
 C = GSL::Blas.dgemm(GSL::Blas::NoTrans, GSL::Blas::NoTrans, 1.0, A, B, 0.0)
 p C
 
-p A*B 
+p A*B

data/examples/bspline.rb

@@ -19,15 +19,15 @@ w = GSL::Vector.alloc(N)
 for i in 0...N do
   xi = (15.0/(N-1)/1)*i
   yi = Math::cos(xi)*Math::exp(-0.1*xi)
-  
+
   sigma = 0.1
   dy = GSL::Ran.gaussian(r, sigma)
   yi += dy
-  
+
   x[i] = xi
   y[i] = yi
   w[i] = sigma
-  
+
 #  printf("%f %f\n", xi, yi)
 end
 

data/examples/complex/functions.rb

@@ -26,8 +26,8 @@ p f.log == z
 p z.pow(f)  # [1.931e-02 1.752e-02], verified with Octave result
 
 # The results below are verified with the Octave results
-p GSL::Complex.sin(z) 
-p GSL::Complex.cos(z) 
+p GSL::Complex.sin(z)
+p GSL::Complex.cos(z)
 p GSL::Complex.tan(z)
 p GSL::Complex.sec(z)
 p GSL::Complex.csc(z)
@@ -43,8 +43,8 @@ p GSL::Complex.arccsc(z)
 p GSL::Complex.arccsc_real(2)
 p GSL::Complex.arccot(z)
 
-p GSL::Complex.sinh(z) 
-p GSL::Complex.cosh(z) 
+p GSL::Complex.sinh(z)
+p GSL::Complex.cosh(z)
 p GSL::Complex.tanh(z)
 p GSL::Complex.sech(z)
 p GSL::Complex.csch(z)

data/examples/complex/mul.rb

@@ -3,7 +3,7 @@ require 'gsl'
 
 a = GSL::Complex.rect(1, 2)
 b = GSL::Complex[3, 4]
- 
+
 p a*b  # -5 + 10i
 p a.mul(b)
 p a

data/examples/const/physconst.rb

@@ -42,7 +42,7 @@ MENU[2] = ["", "ASTRONOMICAL_UNIT", "GRAVITATIONAL_CONSTANT",
 MENU[3] = ["", "ELECTRON_CHARGE", "ELECTRON_VOLT", "UNIFIED_ATOMIC_MASS",
            "MASS_ELECTRON", "MASS_MUON", "MASS_PROTON", "MASS_NEUTRON",
            "FINE_STRUCTURE", "RYDBERG", "BOHR_RADIUS", "ANGSTROM",
-           "BARN", "BOHR_MAGNETON", "NUCLEAR_MAGNETON", 
+           "BARN", "BOHR_MAGNETON", "NUCLEAR_MAGNETON",
            "ELECTRON_MAGNETIC_MOMENT", "PROTON_MAGNETIC_MOMENT",
            "THOMSON_CROSS_SECTION"]
 
@@ -56,10 +56,10 @@ MENU[9] = ["", "POUND_MASS", "OUNCE_MASS", "TON", "METRIC_TON", "UK_TON",
            "TROY_OUNCE", "CARAT", "GRAM_FORCE", "POUND_FORCE",
            "KILOPOUND_FORCE", "POUNDAL"]
 MENU[10] = ["", "CALORIE", "BTU", "THERM", "HORSEPOWER"]
-MENU[11] = ["", "BAR", "STD_ATMOSPHERE", "TORR", "METER_OF_MERCURY", 
+MENU[11] = ["", "BAR", "STD_ATMOSPHERE", "TORR", "METER_OF_MERCURY",
             "INCH_OF_MERCURY", "INCH_OF_WATER", "PSI"]
 MENU[12] = ["", "POISE", "STOKES"]
-MENU[13] = ["", "STILB", "LUMEN", "LUX", "PHOT", "FOOTCANDLE", 
+MENU[13] = ["", "STILB", "LUMEN", "LUX", "PHOT", "FOOTCANDLE",
             "LAMBERT", "FOOTLAMBERT"]
 MENU[14] = ["", "CURIE", "ROENTGEN", "RAD"]
 MENU[15] = ["", "NEWTON", "DYNE", "JOULE", "ERG"]
@@ -96,7 +96,7 @@ def show_submenu(ind, w)
   w.setpos(i+=1, 1)
   str = "    ? ['m' to menu]: "
   w.addstr(str)
-  w.setpos(i, str.size+1)  
+  w.setpos(i, str.size+1)
   w.box('|', '-')
   w.refresh
 end

data/examples/const/travel.rb

@@ -8,8 +8,8 @@ c  = SPEED_OF_LIGHT;
 au = ASTRONOMICAL_UNIT;
 minutes = MINUTE;
 
-# distance stored in meters 
-r_earth = 1.00 * au;  
+# distance stored in meters
+r_earth = 1.00 * au;
 r_mars  = 1.52 * au;
 
 t_min = (r_mars - r_earth) / c;
@@ -30,8 +30,8 @@ c  = SPEED_OF_LIGHT;
 au = ASTRONOMICAL_UNIT;
 minutes = MINUTE;
 
-# distance stored in meters 
-r_earth = 1.00 * au;  
+# distance stored in meters
+r_earth = 1.00 * au;
 r_mars  = 1.52 * au;
 
 t_min = (r_mars - r_earth) / c;

data/examples/eigen/nonsymmv.rb

@@ -13,7 +13,7 @@ p eval
 p evec
 
 =begin
-This can be compared with the corresponding output from GNU OCTAVE, 
+This can be compared with the corresponding output from GNU OCTAVE,
 
   octave> [v,d] = eig(vander([-1 -2 3 4]));
 octave> diag(d)

data/examples/eigen/qhoscillator.rb

@@ -61,13 +61,13 @@ end
 # The differences with Tabelle 5.1 are from NMAX and dx.
 # If we use NMAX=512 and dx=0.02, we obtain the same results (but much slower).
 STDERR.puts("  Eigen values:")
-STDERR.printf("  %2s Exact %5s %10s | %2s Exact %5s %10s\n", 
+STDERR.printf("  %2s Exact %5s %10s | %2s Exact %5s %10s\n",
               "n", "E", "err\(\%\)", "n", "E", "err\(\%\)")
 STDERR.print("  -----------------------------------------------------\n")
 for n1 in 0..6 do
   exact1 = n1 + 0.5
   exact2 = n1 + 7 + 0.5
-  STDERR.printf("  %2d %4.1f %8.5f %+7.5f | %2d %4.1f %8.5f %+7.5f\n", 
+  STDERR.printf("  %2d %4.1f %8.5f %+7.5f | %2d %4.1f %8.5f %+7.5f\n",
                 n1, exact1, eval[n1], (exact1 - eval[n1])/exact1*100,
                 n1+7, exact2, eval[n1+7], (exact2-eval[n1+7])/exact2*100)
 end
@@ -75,7 +75,7 @@ STDERR.flush
 
 # Eigen vectors of n = 0, 1, 2, 10. See p217 "Abb 5.3"
 c = Math::sqrt(1.0/dx)
-vec0 = evec.col(0).scale(c)    
+vec0 = evec.col(0).scale(c)
 vec1 = evec.col(1).scale(c)
 vec2 = evec.col(2).scale(c)
 vec10 = evec.col(10).scale(c)

data/examples/fft/radix2.rb

@@ -1,5 +1,5 @@
 #!/usr/bin/env ruby
-require("gsl") 
+require("gsl")
 
 n = 128
 data = GSL::Vector::Complex.alloc(n)

data/examples/fft/real-halfcomplex.rb

@@ -11,7 +11,7 @@ end
 
 rtable = GSL::FFT::RealWavetable.alloc(n)
 rwork = GSL::FFT::RealWorkspace.alloc(n)
-  
+
 #ffted = data.real_transform(rtable, rwork)
 #ffted = data.real_transform(rtable)
 #ffted = data.real_transform(rwork)
@@ -22,9 +22,9 @@ ffted = data.real_transform()
 for i in 11...n do
   ffted[i] = 0.0
 end
-  
+
 hctable = GSL::FFT::HalfComplexWavetable.alloc(n)
-  
+
 #data2 = ffted.halfcomplex_inverse(hctable, rwork)
 #data2 = ffted.halfcomplex_inverse()
 #data2 = ffted.ifft

data/examples/fft/real-halfcomplex2.rb

@@ -12,15 +12,15 @@ org = data.clone
 
 rtable = GSL::FFT::RealWavetable.alloc(n)
 rwork = GSL::FFT::RealWorkspace.alloc(n)
-  
+
 data.real_transform!(rtable, rwork)
 
 for i in 11...n do
   data[i] = 0.0
 end
-  
+
 hctable = GSL::FFT::HalfComplexWavetable.alloc(n)
-  
+
 #data.halfcomplex_inverse!(hctable, rwork)
 #data.halfcomplex_inverse!()
 #data.ifft!

data/examples/fit/expfit.rb

@@ -22,7 +22,7 @@ p err
 GSL::graph(x, y, y0+amp*GSL::Sf::exp(-b*x))
 
 # This will result in
-# [ 1.019e+00 5.045e+00 1.040e-01 ] 
+# [ 1.019e+00 5.045e+00 1.040e-01 ]
 # [ 3.385e-02 5.395e-02 2.826e-03 ]
 
 # GNUPLOT results:

data/examples/fit/multifit.rb

@@ -2,7 +2,7 @@
 require("gsl")
 
 GSL::Rng.env_setup()
-  
+
 r = GSL::Rng.alloc(GSL::Rng::DEFAULT)
 
 n = 19

data/examples/fit/ndlinear.rb

@@ -2,8 +2,8 @@
 require("gsl")
 
 unless GSL::MultiFit.const_defined?("Ndlinear")
-	puts("The extension library NDLINEAR is not installed.")
-	exit()
+  puts("The extension library NDLINEAR is not installed.")
+  exit()
 end
 
 N_DIM = 3
@@ -13,33 +13,33 @@ N_SUM_PHI = 9
 R_MAX = 3.0
 
 def psi_real_exact(k, l, m, r, theta, phi)
-  rr = GSL::pow(r, l)*Math::exp(-r*r)*GSL::Sf::laguerre_n(k, l + 0.5, 2 * r * r)	 
+  rr = GSL::pow(r, l)*Math::exp(-r*r)*GSL::Sf::laguerre_n(k, l + 0.5, 2 * r * r)
 
-	tt = GSL::Sf::legendre_sphPlm(l, m, Math::cos(theta))
+  tt = GSL::Sf::legendre_sphPlm(l, m, Math::cos(theta))
 
-	pp = Math::cos(m*phi)
+  pp = Math::cos(m*phi)
 
-	rr*tt*pp
+  rr*tt*pp
 end
 
 basis_r = Proc.new { |r, y, params|
-	params.eval(r, y)
+  params.eval(r, y)
 }
 
 basis_theta = Proc.new { |theta, y, params|
-	for i in 0...N_SUM_THETA do
-		y[i] = GSL::Sf::legendre_Pl(i, Math::cos(theta));
-	end
+  for i in 0...N_SUM_THETA do
+    y[i] = GSL::Sf::legendre_Pl(i, Math::cos(theta));
+  end
 }
 
 basis_phi = Proc.new { |phi, y, params|
-	for i in 0...N_SUM_PHI do
-		if i%2 == 0
-			y[i] = Math::cos(i*0.5*phi)
-		else
-			y[i] = Math::sin((i+1.0)*0.5*phi)
-		end
-	end
+  for i in 0...N_SUM_PHI do
+    if i%2 == 0
+      y[i] = Math::cos(i*0.5*phi)
+    else
+      y[i] = Math::sin((i+1.0)*0.5*phi)
+    end
+  end
 }
 
 
@@ -67,19 +67,19 @@ data = GSL::Vector.alloc(NDATA)
 
 for i in 0...NDATA do
 
-	r = rng.uniform()*R_MAX
-	theta = rng.uniform()*Math::PI
-	phi = rng.uniform()*2*Math::PI
+  r = rng.uniform()*R_MAX
+  theta = rng.uniform()*Math::PI
+  phi = rng.uniform()*2*Math::PI
 
-	psi = psi_real_exact(k, l, m, r, theta, phi)
+  psi = psi_real_exact(k, l, m, r, theta, phi)
 
-	dpsi = rng.gaussian(0.05*psi)
+  dpsi = rng.gaussian(0.05*psi)
 
-	vars[i,0] = r
-	vars[i,1] = theta
-	vars[i,2] = phi		
+  vars[i,0] = r
+  vars[i,1] = theta
+  vars[i,2] = phi
 
-	data[i] = psi + dpsi
+  data[i] = psi + dpsi
 end
 
 #GSL::MultiFit::Ndlinear::design(vars, X, ndlinear)
@@ -99,18 +99,18 @@ x = GSL::Vector.alloc(N_DIM)
 
 r = 0.01
 while r < R_MAX do
-	theta = 0.0
-	while theta < Math::PI do
-		phi = 0.0
-		while phi < 2*Math::PI do
-			dV = r*r*Math::sin(theta)*dr*dtheta*dphi
-			x[0] = r
-			x[1] = theta
-			x[2] = phi
-
-			psi_model, err = GSL::MultiFit::Ndlinear.calc(x, coeffs, cov, ndlinear)
-			psi = psi_real_exact(k, l, m, r, theta, phi)
-			err = psi_model - psi
+  theta = 0.0
+  while theta < Math::PI do
+    phi = 0.0
+    while phi < 2*Math::PI do
+      dV = r*r*Math::sin(theta)*dr*dtheta*dphi
+      x[0] = r
+      x[1] = theta
+      x[2] = phi
+
+      psi_model, err = GSL::MultiFit::Ndlinear.calc(x, coeffs, cov, ndlinear)
+      psi = psi_real_exact(k, l, m, r, theta, phi)
+      err = psi_model - psi
       eps_rms += err * err * dV;
       volume += dV;
 
@@ -118,13 +118,13 @@ while r < R_MAX do
          printf("%e %e %e %e\n", r, theta, psi, psi_model)
       end
 
-			phi += dphi
-		end
-		theta += dtheta
-	end
-	printf("\n");
+      phi += dphi
+    end
+    theta += dtheta
+  end
+  printf("\n");
 
-	r += dr
+  r += dr
 end
 
 eps_rms /= volume