rb-gsl 1.16.0.5 → 1.16.0.6

data/Gemfile

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

data/README.md

@@ -1,77 +0,0 @@
-# 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,20 +0,0 @@
-require File.expand_path(%q{../lib/gsl/version}, __FILE__)
-
-require 'bundler/setup'
-require 'rubygems/package_task'
-require 'rake/extensiontask'
-require 'rake/testtask'
-
-Bundler::GemHelper.install_tasks
-
-Rake::TestTask.new do |t|
-  t.libs << 'test'
-  t.libs << 'test/gsl'
-  t.test_files = FileList['test/*.rb', 'test/gsl/*.rb']
-end
-
-spec = eval(IO.read('rb-gsl.gemspec'))
-Gem::PackageTask.new(spec).define
-Rake::ExtensionTask.new(:gsl_native, spec)
-
-task :default => [:compile, :test]

data/THANKS

@@ -1,17 +0,0 @@
-Takeshi Horinouchi
-Masahide Kikkawa
-Seiya Nishizawa (Patches, NArray supports)
-Arno Erpenbeck (ruby-gsl project)
-Austin Newton
-Jordi Burguet Castell (gsl_tensor)
-Cameron McBride (Patches, vector/block manipulations)
-KOSHIRO Tsuyoshi (patches)
-GSL development team
-zunda (patches)
-Christopher Butler
-Andrea Censi <andrea.censi@dis.uniroma1.it> (Bug report on Vector#*)
-Daigo Moriwaki (Bug report)
-Joe Wölfel (Ruby-1.9 support)
-Cyrus Hall (patch for Ruby-1.9)
-Dave MacMahon (patch for FFT and interface refinements)
-Roman Shterenzon (Non-negative Matrix Factorization)

data/examples/alf/alf.gp

@@ -1,15 +0,0 @@
-set multiplot
-set size 0.5, 0.5
-set grid
-
-set origin 0, 0.5
-plot "alf.dat" u 1:2, sin(x)
-
-set origin 0.5, 0.5
-plot "alf.dat" u 1:3, 3*cos(x)*sin(x)
-
-set origin 0, 0
-plot "alf.dat" u 1:4, 3*sin(x)*sin(x)
-
-set origin 0.5, 0
-plot "alf.dat" u 1:5, 1.5*sin(x)*(5*cos(x)*cos(x)-1)

data/examples/alf/alf.rb

@@ -1,32 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-
-exit unless GSL.const_defined?("ALF")
-
-lmax = 3
-#w = GSL::ALF::Workspace.alloc(lmax)
-#w = GSL::ALF.alloc(lmax)
-w = GSL::alf_alloc(lmax)
-w.params(0, 0, GSL::ALF::NORM_NONE)
-
-
-result = GSL::Vector.alloc(GSL::ALF::array_size(lmax))
-deriv = GSL::Vector.alloc(GSL::ALF::array_size(lmax))
-
-ind11 = GSL::ALF::array_index(1, 1)
-ind21 = GSL::ALF::array_index(2, 1)
-ind22 = GSL::ALF::array_index(2, 2)
-ind31 = GSL::ALF::array_index(3, 1)
-
-File.open("alf.dat", "w") do |f|
-  theta = 0.01
-  while theta < 1.99*Math::PI do
-    x = Math.cos(theta)
-#    w.Plm_deriv_array(lmax, x, result, deriv)
-    w.Plm_array(lmax, x, result, deriv)
-    f.printf("%e %e %e %e %e\n", theta, result[ind11], result[ind21], result[ind22], result[ind31])
-    theta += 0.01
-  end
-end
-
-system("gnuplot -persist alf.gp")

data/examples/blas/blas.rb

@@ -1,13 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-
-a = GSL::Vector[0.11, 0.12, 0.13, 0.21, 0.22, 0.23]
-b = GSL::Vector[1011, 1012, 1021, 1022, 1031, 1032]
-
-A = a.matrix_view(2, 3)
-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

data/examples/blas/dnrm2.rb

@@ -1,16 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-include Math
-
-m = GSL::Matrix.alloc(10, 10)
-for i in 0...10
-  for j in 0...10
-    m.set(i, j, sin(i) + cos(j))
-  end
-end
-
-for j in 0...10
-  col = m.col(j)
-  d = GSL::Blas.dnrm2(col)
-  printf("matrix column %d, norm = %g\n", j, d)
-end

data/examples/blas/level1.rb

@@ -1,81 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-include GSL
-include GSL::Blas
-
-x = GSL::Vector[1, 2, 3]
-y = GSL::Vector[4, 5, 6]
-
-p x.blas_ddot(y)
-p GSL::Blas.ddot(x, y)
-
-xz = GSL::Vector::Complex.alloc([1, 0], [2, 0], [3, 0])
-yz = GSL::Vector::Complex.alloc([4, 0], [5, 0], [6, 0])
-p zdotu(xz, yz)
-
-p x.nrm2
-p GSL::Blas.dnrm2(x)
-p xz.nrm2
-p GSL::Blas.dznrm2(xz)
-
-p x.asum
-p GSL::Blas.dasum(x)
-p xz.asum
-p GSL::Blas.dzasum(xz)
-
-x.swap(y)
-p x
-p y
-
-GSL::Blas.dswap(x, y)
-p x
-p y
-
-xz.swap(yz)
-p xz
-p yz
-
-GSL::Blas.zswap(xz, yz)
-p xz
-p yz
-
-p x.axpy(2, y)
-p y
-p GSL::Blas.daxpy(2, x, y)
-p y
-
-p x.axpy!(2, y)
-p y
-p GSL::Blas.daxpy!(2, x, y)
-p y
-
-y = GSL::Vector[4, 5, 6]
-
-az = GSL::Complex.alloc(2, 0)
-p xz.zaxpy(az, yz)
-p yz
-p GSL::Blas.zaxpy(az, xz, yz)
-p yz
-
-az = GSL::Complex[2, 0]
-p xz.axpy!(az, yz)
-p yz
-p GSL::Blas.zaxpy!(az, xz, yz)
-p yz
-
-yz = GSL::Vector::Complex.alloc([4, 0], [5, 0], [6, 0])
-
-p x.scal(2)
-p x
-
-p x.scal!(2)
-p x
-
-x = GSL::Vector[1, 2, 3]
-
-p xz.zscal(az)
-p xz.zdscal(2)
-
-p y
-p GSL::Blas.drot(x, y, 2, 3)
-p drot(x, y, 2, 3)

data/examples/blas/level2.rb

@@ -1,11 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-include GSL::Blas
-
-x = GSL::Vector[1, 2, 3]
-y = GSL::Vector[4, 5, 6]
-
-A = GSL::Matrix.alloc([1, 2, 3], [4, 5, 6], [7, 8, 9])
-# [28, 64, 100]
-p GSL::Blas.dgemv(CblasNoTrans, 2, A, x)
-p dgemv(CblasNoTrans, 2, A, x)

data/examples/blas/level3.rb

@@ -1,12 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-include GSL::Blas
-
-a = GSL::Matrix.alloc([1, 5], [4, 6])
-b = GSL::Matrix.alloc([10, 2], [3, 8])
-c = GSL::Matrix.alloc(2, 2)
-
-p GSL::Blas.dgemm(NoTrans, NoTrans, 1, a, b, 1, c)
-p GSL::Blas.dgemm(a, b)
-p a*b
-

data/examples/bspline.rb

@@ -1,57 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-
-N = 200
-NCOEFFS = 8
-NBREAK = NCOEFFS - 2
-
-GSL::Rng::env_setup()
-r = GSL::Rng.alloc()
-
-bw = GSL::BSpline.alloc(4, NBREAK)
-B = GSL::Vector.alloc(NCOEFFS)
-x = GSL::Vector.alloc(N)
-y = GSL::Vector.alloc(N)
-xx = GSL::Matrix.alloc(N, NCOEFFS)
-w = GSL::Vector.alloc(N)
-
-#printf("#m=0,S=0\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
-
-bw.knots_uniform(0.0, 15.0)
-
-for i in 0...N do
-  xi = x[i]
-  bw.eval(xi, B)
-  for j in 0...NCOEFFS do
-    xx[i,j] = B[j]
-  end
-end
-
-c, cov, chisq = GSL::MultiFit.wlinear(xx, w, y)
-
-#printf("#m=0,S=0\n")
-
-x2 = GSL::Vector.linspace(0, 15, 150)
-y2 = GSL::Vector.alloc(150)
-x2.each_index do |i|
-  bw.eval(x2[i], B)
-  yi, yerr = GSL::MultiFit::linear_est(B, c, cov)
-  y2[i] = yi
-#  printf("%f %f\n", xi, yi)
-end
-
-GSL::graph([x, y], [x2, y2], "-T X -C -X x -Y y -x 0 15 -y -1 1.3")

data/examples/cdf.rb

@@ -1,16 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-
-x = 2.0
-
-P = GSL::Cdf::ugaussian_P(x);
-printf("prob(x < %f) = %f\n", x, P);
-
-Q = GSL::Cdf::ugaussian_Q(x);
-printf("prob(x > %f) = %f\n", x, Q);
-
-x = GSL::Cdf::ugaussian_Pinv(P);
-printf("Pinv(%f) = %f\n", P, x);
-
-x = GSL::Cdf::ugaussian_Qinv(Q);
-printf("Qinv(%f) = %f\n", Q, x);

data/examples/cheb.rb

@@ -1,21 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-
-f = GSL::Function.alloc { |x|
-  if x < 0.5
-    0.25
-  else
-    0.75
-  end
-}
-
-n = 1000
-order = 40
-cs = GSL::Cheb.alloc(order)
-
-x = GSL::Vector.linspace(0, 1, n)
-ff = f.eval(x)
-cs.init(f, 0, 1)
-r10 = cs.eval_n(10, x)
-r40 = cs.eval(x)
-GSL::graph(x, ff, r10, r40)

data/examples/combination.rb

@@ -1,23 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-
-c = 0
-printf("All subsets of {0,1,2,3} by size:\n") ;
-for i in 0...4 do
-  c = GSL::Combination.calloc(4, i);
-  begin
-    printf("{");
-    c.fprintf(STDOUT, " %u");
-    printf(" }\n");
-  end while c.next == GSL::SUCCESS
-end
-
-p c
-c2 = c.clone
-p c2
-
-a = c.data
-p a.class
-p a
-p a.print
-

data/examples/complex/RC-lpf.rb

@@ -1,47 +0,0 @@
-#!/usr/bin/env ruby
-# Transfer function of a RC low-pass filter
-require("gsl")
-include GSL::CONST::NUM
-
-class RC_LPF
-  def initialize(r, c)
-    @r = r
-    @c = c
-    @fp = 1.0/(2*Math::PI*r*c)
-    @omegap = 1.0/(r*c)
-    puts("Pole at #{fp} Hz")
-  end
-
-  def Vout(f)
-    if f.class == GSL::Vector
-      out = GSL::Vector::Complex[f.size]
-      i = 0
-      f.each do |freq|
-        a = 1.0/GSL::Complex[1.0, freq/@fp]
-        out[i] = a
-        i += 1
-      end
-      return out
-    else
-      1.0/GSL::Complex(1.0, f/@fp)
-    end
-  end
-
-  attr_reader :r, :c, :fp
-end
-
-# Create RC filter
-R = 1*KILO    # 1 [kOhm]
-C = 1*MICRO   # 1 [muF]
-lpf = RC_LPF.new(R, C)
-
-# Frequency, 1Hz - 10kHz, 100 divisions
-f = GSL::Vector.logspace2(1, 10*KILO, 100)
-
-# Transfer function
-tf = lpf.Vout(f)
-
-GSL::graph([f, tf.amp.dB], "-C -l x -g 3 -y -100 10 -X 'Frequency [Hz]' -Y 'dB' --toggle-rotate-y-label")
-GSL::graph([f, tf.phase/Math::PI], "-C -l x -g 3 -X 'Frequency [Hz]' -Y 'Phase [Pi]' --toggle-rotate-y-label")
-
-