rb-gsl 1.16.0.5 → 1.16.0.6

data/examples/eigen/vander.rb

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

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

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

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

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

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

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

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

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

@@ -1,256 +0,0 @@
-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