gsl 1.12.109 → 1.14.5

data/examples/linalg/chol_narray.rb

@@ -0,0 +1,15 @@
+#!/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

@@ -0,0 +1,57 @@
+#!/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

@@ -0,0 +1,10 @@
+#!/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

@@ -0,0 +1,67 @@
+#!/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

@@ -0,0 +1,35 @@
+#!/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

@@ -0,0 +1,41 @@
+#!/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

@@ -0,0 +1,34 @@
+#!/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

@@ -0,0 +1,22 @@
+#!/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

@@ -0,0 +1,18 @@
+#!/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

@@ -0,0 +1,31 @@
+#!/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

@@ -0,0 +1,45 @@
+#!/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

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

data/examples/matrix/cat.rb

@@ -0,0 +1,31 @@
+#!/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

@@ -0,0 +1,24 @@
+#!/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