gsl 1.15.3 → 1.16.0.6

data/gsl.gemspec

@@ -0,0 +1,29 @@
+# -*- encoding: utf-8 -*-
+require File.dirname(__FILE__) + '/lib/gsl/version'
+require 'date'
+
+Gem::Specification.new do |s|
+  s.name = 'gsl'
+  s.version = GSL::RUBY_GSL_VERSION
+  s.date = Date.today.to_s
+
+  s.require_paths = %w(lib)
+  s.authors = ['Yoshiki Tsunesada', 'David MacMahon', 'Jens Wille', 'Daniel Mendler']
+  s.summary = 'Ruby interface to the GNU Scientific Library'
+  s.description = 'Ruby/GSL is a Ruby interface to the GNU Scientific Library, for numerical computing with Ruby'
+  s.email = 'mail@daniel-mendler.de'
+  s.extensions  = Dir['ext/**/extconf.rb']
+  s.extra_rdoc_files = Dir['**/*.rdoc']
+  s.files = `git ls-files`.split("\n")
+  s.homepage = 'http://github.com/SciRuby/rb-gsl'
+  s.licenses = ['GPL-2.0']
+  s.rdoc_options = ['--title', "Ruby/GSL (#{GSL::RUBY_GSL_VERSION})", '--charset', 'UTF-8', '--line-numbers', '--all', '--main', 'index.rdoc', '--root', 'rdoc']
+  s.required_ruby_version = '>= 1.9.3'
+  s.requirements = ['GSL (http://www.gnu.org/software/gsl/)']
+
+  s.post_install_message = "#{s.name} can be installed with or without narray support. Please install narray before and reinstall #{s.name} if it is missing."
+
+  s.add_development_dependency 'rake-compiler', '>= 0'
+  s.add_development_dependency 'rake', '>= 0'
+  s.add_development_dependency 'test-unit', '>= 0'
+end

data/lib/gsl.rb

@@ -1,3 +1,8 @@
-require('narray')
-require('rb_gsl')
-require('gsl/oper.rb')
+begin
+  require 'narray'
+rescue LoadError
+end
+
+require 'gsl_native'
+require 'gsl/version'
+require 'gsl/oper'

data/lib/gsl/gnuplot.rb

@@ -1,4 +1,4 @@
-require("gnuplot")
+require 'gnuplot'
 
 class Array
   def to_gplot
@@ -18,7 +18,7 @@ class Array
   end
   def to_gsplot
     f = ""
-    
+
     if ( self[0].kind_of? Array ) then
       x = self[0]
       y = self[1]
@@ -35,7 +35,7 @@ class Array
     else
       self[0].zip( *self[1..-1] ).to_gsplot
     end
-    
+
     f
   end
 end

data/lib/gsl/oper.rb

@@ -1,68 +1,43 @@
-class Fixnum
-  alias :_orig_mul :*
-  alias :_orig_div :/
-  def *(other)
-    if other.kind_of?(GSL::Matrix) or other.kind_of?(GSL::Vector) or other.kind_of?(GSL::Matrix::Int) or other.kind_of?(GSL::Vector::Int) or other.kind_of?(GSL::Vector::Complex)or other.kind_of?(GSL::Matrix::Complex)
-      other.scale(self)
-    else
-      if GSL.have_tensor?
-        if other.kind_of?(GSL::Tensor) or other.kind_of?(GSL::Tensor::Int)
-          other.scale(self)
-        else
-          self._orig_mul(other)
-        end
-      else
-        self._orig_mul(other)
-      end
-    end
-  end
+module GSL::Oper
 
-  def /(other)
-    if other.kind_of?(GSL::Poly) or other.kind_of?(GSL::Poly::Int)
-      a = GSL::Poly[1]; a[0] = self
-      GSL::Rational.new(a, other)
-    elsif other.kind_of?(GSL::Vector::Col) 
-      other.scale(1.0/GSL::pow_2(other.dnrm2))
-    elsif other.kind_of?(GSL::Vector::Int::Col)
-      v = other.to_f
-      v.scale(1.0/GSL::pow_2(v.dnrm2))
-    else
-      self._orig_div(other)
-    end
-  end
-end
+  def self.included(base)
+    base.class_eval {
+      alias_method :_gsl_oper_original_mul, :*
+      alias_method :_gsl_oper_original_div, :/
 
-class Float
-  alias :_orig_mul :*
-  alias :_orig_div :/
+      def *(other)
+        case other
+          when Numeric
+            _gsl_oper_original_mul(other)
+          when GSL::Matrix,          GSL::Vector,
+               GSL::Matrix::Int,     GSL::Vector::Int,
+               GSL::Vector::Complex, GSL::Matrix::Complex,
+               *GSL.have_tensor? ? [GSL::Tensor, GSL::Tensor::Int] : []
+            other.scale(self)
+          else
+            _gsl_oper_original_mul(other)
+        end
+      end
 
-  def *(other)
-    if other.kind_of?(GSL::Matrix) or other.kind_of?(GSL::Vector) or other.kind_of?(GSL::Matrix::Int) or other.kind_of?(GSL::Vector::Int) or other.kind_of?(GSL::Vector::Complex)or other.kind_of?(GSL::Matrix::Complex)
-      other.scale(self)
-    else
-      if GSL.have_tensor?
-        if other.kind_of?(GSL::Tensor) or other.kind_of?(GSL::Tensor::Int)
-          other.scale(self)
-        else
-          self._orig_mul(other)
+      def /(other)
+        case other
+          when Numeric
+            _gsl_oper_original_div(other)
+          when GSL::Poly, GSL::Poly::Int
+            a = GSL::Poly[1]; a[0] = self
+            GSL::Rational.new(a, other)
+          when GSL::Vector::Col
+            other.scale(self / GSL.pow_2(other.dnrm2))
+          when GSL::Vector::Int::Col
+            v = other.to_f
+            v.scale(self / GSL.pow_2(v.dnrm2))
+          else
+            _gsl_oper_original_div(other)
         end
-      else
-        self._orig_mul(other)
       end
-    end
+    }
   end
 
-  def /(other)
-    if other.kind_of?(GSL::Poly) or other.kind_of?(GSL::Poly::Int)
-      a = GSL::Poly[1]; a[0] = self
-      GSL::Rational.new(a, other)
-    elsif other.kind_of?(GSL::Vector::Col) 
-      other.scale(1.0/GSL::pow_2(other.dnrm2))
-    elsif other.kind_of?(GSL::Vector::Int::Col)
-      v = other.to_f
-      v.scale(1.0/GSL::pow_2(v.dnrm2))
-    else
-      self._orig_div(other)
-    end
-  end
 end
+
+[Fixnum, Float].each { |klass| klass.send(:include, GSL::Oper) }

data/lib/gsl/version.rb

@@ -0,0 +1,3 @@
+module GSL
+  RUBY_GSL_VERSION = RB_GSL_VERSION = '1.16.0.6'
+end

data/lib/rbgsl.rb

@@ -1,3 +1 @@
-require('narray')
-require('rb_gsl')
-require('gsl/oper.rb')
+require 'gsl'

data/rdoc/alf.rdoc

@@ -5,18 +5,18 @@
 #
 # The class and method descriptions below are based on references from the document of ALF (alf-1.0/doc/alf.texi) by P.Alken.
 #
-# == {}[link:index.html"name="1] Module structure
+# == Module structure
 # * GSL::ALF (module)
 #   * GSL::ALF::Workspace (Class)
 #
-# == {}[link:index.html"name="2] Creating ALF workspace
+# == Creating ALF workspace
 # ---
 # * GSL::ALF::Workspace.alloc(lmax)
 # * GSL::ALF.alloc(lmax)
 #
 # Creates a workspace for computing associated Legendre polynomials (ALFs). The maximum ALF degree is specified by lmax. The size of this workspace is O(lmax).
 #
-# == {}[link:index.html"name="3] Methods
+# == Methods
 # ---
 # * GSL::ALF::Workspace#params(csphase, cnorm, norm)
 #
@@ -49,7 +49,7 @@
 #   must have enough length to store all the values for the polynomial
 #   P_l^m(x), and the length required can be known using
 #   <tt>ALF::array_size(lmax)</tt>. If a vector is not given, a new vector is
-#   created and returned. 
+#   created and returned.
 #
 #   The indices of <tt>result</tt> (and <tt>deriv</tt> corresponding to the
 #   associated Legendre function of degree <tt>l</tt> and order <tt>m</tt> can
@@ -61,7 +61,7 @@
 # * GSL::ALF::Workspace#Plm_deriv_array(lmax, x, result, deriv)
 #
 #   Compute all associated Legendre polynomials P_l^m(x) and their first
-#   derivatives dP_l^m(x)/dx for 0 <= l <= lmax, 0 <= m <= l. 
+#   derivatives dP_l^m(x)/dx for 0 <= l <= lmax, 0 <= m <= l.
 #
 # ---
 # * GSL::ALF::array_size(lmax)

data/rdoc/blas.rdoc

@@ -1,9 +1,9 @@
 #
 # = BLAS Support
 # The following is the list of the methods defined in <tt>GSL::Blas</tt> module.
-# See {GSL reference}[http://www.gnu.org/software/gsl/manual/gsl-ref_12.html#SEC212"target="_top] for details.
+# See {GSL reference}[https://gnu.org/software/gsl/manual/gsl-ref_12.html#SEC212] for details.
 #
-# == {}[link:index.html"name="1] Level 1
+# == Level 1
 # ---
 # * GSL::Blas::ddot(x, y)
 # * GSL::Vector#blas_ddot(y)
@@ -134,7 +134,7 @@
 # * GSL::Blas::drotm!(x, y, p)
 #
 #
-# == {}[link:index.html"name="2] Level 2
+# == Level 2
 # ---
 # * GSL::Blas::dgemv(trans, a, A, x, b, y)
 # * GSL::Blas::dgemv!(trans, a, A, x, b, y)
@@ -186,7 +186,7 @@
 # * GSL::Blas::zher2!(uplo, a, x, y, A)
 #
 #
-# == {}[link:index.html"name="3] Level 3
+# == Level 3
 # ---
 # * GSL::Blas::dgemm(transA, transB, alpha, A, B, beta, C)
 # * GSL::Blas::dgemm(A, B)
@@ -234,7 +234,7 @@
 # * GSL::Blas::zher2k(uplo, trans, diag, alpha, A, B, beta, C)
 #
 #
-# == {}[link:index.html"name="4] Constants
+# == Constants
 # ---
 # * GSL::Blas::CblasRowMajor
 # * GSL::Blas::RowMajor
@@ -260,10 +260,10 @@
 # * GSL::Blas::Right
 #
 #
-# {prev}[link:files/rdoc/sort_rdoc.html]
-# {next}[link:files/rdoc/linalg_rdoc.html]
+# {prev}[link:rdoc/sort_rdoc.html]
+# {next}[link:rdoc/linalg_rdoc.html]
 #
-# {Reference index}[link:files/rdoc/ref_rdoc.html]
-# {top}[link:files/rdoc/index_rdoc.html]
+# {Reference index}[link:rdoc/ref_rdoc.html]
+# {top}[link:index.html]
 #
 #

data/rdoc/bspline.rdoc

@@ -2,24 +2,24 @@
 # = Basis Splines
 # This chapter describes functions for the computation of smoothing basis splines (B-splines). This is only for GSL-1.9 or later.
 #
-# 1. {Overview}[link:files/rdoc/bspline_rdoc.html#1]
-# 1. {Initializing the B-splines solver}[link:files/rdoc/bspline_rdoc.html#2]
-# 1. {Constructing the knots vector}[link:files/rdoc/bspline_rdoc.html#3]
-# 1. {Evaluation of B-splines}[link:files/rdoc/bspline_rdoc.html#4]
+# 1. {Overview}[link:rdoc/bspline_rdoc.html#label-Overview]
+# 1. {Initializing the B-splines solver}[link:rdoc/bspline_rdoc.html#label-Initializing+the+B-splines+solver]
+# 1. {Constructing the knots vector}[link:rdoc/bspline_rdoc.html#label-Constructing+the+knots+vector]
+# 1. {Evaluation of B-splines}[link:rdoc/bspline_rdoc.html#label-Evaluation+of+B-splines]
 #
-# == {}[link:index.html"name="1] Overview
+# == Overview
 #
-# B-splines are commonly used as basis functions to fit smoothing curves to large data sets. To do this, the abscissa axis is broken up into some number of intervals, where the endpoints of each interval are called breakpoints. These breakpoints are then converted to knots by imposing various continuity and smoothness conditions at each interface. Given a nondecreasing knot vector t = \{t_0, t_1, \dots, t_{n+k-1\, the n basis splines of order k are defined by for i = 0, \dots, n-1. The common case of cubic B-splines is given by k = 4. The above recurrence relation can be evaluated in a numerically stable way by the de Boor algorithm. 
+# B-splines are commonly used as basis functions to fit smoothing curves to large data sets. To do this, the abscissa axis is broken up into some number of intervals, where the endpoints of each interval are called breakpoints. These breakpoints are then converted to knots by imposing various continuity and smoothness conditions at each interface. Given a nondecreasing knot vector t = \{t_0, t_1, \dots, t_{n+k-1\, the n basis splines of order k are defined by for i = 0, \dots, n-1. The common case of cubic B-splines is given by k = 4. The above recurrence relation can be evaluated in a numerically stable way by the de Boor algorithm.
 #
-# If we define appropriate knots on an interval [a,b] then the B-spline basis functions form a complete set on that interval. Therefore we can expand a smoothing function as given enough (x_j, f(x_j)) data pairs. The c_i can be readily obtained from a least-squares fit. 
+# If we define appropriate knots on an interval [a,b] then the B-spline basis functions form a complete set on that interval. Therefore we can expand a smoothing function as given enough (x_j, f(x_j)) data pairs. The c_i can be readily obtained from a least-squares fit.
 #
-# == {}[link:index.html"name="2] Initializing the B-splines solver 
+# == Initializing the B-splines solver
 # ---
 # * GSL::BSpline.alloc(k, nbreak)
 #
-#   This method creates a workspace for computing B-splines of order <tt>k</tt>. The number of breakpoints is given by <tt>nbreak</tt>. This leads to <tt>n = nbreak + k - 2</tt> basis functions. Cubic B-splines are specified by <tt>k = 4</tt>. The size of the workspace is <tt>O(5k + nbreak)</tt>. 
+#   This method creates a workspace for computing B-splines of order <tt>k</tt>. The number of breakpoints is given by <tt>nbreak</tt>. This leads to <tt>n = nbreak + k - 2</tt> basis functions. Cubic B-splines are specified by <tt>k = 4</tt>. The size of the workspace is <tt>O(5k + nbreak)</tt>.
 #
-# == {}[link:index.html"name="3] Constructing the knots vector 
+# == Constructing the knots vector
 # ---
 # * GSL::BSpline#knots(breakpts)
 #
@@ -27,16 +27,16 @@
 # ---
 # * GSL::BSpline#knots_uniform(a, b)
 #
-#   This method assumes uniformly spaced breakpoints on [<tt>a,b</tt>] and constructs the corresponding knot vector using the previously specified <tt>nbreak</tt> parameter. 
-# == {}[link:index.html"name="4] Evaluation of B-splines 
+#   This method assumes uniformly spaced breakpoints on [<tt>a,b</tt>] and constructs the corresponding knot vector using the previously specified <tt>nbreak</tt> parameter.
+# == Evaluation of B-splines
 # ---
 # * GSL::BSpline#eval(x[, B])
 #
-#   This method evaluates all B-spline basis functions at the position <tt>x</tt> and stores them in <tt>B</tt> (if given), so that the ith element of <tt>B</tt> is <tt>B_i(x)</tt>. <tt>B</tt> must be of length <tt>n = nbreak + k - 2</tt>. If <tt>B</tt> is not given, a newly created vector is returned.It is far more efficient to compute all of the basis functions at once than to compute them individually, due to the nature of the defining recurrence relation. 
+#   This method evaluates all B-spline basis functions at the position <tt>x</tt> and stores them in <tt>B</tt> (if given), so that the ith element of <tt>B</tt> is <tt>B_i(x)</tt>. <tt>B</tt> must be of length <tt>n = nbreak + k - 2</tt>. If <tt>B</tt> is not given, a newly created vector is returned.It is far more efficient to compute all of the basis functions at once than to compute them individually, due to the nature of the defining recurrence relation.
 #
-# {prev}[link:files/rdoc/nonlinearfit_rdoc.html]
-# {next}[link:files/rdoc/const_rdoc.html]
+# {prev}[link:rdoc/nonlinearfit_rdoc.html]
+# {next}[link:rdoc/const_rdoc.html]
 #
-# {Reference index}[link:files/rdoc/ref_rdoc.html]
-# {top}[link:files/rdoc/index_rdoc.html]
+# {Reference index}[link:rdoc/ref_rdoc.html]
+# {top}[link:index.html]
 #

data/rdoc/changes.rdoc

@@ -10,12 +10,7 @@
 # The changes that could break old scripts are described below.  They are
 # followed by highlights of the other enhancements and new features.
 #
-# A complete log of changes since Ruby/GSL version 1.10.3 can be obtained by
-# using Subversion's log command to look at the repository...
-#
-#   svn log http://rb-gsl.rubyforge.org/svn/trunk/rb-gsl
-#
-# == {}[link:index.html"name="1] Backwards incompatibilities introduced in Ruby/GSL 1.11.2
+# == Backwards incompatibilities introduced in Ruby/GSL 1.11.2
 #
 # * FFT interface has been extensively changed.
 #
@@ -39,13 +34,13 @@
 #
 # * GSL::Complex objects are now immutable.
 #
-# == {}[link:index.html"name="2] Enhancements and features introduced in Ruby/GSL 1.11.2
+# == Enhancements and features introduced in Ruby/GSL 1.11.2
 #
 # * Now distributed and installable as a Ruby Gem.
 #
 # * Allow nil for real and/or imag in GSL::Vector::Complex#set(i,re,im) to NOT
-#   set that component.  For example, "z=GSL::Vector::Complex[1]; z[0] = [nil,
-#   1]" will set z[0] to GSL::Complex[0,1].
+#   set that component.  For example, <tt>z=GSL::Vector::Complex[1]; z[0] = [nil,
+#   1]</tt> will set \z[0] to \GSL::Complex[0,1].
 #
 # * Added GSL::Vector::Complex#to_s method.
 #

data/rdoc/cheb.rdoc

@@ -1,24 +1,24 @@
 #
 # = Chebyshev Approximations
-# This chapter describes routines for computing Chebyshev approximations to 
-# univariate functions. A Chebyshev approximation is a truncation of the series 
-# f(x) = \sum c_n T_n(x), 
-# where the Chebyshev polynomials T_n(x) = \cos(n \arccos x) 
-# provide an orthogonal basis of polynomials on the interval [-1,1] 
-# with the weight function 1 / \sqrt{1-x^2}. 
-# The first few Chebyshev polynomials are, 
-# T_0(x) = 1, T_1(x) = x, T_2(x) = 2 x^2 - 1. 
-# For further information see Abramowitz & Stegun, Chapter 22. 
-#
-# 1. {GSL::Cheb class}[link:files/rdoc/cheb_rdoc.html#1]
-# 1. {Chebyshev Series Evaluation}[link:files/rdoc/cheb_rdoc.html#2]
-# 1. {Derivatives and Integrals}[link:files/rdoc/cheb_rdoc.html#3]
-# 1. {Examples}[link:files/rdoc/cheb_rdoc.html#4]
-#
-# == {}[link:index.html"name="1] <tt>GSL::Cheb</tt> class
+# This chapter describes routines for computing Chebyshev approximations to
+# univariate functions. A Chebyshev approximation is a truncation of the series
+# f(x) = \sum c_n T_n(x),
+# where the Chebyshev polynomials T_n(x) = \cos(n \arccos x)
+# provide an orthogonal basis of polynomials on the interval [-1,1]
+# with the weight function 1 / \sqrt{1-x^2}.
+# The first few Chebyshev polynomials are,
+# T_0(x) = 1, T_1(x) = x, T_2(x) = 2 x^2 - 1.
+# For further information see Abramowitz & Stegun, Chapter 22.
+#
+# 1. {GSL::Cheb class}[link:rdoc/cheb_rdoc.html#label-Cheb+class]
+# 1. {Chebyshev Series Evaluation}[link:rdoc/cheb_rdoc.html#label-Chebyshev+Series+Evaluation]
+# 1. {Derivatives and Integrals}[link:rdoc/cheb_rdoc.html#label-Derivatives+and+Integrals]
+# 1. {Examples}[link:rdoc/cheb_rdoc.html#label-Example]
+#
+# == <tt>Cheb</tt> class
 #
 # ---
-# * GSL::Cheb.alloc(n) 
+# * GSL::Cheb.alloc(n)
 #
 #   This create an instance of the GSL::Cheb class for a Chebyshev series of order n.
 #
@@ -26,7 +26,7 @@
 # ---
 # * GSL::Cheb#init(f, a, b)
 #
-#   This computes the Chebyshev approximation the function <tt>f</tt> over the range (<tt>a,b</tt>) to the previously specified order. Where <tt>f</tt> is a {GSL::Function}[link:files/rdoc/function_rdoc.html] object. The computation of the Chebyshev approximation is an O(n^2) process, and requires <tt>n</tt> function evaluations.
+#   This computes the Chebyshev approximation the function <tt>f</tt> over the range (<tt>a,b</tt>) to the previously specified order. Where <tt>f</tt> is a {GSL::Function}[link:rdoc/function_rdoc.html] object. The computation of the Chebyshev approximation is an O(n^2) process, and requires <tt>n</tt> function evaluations.
 #
 #   * ex: Approximate a step function defined in (0, 1) by a Chebyshev series of order 40.
 #       f = GSL::Function.alloc { |x|
@@ -40,7 +40,7 @@
 #       cs = GSL::Cheb.alloc(40)
 #       cs.init(f, 0, 1)
 #
-# == {}[link:index.html"name="2] Chebyshev Series Evaluation
+# == Chebyshev Series Evaluation
 # ---
 # * GSL::Cheb#eval(x)
 #
@@ -51,7 +51,7 @@
 #
 #   This evaluates the Chebyshev series at a given point <tt>x</tt>, to (at most) the given order <tt>n</tt>.
 #
-# == {}[link:index.html"name="3] Derivatives and Integrals 
+# == Derivatives and Integrals
 #
 # ---
 # * GSL::Cheb#calc_deriv()
@@ -65,7 +65,7 @@
 #
 #   This computes the integral of the series, and returns a new GSL::Cheb object which contains the computed integral coefficients. The reciever is not changed.
 #
-# == {}[link:index.html"name="4] Example
+# == Example
 #    #!/usr/bin/env ruby
 #    require("gsl")
 #
@@ -91,9 +91,9 @@
 #
 # See also the example scripts in <tt>examples/cheb/</tt>.
 #
-# {prev}[link:files/rdoc/diff_rdoc.html]
-# {next}[link:files/rdoc/sum_rdoc.html]
+# {prev}[link:rdoc/diff_rdoc.html]
+# {next}[link:rdoc/sum_rdoc.html]
 #
-# {Reference index}[link:files/rdoc/ref_rdoc.html]
-# {top}[link:files/rdoc/index_rdoc.html]
+# {Reference index}[link:rdoc/ref_rdoc.html]
+# {top}[link:index.html]
 #