rb-gsl 1.16.0.5 → 1.16.0.6

data/lib/gsl.rb

@@ -1,8 +0,0 @@
-begin
-  require 'narray'
-rescue LoadError
-end
-
-require 'gsl_native'
-require 'gsl/version'
-require 'gsl/oper'

data/lib/gsl/gnuplot.rb

@@ -1,41 +0,0 @@
-require 'gnuplot'
-
-class Array
-  def to_gplot
-    if ( self[0].kind_of? Array ) then
-      tmp = self[0].zip( *self[1..-1] )
-      tmp.collect { |a| a.join(" ") }.join("\n") + "\ne"
-    elsif ( self[0].kind_of? Numeric ) then
-      s = ""
-      self.length.times { |i| s << "#{self[i]}\n" }
-      s
-    elsif ( self[0].kind_of? GSL::Vector ) then
-      tmp = self[0].zip( *self[1..-1] )
-      tmp.collect { |a| a.join(" ") }.join("\n") + "\ne"
-    else
-      self[0].zip( *self[1..-1] ).to_gplot
-    end
-  end
-  def to_gsplot
-    f = ""
-
-    if ( self[0].kind_of? Array ) then
-      x = self[0]
-      y = self[1]
-      d = self[2]
-
-      x.each_with_index do |xv, i|
-        y.each_with_index do |yv, j|
-          f << [ xv, yv, d[i][j] ].join(" ") << "\n"
-        end
-        # f << "\n"
-      end
-    elsif ( self[0].kind_of? Numeric ) then
-      self.length.times do |i| f << "#{self[i]}\n" end
-    else
-      self[0].zip( *self[1..-1] ).to_gsplot
-    end
-
-    f
-  end
-end

data/lib/gsl/oper.rb

@@ -1,43 +0,0 @@
-module GSL::Oper
-
-  def self.included(base)
-    base.class_eval {
-      alias_method :_gsl_oper_original_mul, :*
-      alias_method :_gsl_oper_original_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)
-        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
-      end
-    }
-  end
-
-end
-
-[Fixnum, Float].each { |klass| klass.send(:include, GSL::Oper) }

data/lib/gsl/version.rb

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

data/lib/ool.rb

@@ -1,22 +0,0 @@
-
-module OOL
-  module Conmin
-    class Pgrad_parameters < Array
-      attr_accessor :fmin, :tol, :alpha, :sigma1, :sigma2
-    end
-    class Spg_parameters < Array
-      attr_accessor :fmin, :tol, :M, :alphamin, :alphamax, :gamma
-      attr_accessor :sigma1, :sigma2
-    end
-    class Gencan_parameters < Array
-      attr_accessor :epsgpen, :epsgpsn, :fmin, :udelta0
-      attr_accessor :ucgmia, :ucgmib
-      attr_accessor :cg_scre, :cg_gpnf, :cg_epsi, :cg_epsf
-      attr_accessor :cg_epsnqmp, :cg_maxitnqmp, :nearlyq
-      attr_accessor :nint, :next, :mininterp, :maxextrap
-      attr_accessor :trtype, :eta, :deltamin, :lspgmi, :lspgma
-      attr_accessor :theta, :gamma, :beta, :sigma1, :sigma2
-      attr_accessor :epsrel, :epsabs, :infrel, :infabs
-    end    
-  end
-end

data/lib/ool/conmin.rb

@@ -1,30 +0,0 @@
-
-module OOL
-  module Conmin
-    class Minimizer
-      class Pgrad
-        class Parameters
-          attr_accessor :fmin, :tol, :alpha, :sigma1, :sigma2
-        end
-      end
-      class Spg
-        class Parameters
-          attr_accessor :fmin, :tol, :M, :alphamin, :alphamax, :gamma
-          attr_accessor :sigma1, :sigma2
-        end
-      end
-      class Gencan
-        class Parameters
-          attr_accessor :epsgpen, :epsgpsn, :fmin, :udelta0
-          attr_accessor :ucgmia, :ucgmib
-          attr_accessor :cg_scre, :cg_gpnf, :cg_epsi, :cg_epsf
-          attr_accessor :cg_epsnqmp, :cg_maxitnqmp, :nearlyq
-          attr_accessor :nint, :next, :mininterp, :maxextrap
-          attr_accessor :trtype, :eta, :deltamin, :lspgmi, :lspgma
-          attr_accessor :theta, :gamma, :beta, :sigma1, :sigma2
-          attr_accessor :epsrel, :epsabs, :infrel, :infabs
-        end    
-      end
-    end
-  end
-end

data/lib/rbgsl.rb

@@ -1 +0,0 @@
-require 'gsl'

data/rdoc/alf.rdoc

@@ -1,77 +0,0 @@
-#
-# = ALF: a collection of routines for computing associated Legendre polynomials (ALFs)
-# ALF is an extension library for GSL to compute associated Legendre polynomials developed by Patrick Alken.
-# Ruby/GSL includes interfaces to it if ALF is found during installlation.
-#
-# The class and method descriptions below are based on references from the document of ALF (alf-1.0/doc/alf.texi) by P.Alken.
-#
-# == Module structure
-# * GSL::ALF (module)
-#   * GSL::ALF::Workspace (Class)
-#
-# == 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).
-#
-# == Methods
-# ---
-# * GSL::ALF::Workspace#params(csphase, cnorm, norm)
-#
-#   Sets various parameters for the subsequent computation of ALFs. If
-#   <tt>csphase</tt> is set to a non-zero value, the Condon-Shortley phase of
-#   (-1)^m will be applied to the associated Legendre functions. The
-#   Condon-Shortley phase is included by default. If <tt>cnorm</tt> is set to
-#   zero, the real normalization of the associated Legendre functions will be
-#   used. The default is to use complex normalization. The norm parameter
-#   defines the type of normalization which will be used. The possible values
-#   are as follows.
-#   * ALF::NORM_SCHMIDT: Schmidt semi-normalized associated Legendre polynomials S_l^m(x). (default)
-#   * ALF::NORM_SPHARM: Associated Legendre polynomials Y_l^m(x) suitable for the calculation of spherical harmonics.
-#   * ALF::NORM_ORTHO: Fully orthonormalized associated Legendre polynomials N_l^m(x).
-#   * ALF::NORM_NONE:: Unnormalized associated Legendre polynomials P_l^m(x).
-# ---
-# * GSL::ALF::Workspace#Plm_array(x)
-# * GSL::ALF::Workspace#Plm_array(lmax, x)
-# * GSL::ALF::Workspace#Plm_array(x, result)
-# * GSL::ALF::Workspace#Plm_array(lmax, x, result)
-# * GSL::ALF::Workspace#Plm_array(x, result, deriv)
-# * GSL::ALF::Workspace#Plm_array(lmax, x, result, deriv)
-#
-#   Compute all associated Legendre polynomials P_l^m(x) and optionally their
-#   first derivatives dP_l^m(x)/dx for 0 <= l <= lmax, 0 <= m <= l. The value
-#   of <tt>lmax</tt> cannot exceed the previously specified lmax parameter to
-#   <tt>ALF.alloc</tt>, but may be less. If <tt>lmax</tt> is not given, the
-#   parameter to <tt>ALF.alloc()</tt> is used. The results are stored in
-#   <tt>result</tt>, an instance of <tt>GSL::Vector</tt>. Note that this vector
-#   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.
-#
-#   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
-#   be obtained by calling <tt>ALF::array_index(l, m)</tt>.
-# ---
-# * GSL::ALF::Workspace#Plm_deriv_array(x)
-# * GSL::ALF::Workspace#Plm_deriv_array(lmax, x)
-# * GSL::ALF::Workspace#Plm_deriv_array(x, result, deriv)
-# * 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.
-#
-# ---
-# * GSL::ALF::array_size(lmax)
-#
-#   Returns the size of arrays needed for the array versions of P_l^m(x).
-# ---
-# * GSL::ALF::array_index(l, m)
-#
-#   Returns the array index of results of <tt>Plm_array()</tt> and
-#   <tt>Plm_deriv_array()</tt> corresponding to P_l^m(x) and dP_l^m(x)/dx
-#   respectively. The index is given by l(l + 1)/2 + m.
-#
-#

data/rdoc/blas.rdoc

@@ -1,269 +0,0 @@
-#
-# = BLAS Support
-# The following is the list of the methods defined in <tt>GSL::Blas</tt> module.
-# See {GSL reference}[https://gnu.org/software/gsl/manual/gsl-ref_12.html#SEC212] for details.
-#
-# == Level 1
-# ---
-# * GSL::Blas::ddot(x, y)
-# * GSL::Vector#blas_ddot(y)
-# * GSL::Vector#ddot(y)
-#
-#
-# ---
-# * GSL::Blas::zdotu(x, y)
-# * GSL::Vector::Complex#blas_zdotu(y)
-# * GSL::Vector::Complex#zdotu(y)
-#
-#
-# ---
-# * GSL::Blas::zdotc(x, y)
-# * GSL::Vector::Complex#blas_zdotc(y)
-# * GSL::Vector::Complex#zdotc(y)
-#
-#
-# ---
-# * GSL::Blas::dnrm2(x)
-# * GSL::Vector#blas_dnrm2
-# * GSL::Vector#dnrm2
-#
-#
-# ---
-# * GSL::Blas::dznrm2(x)
-# * GSL::Vector::Complex#blas_dznrm2
-# * GSL::Vector::Complex#dznrm2
-#
-#
-# ---
-# * GSL::Blas::dasum(x)
-# * GSL::Vector#blas_dasum
-# * GSL::Vector#dasum
-#
-#
-# ---
-# * GSL::Blas::dzasum(x)
-# * GSL::Vector::Complex#blas_dzasum
-# * GSL::Vector::Complex#dzasum
-#
-#
-# ---
-# * GSL::Blas::idamax(x)
-# * GSL::Vector#blas_idamax
-# * GSL::Vector#idamax
-#
-#
-# ---
-# * GSL::Blas::izamax(x)
-# * GSL::Vector::Complex#blas_izamax
-# * GSL::Vector::Complex#izamax
-#
-#
-# ---
-# * GSL::Blas::dswap(x, y)
-# * GSL::Vector#blas_dswap(y)
-# * GSL::Vector#dswap(y)
-#
-#
-# ---
-# * GSL::Blas::zswap(x, y)
-# * GSL::Vector::Complex#blas_zswap(y)
-# * GSL::Vector::Complex#zswap(y)
-#
-#
-# ---
-# * GSL::Blas::dcopy(x, y)
-# * GSL::Vector#blas_dcopy(y)
-# * GSL::Vector#dcopy(y)
-#
-#
-# ---
-# * GSL::Blas::zcopy(x, y)
-# * GSL::Vector::Complex#blas_zcopy(y)
-# * GSL::Vector::Complex#zcopy(y)
-#
-#
-# ---
-# * GSL::Blas::daxpy!(a, x, y)
-# * GSL::Vector#blas_daxpy!(a, y)
-# * GSL::Vector#daxpy!(a, y)
-# * GSL::Blas::daxpy(a, x, y)
-# * GSL::Vector#blas_daxpy(a, y)
-# * GSL::Vector#daxpy(a, y)
-#
-#
-# ---
-# * GSL::Blas::zaxpy!(a, x, y)
-# * GSL::Vector::Complex#blas_zaxpy!(a, y)
-# * GSL::Vector::Complex#zaxpy!(a, y)
-# * GSL::Blas::zaxpy(a, x, y)
-# * GSL::Vector::Complex#blas_zaxpy(a, y)
-# * GSL::Vector::Complex#zaxpy(a, y)
-#
-#
-# ---
-# * GSL::Blas::dscal!(a, x)
-# * GSL::Vector#blas_dscal!(a)
-# * GSL::Vector#dscal!(a)
-# * GSL::Blas::dscal(a, x)
-# * GSL::Vector#blas_dscal(a)
-# * GSL::Vector#dscal(a)
-#
-#
-# ---
-# * GSL::Blas::zscal!(a, x)
-# * GSL::Vector::Complex#blas_zscal!(a)
-# * GSL::Vector::Complex#zscal!(a)
-# * GSL::Blas::zscal(a, x)
-# * GSL::Vector::Complex#blas_zscal(a)
-# * GSL::Vector::Complex#zscal(a)
-#
-#
-# ---
-# * GSL::Blas::zdscal!(a, x)
-# * GSL::Vector::Complex#blas_zdscal!(a)
-# * GSL::Vector::Complex#zdscal!(a)
-# * GSL::Blas::zdscal(a, x)
-# * GSL::Vector::Complex#blas_zdscal(a)
-# * GSL::Vector::Complex#zdscal(a)
-#
-#
-# ---
-# * GSL::Blas::drot(x, y, c, s)
-# * GSL::Blas::drot!(x, y, c, s)
-# * GSL::Blas::drotm(x, y, p)
-# * GSL::Blas::drotm!(x, y, p)
-#
-#
-# == Level 2
-# ---
-# * GSL::Blas::dgemv(trans, a, A, x, b, y)
-# * GSL::Blas::dgemv!(trans, a, A, x, b, y)
-# * GSL::Blas::zgemv(trans, a, A, x, b, y)
-# * GSL::Blas::zgemv!(trans, a, A, x, b, y)
-#
-#
-# ---
-# * GSL::Blas::dtrmv(trans, a, A, x, b, y)
-# * GSL::Blas::dtrmv!(trans, a, A, x, b, y)
-# * GSL::Blas::ztrmv(trans, a, A, x, b, y)
-# * GSL::Blas::ztrmv!(trans, a, A, x, b, y)
-#
-#
-# ---
-# * GSL::Blas::dtrsv(uplo, trans, diag, A, x)
-# * GSL::Blas::dtrsv!uplo, trans, diag, A, x)
-# * GSL::Blas::ztrsvuplo, trans, diag, A, x)
-# * GSL::Blas::ztrsv!uplo, trans, diag, A, x)
-#
-#
-# ---
-# * GSL::Blas::dsymv(uplo, a, A, x, b, y)
-# * GSL::Blas::dsymv!(uplo, a, A, x, b, y)
-# * GSL::Blas::zhemv(uplo, a, A, x, b, y)
-# * GSL::Blas::zhemv!(uplo, a, A, x, b, y)
-#
-#
-# ---
-# * GSL::Blas::dger(a, x, y, A)
-# * GSL::Blas::dger!(a, x, y, A)
-# * GSL::Blas::zgeru(a, x, y, A)
-# * GSL::Blas::zgeru!(a, x, y, A)
-# * GSL::Blas::zgerc(a, x, y, A)
-# * GSL::Blas::zgerc!(a, x, y, A)
-#
-#
-# ---
-# * GSL::Blas::dsyr(uplo, a, x, A)
-# * GSL::Blas::dsyr!(uplo, a, x, A)
-# * GSL::Blas::zher(uplo, a, x, A)
-# * GSL::Blas::zher!(uplo, a, x, A)
-#
-#
-# ---
-# * GSL::Blas::dsyr2(uplo, a, x, y, A)
-# * GSL::Blas::dsyr2!(uplo, a, x, y, A)
-# * GSL::Blas::zher2(uplo, a, x, y, A)
-# * GSL::Blas::zher2!(uplo, a, x, y, A)
-#
-#
-# == Level 3
-# ---
-# * GSL::Blas::dgemm(transA, transB, alpha, A, B, beta, C)
-# * GSL::Blas::dgemm(A, B)
-#
-#
-# ---
-# * GSL::Blas::zgemm(transA, transB, alpha, A, B, beta, C)
-# * GSL::Blas::zgemm(A, B)
-#
-#
-# ---
-# * GSL::Blas::dsymm(transA, transB, alpha, A, B, beta, C)
-# * GSL::Blas::dsymm(A, B)
-#
-#
-# ---
-# * GSL::Blas::zsymm(transA, transB, alpha, A, B, beta, C)
-# * GSL::Blas::zsymm(A, B)
-#
-#
-# ---
-# * GSL::Blas::zhemm(transA, transB, alpha, A, B, beta, C)
-# * GSL::Blas::zhemm(A, B)
-#
-#
-# ---
-# * GSL::Blas::dtrmm(side, uplo, trans, diag, alpha, A, B)
-# * GSL::Blas::ztrmm(side, uplo, trans, diag, alpha, A, B)
-#
-#
-# ---
-# * GSL::Blas::dtrsm(side, uplo, trans, diag, alpha, A, B)
-# * GSL::Blas::ztrsm(side, uplo, trans, diag, alpha, A, B)
-#
-#
-# ---
-# * GSL::Blas::dsyrk(uplo, trans, diag, alpha, A, beta, C)
-# * GSL::Blas::zsyrk(uplo, trans, diag, alpha, A, beta, C)
-# * GSL::Blas::zherk(uplo, trans, diag, alpha, A, beta, C)
-#
-#
-# ---
-# * GSL::Blas::dsyr2k(uplo, trans, diag, alpha, A, B, beta, C)
-# * GSL::Blas::zsyr2k(uplo, trans, diag, alpha, A, B, beta, C)
-# * GSL::Blas::zher2k(uplo, trans, diag, alpha, A, B, beta, C)
-#
-#
-# == Constants
-# ---
-# * GSL::Blas::CblasRowMajor
-# * GSL::Blas::RowMajor
-# * GSL::Blas::CblasColMajor
-# * GSL::Blas::ColMajor
-# * GSL::Blas::CblasNoTrans
-# * GSL::Blas::NoTrans
-# * GSL::Blas::CblasTrans
-# * GSL::Blas::Trans
-# * GSL::Blas::CblasConjTrans
-# * GSL::Blas::ConjTrans
-# * GSL::Blas::CblasUpper
-# * GSL::Blas::Upper
-# * GSL::Blas::CblasTrans
-# * GSL::Blas::Lower
-# * GSL::Blas::CblasNonUnit
-# * GSL::Blas::NonUnit
-# * GSL::Blas::CblasUnit
-# * GSL::Blas::Unit
-# * GSL::Blas::CblasLeft
-# * GSL::Blas::Left
-# * GSL::Blas::CblasRight
-# * GSL::Blas::Right
-#
-#
-# {prev}[link:rdoc/sort_rdoc.html]
-# {next}[link:rdoc/linalg_rdoc.html]
-#
-# {Reference index}[link:rdoc/ref_rdoc.html]
-# {top}[link:index.html]
-#
-#