rb-gsl 1.15.3.1

data/examples/eigen/vander.rb

@@ -0,0 +1,41 @@
+#!/usr/bin/env ruby
+require("gsl")
+
+m = GSL::Matrix.vandermonde([1, 2, 3, 4])
+p m
+eval, evec = m.eigen_nonsymmv
+
+p eval.real
+p evec.real
+
+GSL::Eigen::nonsymmv_sort(eval, evec, GSL::Eigen::SORT_ABS_ASC)
+p eval.real
+p evec.real
+
+
+=begin
+Octave result:
+
+octave:1> m = vander([1 2 3 4])
+m =
+
+   1   1   1   1
+   8   4   2   1
+  27   9   3   1
+  64  16   4   1
+
+octave:2> [v, d] = eig(m)
+v =
+
+   0.106363   0.137820  -0.129196   0.052254
+   0.228569   0.036136   0.407586  -0.375568
+   0.462688  -0.289537   0.376326   0.794157
+   0.849919  -0.946503  -0.821925  -0.474902
+
+d =
+
+   15.48974    0.00000    0.00000    0.00000
+    0.00000   -7.70629    0.00000    0.00000
+    0.00000    0.00000    1.29422    0.00000
+    0.00000    0.00000    0.00000   -0.07768
+=end

data/examples/fft/fft.rb

@@ -0,0 +1,17 @@
+#!/usr/bin/env ruby
+require("gsl")
+
+N = 2048
+SAMPLING = 1000   # 1 kHz
+TMAX = 1.0/SAMPLING*N
+FREQ1 = 50
+FREQ2 = 120
+t = GSL::Vector.linspace(0, TMAX, N)
+x = GSL::Sf::sin(2*Math::PI*FREQ1*t) + GSL::Sf::sin(2*Math::PI*FREQ2*t)
+y = x.fft
+
+y2 = y.subvector(1, N-2).to_complex2
+mag = y2.abs
+phase = y2.arg
+f = GSL::Vector.linspace(0, SAMPLING/2, mag.size)
+GSL::graph(f, mag, "-C -g 3 -x 0 200 -X 'Frequency [Hz]'")

data/examples/fft/fft2.rb

@@ -0,0 +1,17 @@
+#!/usr/bin/env ruby
+require("gsl")
+
+N = 2048
+SAMPLING = 1000   # 1 kHz
+TMAX = 1.0/SAMPLING*N
+FREQ1 = 50
+FREQ2 = 120
+t = GSL::Vector.linspace(0, TMAX, N)
+x = GSL::Sf::sin(2*Math::PI*FREQ1*t) + GSL::Sf::sin(2*Math::PI*FREQ2*t)
+x.fft!
+
+y2 = x.subvector(1, N-2).to_complex2
+mag = y2.abs
+phase = y2.arg
+f = GSL::Vector.linspace(0, SAMPLING/2, mag.size)
+GSL::graph(f, mag, "-C -g 3 -x 0 200 -X 'Frequency [Hz]'")

data/examples/fft/forward.rb

@@ -0,0 +1,25 @@
+#!/usr/bin/env ruby
+require("gsl")
+
+n = 630
+data = GSL::Vector::Complex::alloc(n)
+
+table = GSL::FFT::ComplexWavetable.alloc(n)
+space = GSL::FFT::ComplexWorkspace.alloc(n)
+
+data[0] = 1.0
+for i in 1..10 do
+  data[i]   = GSL::Complex[1.0, 0.0]
+  data[n-i] = GSL::Complex[1.0, 0.0]
+end
+
+# Select whichever you like
+#ffted = data.forward(table, space)
+#ffted = data.forward(table)
+#ffted = data.forward(space)
+#ffted = data.forward()
+#ffted = data.transform(table, space, GSL::FFT::Forward)
+#ffted = data.transform(GSL::FFT::Forward)
+ffted = data.transform(GSL::FFT::Forward)
+ffted /= Math::sqrt(n)
+GSL::graph(nil, data.re, ffted.re, "-C -g 3")

data/examples/fft/forward2.rb

@@ -0,0 +1,26 @@
+#!/usr/bin/env ruby
+require("gsl")
+
+n = 630
+data = GSL::Vector::Complex.alloc(n)
+
+table = GSL::FFT::ComplexWavetable.alloc(n)
+space = GSL::FFT::ComplexWorkspace.alloc(n)
+
+data[0] = 1.0
+for i in 1..10 do
+  data[i] = GSL::Complex[1.0, 0.0]
+  data[n-i] = GSL::Complex[1.0, 0.0]
+end
+org = data.clone
+
+# Select whichever you like
+#data.forward!(table, space)
+#data.forward!(table)
+#data.forward!(space)
+#data.forward!()
+#data.transform!(table, space, GSL::FFT::Forward)
+#data.transform!(GSL::FFT::Forward)
+data.forward!(table, space)
+data /= Math::sqrt(n)
+GSL::graph(nil, org.re, data.re, "-C -g 3")

data/examples/fft/radix2.rb

@@ -0,0 +1,18 @@
+#!/usr/bin/env ruby
+require("gsl") 
+
+n = 128
+data = GSL::Vector::Complex.alloc(n)
+
+data[0] = 1.0
+c = GSL::Complex[1, 0]
+for i in 1..10 do
+  data[i] = c
+  data[n-i] = c
+end
+
+ffted = data.radix2_forward
+ffted /= Math::sqrt(n)
+
+GSL::graph(nil, data.re, ffted.re, "-T X -C -g 3 -L 'Radix-2' -x 0 #{data.size/2}")
+

data/examples/fft/real-halfcomplex.rb

@@ -0,0 +1,33 @@
+#!/usr/bin/env ruby
+require("gsl")
+
+n = 100
+data = GSL::Vector.alloc(n)
+#data = NArray.float(n)
+
+for i in (n/3)...(2*n/3) do
+  data[i] = 1.0
+end
+
+rtable = GSL::FFT::RealWavetable.alloc(n)
+rwork = GSL::FFT::RealWorkspace.alloc(n)
+  
+#ffted = data.real_transform(rtable, rwork)
+#ffted = data.real_transform(rtable)
+#ffted = data.real_transform(rwork)
+#ffted = data.real_transform()
+#ffted = data.fft
+ffted = data.real_transform()
+
+for i in 11...n do
+  ffted[i] = 0.0
+end
+  
+hctable = GSL::FFT::HalfComplexWavetable.alloc(n)
+  
+#data2 = ffted.halfcomplex_inverse(hctable, rwork)
+#data2 = ffted.halfcomplex_inverse()
+#data2 = ffted.ifft
+data2 = ffted.halfcomplex_inverse()
+
+GSL::graph(nil, data, data2, "-T X -C -g 3 -L 'Real-halfcomplex' -x 0 #{data.size}")

data/examples/fft/real-halfcomplex2.rb

@@ -0,0 +1,30 @@
+#!/usr/bin/env ruby
+require("gsl")
+
+n = 100
+data = GSL::Vector.alloc(n)
+#data = NArray.float(n)
+
+for i in (n/3)...(2*n/3) do
+  data[i] = 1.0
+end
+org = data.clone
+
+rtable = GSL::FFT::RealWavetable.alloc(n)
+rwork = GSL::FFT::RealWorkspace.alloc(n)
+  
+data.real_transform!(rtable, rwork)
+
+for i in 11...n do
+  data[i] = 0.0
+end
+  
+hctable = GSL::FFT::HalfComplexWavetable.alloc(n)
+  
+#data.halfcomplex_inverse!(hctable, rwork)
+#data.halfcomplex_inverse!()
+#data.ifft!
+#data.halfcomplex_inverse!(hctable, rwork)
+data.ifft!(hctable, rwork)
+
+GSL::graph(nil, org, data, "-T X -C -g 3 -L 'Real-halfcomplex' -x 0 #{data.size}")

data/examples/fft/realradix2.rb

@@ -0,0 +1,19 @@
+#!/usr/bin/env ruby
+require("gsl")
+
+n = 128
+data = GSL::Vector.alloc(n)
+#data = NArray.float(n)
+
+data[0] = 1.0
+for i in 1..10 do
+  data[i] = 1.0
+  data[n-i] = 1.0
+end
+
+#ffted = data.radix2_transform(GSL::FFT:FORWARD)
+#ffted = data.real_radix2_transform(GSL::FFT::FORWARD)
+# etc...
+ffted = data.radix2_transform(GSL::FFT::FORWARD)
+ffted /= Math::sqrt(n)
+GSL::graph(nil, data, ffted, "-T X -C -g 3 -L 'Real Radix-2' -x 0 #{data.size}")

data/examples/fft/sunspot.dat

@@ -0,0 +1,256 @@
+1716   47.0
+1717   63.0
+1718   60.0
+1719   39.0
+1720   28.0
+1721   26.0
+1722   22.0
+1723   11.0
+1724   21.0
+1725   40.0
+1726   78.0
+1727  122.0
+1728  103.0
+1729   73.0
+1730   47.0
+1731   35.0
+1732   11.0
+1733    5.0
+1734   16.0
+1735   34.0
+1736   70.0
+1737   81.0
+1738  111.0
+1739  101.0
+1740   73.0
+1741   40.0
+1742   20.0
+1743   16.0
+1744    5.0
+1745   11.0
+1746   22.0
+1747   40.0
+1748   60.0
+1749   80.9
+1750   83.4
+1751   47.7
+1752   47.8
+1753   30.7
+1754   12.2
+1755    9.6
+1756   10.2
+1757   32.4
+1758   47.6
+1759   54.0
+1760   62.9
+1761   85.9
+1762   61.2
+1763   45.1
+1764   36.4
+1765   20.9
+1766   11.4
+1767   37.8
+1768   69.8
+1769  106.1
+1770  100.8
+1771   81.6
+1772   66.5
+1773   34.8
+1774   30.6
+1775    7.0
+1776   19.8
+1777   92.5
+1778  154.4
+1779  125.9
+1780   84.8
+1781   68.1
+1782   38.5
+1783   22.8
+1784   10.2
+1785   24.1
+1786   82.9
+1787  132.0
+1788  130.9
+1789  118.1
+1790   89.9
+1791   66.6
+1792   60.0
+1793   46.9
+1794   41.0
+1795   21.3
+1796   16.0
+1797    6.4
+1798    4.1
+1799    6.8
+1800   14.5
+1801   34.0
+1802   45.0
+1803   43.1
+1804   47.5
+1805   42.2
+1806   28.1
+1807   10.1
+1808    8.1
+1809    2.5
+1810    0.0
+1811    1.4
+1812    5.0
+1813   12.2
+1814   13.9
+1815   35.4
+1816   45.8
+1817   41.1
+1818   30.1
+1819   23.9
+1820   15.6
+1821    6.6
+1822    4.0
+1823    1.8
+1824    8.5
+1825   16.6
+1826   36.3
+1827   49.6
+1828   64.2
+1829   67.0
+1830   70.9
+1831   47.8
+1832   27.5
+1833    8.5
+1834   13.2
+1835   56.9
+1836  121.5
+1837  138.3
+1838  103.2
+1839   85.7
+1840   64.6
+1841   36.7
+1842   24.2
+1843   10.7
+1844   15.0
+1845   40.1
+1846   61.5
+1847   98.5
+1848  124.7
+1849   96.3
+1850   66.6
+1851   64.5
+1852   54.1
+1853   39.0
+1854   20.6
+1855    6.7
+1856    4.3
+1857   22.7
+1858   54.8
+1859   93.8
+1860   95.8
+1861   77.2
+1862   59.1
+1863   44.0
+1864   47.0
+1865   30.5
+1866   16.3
+1867    7.3
+1868   37.6
+1869   74.0
+1870  139.0
+1871  111.2
+1872  101.6
+1873   66.2
+1874   44.7
+1875   17.0
+1876   11.3
+1877   12.4
+1878    3.4
+1879    6.0
+1880   32.3
+1881   54.3
+1882   59.7
+1883   63.7
+1884   63.5
+1885   52.2
+1886   25.4
+1887   13.1
+1888    6.8
+1889    6.3
+1890    7.1
+1891   35.6
+1892   73.0
+1893   85.1
+1894   78.0
+1895   64.0
+1896   41.8
+1897   26.2
+1898   26.7
+1899   12.1
+1900    9.5
+1901    2.7
+1902    5.0
+1903   24.4
+1904   42.0
+1905   63.5
+1906   53.8
+1907   62.0
+1908   48.5
+1909   43.9
+1910   18.6
+1911    5.7
+1912    3.6
+1913    1.4
+1914    9.6
+1915   47.4
+1916   57.1
+1917  103.9
+1918   80.6
+1919   63.6
+1920   37.6
+1921   26.1
+1922   14.2
+1923    5.8
+1924   16.7
+1925   44.3
+1926   63.9
+1927   69.0
+1928   77.8
+1929   64.9
+1930   35.7
+1931   21.2
+1932   11.1
+1933    5.7
+1934    8.7
+1935   36.1
+1936   79.7
+1937  114.4
+1938  109.6
+1939   88.8
+1940   67.8
+1941   47.5
+1942   30.6
+1943   16.3
+1944    9.6
+1945   33.2
+1946   92.6
+1947  151.6
+1948  136.3
+1949  134.7
+1950   83.9
+1951   69.4
+1952   31.5
+1953   13.9
+1954    4.4
+1955   38.0
+1956  141.7
+1957  190.2
+1958  184.8
+1959  159.0
+1960  112.3
+1961   53.9
+1962   37.5
+1963   27.9
+1964   10.2
+1965   15.1
+1966   47.0
+1967   93.8
+1968  105.9
+1969  105.5
+1970  104.5
+1971   66.6