rb-gsl 1.16.0.5 → 1.16.0.6

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (612) hide show
  1. checksums.yaml +4 -4
  2. data/rb-gsl.gemspec +5 -10
  3. metadata +10 -713
  4. data/.gitignore +0 -12
  5. data/.travis.yml +0 -24
  6. data/AUTHORS +0 -12
  7. data/COPYING +0 -341
  8. data/ChangeLog +0 -621
  9. data/Gemfile +0 -4
  10. data/README.md +0 -77
  11. data/Rakefile +0 -20
  12. data/THANKS +0 -17
  13. data/examples/alf/alf.gp +0 -15
  14. data/examples/alf/alf.rb +0 -32
  15. data/examples/blas/blas.rb +0 -13
  16. data/examples/blas/dnrm2.rb +0 -16
  17. data/examples/blas/level1.rb +0 -81
  18. data/examples/blas/level2.rb +0 -11
  19. data/examples/blas/level3.rb +0 -12
  20. data/examples/bspline.rb +0 -57
  21. data/examples/cdf.rb +0 -16
  22. data/examples/cheb.rb +0 -21
  23. data/examples/combination.rb +0 -23
  24. data/examples/complex/RC-lpf.rb +0 -47
  25. data/examples/complex/add.rb +0 -36
  26. data/examples/complex/coerce.rb +0 -14
  27. data/examples/complex/complex.rb +0 -25
  28. data/examples/complex/fpmi.rb +0 -70
  29. data/examples/complex/functions.rb +0 -77
  30. data/examples/complex/michelson.rb +0 -36
  31. data/examples/complex/mul.rb +0 -28
  32. data/examples/complex/oscillator.rb +0 -17
  33. data/examples/complex/set.rb +0 -37
  34. data/examples/const/physconst.rb +0 -151
  35. data/examples/const/travel.rb +0 -45
  36. data/examples/deriv/demo.rb +0 -13
  37. data/examples/deriv/deriv.rb +0 -36
  38. data/examples/deriv/diff.rb +0 -35
  39. data/examples/dht.rb +0 -42
  40. data/examples/dirac.rb +0 -56
  41. data/examples/eigen/eigen.rb +0 -34
  42. data/examples/eigen/herm.rb +0 -22
  43. data/examples/eigen/narray.rb +0 -9
  44. data/examples/eigen/nonsymm.rb +0 -37
  45. data/examples/eigen/nonsymmv.rb +0 -43
  46. data/examples/eigen/qhoscillator.gp +0 -35
  47. data/examples/eigen/qhoscillator.rb +0 -90
  48. data/examples/eigen/vander.rb +0 -41
  49. data/examples/fft/fft.rb +0 -17
  50. data/examples/fft/fft2.rb +0 -17
  51. data/examples/fft/forward.rb +0 -25
  52. data/examples/fft/forward2.rb +0 -26
  53. data/examples/fft/radix2.rb +0 -18
  54. data/examples/fft/real-halfcomplex.rb +0 -33
  55. data/examples/fft/real-halfcomplex2.rb +0 -30
  56. data/examples/fft/realradix2.rb +0 -19
  57. data/examples/fft/sunspot.dat +0 -256
  58. data/examples/fft/sunspot.rb +0 -16
  59. data/examples/fit/expdata.dat +0 -20
  60. data/examples/fit/expfit.rb +0 -31
  61. data/examples/fit/gaussfit.rb +0 -29
  62. data/examples/fit/gaussian_2peaks.rb +0 -34
  63. data/examples/fit/hillfit.rb +0 -40
  64. data/examples/fit/lognormal.rb +0 -26
  65. data/examples/fit/lorentzfit.rb +0 -22
  66. data/examples/fit/multifit.rb +0 -72
  67. data/examples/fit/ndlinear.rb +0 -133
  68. data/examples/fit/nonlinearfit.rb +0 -89
  69. data/examples/fit/plot.gp +0 -36
  70. data/examples/fit/polyfit.rb +0 -9
  71. data/examples/fit/powerfit.rb +0 -21
  72. data/examples/fit/sigmoidfit.rb +0 -40
  73. data/examples/fit/sinfit.rb +0 -22
  74. data/examples/fit/wlinear.rb +0 -46
  75. data/examples/fresnel.rb +0 -11
  76. data/examples/function/function.rb +0 -36
  77. data/examples/function/log.rb +0 -7
  78. data/examples/function/min.rb +0 -33
  79. data/examples/function/sin.rb +0 -10
  80. data/examples/function/synchrotron.rb +0 -18
  81. data/examples/gallery/butterfly.rb +0 -7
  82. data/examples/gallery/cayley.rb +0 -12
  83. data/examples/gallery/cornu.rb +0 -23
  84. data/examples/gallery/eight.rb +0 -11
  85. data/examples/gallery/koch.rb +0 -40
  86. data/examples/gallery/lemniscate.rb +0 -11
  87. data/examples/gallery/polar.rb +0 -11
  88. data/examples/gallery/rgplot/cossin.rb +0 -35
  89. data/examples/gallery/rgplot/rgplot.replaced +0 -0
  90. data/examples/gallery/rgplot/roesller.rb +0 -55
  91. data/examples/gallery/roesller.rb +0 -39
  92. data/examples/gallery/scarabaeus.rb +0 -14
  93. data/examples/histogram/cauchy.rb +0 -27
  94. data/examples/histogram/cauchy.sh +0 -2
  95. data/examples/histogram/exponential.rb +0 -19
  96. data/examples/histogram/gauss.rb +0 -16
  97. data/examples/histogram/gsl-histogram.rb +0 -40
  98. data/examples/histogram/histo2d.rb +0 -31
  99. data/examples/histogram/histo3d.rb +0 -34
  100. data/examples/histogram/histogram-pdf.rb +0 -27
  101. data/examples/histogram/histogram.rb +0 -26
  102. data/examples/histogram/integral.rb +0 -28
  103. data/examples/histogram/poisson.rb +0 -27
  104. data/examples/histogram/power.rb +0 -25
  105. data/examples/histogram/rebin.rb +0 -17
  106. data/examples/histogram/smp.dat +0 -5
  107. data/examples/histogram/xexp.rb +0 -21
  108. data/examples/integration/ahmed.rb +0 -21
  109. data/examples/integration/cosmology.rb +0 -75
  110. data/examples/integration/friedmann.gp +0 -16
  111. data/examples/integration/friedmann.rb +0 -35
  112. data/examples/integration/gamma-zeta.rb +0 -35
  113. data/examples/integration/integration.rb +0 -22
  114. data/examples/integration/qag.rb +0 -8
  115. data/examples/integration/qag2.rb +0 -14
  116. data/examples/integration/qag3.rb +0 -8
  117. data/examples/integration/qagi.rb +0 -28
  118. data/examples/integration/qagi2.rb +0 -49
  119. data/examples/integration/qagiu.rb +0 -29
  120. data/examples/integration/qagp.rb +0 -20
  121. data/examples/integration/qags.rb +0 -14
  122. data/examples/integration/qawc.rb +0 -18
  123. data/examples/integration/qawf.rb +0 -41
  124. data/examples/integration/qawo.rb +0 -29
  125. data/examples/integration/qaws.rb +0 -30
  126. data/examples/integration/qng.rb +0 -17
  127. data/examples/interp/demo.gp +0 -20
  128. data/examples/interp/demo.rb +0 -45
  129. data/examples/interp/interp.rb +0 -37
  130. data/examples/interp/points +0 -10
  131. data/examples/interp/spline.rb +0 -20
  132. data/examples/jacobi/deriv.rb +0 -40
  133. data/examples/jacobi/integrate.rb +0 -34
  134. data/examples/jacobi/interp.rb +0 -43
  135. data/examples/jacobi/jacobi.rb +0 -11
  136. data/examples/linalg/HH.rb +0 -15
  137. data/examples/linalg/HH_narray.rb +0 -13
  138. data/examples/linalg/LQ_solve.rb +0 -73
  139. data/examples/linalg/LU.rb +0 -84
  140. data/examples/linalg/LU2.rb +0 -31
  141. data/examples/linalg/LU_narray.rb +0 -24
  142. data/examples/linalg/PTLQ.rb +0 -47
  143. data/examples/linalg/QR.rb +0 -18
  144. data/examples/linalg/QRPT.rb +0 -47
  145. data/examples/linalg/QR_solve.rb +0 -78
  146. data/examples/linalg/QR_solve_narray.rb +0 -13
  147. data/examples/linalg/SV.rb +0 -16
  148. data/examples/linalg/SV_narray.rb +0 -12
  149. data/examples/linalg/SV_solve.rb +0 -49
  150. data/examples/linalg/chol.rb +0 -29
  151. data/examples/linalg/chol_narray.rb +0 -15
  152. data/examples/linalg/complex.rb +0 -57
  153. data/examples/linalg/invert_narray.rb +0 -10
  154. data/examples/math/const.rb +0 -67
  155. data/examples/math/elementary.rb +0 -35
  156. data/examples/math/functions.rb +0 -41
  157. data/examples/math/inf_nan.rb +0 -34
  158. data/examples/math/minmax.rb +0 -22
  159. data/examples/math/power.rb +0 -18
  160. data/examples/math/test.rb +0 -31
  161. data/examples/matrix/a.dat +0 -0
  162. data/examples/matrix/add.rb +0 -45
  163. data/examples/matrix/b.dat +0 -4
  164. data/examples/matrix/cat.rb +0 -31
  165. data/examples/matrix/colvectors.rb +0 -24
  166. data/examples/matrix/complex.rb +0 -41
  167. data/examples/matrix/det.rb +0 -29
  168. data/examples/matrix/diagonal.rb +0 -23
  169. data/examples/matrix/get_all.rb +0 -159
  170. data/examples/matrix/hilbert.rb +0 -31
  171. data/examples/matrix/iterator.rb +0 -19
  172. data/examples/matrix/matrix.rb +0 -57
  173. data/examples/matrix/minmax.rb +0 -53
  174. data/examples/matrix/mul.rb +0 -39
  175. data/examples/matrix/rand.rb +0 -20
  176. data/examples/matrix/read.rb +0 -29
  177. data/examples/matrix/rowcol.rb +0 -47
  178. data/examples/matrix/set.rb +0 -41
  179. data/examples/matrix/set_all.rb +0 -100
  180. data/examples/matrix/view.rb +0 -32
  181. data/examples/matrix/view_all.rb +0 -148
  182. data/examples/matrix/write.rb +0 -23
  183. data/examples/min.rb +0 -29
  184. data/examples/monte/miser.rb +0 -47
  185. data/examples/monte/monte.rb +0 -47
  186. data/examples/monte/plain.rb +0 -47
  187. data/examples/monte/vegas.rb +0 -46
  188. data/examples/multimin/bundle.rb +0 -66
  189. data/examples/multimin/cqp.rb +0 -109
  190. data/examples/multimin/fdfminimizer.rb +0 -40
  191. data/examples/multimin/fminimizer.rb +0 -41
  192. data/examples/multiroot/demo.rb +0 -36
  193. data/examples/multiroot/fdfsolver.rb +0 -50
  194. data/examples/multiroot/fsolver.rb +0 -33
  195. data/examples/multiroot/fsolver2.rb +0 -32
  196. data/examples/multiroot/fsolver3.rb +0 -26
  197. data/examples/narray/histogram.rb +0 -14
  198. data/examples/narray/mandel.rb +0 -27
  199. data/examples/narray/narray.rb +0 -28
  200. data/examples/narray/narray2.rb +0 -44
  201. data/examples/narray/sf.rb +0 -26
  202. data/examples/ntuple/create.rb +0 -17
  203. data/examples/ntuple/project.rb +0 -31
  204. data/examples/odeiv/binarysystem.gp +0 -23
  205. data/examples/odeiv/binarysystem.rb +0 -104
  206. data/examples/odeiv/demo.gp +0 -24
  207. data/examples/odeiv/demo.rb +0 -69
  208. data/examples/odeiv/demo2.gp +0 -26
  209. data/examples/odeiv/duffing.rb +0 -45
  210. data/examples/odeiv/frei1.rb +0 -109
  211. data/examples/odeiv/frei2.rb +0 -76
  212. data/examples/odeiv/legendre.rb +0 -52
  213. data/examples/odeiv/odeiv.rb +0 -32
  214. data/examples/odeiv/odeiv2.rb +0 -45
  215. data/examples/odeiv/oscillator.rb +0 -42
  216. data/examples/odeiv/sedov.rb +0 -97
  217. data/examples/odeiv/whitedwarf.gp +0 -40
  218. data/examples/odeiv/whitedwarf.rb +0 -158
  219. data/examples/ool/conmin.rb +0 -100
  220. data/examples/ool/gencan.rb +0 -99
  221. data/examples/ool/pgrad.rb +0 -100
  222. data/examples/ool/spg.rb +0 -100
  223. data/examples/pdf/bernoulli.rb +0 -5
  224. data/examples/pdf/beta.rb +0 -7
  225. data/examples/pdf/binomiral.rb +0 -10
  226. data/examples/pdf/cauchy.rb +0 -6
  227. data/examples/pdf/chisq.rb +0 -8
  228. data/examples/pdf/exponential.rb +0 -7
  229. data/examples/pdf/exppow.rb +0 -6
  230. data/examples/pdf/fdist.rb +0 -7
  231. data/examples/pdf/flat.rb +0 -7
  232. data/examples/pdf/gamma.rb +0 -8
  233. data/examples/pdf/gauss-tail.rb +0 -5
  234. data/examples/pdf/gauss.rb +0 -6
  235. data/examples/pdf/geometric.rb +0 -5
  236. data/examples/pdf/gumbel.rb +0 -6
  237. data/examples/pdf/hypergeometric.rb +0 -11
  238. data/examples/pdf/landau.rb +0 -5
  239. data/examples/pdf/laplace.rb +0 -7
  240. data/examples/pdf/logarithmic.rb +0 -5
  241. data/examples/pdf/logistic.rb +0 -6
  242. data/examples/pdf/lognormal.rb +0 -6
  243. data/examples/pdf/neg-binomiral.rb +0 -10
  244. data/examples/pdf/pareto.rb +0 -7
  245. data/examples/pdf/pascal.rb +0 -10
  246. data/examples/pdf/poisson.rb +0 -5
  247. data/examples/pdf/rayleigh-tail.rb +0 -6
  248. data/examples/pdf/rayleigh.rb +0 -6
  249. data/examples/pdf/tdist.rb +0 -6
  250. data/examples/pdf/weibull.rb +0 -8
  251. data/examples/permutation/ex1.rb +0 -22
  252. data/examples/permutation/permutation.rb +0 -16
  253. data/examples/poly/bell.rb +0 -6
  254. data/examples/poly/bessel.rb +0 -6
  255. data/examples/poly/cheb.rb +0 -6
  256. data/examples/poly/cheb_II.rb +0 -6
  257. data/examples/poly/cubic.rb +0 -9
  258. data/examples/poly/demo.rb +0 -20
  259. data/examples/poly/eval.rb +0 -28
  260. data/examples/poly/eval_derivs.rb +0 -14
  261. data/examples/poly/fit.rb +0 -21
  262. data/examples/poly/hermite.rb +0 -6
  263. data/examples/poly/poly.rb +0 -13
  264. data/examples/poly/quadratic.rb +0 -25
  265. data/examples/random/diffusion.rb +0 -34
  266. data/examples/random/gaussian.rb +0 -9
  267. data/examples/random/generator.rb +0 -27
  268. data/examples/random/hdsobol.rb +0 -21
  269. data/examples/random/poisson.rb +0 -9
  270. data/examples/random/qrng.rb +0 -19
  271. data/examples/random/randomwalk.rb +0 -37
  272. data/examples/random/randomwalk2d.rb +0 -19
  273. data/examples/random/rayleigh.rb +0 -36
  274. data/examples/random/rng.rb +0 -33
  275. data/examples/random/rngextra.rb +0 -14
  276. data/examples/roots/bisection.rb +0 -25
  277. data/examples/roots/brent.rb +0 -43
  278. data/examples/roots/demo.rb +0 -30
  279. data/examples/roots/newton.rb +0 -46
  280. data/examples/roots/recombination.gp +0 -11
  281. data/examples/roots/recombination.rb +0 -61
  282. data/examples/roots/steffenson.rb +0 -48
  283. data/examples/sf/ShiChi.rb +0 -6
  284. data/examples/sf/SiCi.rb +0 -6
  285. data/examples/sf/airy_Ai.rb +0 -8
  286. data/examples/sf/airy_Bi.rb +0 -8
  287. data/examples/sf/bessel_IK.rb +0 -12
  288. data/examples/sf/bessel_JY.rb +0 -13
  289. data/examples/sf/beta_inc.rb +0 -9
  290. data/examples/sf/clausen.rb +0 -6
  291. data/examples/sf/dawson.rb +0 -5
  292. data/examples/sf/debye.rb +0 -9
  293. data/examples/sf/dilog.rb +0 -6
  294. data/examples/sf/ellint.rb +0 -6
  295. data/examples/sf/expint.rb +0 -8
  296. data/examples/sf/fermi.rb +0 -10
  297. data/examples/sf/gamma_inc_P.rb +0 -9
  298. data/examples/sf/gegenbauer.rb +0 -8
  299. data/examples/sf/hyperg.rb +0 -7
  300. data/examples/sf/laguerre.rb +0 -19
  301. data/examples/sf/lambertW.rb +0 -5
  302. data/examples/sf/legendre_P.rb +0 -10
  303. data/examples/sf/lngamma.rb +0 -5
  304. data/examples/sf/psi.rb +0 -54
  305. data/examples/sf/sphbessel.gp +0 -27
  306. data/examples/sf/sphbessel.rb +0 -30
  307. data/examples/sf/synchrotron.rb +0 -5
  308. data/examples/sf/transport.rb +0 -10
  309. data/examples/sf/zetam1.rb +0 -5
  310. data/examples/siman.rb +0 -44
  311. data/examples/sort/heapsort.rb +0 -23
  312. data/examples/sort/heapsort_vector_complex.rb +0 -21
  313. data/examples/sort/sort.rb +0 -23
  314. data/examples/sort/sort2.rb +0 -16
  315. data/examples/stats/mean.rb +0 -17
  316. data/examples/stats/statistics.rb +0 -18
  317. data/examples/stats/test.rb +0 -9
  318. data/examples/sum.rb +0 -34
  319. data/examples/tamu_anova.rb +0 -18
  320. data/examples/vector/a.dat +0 -0
  321. data/examples/vector/add.rb +0 -56
  322. data/examples/vector/b.dat +0 -4
  323. data/examples/vector/c.dat +0 -3
  324. data/examples/vector/collect.rb +0 -26
  325. data/examples/vector/compare.rb +0 -28
  326. data/examples/vector/complex.rb +0 -51
  327. data/examples/vector/complex_get_all.rb +0 -85
  328. data/examples/vector/complex_set_all.rb +0 -131
  329. data/examples/vector/complex_view_all.rb +0 -77
  330. data/examples/vector/connect.rb +0 -22
  331. data/examples/vector/decimate.rb +0 -38
  332. data/examples/vector/diff.rb +0 -31
  333. data/examples/vector/filescan.rb +0 -17
  334. data/examples/vector/floor.rb +0 -23
  335. data/examples/vector/get_all.rb +0 -82
  336. data/examples/vector/gnuplot.rb +0 -38
  337. data/examples/vector/graph.rb +0 -28
  338. data/examples/vector/histogram.rb +0 -22
  339. data/examples/vector/linspace.rb +0 -24
  340. data/examples/vector/log.rb +0 -17
  341. data/examples/vector/logic.rb +0 -33
  342. data/examples/vector/logspace.rb +0 -25
  343. data/examples/vector/minmax.rb +0 -47
  344. data/examples/vector/mul.rb +0 -49
  345. data/examples/vector/narray.rb +0 -46
  346. data/examples/vector/read.rb +0 -29
  347. data/examples/vector/set.rb +0 -35
  348. data/examples/vector/set_all.rb +0 -121
  349. data/examples/vector/smpv.dat +0 -15
  350. data/examples/vector/test.rb +0 -43
  351. data/examples/vector/test_gslblock.rb +0 -58
  352. data/examples/vector/vector.rb +0 -110
  353. data/examples/vector/view.rb +0 -35
  354. data/examples/vector/view_all.rb +0 -73
  355. data/examples/vector/where.rb +0 -29
  356. data/examples/vector/write.rb +0 -24
  357. data/examples/vector/zip.rb +0 -34
  358. data/examples/wavelet/ecg.dat +0 -256
  359. data/examples/wavelet/wavelet1.rb +0 -50
  360. data/ext/gsl_native/alf.c +0 -206
  361. data/ext/gsl_native/array.c +0 -553
  362. data/ext/gsl_native/array_complex.c +0 -245
  363. data/ext/gsl_native/blas.c +0 -28
  364. data/ext/gsl_native/blas1.c +0 -733
  365. data/ext/gsl_native/blas2.c +0 -1088
  366. data/ext/gsl_native/blas3.c +0 -880
  367. data/ext/gsl_native/block.c +0 -40
  368. data/ext/gsl_native/block_source.h +0 -885
  369. data/ext/gsl_native/bspline.c +0 -122
  370. data/ext/gsl_native/bundle.c +0 -3
  371. data/ext/gsl_native/cdf.c +0 -740
  372. data/ext/gsl_native/cheb.c +0 -531
  373. data/ext/gsl_native/combination.c +0 -275
  374. data/ext/gsl_native/common.c +0 -293
  375. data/ext/gsl_native/complex.c +0 -1002
  376. data/ext/gsl_native/const.c +0 -331
  377. data/ext/gsl_native/const_additional.c +0 -99
  378. data/ext/gsl_native/cqp.c +0 -283
  379. data/ext/gsl_native/deriv.c +0 -187
  380. data/ext/gsl_native/dht.c +0 -353
  381. data/ext/gsl_native/diff.c +0 -164
  382. data/ext/gsl_native/dirac.c +0 -388
  383. data/ext/gsl_native/eigen.c +0 -2322
  384. data/ext/gsl_native/error.c +0 -193
  385. data/ext/gsl_native/extconf.rb +0 -118
  386. data/ext/gsl_native/fft.c +0 -1095
  387. data/ext/gsl_native/fit.c +0 -204
  388. data/ext/gsl_native/fresnel.c +0 -312
  389. data/ext/gsl_native/function.c +0 -518
  390. data/ext/gsl_native/geometry.c +0 -139
  391. data/ext/gsl_native/graph.c +0 -1590
  392. data/ext/gsl_native/gsl.c +0 -259
  393. data/ext/gsl_native/gsl_narray.c +0 -794
  394. data/ext/gsl_native/histogram.c +0 -1964
  395. data/ext/gsl_native/histogram2d.c +0 -1042
  396. data/ext/gsl_native/histogram3d.c +0 -884
  397. data/ext/gsl_native/histogram3d_source.c +0 -749
  398. data/ext/gsl_native/histogram_find.c +0 -99
  399. data/ext/gsl_native/histogram_oper.c +0 -150
  400. data/ext/gsl_native/ieee.c +0 -88
  401. data/ext/gsl_native/include/rb_gsl.h +0 -136
  402. data/ext/gsl_native/include/rb_gsl_array.h +0 -214
  403. data/ext/gsl_native/include/rb_gsl_cheb.h +0 -19
  404. data/ext/gsl_native/include/rb_gsl_common.h +0 -348
  405. data/ext/gsl_native/include/rb_gsl_complex.h +0 -25
  406. data/ext/gsl_native/include/rb_gsl_const.h +0 -23
  407. data/ext/gsl_native/include/rb_gsl_dirac.h +0 -6
  408. data/ext/gsl_native/include/rb_gsl_eigen.h +0 -17
  409. data/ext/gsl_native/include/rb_gsl_fft.h +0 -49
  410. data/ext/gsl_native/include/rb_gsl_fit.h +0 -23
  411. data/ext/gsl_native/include/rb_gsl_function.h +0 -23
  412. data/ext/gsl_native/include/rb_gsl_graph.h +0 -68
  413. data/ext/gsl_native/include/rb_gsl_histogram.h +0 -63
  414. data/ext/gsl_native/include/rb_gsl_histogram3d.h +0 -97
  415. data/ext/gsl_native/include/rb_gsl_integration.h +0 -17
  416. data/ext/gsl_native/include/rb_gsl_interp.h +0 -41
  417. data/ext/gsl_native/include/rb_gsl_linalg.h +0 -21
  418. data/ext/gsl_native/include/rb_gsl_math.h +0 -20
  419. data/ext/gsl_native/include/rb_gsl_odeiv.h +0 -18
  420. data/ext/gsl_native/include/rb_gsl_poly.h +0 -67
  421. data/ext/gsl_native/include/rb_gsl_rational.h +0 -30
  422. data/ext/gsl_native/include/rb_gsl_rng.h +0 -20
  423. data/ext/gsl_native/include/rb_gsl_root.h +0 -22
  424. data/ext/gsl_native/include/rb_gsl_sf.h +0 -110
  425. data/ext/gsl_native/include/rb_gsl_statistics.h +0 -17
  426. data/ext/gsl_native/include/rb_gsl_tensor.h +0 -43
  427. data/ext/gsl_native/include/rb_gsl_with_narray.h +0 -31
  428. data/ext/gsl_native/include/templates_off.h +0 -87
  429. data/ext/gsl_native/include/templates_on.h +0 -241
  430. data/ext/gsl_native/integration.c +0 -1154
  431. data/ext/gsl_native/interp.c +0 -499
  432. data/ext/gsl_native/jacobi.c +0 -733
  433. data/ext/gsl_native/linalg.c +0 -3915
  434. data/ext/gsl_native/linalg_complex.c +0 -726
  435. data/ext/gsl_native/math.c +0 -706
  436. data/ext/gsl_native/matrix.c +0 -36
  437. data/ext/gsl_native/matrix_complex.c +0 -1733
  438. data/ext/gsl_native/matrix_double.c +0 -557
  439. data/ext/gsl_native/matrix_int.c +0 -255
  440. data/ext/gsl_native/matrix_source.h +0 -2708
  441. data/ext/gsl_native/min.c +0 -219
  442. data/ext/gsl_native/monte.c +0 -978
  443. data/ext/gsl_native/multifit.c +0 -1862
  444. data/ext/gsl_native/multimin.c +0 -778
  445. data/ext/gsl_native/multimin_fsdf.c +0 -156
  446. data/ext/gsl_native/multiroots.c +0 -952
  447. data/ext/gsl_native/multiset.c +0 -210
  448. data/ext/gsl_native/ndlinear.c +0 -320
  449. data/ext/gsl_native/nmf.c +0 -171
  450. data/ext/gsl_native/nmf_wrap.c +0 -75
  451. data/ext/gsl_native/ntuple.c +0 -469
  452. data/ext/gsl_native/odeiv.c +0 -947
  453. data/ext/gsl_native/ool.c +0 -879
  454. data/ext/gsl_native/permutation.c +0 -598
  455. data/ext/gsl_native/poly.c +0 -39
  456. data/ext/gsl_native/poly2.c +0 -265
  457. data/ext/gsl_native/poly_source.h +0 -1871
  458. data/ext/gsl_native/qrng.c +0 -160
  459. data/ext/gsl_native/randist.c +0 -1848
  460. data/ext/gsl_native/rational.c +0 -480
  461. data/ext/gsl_native/rng.c +0 -595
  462. data/ext/gsl_native/root.c +0 -407
  463. data/ext/gsl_native/sf.c +0 -1446
  464. data/ext/gsl_native/sf_airy.c +0 -200
  465. data/ext/gsl_native/sf_bessel.c +0 -871
  466. data/ext/gsl_native/sf_clausen.c +0 -28
  467. data/ext/gsl_native/sf_coulomb.c +0 -206
  468. data/ext/gsl_native/sf_coupling.c +0 -121
  469. data/ext/gsl_native/sf_dawson.c +0 -29
  470. data/ext/gsl_native/sf_debye.c +0 -148
  471. data/ext/gsl_native/sf_dilog.c +0 -43
  472. data/ext/gsl_native/sf_elementary.c +0 -46
  473. data/ext/gsl_native/sf_ellint.c +0 -206
  474. data/ext/gsl_native/sf_elljac.c +0 -30
  475. data/ext/gsl_native/sf_erfc.c +0 -89
  476. data/ext/gsl_native/sf_exp.c +0 -169
  477. data/ext/gsl_native/sf_expint.c +0 -201
  478. data/ext/gsl_native/sf_fermi_dirac.c +0 -148
  479. data/ext/gsl_native/sf_gamma.c +0 -343
  480. data/ext/gsl_native/sf_gegenbauer.c +0 -97
  481. data/ext/gsl_native/sf_hyperg.c +0 -203
  482. data/ext/gsl_native/sf_laguerre.c +0 -113
  483. data/ext/gsl_native/sf_lambert.c +0 -47
  484. data/ext/gsl_native/sf_legendre.c +0 -368
  485. data/ext/gsl_native/sf_log.c +0 -105
  486. data/ext/gsl_native/sf_mathieu.c +0 -235
  487. data/ext/gsl_native/sf_power.c +0 -47
  488. data/ext/gsl_native/sf_psi.c +0 -89
  489. data/ext/gsl_native/sf_synchrotron.c +0 -48
  490. data/ext/gsl_native/sf_transport.c +0 -76
  491. data/ext/gsl_native/sf_trigonometric.c +0 -210
  492. data/ext/gsl_native/sf_zeta.c +0 -115
  493. data/ext/gsl_native/signal.c +0 -303
  494. data/ext/gsl_native/siman.c +0 -713
  495. data/ext/gsl_native/sort.c +0 -207
  496. data/ext/gsl_native/spline.c +0 -377
  497. data/ext/gsl_native/stats.c +0 -787
  498. data/ext/gsl_native/sum.c +0 -168
  499. data/ext/gsl_native/tamu_anova.c +0 -56
  500. data/ext/gsl_native/tensor.c +0 -35
  501. data/ext/gsl_native/tensor_source.h +0 -1122
  502. data/ext/gsl_native/vector.c +0 -35
  503. data/ext/gsl_native/vector_complex.c +0 -2241
  504. data/ext/gsl_native/vector_double.c +0 -1433
  505. data/ext/gsl_native/vector_int.c +0 -202
  506. data/ext/gsl_native/vector_source.h +0 -3321
  507. data/ext/gsl_native/wavelet.c +0 -923
  508. data/lib/gsl.rb +0 -8
  509. data/lib/gsl/gnuplot.rb +0 -41
  510. data/lib/gsl/oper.rb +0 -43
  511. data/lib/gsl/version.rb +0 -3
  512. data/lib/ool.rb +0 -22
  513. data/lib/ool/conmin.rb +0 -30
  514. data/lib/rbgsl.rb +0 -1
  515. data/rdoc/alf.rdoc +0 -77
  516. data/rdoc/blas.rdoc +0 -269
  517. data/rdoc/bspline.rdoc +0 -42
  518. data/rdoc/changes.rdoc +0 -159
  519. data/rdoc/cheb.rdoc +0 -99
  520. data/rdoc/cholesky_complex.rdoc +0 -46
  521. data/rdoc/combi.rdoc +0 -125
  522. data/rdoc/complex.rdoc +0 -210
  523. data/rdoc/const.rdoc +0 -546
  524. data/rdoc/dht.rdoc +0 -122
  525. data/rdoc/diff.rdoc +0 -133
  526. data/rdoc/ehandling.rdoc +0 -50
  527. data/rdoc/eigen.rdoc +0 -401
  528. data/rdoc/fft.rdoc +0 -535
  529. data/rdoc/fit.rdoc +0 -284
  530. data/rdoc/function.rdoc +0 -94
  531. data/rdoc/graph.rdoc +0 -137
  532. data/rdoc/hist.rdoc +0 -409
  533. data/rdoc/hist2d.rdoc +0 -279
  534. data/rdoc/hist3d.rdoc +0 -112
  535. data/rdoc/integration.rdoc +0 -398
  536. data/rdoc/interp.rdoc +0 -231
  537. data/rdoc/intro.rdoc +0 -27
  538. data/rdoc/linalg.rdoc +0 -681
  539. data/rdoc/linalg_complex.rdoc +0 -88
  540. data/rdoc/math.rdoc +0 -276
  541. data/rdoc/matrix.rdoc +0 -1093
  542. data/rdoc/min.rdoc +0 -189
  543. data/rdoc/monte.rdoc +0 -234
  544. data/rdoc/multimin.rdoc +0 -312
  545. data/rdoc/multiroot.rdoc +0 -293
  546. data/rdoc/narray.rdoc +0 -177
  547. data/rdoc/ndlinear.rdoc +0 -250
  548. data/rdoc/nonlinearfit.rdoc +0 -348
  549. data/rdoc/ntuple.rdoc +0 -88
  550. data/rdoc/odeiv.rdoc +0 -378
  551. data/rdoc/perm.rdoc +0 -221
  552. data/rdoc/poly.rdoc +0 -335
  553. data/rdoc/qrng.rdoc +0 -90
  554. data/rdoc/randist.rdoc +0 -233
  555. data/rdoc/ref.rdoc +0 -93
  556. data/rdoc/rng.rdoc +0 -203
  557. data/rdoc/roots.rdoc +0 -305
  558. data/rdoc/sf.rdoc +0 -1622
  559. data/rdoc/siman.rdoc +0 -89
  560. data/rdoc/sort.rdoc +0 -94
  561. data/rdoc/start.rdoc +0 -16
  562. data/rdoc/stats.rdoc +0 -219
  563. data/rdoc/sum.rdoc +0 -65
  564. data/rdoc/tensor.rdoc +0 -251
  565. data/rdoc/tut.rdoc +0 -5
  566. data/rdoc/use.rdoc +0 -177
  567. data/rdoc/vector.rdoc +0 -1243
  568. data/rdoc/vector_complex.rdoc +0 -347
  569. data/rdoc/wavelet.rdoc +0 -218
  570. data/test/gsl/blas_test.rb +0 -79
  571. data/test/gsl/bspline_test.rb +0 -63
  572. data/test/gsl/cdf_test.rb +0 -1512
  573. data/test/gsl/cheb_test.rb +0 -80
  574. data/test/gsl/combination_test.rb +0 -100
  575. data/test/gsl/complex_test.rb +0 -20
  576. data/test/gsl/const_test.rb +0 -29
  577. data/test/gsl/deriv_test.rb +0 -62
  578. data/test/gsl/dht_test.rb +0 -79
  579. data/test/gsl/diff_test.rb +0 -53
  580. data/test/gsl/eigen_test.rb +0 -563
  581. data/test/gsl/err_test.rb +0 -23
  582. data/test/gsl/fit_test.rb +0 -101
  583. data/test/gsl/histo_test.rb +0 -14
  584. data/test/gsl/index_test.rb +0 -61
  585. data/test/gsl/integration_test.rb +0 -274
  586. data/test/gsl/interp_test.rb +0 -27
  587. data/test/gsl/linalg_test.rb +0 -463
  588. data/test/gsl/matrix_nmf_test.rb +0 -37
  589. data/test/gsl/matrix_test.rb +0 -98
  590. data/test/gsl/min_test.rb +0 -89
  591. data/test/gsl/monte_test.rb +0 -77
  592. data/test/gsl/multifit_test.rb +0 -753
  593. data/test/gsl/multimin_test.rb +0 -157
  594. data/test/gsl/multiroot_test.rb +0 -135
  595. data/test/gsl/multiset_test.rb +0 -52
  596. data/test/gsl/odeiv_test.rb +0 -275
  597. data/test/gsl/oper_test.rb +0 -98
  598. data/test/gsl/poly_test.rb +0 -338
  599. data/test/gsl/qrng_test.rb +0 -94
  600. data/test/gsl/quartic_test.rb +0 -28
  601. data/test/gsl/randist_test.rb +0 -122
  602. data/test/gsl/rng_test.rb +0 -303
  603. data/test/gsl/roots_test.rb +0 -78
  604. data/test/gsl/sf_test.rb +0 -2079
  605. data/test/gsl/stats_test.rb +0 -122
  606. data/test/gsl/sum_test.rb +0 -69
  607. data/test/gsl/tensor_test.rb +0 -396
  608. data/test/gsl/vector_test.rb +0 -223
  609. data/test/gsl/wavelet_test.rb +0 -130
  610. data/test/gsl_test.rb +0 -321
  611. data/test/test_helper.rb +0 -42
  612. data/uncrustify.cfg +0 -1693
@@ -1,221 +0,0 @@
1
- #
2
- # = Permutations
3
- # Contents:
4
- # 1. {Permuation allocations}[link:rdoc/perm_rdoc.html#label-Permuation+allocations]
5
- # 1. {Methods}[link:rdoc/perm_rdoc.html#label-Methods]
6
- # 1. {Accessing permutation elements}[link:rdoc/perm_rdoc.html#label-Accessing+permutation+elements]
7
- # 1. {Permuation properties}[link:rdoc/perm_rdoc.html#label-Permutation+properties]
8
- # 1. {Permuation functions}[link:rdoc/perm_rdoc.html#label-Permutation+functions]
9
- # 1. {Reading and writing permutations}[link:rdoc/perm_rdoc.html#label-Reading+and+writing+permutations]
10
- # 1. {Permutations in cyclic form}[link:rdoc/perm_rdoc.html#label-Permutations+in+cyclic+Form]
11
- # 1. {Applying Permutations}[link:rdoc/perm_rdoc.html#label-Applying+Permutations]
12
- #
13
- # == Permuation allocations
14
- # ---
15
- # * GSL::Permutation.alloc(n)
16
- #
17
- # These functions create a new permutation of size <tt>n</tt>.
18
- # The permutation is not initialized and its elements are undefined.
19
- # Use <tt>GSL::Permutation.calloc</tt> if you want to create a permutation
20
- # which is initialized to the identity.
21
- #
22
- # ---
23
- # * GSL::Permutation.calloc(n)
24
- #
25
- # This creates a new permutation of size <tt>n</tt> and initializes it to the identity.
26
- #
27
- # == Methods
28
- # ---
29
- # * GSL::Permutation#init()
30
- #
31
- # This initializes the permutation to the identity, i.e. (0,1,2,...,n-1).
32
- #
33
- # ---
34
- # * GSL::Permutation.memcpy(dest, src)
35
- #
36
- # This method copies the elements of the permutation <tt>src</tt>
37
- # into the permutation <tt>dest</tt>. The two permutations must have the same size.
38
- #
39
- # ---
40
- # * GSL::Permutation#clone
41
- #
42
- # This creates a new permutation with the same elements of <tt>self</tt>.
43
- #
44
- # === Accessing permutation elements
45
- #
46
- # ---
47
- # * GSL::Permutation#get(i)
48
- #
49
- # Returns the value of the <tt>i</tt>-th element of the permutation.
50
- #
51
- # ---
52
- # * GSL::Permutation#swap(i, j)
53
- #
54
- # This exchanges the <tt>i</tt>-th and <tt>j</tt>-th elements of the permutation.
55
- #
56
- # === Permutation properties
57
- # ---
58
- # * GSL::Permutation#size
59
- #
60
- # Returns the size of the permutation.
61
- # ---
62
- # * GSL::Permutation#valid
63
- #
64
- # This checks that the permutation <tt>self</tt> is valid.
65
- # The n elements should contain each of the numbers 0 .. n-1 once and only once.
66
- #
67
- # ---
68
- # * GSL::Permutation#valid?
69
- #
70
- # This returns true if the permutation <tt>self</tt> is valid, and false otherwise.
71
- #
72
- # === Permutation functions
73
- #
74
- # ---
75
- # * GSL::Permutation#reverse
76
- #
77
- # This reverses the elements of the permutation <tt>self</tt>.
78
- # ---
79
- # * GSL::Permutation#inverse
80
- #
81
- # This computes the inverse of the permutation <tt>self</tt>, and returns
82
- # as a new permutation.
83
- #
84
- # ---
85
- # * GSL::Permutation#next
86
- #
87
- # This method advances the permutation <tt>self</tt> to the next permutation in
88
- # lexicographic order and returns <tt>GSL::SUCCESS</tt>. If no further permutations
89
- # are available it returns <tt>GSL::FAILURE</tt> and leaves <tt>self</tt> unmodified.
90
- # Starting with the identity permutation and repeatedly applying this function
91
- # will iterate through all possible permutations of a given order.
92
- # ---
93
- # * GSL::Permutation#prev
94
- #
95
- # This method steps backwards from the permutation <tt>self</tt> to the previous
96
- # permutation in lexicographic order, returning <tt>GSL_SUCCESS</tt>.
97
- # If no previous permutation is available it returns <tt>GSL_FAILURE</tt>
98
- # and leaves <tt>self</tt> unmodified.
99
- #
100
- # === Reading and writing permutations
101
- # ---
102
- # * GSL::Permutation#fwrite(io)
103
- # * GSL::Permutation#fwrite(filename)
104
- # * GSL::Permutation#fread(io)
105
- # * GSL::Permutation#fread(filename)
106
- # * GSL::Permutation#fprintf(io, format = "%u\n")
107
- # * GSL::Permutation#fprintf(filename, format = "%u\n")
108
- # * GSL::Permutation#fscanf(io)
109
- # * GSL::Permutation#fscanf(filename)
110
- #
111
- #
112
- # === Permutations in cyclic Form
113
- # A permutation can be represented in both <tt>linear</tt> and
114
- # <tt>cyclic</tt> notations. The functions described in this section convert
115
- # between the two forms. The linear notation is an index mapping, and has
116
- # already been described above. The cyclic notation expresses a
117
- # permutation as a series of circular rearrangements of groups
118
- # of elements, or <tt>cycles</tt>.
119
- #
120
- # For example, under the cycle (1 2 3), 1 is replaced by 2, 2 is replaced
121
- # by 3 and 3 is replaced by 1 in a circular fashion. Cycles of different
122
- # sets of elements can be combined independently, for example (1 2 3) (4 5)
123
- # combines the cycle (1 2 3) with the cycle (4 5), which is an exchange of
124
- # elements 4 and 5. A cycle of length one represents an element which is
125
- # unchanged by the permutation and is referred to as a <tt>singleton</tt>.
126
- #
127
- # It can be shown that every permutation can be decomposed into combinations
128
- # of cycles. The decomposition is not unique, but can always be rearranged
129
- # into a standard <tt>canonical form</tt> by a reordering of elements.
130
- # The library uses the canonical form defined in Knuth's
131
- # <tt>Art of Computer Programming</tt> (Vol 1, 3rd Ed, 1997) Section 1.3.3, p.178.
132
- #
133
- # The procedure for obtaining the canonical form given by Knuth is,
134
- #
135
- #
136
- # 1. Write all singleton cycles explicitly
137
- # 1. Within each cycle, put the smallest number first
138
- # 1. Order the cycles in decreasing order of the first number in the cycle.
139
- #
140
- # For example, the linear representation (2 4 3 0 1) is represented as
141
- # (1 4) (0 2 3) in canonical form. The permutation corresponds to an
142
- # exchange of elements 1 and 4, and rotation of elements 0, 2 and 3.
143
- #
144
- # The important property of the canonical form is that it can be reconstructed
145
- # from the contents of each cycle without the brackets. In addition, by removing
146
- # the brackets it can be considered as a linear representation of a different
147
- # permutation. In the example given above the permutation (2 4 3 0 1) would
148
- # become (1 4 0 2 3). This mapping has many applications in the theory of
149
- # permutations.
150
- #
151
- # ---
152
- # * GSL::Permutation#linear_to_canonical
153
- # * GSL::Permutation#to_canonical
154
- #
155
- # Computes the canonical form of the permutation <tt>self</tt> and
156
- # returns it as a new <tt>GSL::Permutation</tt>.
157
- #
158
- # ---
159
- # * GSL::Permutation#canonical_to_linear
160
- # * GSL::Permutation#to_linear
161
- #
162
- # Converts a permutation <tt>self</tt> in canonical form back into linear
163
- # form and returns it as a new <tt>GSL::Permutation</tt>.
164
- #
165
- #
166
- # ---
167
- # * GSL::Permutation#inversions
168
- #
169
- # Counts the number of inversions in the permutation <tt>self</tt>.
170
- # An inversion is any pair of elements that are not in order.
171
- # For example, the permutation 2031 has three inversions, corresponding
172
- # to the pairs (2,0) (2,1) and (3,1).
173
- # The identity permutation has no inversions.
174
- #
175
- # ---
176
- # * GSL::Permutation#linear_cycles
177
- #
178
- # Counts the number of cycles in the permutation <tt>self</tt>,
179
- # given in linear form.
180
- #
181
- # ---
182
- # * GSL::Permutation#canonical_cycles
183
- #
184
- # Counts the number of cycles in the permutation <tt>self</tt>,
185
- # given in canonical form.
186
- #
187
- # == Applying Permutations
188
- # ---
189
- # * GSL::Permutation::permute(v)
190
- #
191
- # Applies the permutation <tt>self</tt> to the elements of the vector <tt>v</tt>,
192
- # considered as a row-vector acted on by a permutation matrix from the
193
- # right, v' = v P. The j-th column of the permutation matrix P is
194
- # given by the p_j-th column of the identity matrix.
195
- # The permutation <tt>self</tt> and the vector <tt>v</tt> must have the same length.
196
- # ---
197
- # * GSL::Permutation::permute_inverse(v)
198
- #
199
- # Applies the inverse of the permutation <tt>self</tt> to the elements of
200
- # the vector <tt>v</tt>, considered as a row-vector acted on by an inverse
201
- # permutation matrix from the right, v' = v P^T.
202
- # Note that for permutation matrices the inverse is the same as the
203
- # transpose. The j-th column of the permutation matrix P is given by
204
- # the p_j-th column of the identity matrix.
205
- # The permutation <tt>self</tt> and the vector <tt>v</tt> must have the same length.
206
- # ---
207
- # * GSL::Permutation.mul(pa, pb)
208
- #
209
- # Combines the two permutations <tt>pa</tt> and <tt>pb</tt> into a single
210
- # permutation <tt>p</tt> and returns it.
211
- # The permutation <tt>p</tt> is equivalent to applying <tt>pb</tt> first
212
- # and then <tt>pa</tt>.
213
- #
214
- #
215
- # {prev}[link:rdoc/matrix_rdoc.html]
216
- # {next}[link:rdoc/combi_rdoc.html]
217
- #
218
- # {Reference index}[link:rdoc/ref_rdoc.html]
219
- # {top}[link:index.html]
220
- #
221
- #
@@ -1,335 +0,0 @@
1
- #
2
- # = Polynomials
3
- # Contents:
4
- # 1. {Polynomial Evaluation}[link:rdoc/poly_rdoc.html#label-Polynomial+Evaluation]
5
- # 1. {Solving polynomial equations}[link:rdoc/poly_rdoc.html#label-Solving+polynomial+equations]
6
- # 1. {Quadratic Equations}[link:rdoc/poly_rdoc.html#label-Quadratic+Equations]
7
- # 1. {Cubic Equations}[link:rdoc/poly_rdoc.html#label-Cubic+Equations]
8
- # 1. {General Polynomial Equations}[link:rdoc/poly_rdoc.html#label-General+Polynomial+Equations]
9
- # 1. {GSL::Poly class}[link:rdoc/poly_rdoc.html#label-Poly+class]
10
- # 1. {Constructors}[link:rdoc/poly_rdoc.html#label-Constructors]
11
- # 1. {Methods}[link:rdoc/poly_rdoc.html#label-Instance+Methods]
12
- # 1. {Polynomial Fitting}[link:rdoc/poly_rdoc.html#label-Polynomial+fitting]
13
- # 1. {Divided-difference representations}[link:rdoc/poly_rdoc.html#label-Divided-difference+representations]
14
- # 1. {Extensions}[link:rdoc/poly_rdoc.html#label-Extensions]
15
- # 1. {Special Polynomials}[link:rdoc/poly_rdoc.html#label-Special+Polynomials]
16
- # 1. {Polynomial Operations}[link:rdoc/poly_rdoc.html#label-Polynomial+Operations]
17
- #
18
- # == Polynomial Evaluation
19
- # ---
20
- # * GSL::Poly.eval(c, x)
21
- #
22
- # Evaluates the polynomial <tt>c[0] + c[1]x + c[2]x^2 + ...</tt>.
23
- # The polynomial coefficients <tt>c</tt> can be an <tt>Array</tt>,
24
- # a <tt>GSL::Vector</tt>, or an <tt>NArray</tt>. The evaluation point <tt>x</tt>
25
- # is a <tt>Numeric</tt>, <tt>Array</tt>, <tt>GSL::Vector</tt> or <tt>NArray</tt>.
26
- # From GSL 1.11, <tt>x</tt> can be a complex number, and <tt>c</tt> can be a complex polynomial given by a <tt>GSL::Vector::Complex</tt> or an <tt>Array</tt>.
27
- #
28
- # Ex)
29
- # >> require("gsl")
30
- # => true
31
- # >> GSL::Poly.eval([1, 2, 3], 2)
32
- # => 17.0
33
- # >> GSL::Poly.eval(GSL::Vector[1, 2, 3], 2)
34
- # => 17.0
35
- # >> GSL::Poly.eval(NArray[1.0, 2, 3], 2)
36
- # => 17.0
37
- # >> GSL::Poly.eval([1, 2, 3], [1, 2, 3])
38
- # => [6.0, 17.0, 34.0]
39
- # >> GSL::Poly.eval([1, 2, 3], GSL::Vector[1, 2, 3])
40
- # => GSL::Vector
41
- # [ 6.000e+00 1.700e+01 3.400e+01 ]
42
- # >> GSL::Poly.eval([1, 2, 3], NArray[1.0, 2, 3])
43
- # => NArray.float(3):
44
- # [ 6.0, 17.0, 34.0 ]
45
- #
46
- # ---
47
- # * GSL::Poly.eval_derivs(c, x)
48
- # * GSL::Poly.eval_derivs(c, x, lenres)
49
- #
50
- # (GSL-1.13) Evaluate and return a polynomial and its derivatives. The output contains the values of d^k P/d x^k for the specified value of x starting with k = 0. The input polynomial <tt>c</tt> can be an <tt>Array</tt>, <tt>GSL::Poly</tt> or an <tt>NArray</tt>. If <tt>lenres</tt> is not given, <tt>lenres = LENGTH(c) + 1</tt> is used, therefore the last element of the output is 0.
51
- #
52
- # ---
53
- # * GSL::Poly#eval_derivs(x)
54
- # * GSL::Poly#eval_derivs(x, lenres)
55
- #
56
- # (GSL-1.13) Evaluate and return a polynomial and its derivatives. The output contains the values of d^k P/d x^k for the specified value of x starting with k = 0. If <tt>lenres</tt> is not given, <tt>lenres = LENGTH(self) + 1</tt> is used, therefore the last element of the output is 0.
57
- #
58
- # Ex.)
59
- # >> ary = [1, 2, 3]
60
- # => [1, 2, 3]
61
- # >> GSL::Poly.eval_derivs(ary, 1)
62
- # => [6.0, 8.0, 6.0, 0.0]
63
- # >> na = NArray[1.0, 2, 3]
64
- # => NArray.float(3):
65
- # [ 1.0, 2.0, 3.0 ]
66
- # >> GSL::Poly.eval_derivs(na, 1)
67
- # => NArray.float(4):
68
- # [ 6.0, 8.0, 6.0, 0.0 ]
69
- # >> poly = GSL::Poly[1.0, 2, 3]
70
- # => GSL::Poly
71
- # [ 1.000e+00 2.000e+00 3.000e+00 ]
72
- # >> GSL::Poly.eval_derivs(poly, 1)
73
- # => GSL::Poly
74
- # [ 6.000e+00 8.000e+00 6.000e+00 0.000e+00 ]
75
- # >> poly.eval_derivs(1)
76
- # => GSL::Poly
77
- # [ 6.000e+00 8.000e+00 6.000e+00 0.000e+00 ]
78
- # >> poly.eval_derivs(1, 3)
79
- # => GSL::Poly
80
- # [ 6.000e+00 8.000e+00 6.000e+00 ]
81
- #
82
- # == Solving polynomial equations
83
- # === Quadratic Equations
84
- # ---
85
- # * GSL::Poly::solve_quadratic(a, b, c)
86
- # * GSL::Poly::solve_quadratic([a, b, c])
87
- #
88
- # Find the real roots of the quadratic equation,
89
- # a x^2 + b x + c = 0
90
- # The coefficients are given by 3 numbers, or a Ruby array,
91
- # or a <tt>GSL::Vector</tt> object. The roots are returned as a <tt>GSL::Vector</tt>.
92
- #
93
- # * Ex: z^2 - 3z + 2 = 0
94
- # >> GSL::Poly::solve_quadratic(1, -3, 2)
95
- # => GSL::Vector:
96
- # [ 1.000e+00 2.000e+00 ]
97
- #
98
- #
99
- # ---
100
- # * GSL::Poly::complex_solve_quadratic(a, b, c)
101
- # * GSL::Poly::complex_solve_quadratic([a, b, c])
102
- #
103
- # Find the complex roots of the quadratic equation,
104
- # a z^2 + b z + z = 0
105
- # The coefficients are given by 3 numbers or a Ruby array, or a
106
- # <tt>GSL::Vector</tt>.
107
- # The roots are returned as a <tt>GSL::Vector::Complex</tt> of two elements.
108
- #
109
- # * Ex: z^2 - 3z + 2 = 0
110
- # >> require("gsl")
111
- # => true
112
- # >> GSL::Poly::complex_solve_quadratic(1, -3, 2)
113
- # [ [1.000e+00 0.000e+00] [2.000e+00 0.000e+00] ]
114
- # => #<GSL::Vector::Complex:0x764014>
115
- # >> GSL::Poly::complex_solve_quadratic(1, -3, 2).real <--- Real part
116
- # => GSL::Vector::View:
117
- # [ 1.000e+00 2.000e+00 ]
118
- #
119
- # === Cubic Equations
120
- # ---
121
- # * GSL::Poly::solve_cubic(same as solve_quadratic)
122
- #
123
- # This method finds the real roots of the cubic equation,
124
- # x^3 + a x^2 + b x + c = 0
125
- #
126
- # ---
127
- # * GSL::Poly::complex_solve_cubic(same as solve_cubic)
128
- #
129
- # This method finds the complex roots of the cubic equation,
130
- # z^3 + a z^2 + b z + c = 0
131
- #
132
- # === General Polynomial Equations
133
- # ---
134
- # * GSL::Poly::complex_solve(c0, c1, c2,,, )
135
- # * GSL::Poly::solve(c0, c1, c2,,, )
136
- #
137
- # Find the complex roots of the polynomial equation. Note that
138
- # the coefficients are given by "ascending" order.
139
- #
140
- # * Ex: x^2 - 3 x + 2 == 0
141
- # >> GSL::Poly::complex_solve(2, -3, 1) <--- different from Poly::quadratic_solve
142
- # [ [1.000e+00 0.000e+00] [2.000e+00 0.000e+00] ]
143
- # => #<GSL::Vector::Complex:0x75e614>
144
- #
145
- # == Poly class
146
- # This class expresses polynomials of arbitrary orders.
147
- #
148
- # === Constructors
149
- # ---
150
- # * GSL::Poly.alloc(c0, c1, c2, ....)
151
- # * GSL::Poly[c0, c1, c2, ....]
152
- #
153
- # This creates an instance of the <tt>GSL::Poly</tt> class,
154
- # which represents a polynomial
155
- # c0 + c1 x + c2 x^2 + ....
156
- # This class is derived from <tt>GSL::Vector</tt>.
157
- #
158
- # * Ex: x^2 - 3 x + 2
159
- # poly = GSL::Poly.alloc([2, -3, 1])
160
- #
161
- # === Instance Methods
162
- # ---
163
- # * GSL::Poly#eval(x)
164
- # * GSL::Poly#at(x)
165
- #
166
- # Evaluates the polynomial
167
- # <tt>c[0] + c[1] x + c[2] x^2 + ... + c[len-1] x^{len-1}</tt>
168
- # using Horner's method for stability. The argument <tt>x</tt> is a
169
- # <tt>Numeric</tt>, <tt>GSL::Vector, Matrix</tt> or an <tt>Array</tt>.
170
- #
171
- # ---
172
- # * GSL::Poly#solve_quadratic
173
- #
174
- # Solve the quadratic equation.
175
- #
176
- # * Ex: z^2 - 3 z + 2 = 0:
177
- # >> a = GSL::Poly[2, -3, 1]
178
- # => GSL::Poly:
179
- # [ 2.000e+00 -3.000e+00 1.000e+00 ]
180
- # >> a.solve_quadratic
181
- # => GSL::Vector:
182
- # [ 1.000e+00 2.000e+00 ]
183
- #
184
- # ---
185
- # * GSL::Poly#solve_cubic
186
- #
187
- # Solve the cubic equation.
188
- #
189
- # ---
190
- # * GSL::Poly#complex_solve
191
- # * GSL::Poly#solve
192
- # * GSL::Poly#roots
193
- #
194
- # These methods find the complex roots of the quadratic equation,
195
- # c0 + c1 z + c2 z^2 + .... = 0
196
- #
197
- # * Ex: z^2 - 3 z + 2 = 0:
198
- # >> a = GSL::Poly[2, -3, 1]
199
- # => GSL::Poly:
200
- # [ 2.000e+00 -3.000e+00 1.000e+00 ]
201
- # >> a.solve
202
- # [ [1.000e+00 0.000e+00] [2.000e+00 0.000e+00] ]
203
- # => #<GSL::Vector::Complex:0x35db28>
204
- #
205
- # == Polynomial fitting
206
- # ---
207
- # * GSL::Poly.fit(x, y, order)
208
- # * GSL::Poly.wfit(x, w, y, order)
209
- #
210
- # Finds the coefficient of a polynomial of order <tt>order</tt>
211
- # that fits the vector data (<tt>x, y</tt>) in a least-square sense.
212
- # This provides a higher-level interface to the method
213
- # {GSL::Multifit#linear}[link:rdoc/fit_rdoc.html] in a case of polynomial fitting.
214
- #
215
- # Example:
216
- # #!/usr/bin/env ruby
217
- # require("gsl")
218
- #
219
- # x = GSL::Vector[1, 2, 3, 4, 5]
220
- # y = GSL::Vector[5.5, 43.1, 128, 290.7, 498.4]
221
- # # The results are stored in a polynomial "coef"
222
- # coef, cov, chisq, status = Poly.fit(x, y, 3)
223
- #
224
- # x2 = GSL::Vector.linspace(1, 5, 20)
225
- # graph([x, y], [x2, coef.eval(x2)], "-C -g 3 -S 4")
226
- #
227
- # == Divided-difference representations
228
- #
229
- # ---
230
- # * GSL::Poly::dd_init(xa, ya)
231
- #
232
- # This method computes a divided-difference representation of the
233
- # interpolating polynomial for the points <tt>(xa, ya)</tt>.
234
- #
235
- # ---
236
- # * GSL::Poly::DividedDifference#eval(x)
237
- #
238
- # This method evaluates the polynomial stored in divided-difference form
239
- # <tt>self</tt> at the point <tt>x</tt>.
240
- #
241
- # ---
242
- # * GSL::Poly::DividedDifference#taylor(xp)
243
- #
244
- # This method converts the divided-difference representation of a polynomial
245
- # to a Taylor expansion. On output the Taylor coefficients of the polynomial
246
- # expanded about the point <tt>xp</tt> are returned.
247
- #
248
- # == Extensions
249
- # === Special Polynomials
250
- # ---
251
- # * GSL::Poly.hermite(n)
252
- #
253
- # This returns coefficients of the <tt>n</tt>-th order Hermite polynomial, <tt>H(x; n)</tt>.
254
- # For order of <tt>n</tt> >= 3, this method uses the recurrence relation
255
- # H(x; n+1) = 2 x H(x; n) - 2 n H(x; n-1)
256
- # * Ex:
257
- # >> GSL::Poly.hermite(2)
258
- # => GSL::Poly::Int:
259
- # [ -2 0 4 ] <----- 4x^2 - 2
260
- # >> GSL::Poly.hermite(5)
261
- # => GSL::Poly::Int:
262
- # [ 0 120 0 -160 0 32 ] <----- 32x^5 - 160x^3 + 120x
263
- # >> GSL::Poly.hermite(7)
264
- # => GSL::Poly::Int:
265
- # [ 0 -1680 0 3360 0 -1344 0 128 ]
266
- #
267
- # ---
268
- # * GSL::Poly.cheb(n)
269
- # * GSL::Poly.chebyshev(n)
270
- #
271
- # Return the coefficients of the <tt>n</tt>-th order Chebyshev polynomial, <tt>T(x; n</tt>.
272
- # For order of <tt>n</tt> >= 3, this method uses the recurrence relation
273
- # T(x; n+1) = 2 x T(x; n) - T(x; n-1)
274
- #
275
- # ---
276
- # * GSL::Poly.cheb_II(n)
277
- # * GSL::Poly.chebyshev_II(n)
278
- #
279
- # Return the coefficients of the <tt>n</tt>-th order Chebyshev polynomial of type II,
280
- # <tt>U(x; n</tt>.
281
- # U(x; n+1) = 2 x U(x; n) - U(x; n-1)
282
- #
283
- # ---
284
- # * GSL::Poly.bell(n)
285
- #
286
- # Bell polynomial
287
- #
288
- # ---
289
- # * GSL::Poly.bessel(n)
290
- #
291
- # Bessel polynomial
292
- #
293
- # ---
294
- # * GSL::Poly.laguerre(n)
295
- #
296
- # Retunrs the coefficients of the <tt>n</tt>-th order Laguerre polynomial
297
- # multiplied by n!.
298
- #
299
- # Ex:
300
- # rb(main):001:0> require("gsl")
301
- # => true
302
- # >> GSL::Poly.laguerre(0)
303
- # => GSL::Poly::Int:
304
- # [ 1 ] <--- 1
305
- # >> GSL::Poly.laguerre(1)
306
- # => GSL::Poly::Int:
307
- # [ 1 -1 ] <--- -x + 1
308
- # >> GSL::Poly.laguerre(2)
309
- # => GSL::Poly::Int:
310
- # [ 2 -4 1 ] <--- (x^2 - 4x + 2)/2!
311
- # >> GSL::Poly.laguerre(3)
312
- # => GSL::Poly::Int:
313
- # [ 6 -18 9 -1 ] <--- (-x^3 + 9x^2 - 18x + 6)/3!
314
- # >> GSL::Poly.laguerre(4)
315
- # => GSL::Poly::Int:
316
- # [ 24 -96 72 -16 1 ] <--- (x^4 - 16x^3 + 72x^2 - 96x + 24)/4!
317
- #
318
- # === Polynomial Operations
319
- # ---
320
- # * GSL::Poly#conv
321
- # * GSL::Poly#deconv
322
- # * GSL::Poly#reduce
323
- # * GSL::Poly#deriv
324
- # * GSL::Poly#integ
325
- # * GSL::Poly#compan
326
- #
327
- #
328
- # {prev}[link:rdoc/complex_rdoc.html]
329
- # {next}[link:rdoc/sf_rdoc.html]
330
- #
331
- # {Reference index}[link:rdoc/ref_rdoc.html]
332
- # {top}[link:index.html]
333
- #
334
- #
335
- #