gsl 1.12.109 → 1.14.5

data/examples/vector/view.rb

@@ -0,0 +1,35 @@
+#!/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 test vector v
+v = GSL::Vector.alloc(1, 2, 3, 4, 5, 6, 7, 8, 9)
+
+# Create Vector::View of Vector v starting at element 2 and spanning 3 elements
+vv = v.subvector(2, 3)
+
+# Set element 1 of View to 9.0
+vv.set(1, 9)
+
+# Show that Vector is also modified
+v
+
+# Create Matrix::View of Vector
+m = v.matrix_view(3, 3)
+
+# Create Vector::View spanning the entire Vector
+v2 = v.view
+
+# Create Vector::View of Vector starting at element 0 and spanning 3 elements
+v2 = v.view(3)
+
+# Show size of Vector::View
+v2.size
+END

data/examples/vector/view_all.rb

@@ -0,0 +1,73 @@
+#!/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
+# These examples show all(?) the ways that a Vector::View can be created using
+# Vector#subvector or its alias Vector#view.  Note that Vector#get or, more
+# commonly, its alias Vector#[] can also be used to create a Vector::View.  See
+# examples/vector/get_all.rb for more examples.
+
+# Create test vector v
+v = GSL::Vector.indgen(9)
+
+# Vector#view with zero args returns a Vector::View of entire Vector
+v.view
+
+# Vector#view with one Fixnum argument, i, returns a Vector::View of the first
+# i'th elements if i is positive or the last i'th elements if i is negative.
+v.view(3)
+v.view(-3)
+
+# Vector#view with one Range argument returns a Vector::View of the specified
+# elements.  If the begin value is greater than the end value, the View will
+# have the elements in reverse order.  If begin and/or end value is negative,
+# the value is taken to be "size-n".
+v.view(1..4)
+v.view(4..1)
+v.view(1...4)
+v.view(4...1)
+
+v.view(4..-2)
+v.view(-2..4)
+v.view(4...-2)
+v.view(-2...4)
+
+v.view(-4..8)
+v.view(8..-4)
+v.view(-4...8)
+v.view(8...-4)
+
+v.view(-5..-2)
+v.view(-2..-5)
+v.view(-5...-2)
+v.view(-2...-5)
+
+# Vector#view with a Range argument and a Fixnum argument is like a single
+# Range argument, but with a stride given by the Fixnum argument.
+v.view(1..7, 3)
+v.view(7..1, 3)
+v.view(1...7, 3)
+v.view(7...1, 3)
+
+# Vector#view with two Fixnum arguments is offset, length.  If offset is
+# negative, it means size+offset.  If length is negative, it means step is -1.
+v.view(2, 4)
+v.view(4, 2)
+v.view(-4, 2)
+v.view(-2, -4)
+
+# Vector#view with three Fixnum arguments is offset, stride, length.  If offset
+# is negative, it means size+offset.  If length is negative, the sign of both
+# stride and length is inverted.
+v.view(1, 2, 3)
+v.view(1, -2, -3)
+v.view(-1, -2, 3)
+v.view(-1, 2, -3)
+END

data/examples/vector/where.rb

@@ -0,0 +1,29 @@
+#!/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 test Vector::Int
+z = GSL::Vector::Int[1, 0, 0, 5, 2, 0, 9]
+
+
+# #where without a block indices of non-zero elements.
+z.where
+
+# #where with a block returns indices of elements for which block returns true.
+z.where {|e| e >= 2}
+
+# #where2 without a block returns ["indices of non-zero elemens", "indices of
+# zero elements"]
+z.where2
+
+# #where2 with a block returns ["indices of elements for which block returned
+# true", "indices of elements for which block returned false"]
+z.where2 {|e| e >= 2}
+END

data/examples/vector/write.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 test Vector
+v = GSL::Vector.alloc([1, 2, 3, 4])
+
+# Write Vector in binary format to file a.dat
+v.fwrite("a.dat")
+
+# Create another test Vector
+v2 = GSL::Vector.alloc([5, 6, 7, 8])
+
+# Write Vector in text format to file b.dat
+v2.fprintf("b.dat")
+END
+

data/examples/vector/zip.rb

@@ -0,0 +1,34 @@
+#!/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 three test vectors of varying lengths
+a = GSL::Vector[0..4]
+b = GSL::Vector[2, 3, 4]
+c = GSL::Vector[5, 7, 4, 8, 9, 2]
+
+# Zip them together using Vector#zip
+a.zip(b, c)
+
+# Zip them together using GSL::Vector.zip
+GSL::Vector.zip(a, b, c)
+
+# Convert test Vectors to Vector::Complex
+aa = a.to_complex
+bb = b.to_complex
+cc = c.to_complex
+
+# Zip Vector::Complex objects together using Vector::Complex#zip
+aa.zip(bb, cc)
+
+
+# Zip Vector::Complex objects together using GSL::Vector::Complex.zip
+GSL::Vector::Complex.zip(aa, bb, cc)
+END

data/examples/wavelet/ecg.dat

@@ -0,0 +1,256 @@
+0.0462458
+0.0462458
+0.0512458
+0.0712458
+0.0712458
+0.0662458
+0.0962458
+0.1012458
+0.1162458
+0.1212458
+0.1162458
+0.1062458
+0.0912458
+0.1012458
+0.0962458
+0.0962458
+0.0962458
+0.0912458
+0.0862458
+0.0812458
+0.0862458
+0.1012458
+0.1112458
+0.1162458
+0.0762458
+0.0362458
+0.0362458
+0.0212458
+0.0112458
+-0.0087542
+-0.0087542
+-0.0037542
+0.0062458
+0.0012458
+0.0062458
+-0.0037542
+-0.0187542
+-0.0237542
+-0.0187542
+-0.0187542
+-0.0287542
+-0.0237542
+-0.0337542
+-0.0087542
+-0.0137542
+-0.0087542
+0.0012458
+-0.0237542
+-0.0337542
+-0.0187542
+-0.0087542
+-0.0037542
+-0.0087542
+-0.0287542
+-0.0437542
+-0.0387542
+-0.0587542
+-0.1037542
+-0.1237542
+-0.1537542
+-0.1887542
+-0.2137542
+-0.1837542
+-0.0937542
+0.0212458
+0.1612458
+0.3062458
+0.5562458
+0.8112458
+1.0412458
+1.1962458
+1.2612458
+1.2262458
+1.0762458
+0.8112458
+0.4862458
+0.2112458
+0.0512458
+-0.0687542
+-0.1287542
+-0.1537542
+-0.1337542
+-0.1037542
+-0.0687542
+-0.0687542
+-0.0637542
+-0.0687542
+-0.0587542
+-0.0587542
+-0.0587542
+-0.0737542
+-0.0637542
+-0.0637542
+-0.0637542
+-0.0537542
+-0.0737542
+-0.0887542
+-0.0887542
+-0.0787542
+-0.0737542
+-0.0687542
+-0.0837542
+-0.0737542
+-0.0637542
+-0.0537542
+-0.0687542
+-0.0687542
+-0.0837542
+-0.0887542
+-0.0887542
+-0.0687542
+-0.0687542
+-0.0737542
+-0.0837542
+-0.0937542
+-0.0787542
+-0.0887542
+-0.0837542
+-0.0887542
+-0.0937542
+-0.0887542
+-0.0787542
+-0.0787542
+-0.0737542
+-0.0687542
+-0.0837542
+-0.0887542
+-0.0687542
+-0.0687542
+-0.0637542
+-0.0637542
+-0.0887542
+-0.0837542
+-0.0737542
+-0.0687542
+-0.0537542
+-0.0687542
+-0.0737542
+-0.0887542
+-0.0787542
+-0.0687542
+-0.0687542
+-0.0637542
+-0.0837542
+-0.0937542
+-0.0937542
+-0.0787542
+-0.0737542
+-0.0837542
+-0.0937542
+-0.0987542
+-0.0987542
+-0.0887542
+-0.0937542
+-0.1037542
+-0.0987542
+-0.1137542
+-0.1087542
+-0.1087542
+-0.0987542
+-0.1087542
+-0.1287542
+-0.1337542
+-0.1287542
+-0.1137542
+-0.1237542
+-0.1287542
+-0.1337542
+-0.1487542
+-0.1387542
+-0.1337542
+-0.1287542
+-0.1337542
+-0.1487542
+-0.1537542
+-0.1387542
+-0.1287542
+-0.1237542
+-0.1187542
+-0.1137542
+-0.1187542
+-0.0887542
+-0.0737542
+-0.0487542
+-0.0437542
+-0.0387542
+-0.0437542
+-0.0187542
+-0.0037542
+0.0062458
+0.0012458
+-0.0087542
+-0.0087542
+0.0012458
+0.0112458
+0.0212458
+0.0212458
+0.0012458
+0.0012458
+0.0062458
+0.0162458
+0.0162458
+0.0262458
+0.0012458
+-0.0087542
+0.0162458
+0.0112458
+0.0212458
+0.0212458
+0.0012458
+-0.0037542
+0.0112458
+0.0162458
+0.0062458
+0.0162458
+0.0062458
+0.0062458
+0.0112458
+0.0262458
+0.0312458
+0.0162458
+0.0112458
+0.0012458
+0.0062458
+0.0212458
+0.0062458
+0.0062458
+0.0062458
+-0.0087542
+0.0062458
+0.0012458
+0.0062458
+-0.0037542
+-0.0137542
+-0.0187542
+-0.0137542
+-0.0137542
+-0.0087542
+-0.0037542
+-0.0237542
+-0.0287542
+-0.0237542
+-0.0137542
+-0.0087542
+-0.0087542
+-0.0237542
+-0.0237542
+-0.0237542
+0.0012458
+-0.0087542
+-0.0137542
+-0.0187542
+-0.0337542
+-0.0137542
+-0.0087542
+-0.0087542