rb-gsl 1.16.0.5 → 1.16.0.6

data/examples/linalg/chol_narray.rb

@@ -1,15 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-include GSL
-include Linalg
-
-m = NArray[[4.0, 2], [2, 3]]
-c = Cholesky.decomp(m)
-p c
-
-b = NArray[1.0, 2]
-p Cholesky.solve(c, b)    # Expected [-0.125, 0.75]
-
-b = NArray[1.0, 2]
-Cholesky.svx(c, b)
-p b

data/examples/linalg/complex.rb

@@ -1,57 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-include GSL
-include Linalg
-
-m = Matrix::alloc([0.18, 0.60, 0.57, 0.96], [0.41, 0.24, 0.99, 0.58],
-                  [0.14, 0.30, 0.97, 0.66], [0.51, 0.13, 0.19, 0.85])
-m.inv.print
-p m.det
-
-#lu, perm = m.LU_decomp
-#LU.invert(lu, perm).print
-#p lu.det
-
-zm = Matrix::Complex.alloc(4, 4)
-
-zm.set_row(0, [0.18, 0], [0.60, 0], [0.57, 0], [0.96, 0])
-zm.set_row(1, [0.41, 0], [0.24, 0], [0.99, 0], [0.58, 0])
-zm.set_row(2, [0.14, 0], [0.30, 0], [0.97, 0], [0.66, 0])
-zm.set_row(3, [0.51, 0], [0.13, 0], [0.19, 0], [0.85, 0])
-
-p zm.inv
-p zm.det
-
-b = Vector::Complex.alloc(4)
-b.set(1.0, 2, 3, 4)
-
-#p zm.LU_solve(b)
-
-#lu2, perm2, signum = zm.LU_decomp
-lu2, perm2, signum = LU::decomp(zm)
-p lu2
-p lu2.det(signum)
-
-p x = lu2.solve(perm2, b)
-
-p zm*x
-
-p lu2.LU_invert(perm2)
-p lu2.LU_det(signum)
-
-#p zm.LU_svx(b)
-
-p Linalg::Complex::LU_invert(zm)
-p Linalg::Complex::LU::invert(zm)
-p Linalg::Complex::LU_invert(lu2, perm2)
-p Linalg::Complex::LU::invert(lu2, perm2)
-p zm.invert
-p lu2.invert(perm2)
-
-p Linalg::Complex::LU_det(zm)
-p Linalg::Complex::LU::det(zm)
-p Linalg::Complex::LU_det(lu2, signum)
-p Linalg::Complex::LU::det(lu2, signum)
-p zm.det
-p lu2.det(signum)
-

data/examples/linalg/invert_narray.rb

@@ -1,10 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-include GSL
-
-na = NArray[[1.0, 4], [2, 3]]
-lu, perm, signum = Linalg::LU.decomp(na)
-
-inv = Linalg::LU.invert(lu, perm)
-p inv
-

data/examples/math/const.rb

@@ -1,67 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-include GSL
-
-puts("e")
-p GSL::M_E
-p Math::exp(1)
-
-
-puts("log10(e)")
-p GSL::M_LOG10E
-p Math::log10(M_E)
-
-puts("sqrt(2)")
-p GSL::M_SQRT2
-p Math::sqrt(2)
-
-puts("sqrt(1/2)")
-p GSL::M_SQRT1_2
-p Math::sqrt(1.0/2)
-
-puts("sqrt(3)")
-p GSL::M_SQRT3
-p Math::sqrt(3)
-
-puts("pi")
-p GSL::M_PI
-p Math::PI
-
-puts("pi/2")
-p GSL::M_PI_2
-p Math::PI/2.0
-
-puts("pi/4")
-p GSL::M_PI_4
-p Math::PI/4.0
-
-puts("sqrt(pi)")
-p GSL::M_SQRTPI
-p Math::sqrt(Math::PI)
-
-puts("2/sqrt(pi)")
-p GSL::M_2_SQRTPI
-p 2.0/Math::sqrt(Math::PI)
-
-puts("1/pi")
-p GSL::M_1_PI
-p 1.0/Math::PI
-
-puts("2/pi")
-p GSL::M_2_PI
-p 2.0/Math::PI
-
-puts("ln10")
-p GSL::M_LN10
-p Math::log(10)
-
-puts("ln2")
-p GSL::M_LN2
-p Math::log(2)
-
-puts("ln(pi)")
-p GSL::M_LNPI
-p Math::log(Math::PI)
-
-puts("euler constant 0.5772156...")
-p GSL::M_EULER

data/examples/math/elementary.rb

@@ -1,35 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-
-puts("gsl_log1p")
-p GSL::log1p(0.1)
-p Math::log(1 + 0.1)
-
-puts("gsl_expm1")
-p GSL::expm1(0.1)
-p Math::exp(0.1) - 1.0
-
-puts("gsl_hypot")
-p GSL::hypot(20000, 30000)
-p Math::sqrt(20000*20000 + 30000*30000)
-
-puts("gsl_acosh")
-p GSL::acosh(1.5)
-p Math::acosh(1.5)
-
-puts("gsl_asinh")
-p GSL::asinh(1.5)
-p Math::asinh(1.5)
-
-puts("gsl_atanh")
-p GSL::atanh(0.5)
-p Math::atanh(0.5)
-
-puts("gsl_ldexp")
-p GSL::ldexp(1.5, 8)
-p Math::ldexp(1.5, 8)
-
-puts("gsl_frexp")
-p GSL::frexp(100)
-p Math::frexp(100)
-

data/examples/math/functions.rb

@@ -1,41 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-
-puts("gsl_log1p")
-p GSL::log1p(0.1)
-
-a = GSL::Vector.alloc(0.1, 1, 2, 3, 4)
-p GSL::log1p(a)
-
-b = GSL::Matrix.alloc([1, 2], [3, 4])
-p GSL::log1p(b)
-exit
-
-puts("gsl_expm1")
-p GSL::expm1(0.1)
-p Math::exp(0.1) - 1.0
-
-puts("gsl_hypot")
-p GSL::hypot(20000, 30000)
-p Math::sqrt(20000*20000 + 30000*30000)
-
-puts("gsl_acosh")
-p GSL::acosh(1.5)
-p Math::acosh(1.5)
-
-puts("gsl_asinh")
-p GSL::asinh(1.5)
-p Math::asinh(1.5)
-
-puts("gsl_atanh")
-p GSL::atanh(0.5)
-p Math::atanh(0.5)
-
-puts("gsl_ldexp")
-p GSL::ldexp(1.5, 8)
-p Math::ldexp(1.5, 8)
-
-puts("gsl_frexp")
-p GSL::frexp(100)
-p Math::frexp(100)
-

data/examples/math/inf_nan.rb

@@ -1,34 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-include GSL
-
-puts("GSL_POSINF")
-p GSL::POSINF
-p GSL::POSINF.class
-p isinf(GSL::POSINF)
-p isinf?(GSL::POSINF)
-p isnan(GSL::POSINF)
-p isnan?(GSL::POSINF)
-p finite(GSL::POSINF)
-p finite?(GSL::POSINF)
-
-puts("GSL_NEGINF")
-p GSL::NEGINF
-p GSL::NEGINF.class
-p isinf(GSL::NEGINF)
-p isinf?(GSL::NEGINF)
-p isnan(GSL::NEGINF)
-p isnan?(GSL::NEGINF)
-p finite(GSL::NEGINF)
-p finite?(GSL::NEGINF)
-
-puts("GSL_NAN")
-p GSL::NAN
-p GSL::NAN.class
-p isinf(GSL::NAN)
-p isinf?(GSL::NAN)
-p isnan(GSL::NAN)
-p isnan?(GSL::NAN)
-p finite(GSL::NAN)
-p finite?(GSL::NAN
-)

data/examples/math/minmax.rb

@@ -1,22 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-
-p GSL::MAX(1, 2)
-p GSL::MAX(1, 2.0)
-p GSL::MAX(3, 1)
-p GSL::MAX(3.0, 1)
-
-p GSL::MIN(1, 2)
-p GSL::MIN(1.0, 2.0)
-p GSL::MIN(3, 1.0)
-p GSL::MIN(3.0, 1)
-
-p GSL::MAX_INT(5, 6)
-p GSL::MAX_INT(5, 6.0)
-p GSL::MAX_INT(6, 5)
-p GSL::MAX_INT(6.0, 5)
-
-p GSL::MIN_INT(2, 2)
-p GSL::MIN_INT(2.0, 2.0)
-p GSL::MIN_INT(7, 2.0)
-p GSL::MIN_INT(7.0, 2)

data/examples/math/power.rb

@@ -1,18 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-
-p GSL::pow_int(8, 3)
-p 8*8*8.0
-
-p GSL::pow_2(2)
-p GSL::pow_3(2)
-p GSL::pow_4(2)
-p GSL::pow_5(2)
-p GSL::pow_6(2)
-p GSL::pow_7(2)
-p GSL::pow_8(2)
-p GSL::pow_9(2)
-
-p GSL::pow(3.2, 4.5)
-p GSL::pow([1, 2, 3, 4, 5], 3)
-p GSL::pow(GSL::Vector[1, 2, 3, 4, 5], 3)

data/examples/math/test.rb

@@ -1,31 +0,0 @@
-#!/usr/bin/env ruby
-require("gsl")
-
-puts("sign of 5.0")
-p GSL::SIGN(5.0)
-
-puts("sign of -2.0")
-p GSL::SIGN(-2.0)
-
-puts("sign of 0: positive")
-p GSL::SIGN(0)
-
-puts("Is 1 odd?")
-p GSL::IS_ODD(1)
-p GSL::IS_ODD?(1)
-puts("Is 1 even?")
-p GSL::IS_EVEN(1)
-p GSL::IS_EVEN?(1)
-
-puts("Is 4 odd?")
-p GSL::IS_ODD(4)
-p GSL::IS_ODD?(4)
-puts("Is 4 even?")
-p GSL::IS_EVEN(4)
-p GSL::IS_EVEN?(4)
-
-puts("frcmp")
-
-p GSL::fcmp(1.0, 1.0, 1e-10)
-p GSL::fcmp(1.1, 1.0, 1e-10)
-p GSL::fcmp(1.0, 1.1, 1e-10)

data/examples/matrix/add.rb

@@ -1,45 +0,0 @@
-#!/usr/bin/env ruby
-# Turn on warnings
-$-w = true
-
-require 'irb/xmp'
-require 'gsl'
-
-# Apparently, IRB::Frame has a bug that prevents the defaults from working, so
-# an XMP instance must be created explicitly this way instead of using the
-# otherwise convenient xmp method.
-XMP.new(IRB::Frame.top(-1)).puts <<END
-# Create to 3x3 matrices
-a = GSL::Matrix.alloc([1, 2, 3, 4, 5, 6, 7, 8, 9], 3, 3)
-b = GSL::Matrix.alloc([6, 7, 8], [2, 3, 4], [3, 4, 5])
-
-# Add together to create new matrix
-a + b
-
-# Add b to a, modifying a
-a += b
-
-# Show that a is modified
-a
-
-# Subtract b from a, modifying a
-a -= b
-
-# Add 2 to all elements of a, creating new matrix
-a + 2
-
-# Another way to add 2 to all elements of a, creating new matrix
-2 + a
-
-# Add 2 to all elements of a, modifying a
-a += 2
-
-# Show that a is modified
-a
-
-# Subtract 2 from all elements of a, modifying a
-a -= 2
-
-# Show that a is modified
-a
-END

data/examples/matrix/b.dat

@@ -1,4 +0,0 @@
-5
-6
-7
-8

data/examples/matrix/cat.rb

@@ -1,31 +0,0 @@
-#!/usr/bin/env ruby
-# Turn on warnings
-$-w = true
-
-require 'irb/xmp'
-require 'gsl'
-
-# Apparently, IRB::Frame has a bug that prevents the defaults from working, so
-# an XMP instance must be created explicitly this way instead of using the
-# otherwise convenient xmp method.
-XMP.new(IRB::Frame.top(-1)).puts <<END
-# Create a 2x2 Matrix::Int and a 2x3 Matrix::Int
-a = GSL::Matrix::Int[1..4, 2, 2]
-b = GSL::Matrix::Int[5..10, 2, 3]
-
-# Concatenate them horizontally using Matrix::Int#horzcat
-a.horzcat(b)
-
-# Concatenate them horizontally using Matrix::Int.horzcat
-GSL::Matrix::Int.horzcat(a, b)
-
-# Create a 2x2 Matrix::Int and a 3x2 Matrix::Int
-a = GSL::Matrix::Int[1..4, 2, 2]
-b = GSL::Matrix::Int[5..10, 3, 2]
-
-# Concatenate them vertically using Matrix::Int#vertcat
-a.vertcat(b)
-
-# Concatenate them vertically using Matrix::Int.vertcat
-GSL::Matrix::Int.vertcat(a, b)
-END

data/examples/matrix/colvectors.rb

@@ -1,24 +0,0 @@
-#!/usr/bin/env ruby
-# Turn on warnings
-$-w = true
-
-require 'irb/xmp'
-require 'gsl'
-
-# Apparently, IRB::Frame has a bug that prevents the defaults from working, so
-# an XMP instance must be created explicitly this way instead of using the
-# otherwise convenient xmp method.
-XMP.new(IRB::Frame.top(-1)).puts <<END
-# Create two 3x3 matrices
-m = GSL::Matrix.alloc(1..9, 3, 3)
-n = GSL::Matrix.alloc(11..19, 3, 3)
-
-# Get Vector::Col::View for column 1 of Matrix a
-a = m.col(1)
-
-# Get Vector::Col::View for column 2 of Matrix b
-b = n.col(2)
-
-# Create new Matrix from Vector::Col::Views
-c = GSL::Matrix[a, b]
-END