gsl 1.15.3 → 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 (446) hide show
  1. checksums.yaml +7 -0
  2. data/.gitignore +12 -0
  3. data/.travis.yml +24 -0
  4. data/AUTHORS +10 -2
  5. data/COPYING +341 -339
  6. data/ChangeLog +612 -554
  7. data/Gemfile +4 -0
  8. data/README.md +77 -0
  9. data/Rakefile +14 -93
  10. data/THANKS +1 -1
  11. data/examples/blas/blas.rb +1 -1
  12. data/examples/bspline.rb +3 -3
  13. data/examples/complex/functions.rb +4 -4
  14. data/examples/complex/mul.rb +1 -1
  15. data/examples/const/physconst.rb +4 -4
  16. data/examples/const/travel.rb +4 -4
  17. data/examples/eigen/nonsymmv.rb +1 -1
  18. data/examples/eigen/qhoscillator.rb +3 -3
  19. data/examples/fft/radix2.rb +1 -1
  20. data/examples/fft/real-halfcomplex.rb +3 -3
  21. data/examples/fft/real-halfcomplex2.rb +3 -3
  22. data/examples/fit/expfit.rb +1 -1
  23. data/examples/fit/multifit.rb +1 -1
  24. data/examples/fit/ndlinear.rb +44 -44
  25. data/examples/fit/nonlinearfit.rb +1 -1
  26. data/examples/fit/wlinear.rb +3 -3
  27. data/examples/function/function.rb +1 -1
  28. data/examples/function/min.rb +1 -1
  29. data/examples/function/synchrotron.rb +2 -2
  30. data/examples/gallery/koch.rb +1 -1
  31. data/examples/histogram/cauchy.rb +2 -2
  32. data/examples/histogram/exponential.rb +1 -1
  33. data/examples/histogram/histo3d.rb +1 -1
  34. data/examples/histogram/histogram-pdf.rb +2 -2
  35. data/examples/histogram/xexp.rb +1 -1
  36. data/examples/integration/ahmed.rb +2 -2
  37. data/examples/integration/cosmology.rb +7 -7
  38. data/examples/integration/friedmann.rb +4 -4
  39. data/examples/integration/qng.rb +1 -1
  40. data/examples/interp/demo.rb +2 -2
  41. data/examples/linalg/LQ_solve.rb +1 -1
  42. data/examples/linalg/LU.rb +1 -1
  43. data/examples/linalg/LU2.rb +1 -1
  44. data/examples/linalg/LU_narray.rb +1 -1
  45. data/examples/linalg/PTLQ.rb +1 -1
  46. data/examples/linalg/QRPT.rb +1 -1
  47. data/examples/linalg/chol.rb +1 -1
  48. data/examples/linalg/chol_narray.rb +1 -1
  49. data/examples/linalg/complex.rb +1 -1
  50. data/examples/math/elementary.rb +1 -1
  51. data/examples/math/functions.rb +1 -1
  52. data/examples/math/inf_nan.rb +1 -1
  53. data/examples/math/minmax.rb +1 -1
  54. data/examples/math/power.rb +1 -1
  55. data/examples/math/test.rb +1 -1
  56. data/examples/min.rb +1 -1
  57. data/examples/multimin/bundle.rb +1 -1
  58. data/examples/multimin/cqp.rb +17 -17
  59. data/examples/multiroot/fsolver3.rb +1 -1
  60. data/examples/odeiv/binarysystem.rb +12 -12
  61. data/examples/odeiv/demo.rb +3 -3
  62. data/examples/odeiv/frei1.rb +7 -7
  63. data/examples/odeiv/frei2.rb +4 -4
  64. data/examples/odeiv/oscillator.rb +1 -1
  65. data/examples/odeiv/sedov.rb +3 -3
  66. data/examples/odeiv/whitedwarf.rb +11 -11
  67. data/examples/permutation/ex1.rb +2 -2
  68. data/examples/permutation/permutation.rb +1 -1
  69. data/examples/poly/demo.rb +1 -1
  70. data/examples/random/diffusion.rb +1 -1
  71. data/examples/random/generator.rb +2 -2
  72. data/examples/random/randomwalk.rb +3 -3
  73. data/examples/random/rng.rb +1 -1
  74. data/examples/roots/bisection.rb +1 -1
  75. data/examples/roots/brent.rb +1 -1
  76. data/examples/roots/demo.rb +1 -1
  77. data/examples/roots/newton.rb +2 -2
  78. data/examples/roots/recombination.gp +0 -1
  79. data/examples/sf/hyperg.rb +1 -1
  80. data/examples/sf/sphbessel.rb +1 -1
  81. data/examples/sort/sort.rb +1 -1
  82. data/examples/tamu_anova.rb +4 -4
  83. data/examples/vector/add.rb +1 -1
  84. data/examples/vector/decimate.rb +1 -1
  85. data/examples/vector/gnuplot.rb +8 -8
  86. data/examples/vector/vector.rb +2 -2
  87. data/examples/wavelet/wavelet1.rb +1 -1
  88. data/ext/{alf.c → gsl_native/alf.c} +10 -10
  89. data/ext/{array.c → gsl_native/array.c} +70 -159
  90. data/ext/{array_complex.c → gsl_native/array_complex.c} +63 -66
  91. data/ext/{blas.c → gsl_native/blas.c} +2 -3
  92. data/ext/{blas1.c → gsl_native/blas1.c} +35 -36
  93. data/ext/{blas2.c → gsl_native/blas2.c} +57 -62
  94. data/ext/{blas3.c → gsl_native/blas3.c} +57 -58
  95. data/ext/{block.c → gsl_native/block.c} +14 -18
  96. data/ext/{block_source.c → gsl_native/block_source.h} +110 -112
  97. data/ext/gsl_native/bspline.c +122 -0
  98. data/ext/{bundle.c → gsl_native/bundle.c} +0 -0
  99. data/ext/{cdf.c → gsl_native/cdf.c} +79 -93
  100. data/ext/{cheb.c → gsl_native/cheb.c} +78 -89
  101. data/ext/{combination.c → gsl_native/combination.c} +11 -19
  102. data/ext/{common.c → gsl_native/common.c} +9 -41
  103. data/ext/{complex.c → gsl_native/complex.c} +116 -118
  104. data/ext/gsl_native/const.c +331 -0
  105. data/ext/{const_additional.c → gsl_native/const_additional.c} +13 -34
  106. data/ext/gsl_native/cqp.c +283 -0
  107. data/ext/{deriv.c → gsl_native/deriv.c} +25 -33
  108. data/ext/{dht.c → gsl_native/dht.c} +23 -31
  109. data/ext/{diff.c → gsl_native/diff.c} +26 -28
  110. data/ext/{dirac.c → gsl_native/dirac.c} +45 -46
  111. data/ext/{eigen.c → gsl_native/eigen.c} +1044 -1095
  112. data/ext/{error.c → gsl_native/error.c} +18 -18
  113. data/ext/gsl_native/extconf.rb +118 -0
  114. data/ext/{fft.c → gsl_native/fft.c} +197 -204
  115. data/ext/{fit.c → gsl_native/fit.c} +17 -18
  116. data/ext/gsl_native/fresnel.c +312 -0
  117. data/ext/{function.c → gsl_native/function.c} +37 -43
  118. data/ext/{geometry.c → gsl_native/geometry.c} +16 -16
  119. data/ext/{graph.c → gsl_native/graph.c} +39 -89
  120. data/ext/{gsl.c → gsl_native/gsl.c} +12 -33
  121. data/ext/{gsl_narray.c → gsl_native/gsl_narray.c} +20 -30
  122. data/ext/{histogram.c → gsl_native/histogram.c} +133 -160
  123. data/ext/{histogram2d.c → gsl_native/histogram2d.c} +78 -104
  124. data/ext/{histogram3d.c → gsl_native/histogram3d.c} +76 -76
  125. data/ext/{histogram3d_source.c → gsl_native/histogram3d_source.c} +196 -197
  126. data/ext/{histogram_find.c → gsl_native/histogram_find.c} +32 -34
  127. data/ext/{histogram_oper.c → gsl_native/histogram_oper.c} +43 -52
  128. data/ext/{ieee.c → gsl_native/ieee.c} +9 -21
  129. data/{include → ext/gsl_native/include}/rb_gsl.h +4 -26
  130. data/{include → ext/gsl_native/include}/rb_gsl_array.h +15 -39
  131. data/{include → ext/gsl_native/include}/rb_gsl_cheb.h +0 -2
  132. data/{include → ext/gsl_native/include}/rb_gsl_common.h +61 -61
  133. data/{include → ext/gsl_native/include}/rb_gsl_complex.h +1 -1
  134. data/{include → ext/gsl_native/include}/rb_gsl_const.h +0 -6
  135. data/ext/gsl_native/include/rb_gsl_dirac.h +6 -0
  136. data/{include → ext/gsl_native/include}/rb_gsl_eigen.h +1 -1
  137. data/{include → ext/gsl_native/include}/rb_gsl_fft.h +0 -13
  138. data/{include → ext/gsl_native/include}/rb_gsl_fit.h +0 -2
  139. data/{include → ext/gsl_native/include}/rb_gsl_function.h +0 -4
  140. data/{include → ext/gsl_native/include}/rb_gsl_graph.h +2 -4
  141. data/{include → ext/gsl_native/include}/rb_gsl_histogram.h +8 -8
  142. data/{include → ext/gsl_native/include}/rb_gsl_histogram3d.h +50 -50
  143. data/{include → ext/gsl_native/include}/rb_gsl_integration.h +1 -1
  144. data/{include → ext/gsl_native/include}/rb_gsl_interp.h +0 -5
  145. data/{include → ext/gsl_native/include}/rb_gsl_linalg.h +2 -6
  146. data/{include → ext/gsl_native/include}/rb_gsl_math.h +0 -6
  147. data/{include → ext/gsl_native/include}/rb_gsl_odeiv.h +0 -3
  148. data/{include → ext/gsl_native/include}/rb_gsl_poly.h +3 -7
  149. data/{include → ext/gsl_native/include}/rb_gsl_rational.h +1 -8
  150. data/{include → ext/gsl_native/include}/rb_gsl_rng.h +0 -1
  151. data/{include → ext/gsl_native/include}/rb_gsl_root.h +1 -1
  152. data/{include → ext/gsl_native/include}/rb_gsl_sf.h +39 -48
  153. data/{include → ext/gsl_native/include}/rb_gsl_statistics.h +1 -1
  154. data/{include → ext/gsl_native/include}/rb_gsl_tensor.h +0 -2
  155. data/{include → ext/gsl_native/include}/rb_gsl_with_narray.h +3 -1
  156. data/{include → ext/gsl_native/include}/templates_off.h +0 -0
  157. data/{include → ext/gsl_native/include}/templates_on.h +1 -1
  158. data/ext/{integration.c → gsl_native/integration.c} +164 -189
  159. data/ext/{interp.c → gsl_native/interp.c} +25 -38
  160. data/ext/gsl_native/jacobi.c +733 -0
  161. data/ext/{linalg.c → gsl_native/linalg.c} +462 -589
  162. data/ext/{linalg_complex.c → gsl_native/linalg_complex.c} +93 -106
  163. data/ext/{math.c → gsl_native/math.c} +48 -67
  164. data/ext/{matrix.c → gsl_native/matrix.c} +13 -16
  165. data/ext/{matrix_complex.c → gsl_native/matrix_complex.c} +119 -123
  166. data/ext/{matrix_double.c → gsl_native/matrix_double.c} +79 -82
  167. data/ext/{matrix_int.c → gsl_native/matrix_int.c} +53 -54
  168. data/ext/{matrix_source.c → gsl_native/matrix_source.h} +292 -318
  169. data/ext/{min.c → gsl_native/min.c} +45 -76
  170. data/ext/{monte.c → gsl_native/monte.c} +50 -64
  171. data/ext/{multifit.c → gsl_native/multifit.c} +142 -151
  172. data/ext/{multimin.c → gsl_native/multimin.c} +64 -92
  173. data/ext/{multimin_fsdf.c → gsl_native/multimin_fsdf.c} +16 -16
  174. data/ext/{multiroots.c → gsl_native/multiroots.c} +73 -76
  175. data/ext/{multiset.c → gsl_native/multiset.c} +4 -8
  176. data/ext/{ndlinear.c → gsl_native/ndlinear.c} +320 -321
  177. data/ext/{nmf.c → gsl_native/nmf.c} +11 -11
  178. data/ext/{nmf_wrap.c → gsl_native/nmf_wrap.c} +1 -1
  179. data/ext/{ntuple.c → gsl_native/ntuple.c} +23 -23
  180. data/ext/{odeiv.c → gsl_native/odeiv.c} +101 -116
  181. data/ext/gsl_native/ool.c +879 -0
  182. data/ext/{permutation.c → gsl_native/permutation.c} +39 -37
  183. data/ext/{poly.c → gsl_native/poly.c} +10 -13
  184. data/ext/{poly2.c → gsl_native/poly2.c} +16 -16
  185. data/ext/{poly_source.c → gsl_native/poly_source.h} +249 -293
  186. data/ext/{qrng.c → gsl_native/qrng.c} +9 -20
  187. data/ext/{randist.c → gsl_native/randist.c} +222 -247
  188. data/ext/{rational.c → gsl_native/rational.c} +12 -12
  189. data/ext/{rng.c → gsl_native/rng.c} +30 -47
  190. data/ext/{root.c → gsl_native/root.c} +47 -48
  191. data/ext/{sf.c → gsl_native/sf.c} +196 -244
  192. data/ext/{sf_airy.c → gsl_native/sf_airy.c} +2 -2
  193. data/ext/{sf_bessel.c → gsl_native/sf_bessel.c} +7 -7
  194. data/ext/{sf_clausen.c → gsl_native/sf_clausen.c} +1 -1
  195. data/ext/{sf_coulomb.c → gsl_native/sf_coulomb.c} +40 -40
  196. data/ext/{sf_coupling.c → gsl_native/sf_coupling.c} +30 -30
  197. data/ext/{sf_dawson.c → gsl_native/sf_dawson.c} +1 -1
  198. data/ext/{sf_debye.c → gsl_native/sf_debye.c} +1 -10
  199. data/ext/{sf_dilog.c → gsl_native/sf_dilog.c} +1 -1
  200. data/ext/{sf_elementary.c → gsl_native/sf_elementary.c} +3 -3
  201. data/ext/{sf_ellint.c → gsl_native/sf_ellint.c} +43 -43
  202. data/ext/{sf_elljac.c → gsl_native/sf_elljac.c} +3 -3
  203. data/ext/{sf_erfc.c → gsl_native/sf_erfc.c} +1 -5
  204. data/ext/{sf_exp.c → gsl_native/sf_exp.c} +3 -3
  205. data/ext/{sf_expint.c → gsl_native/sf_expint.c} +2 -12
  206. data/ext/{sf_fermi_dirac.c → gsl_native/sf_fermi_dirac.c} +1 -1
  207. data/ext/{sf_gamma.c → gsl_native/sf_gamma.c} +2 -6
  208. data/ext/{sf_gegenbauer.c → gsl_native/sf_gegenbauer.c} +1 -1
  209. data/ext/{sf_hyperg.c → gsl_native/sf_hyperg.c} +1 -1
  210. data/ext/{sf_laguerre.c → gsl_native/sf_laguerre.c} +4 -4
  211. data/ext/{sf_lambert.c → gsl_native/sf_lambert.c} +1 -1
  212. data/ext/{sf_legendre.c → gsl_native/sf_legendre.c} +1 -1
  213. data/ext/{sf_log.c → gsl_native/sf_log.c} +4 -4
  214. data/ext/gsl_native/sf_mathieu.c +235 -0
  215. data/ext/{sf_power.c → gsl_native/sf_power.c} +1 -1
  216. data/ext/{sf_psi.c → gsl_native/sf_psi.c} +3 -12
  217. data/ext/{sf_synchrotron.c → gsl_native/sf_synchrotron.c} +1 -1
  218. data/ext/{sf_transport.c → gsl_native/sf_transport.c} +1 -1
  219. data/ext/{sf_trigonometric.c → gsl_native/sf_trigonometric.c} +4 -4
  220. data/ext/{sf_zeta.c → gsl_native/sf_zeta.c} +1 -5
  221. data/ext/{signal.c → gsl_native/signal.c} +63 -68
  222. data/ext/{siman.c → gsl_native/siman.c} +45 -49
  223. data/ext/{sort.c → gsl_native/sort.c} +6 -7
  224. data/ext/{spline.c → gsl_native/spline.c} +28 -46
  225. data/ext/{stats.c → gsl_native/stats.c} +105 -118
  226. data/ext/{sum.c → gsl_native/sum.c} +34 -34
  227. data/ext/{tamu_anova.c → gsl_native/tamu_anova.c} +1 -1
  228. data/ext/{tensor.c → gsl_native/tensor.c} +8 -11
  229. data/ext/{tensor_source.c → gsl_native/tensor_source.h} +147 -148
  230. data/ext/{vector.c → gsl_native/vector.c} +11 -14
  231. data/ext/{vector_complex.c → gsl_native/vector_complex.c} +179 -184
  232. data/ext/{vector_double.c → gsl_native/vector_double.c} +178 -183
  233. data/ext/{vector_int.c → gsl_native/vector_int.c} +27 -29
  234. data/ext/{vector_source.c → gsl_native/vector_source.h} +428 -443
  235. data/ext/{wavelet.c → gsl_native/wavelet.c} +224 -246
  236. data/gsl.gemspec +29 -0
  237. data/lib/gsl.rb +8 -3
  238. data/lib/gsl/gnuplot.rb +3 -3
  239. data/lib/gsl/oper.rb +35 -60
  240. data/lib/gsl/version.rb +3 -0
  241. data/lib/rbgsl.rb +1 -3
  242. data/rdoc/alf.rdoc +5 -5
  243. data/rdoc/blas.rdoc +9 -9
  244. data/rdoc/bspline.rdoc +17 -17
  245. data/rdoc/changes.rdoc +4 -9
  246. data/rdoc/cheb.rdoc +25 -25
  247. data/rdoc/cholesky_complex.rdoc +21 -21
  248. data/rdoc/combi.rdoc +37 -37
  249. data/rdoc/complex.rdoc +22 -22
  250. data/rdoc/const.rdoc +47 -47
  251. data/rdoc/dht.rdoc +49 -49
  252. data/rdoc/diff.rdoc +42 -42
  253. data/rdoc/ehandling.rdoc +6 -6
  254. data/rdoc/eigen.rdoc +153 -153
  255. data/rdoc/fft.rdoc +146 -146
  256. data/rdoc/fit.rdoc +109 -109
  257. data/rdoc/function.rdoc +11 -11
  258. data/rdoc/graph.rdoc +17 -17
  259. data/rdoc/hist.rdoc +103 -103
  260. data/rdoc/hist2d.rdoc +42 -42
  261. data/rdoc/hist3d.rdoc +9 -9
  262. data/rdoc/integration.rdoc +110 -110
  263. data/rdoc/interp.rdoc +71 -71
  264. data/rdoc/intro.rdoc +8 -8
  265. data/rdoc/linalg.rdoc +188 -188
  266. data/rdoc/linalg_complex.rdoc +1 -1
  267. data/rdoc/math.rdoc +58 -58
  268. data/rdoc/matrix.rdoc +275 -275
  269. data/rdoc/min.rdoc +57 -57
  270. data/rdoc/monte.rdoc +22 -22
  271. data/rdoc/multimin.rdoc +95 -95
  272. data/rdoc/multiroot.rdoc +80 -80
  273. data/rdoc/narray.rdoc +32 -32
  274. data/rdoc/ndlinear.rdoc +54 -54
  275. data/rdoc/nonlinearfit.rdoc +100 -100
  276. data/rdoc/ntuple.rdoc +31 -31
  277. data/rdoc/odeiv.rdoc +88 -88
  278. data/rdoc/perm.rdoc +90 -90
  279. data/rdoc/poly.rdoc +66 -66
  280. data/rdoc/qrng.rdoc +21 -21
  281. data/rdoc/randist.rdoc +82 -82
  282. data/rdoc/ref.rdoc +57 -57
  283. data/rdoc/rng.rdoc +85 -85
  284. data/rdoc/roots.rdoc +57 -57
  285. data/rdoc/sf.rdoc +428 -428
  286. data/rdoc/siman.rdoc +19 -19
  287. data/rdoc/sort.rdoc +30 -30
  288. data/rdoc/start.rdoc +8 -8
  289. data/rdoc/stats.rdoc +52 -52
  290. data/rdoc/sum.rdoc +12 -12
  291. data/rdoc/tensor.rdoc +31 -31
  292. data/rdoc/tut.rdoc +1 -1
  293. data/rdoc/use.rdoc +39 -39
  294. data/rdoc/vector.rdoc +188 -188
  295. data/rdoc/vector_complex.rdoc +24 -24
  296. data/rdoc/wavelet.rdoc +46 -46
  297. data/test/gsl/blas_test.rb +79 -0
  298. data/test/gsl/bspline_test.rb +63 -0
  299. data/test/gsl/cdf_test.rb +1512 -0
  300. data/test/gsl/cheb_test.rb +80 -0
  301. data/test/gsl/combination_test.rb +100 -0
  302. data/test/gsl/complex_test.rb +20 -0
  303. data/test/gsl/const_test.rb +29 -0
  304. data/test/gsl/deriv_test.rb +62 -0
  305. data/test/gsl/dht_test.rb +79 -0
  306. data/test/gsl/diff_test.rb +53 -0
  307. data/test/gsl/eigen_test.rb +563 -0
  308. data/test/gsl/err_test.rb +23 -0
  309. data/test/gsl/fit_test.rb +101 -0
  310. data/test/gsl/histo_test.rb +14 -0
  311. data/test/gsl/index_test.rb +61 -0
  312. data/test/gsl/integration_test.rb +274 -0
  313. data/test/gsl/interp_test.rb +27 -0
  314. data/test/gsl/linalg_test.rb +463 -0
  315. data/test/gsl/matrix_nmf_test.rb +37 -0
  316. data/test/gsl/matrix_test.rb +98 -0
  317. data/test/gsl/min_test.rb +89 -0
  318. data/test/gsl/monte_test.rb +77 -0
  319. data/test/gsl/multifit_test.rb +753 -0
  320. data/test/gsl/multimin_test.rb +157 -0
  321. data/test/gsl/multiroot_test.rb +135 -0
  322. data/test/gsl/multiset_test.rb +52 -0
  323. data/test/gsl/odeiv_test.rb +275 -0
  324. data/test/gsl/oper_test.rb +98 -0
  325. data/test/gsl/poly_test.rb +338 -0
  326. data/test/gsl/qrng_test.rb +94 -0
  327. data/test/gsl/quartic_test.rb +28 -0
  328. data/test/gsl/randist_test.rb +122 -0
  329. data/test/gsl/rng_test.rb +303 -0
  330. data/test/gsl/roots_test.rb +78 -0
  331. data/test/gsl/sf_test.rb +2079 -0
  332. data/test/gsl/stats_test.rb +122 -0
  333. data/test/gsl/sum_test.rb +69 -0
  334. data/test/gsl/tensor_test.rb +396 -0
  335. data/test/gsl/vector_test.rb +223 -0
  336. data/test/gsl/wavelet_test.rb +130 -0
  337. data/test/gsl_test.rb +321 -0
  338. data/test/test_helper.rb +42 -0
  339. data/uncrustify.cfg +1693 -0
  340. metadata +337 -378
  341. data/README +0 -32
  342. data/VERSION +0 -1
  343. data/ext/bspline.c +0 -130
  344. data/ext/const.c +0 -673
  345. data/ext/cqp.c +0 -283
  346. data/ext/extconf.rb +0 -295
  347. data/ext/fcmp.c +0 -66
  348. data/ext/fresnel.c +0 -312
  349. data/ext/jacobi.c +0 -739
  350. data/ext/ool.c +0 -879
  351. data/ext/oper_complex_source.c +0 -253
  352. data/ext/sf_mathieu.c +0 -238
  353. data/include/rb_gsl_config.h +0 -62
  354. data/include/rb_gsl_dirac.h +0 -13
  355. data/rdoc/index.rdoc +0 -62
  356. data/rdoc/rngextra.rdoc +0 -11
  357. data/rdoc/screenshot.rdoc +0 -40
  358. data/setup.rb +0 -1585
  359. data/tests/blas/amax.rb +0 -14
  360. data/tests/blas/asum.rb +0 -16
  361. data/tests/blas/axpy.rb +0 -25
  362. data/tests/blas/copy.rb +0 -23
  363. data/tests/blas/dot.rb +0 -23
  364. data/tests/bspline.rb +0 -53
  365. data/tests/cdf.rb +0 -1388
  366. data/tests/cheb.rb +0 -112
  367. data/tests/combination.rb +0 -123
  368. data/tests/complex.rb +0 -17
  369. data/tests/const.rb +0 -24
  370. data/tests/deriv.rb +0 -85
  371. data/tests/dht/dht1.rb +0 -17
  372. data/tests/dht/dht2.rb +0 -23
  373. data/tests/dht/dht3.rb +0 -23
  374. data/tests/dht/dht4.rb +0 -23
  375. data/tests/diff.rb +0 -78
  376. data/tests/eigen/eigen.rb +0 -220
  377. data/tests/eigen/gen.rb +0 -105
  378. data/tests/eigen/genherm.rb +0 -66
  379. data/tests/eigen/gensymm.rb +0 -68
  380. data/tests/eigen/nonsymm.rb +0 -53
  381. data/tests/eigen/nonsymmv.rb +0 -53
  382. data/tests/eigen/symm-herm.rb +0 -74
  383. data/tests/err.rb +0 -58
  384. data/tests/fit.rb +0 -124
  385. data/tests/gsl_test.rb +0 -118
  386. data/tests/gsl_test2.rb +0 -110
  387. data/tests/histo.rb +0 -12
  388. data/tests/integration/integration1.rb +0 -72
  389. data/tests/integration/integration2.rb +0 -71
  390. data/tests/integration/integration3.rb +0 -71
  391. data/tests/integration/integration4.rb +0 -71
  392. data/tests/interp.rb +0 -45
  393. data/tests/linalg/HH.rb +0 -64
  394. data/tests/linalg/LU.rb +0 -47
  395. data/tests/linalg/QR.rb +0 -77
  396. data/tests/linalg/SV.rb +0 -24
  397. data/tests/linalg/TDN.rb +0 -116
  398. data/tests/linalg/TDS.rb +0 -122
  399. data/tests/linalg/bidiag.rb +0 -73
  400. data/tests/linalg/cholesky.rb +0 -20
  401. data/tests/linalg/linalg.rb +0 -158
  402. data/tests/matrix/matrix_complex_test.rb +0 -36
  403. data/tests/matrix/matrix_nmf_test.rb +0 -39
  404. data/tests/matrix/matrix_test.rb +0 -48
  405. data/tests/min.rb +0 -99
  406. data/tests/monte/miser.rb +0 -31
  407. data/tests/monte/vegas.rb +0 -45
  408. data/tests/multifit/test_2dgauss.rb +0 -112
  409. data/tests/multifit/test_brown.rb +0 -90
  410. data/tests/multifit/test_enso.rb +0 -246
  411. data/tests/multifit/test_filip.rb +0 -155
  412. data/tests/multifit/test_gauss.rb +0 -97
  413. data/tests/multifit/test_longley.rb +0 -110
  414. data/tests/multifit/test_multifit.rb +0 -52
  415. data/tests/multimin.rb +0 -139
  416. data/tests/multiroot.rb +0 -131
  417. data/tests/multiset.rb +0 -52
  418. data/tests/narray/blas_dnrm2.rb +0 -20
  419. data/tests/odeiv.rb +0 -353
  420. data/tests/poly/poly.rb +0 -290
  421. data/tests/poly/special.rb +0 -65
  422. data/tests/qrng.rb +0 -131
  423. data/tests/quartic.rb +0 -29
  424. data/tests/randist.rb +0 -134
  425. data/tests/rng.rb +0 -305
  426. data/tests/roots.rb +0 -76
  427. data/tests/run-test.sh +0 -17
  428. data/tests/sf/gsl_test_sf.rb +0 -249
  429. data/tests/sf/test_airy.rb +0 -83
  430. data/tests/sf/test_bessel.rb +0 -306
  431. data/tests/sf/test_coulomb.rb +0 -17
  432. data/tests/sf/test_dilog.rb +0 -25
  433. data/tests/sf/test_gamma.rb +0 -209
  434. data/tests/sf/test_hyperg.rb +0 -356
  435. data/tests/sf/test_legendre.rb +0 -227
  436. data/tests/sf/test_mathieu.rb +0 -59
  437. data/tests/sf/test_mode.rb +0 -19
  438. data/tests/sf/test_sf.rb +0 -839
  439. data/tests/stats.rb +0 -174
  440. data/tests/stats_mt.rb +0 -16
  441. data/tests/sum.rb +0 -98
  442. data/tests/sys.rb +0 -323
  443. data/tests/tensor.rb +0 -419
  444. data/tests/vector/vector_complex_test.rb +0 -101
  445. data/tests/vector/vector_test.rb +0 -141
  446. data/tests/wavelet.rb +0 -142
@@ -9,7 +9,7 @@ func = GSL::MultiRoot::Function.alloc(2) { |x, f|
9
9
  f[0] = -2.0*x0*x0 + 3.0*x0*x1 + 4.0*sin(x1) - 6.0
10
10
  f[1] = 3.0*x0*x0 - 2.0*x0*x1*x1 + 3.0*cos(x0) + 4.0
11
11
  }
12
-
12
+
13
13
  p func.solve([1.0, 2.0].to_gv, 1000, 1e-7, "hybrids")
14
14
  p func.solve([1.0, 2.0].to_gv, 1000, "broyden")
15
15
  p func.solve([1.0, 2.0], "hybrid")
@@ -2,13 +2,13 @@
2
2
  # 19/Apr/2004 by Yoshiki Tsunesada
3
3
  #
4
4
  # This is an example to calculate the orbital evolution of
5
- # a double neutron star (binary) system. General relativity predicts
6
- # that the binary orbital decays by radiating gravitational waves,
7
- # and the two stars will coalesce in time scale of 100-1000 Mega-years.
8
- # The values used here are of the binary system J0730-3039 discovered
9
- # in 2003 (Burgay et al., Nature 2003). The result shows that the two
10
- # neutron stars will merge after about 85 Mega-years. From the age of
11
- # the system 100 Mega-year, the lifetime of the system is estimated
5
+ # a double neutron star (binary) system. General relativity predicts
6
+ # that the binary orbital decays by radiating gravitational waves,
7
+ # and the two stars will coalesce in time scale of 100-1000 Mega-years.
8
+ # The values used here are of the binary system J0730-3039 discovered
9
+ # in 2003 (Burgay et al., Nature 2003). The result shows that the two
10
+ # neutron stars will merge after about 85 Mega-years. From the age of
11
+ # the system 100 Mega-year, the lifetime of the system is estimated
12
12
  # about 185 Mega-years.
13
13
  #
14
14
  # References:
@@ -16,7 +16,7 @@
16
16
  # 2. Shapiro & Teukolsky, "Black holes, white dwarfs and neutron stars"
17
17
  # John Wiley and Sans (1983)
18
18
  #
19
-
19
+
20
20
  require("gsl")
21
21
  include Math
22
22
 
@@ -33,7 +33,7 @@ GMsolarC3 = 4.925490947e-6
33
33
  MegaYear = 3600*24*365*1e6
34
34
 
35
35
  # Time evolution of the binary orbital period and the eccentricity
36
- # due to gravitational radiation.
36
+ # due to gravitational radiation.
37
37
  # The calculation is based on general relativity (See e.g. Ref.2).
38
38
  # y[0]: orbital period (pb)
39
39
  # y[1]: eccentricity (e)
@@ -45,7 +45,7 @@ deriv = Proc.new { |t, y, dydt, binary|
45
45
  e = y[1] # eccentricity
46
46
  m1 = binary.m1 # neutron star masses
47
47
  m2 = binary.m2
48
- totalM = m1 + m2 # total mass
48
+ totalM = m1 + m2 # total mass
49
49
  mu = m1*m2/totalM # reduced mass
50
50
  mm = mu*GSL::pow(totalM, 2.0/3.0)
51
51
  f_e = GSL::pow(1.0 - e*e, -3.5)*(1.0 + (73.0/24.0 + 37.0/96.0*e*e)*e*e);
@@ -85,8 +85,8 @@ tend = 2500*MegaYear
85
85
  h = 1.0*MegaYear
86
86
 
87
87
  begin
88
- file = File.open("binarysystem.dat", "w")
89
- while t < tend
88
+ file = File.open("binarysystem.dat", "w")
89
+ while t < tend
90
90
  t, h, status = solver.apply(t, tend, h, y)
91
91
  break if status != GSL::SUCCESS
92
92
  break if GSL::isnan?(y[0])
@@ -1,5 +1,5 @@
1
1
  #!/usr/bin/env ruby
2
- # Solve
2
+ # Solve
3
3
  # dydt = -2y ---> y(t) = exp(-2t)
4
4
 
5
5
  require("gsl")
@@ -40,11 +40,11 @@ def odeiv_compare_algorithm(solver, steptype, t0, tend, h0, y0, outfile)
40
40
  end
41
41
 
42
42
  if GSL::VERSION >= "1.5.90"
43
- ALGORITHMS = ["rk2", "rk4", "rkf45", "rkck", "rk8pd", "rk2imp", "rk4imp",
43
+ ALGORITHMS = ["rk2", "rk4", "rkf45", "rkck", "rk8pd", "rk2imp", "rk4imp",
44
44
  "bsimp", "gear1", "gear2", "rk2simp"]
45
45
  gpfile = "demo2.gp"
46
46
  else
47
- ALGORITHMS = ["rk2", "rk4", "rkf45", "rkck", "rk8pd", "rk2imp", "rk4imp",
47
+ ALGORITHMS = ["rk2", "rk4", "rkf45", "rkck", "rk8pd", "rk2imp", "rk4imp",
48
48
  "bsimp", "gear1", "gear2"]
49
49
  gpfile = "demo.gp"
50
50
  end
@@ -1,13 +1,13 @@
1
1
  #!/usr/bin/env ruby
2
2
  # Solve Schroedinger equation
3
3
  #
4
- # This example is taken from frei1.cpp
4
+ # This example is taken from frei1.cpp
5
5
  # in "Numerische Physik" p201-204 (Springer),
6
6
  # which simulates the time evolution of a probability density.
7
7
  #
8
8
  # Name: frei1.cpp
9
- # Zweck: Simuliert ein quantenmechanisches freies Teilchen
10
- # Gleichung: Schroedingergleichung ohne Potential verwendete
9
+ # Zweck: Simuliert ein quantenmechanisches freies Teilchen
10
+ # Gleichung: Schroedingergleichung ohne Potential verwendete
11
11
  # Bibiliothek: GSL
12
12
  #
13
13
  # Reference:
@@ -20,7 +20,7 @@ require("gsl")
20
20
  #NMAX = 8192
21
21
  NMAX = 256
22
22
 
23
- # The wave equation:
23
+ # The wave equation:
24
24
  # calculate time derivative of the wave function.
25
25
  # The second spatial derivative is approximated by
26
26
  # d2_psi/dx2 ~ (psi[n+1] - 2*psi[n] + pxi[n-1])/(dx*dx)
@@ -68,7 +68,7 @@ sum = 0.0
68
68
  for n in 0...NMAX do
69
69
  x = (n-NMAX/2) * $dx
70
70
  psi[n] = Math::exp(-GSL::pow_2(x/alpha)/2)
71
- sum += GSL::pow_2(psi[n])
71
+ sum += GSL::pow_2(psi[n])
72
72
  end
73
73
  sum = 1.0/Math::sqrt(sum)
74
74
 
@@ -89,12 +89,12 @@ IO.popen("graph -T X -C -g 3", "w") do |io|
89
89
  c = GSL::Odeiv::Control.y_new(atol, rtol)
90
90
  evolve = GSL::Odeiv::Evolve.alloc(2*NMAX)
91
91
  sys = GSL::Odeiv::System.alloc(f, 2*NMAX)
92
-
92
+
93
93
  t = 0.0
94
94
  for n in 1..n_out do
95
95
  t1 = n*dt
96
96
  STDOUT.printf("t = %2.1f (%2d/%2d)\n", t1-dt, n, n_out)
97
- while t < t1
97
+ while t < t1
98
98
  t, h, status = evolve.apply(c, step, sys, t, t1, h, psi)
99
99
  break if status != GSL::SUCCESS
100
100
  end
@@ -1,5 +1,5 @@
1
1
  #!/usr/bin/env ruby
2
- # This example is taken from frei2.cpp
2
+ # This example is taken from frei2.cpp
3
3
  # in "Numerische Physik" p205-206 (Springer).
4
4
  #
5
5
  # Reference:
@@ -45,10 +45,10 @@ IO.popen("graph -T X -C -g 3", "w") do |io|
45
45
  io.printf("%e %e\n", x, psi[n1].abs)
46
46
  end
47
47
  io.printf("\n")
48
-
48
+
49
49
  psi_p.radix2_forward!
50
-
51
-
50
+
51
+
52
52
  t = 0.0
53
53
  for n in 1..n_out do
54
54
  t1 = n*dt
@@ -5,7 +5,7 @@
5
5
  # b: resist
6
6
  # f: external force
7
7
 
8
- require("gsl")
8
+ require("gsl")
9
9
 
10
10
  dim = 2
11
11
 
@@ -18,16 +18,16 @@ include Math
18
18
  Sedov = Proc.new { |logx, y, dydlogx, sh|
19
19
  a = -5.0*(6.0*y[2] - 15.0*y[2]*y[0]*sh + 2.0*y[0]*y[1] - 7.0*y[0]*y[0]*y[1] + 5.0*GSL::pow_3(y[0])*y[1])
20
20
  b = -25.0*y[2]*sh + 4.0*y[1] - 20.0*y[0]*y[1] + 25.0*y[0]*y[0]*y[1]
21
- dydlogx[0] = a/b
21
+ dydlogx[0] = a/b
22
22
 
23
23
  bb = -5.0*(6.0*y[2] - 15.0*y[2]*y[0]*sh + 2.0*y[0]*y[1] - 7.0*y[0]*y[0]*y[1])
24
24
 
25
25
  a = -5.0*(-30.0*y[2]*y[1] + 2.0*y[0]*y[1]*y[1] - 25.0*GSL::pow_2(y[0]*y[1]) + 50.0*GSL::pow_3(y[0])*GSL::pow_2(y[1]))
26
26
  bb = (-2.0 + 5.0*y[0])*(-25.0*y[2]*sh + 4.0*y[1] - 20.0*y[0]*y[1] + 25.0*y[0]*y[0]*y[1])
27
- dydlogx[1] = a/bb
27
+ dydlogx[1] = a/bb
28
28
 
29
29
  a = -5.0*y[2]*(-10.0*y[2]*sh + 4.0*y[1] - 14.0*y[0]*y[1] + 10.0*y[0]*y[0]*y[1] - y[0]*sh*y[1] + 10.0*y[0]*y[0]*sh*y[1])
30
- dydlogx[2] = a/b
30
+ dydlogx[2] = a/b
31
31
  }
32
32
 
33
33
  DIM = 3
@@ -15,8 +15,8 @@ module Degenerate
15
15
  hbar = PLANCKS_CONSTANT_HBAR # Planck's constant
16
16
  me = MASS_ELECTRON # Electron mass
17
17
  mn = MASS_NEUTRON # Neutron mass
18
- mu = UNIFIED_ATOMIC_MASS
19
- c = SPEED_OF_LIGHT
18
+ mu = UNIFIED_ATOMIC_MASS
19
+ c = SPEED_OF_LIGHT
20
20
  ELambda = hbar/me/c # Compton length of electron
21
21
  NLambda = hbar/mn/c
22
22
  MeC2 = me*c*c # Electron rest mass energy
@@ -29,24 +29,24 @@ module Degenerate
29
29
  tmp = sqrt(1.0 + x*x)
30
30
  (x*tmp*(2.0*x*x/3.0 - 1.0) + log(x + tmp))/8/PI/PI
31
31
  end
32
-
32
+
33
33
  def chi(x)
34
34
  tmp = sqrt(1.0 + x*x)
35
35
  (x*tmp*(1.0 + 2*x*x) - log(x + tmp))/8/PI/PI
36
36
  end
37
-
37
+
38
38
  def xe_rho_mue(rho, mue)
39
39
  Factor_xe*GSL::pow(rho/mue, 1.0/3.0)
40
40
  end
41
-
41
+
42
42
  def xn_rho(rho)
43
43
  Factor_xn*GSL::pow(rho, 1.0/3.0)
44
44
  end
45
-
45
+
46
46
  end
47
47
 
48
48
  # Polytrope gas sphere
49
- module Polytrope
49
+ module Polytrope
50
50
 
51
51
  # Lane-Emden equation
52
52
  # n: polytrope index
@@ -59,7 +59,7 @@ module Polytrope
59
59
  dim = 2
60
60
  y = GSL::Vector[1.0, 0.0]
61
61
  dydx = GSL::Vector.alloc(dim)
62
-
62
+
63
63
  solver = Solver.alloc(Step::RKF45, [1e-6, 0], EmdenEq, dim)
64
64
  solver.set_params(n)
65
65
  solver.reset
@@ -71,7 +71,7 @@ module Polytrope
71
71
  x = 0.0001
72
72
  xend = 10.0
73
73
  h = 1e-6
74
-
74
+
75
75
  file = File.open("polytrope.dat", "w")
76
76
  i = 0
77
77
  while x < xend
@@ -94,10 +94,10 @@ module Polytrope
94
94
  vdy2 = vdy.subvector(0, i)
95
95
  spline = GSL::Spline.alloc(GSL::Interp::AKIMA, i)
96
96
  spline.init(vy2.reverse, vx2.reverse)
97
-
97
+
98
98
  # Find star surface:
99
99
  # Star sufrace is defined as the zero point of density structure function
100
- x1 = spline.eval(0.0)
100
+ x1 = spline.eval(0.0)
101
101
  spline.init(vx2, vdy2)
102
102
  yx2 = spline.eval(x1).abs
103
103
  return [x1, yx2*x1*x1]
@@ -11,12 +11,12 @@ pp.init
11
11
  pp.fprintf(STDOUT, " %u")
12
12
  printf("\n")
13
13
 
14
- puts(" random permutation:");
14
+ puts(" random permutation:");
15
15
  r.shuffle(pp)
16
16
  pp.fprintf(STDOUT, " %u")
17
17
  printf("\n");
18
18
 
19
- puts("inverse permutation:");
19
+ puts("inverse permutation:");
20
20
  q = pp.inverse
21
21
  q.fprintf(STDOUT, " %u")
22
22
  printf ("\n");
@@ -1,5 +1,5 @@
1
1
  #!/usr/bin/env ruby
2
- require("gsl")
2
+ require("gsl")
3
3
 
4
4
  p = GSL::Permutation::alloc(10)
5
5
  p.init
@@ -6,7 +6,7 @@ include GSL
6
6
  poly = Poly[1.5, -1.25, -3.75, 0, 1]
7
7
  # Solve the equation p(x) == 0
8
8
  root = poly.solve # Vector::Complex
9
- # Extract only the real parts
9
+ # Extract only the real parts
10
10
  # (imaginary parts are zero for this case)
11
11
  re = root.real # Vector::View
12
12
 
@@ -19,7 +19,7 @@ h[2] = GSL::Histogram.alloc(61, -30, 30)
19
19
 
20
20
  i = 0
21
21
  for n in [6, 50, 100] do
22
- M.times do
22
+ M.times do
23
23
  s = 0
24
24
  n.times do
25
25
  ds = rng.get%2 == 0 ? 1 : -1
@@ -7,8 +7,8 @@ names = ["default", "mt19937", "mt19937_1999", "mt19937_1998", "ranlxs0", "ranlx
7
7
  "gfsr4", "rand", "random_bsd", "random8_bsd", "random32_bsd", "random64_bsd",
8
8
  "random128_bsd", "random256_bsd", "random_libc5", "random_glibc2", "rand48", "ran0",
9
9
  "ran1", "ran2", "ran3", "ranf", "ranmar", "r250", "tt800", "vax", "transputer",
10
- "randu", "minstd", "uni", "uni32", "slatec", "zuf", "borosh13", "coveyou",
11
- "fishman18", "fishman20", "fishman2x", "knuthran2", "knuthran", "lecuyer21",
10
+ "randu", "minstd", "uni", "uni32", "slatec", "zuf", "borosh13", "coveyou",
11
+ "fishman18", "fishman20", "fishman2x", "knuthran2", "knuthran", "lecuyer21",
12
12
  "waterman14"]
13
13
 
14
14
  names.each do |name|
@@ -17,7 +17,7 @@ sigma = Math::sqrt(N).to_i
17
17
 
18
18
  h = GSL::Histogram.alloc(8*sigma+1, [-4*sigma-0.5, 4*sigma+0.5])
19
19
 
20
- M.times do
20
+ M.times do
21
21
  s = 0
22
22
  N.times do
23
23
  ds = rng.get%2 == 0 ? 1 : -1
@@ -27,8 +27,8 @@ M.times do
27
27
  end
28
28
 
29
29
  x = GSL::Vector.linspace(-40, 40, 80)
30
- y = GSL::Ran::gaussian_pdf(x, sigma)*M*2
31
- # Factor 2 is not important, but necessary
30
+ y = GSL::Ran::gaussian_pdf(x, sigma)*M*2
31
+ # Factor 2 is not important, but necessary
32
32
  # because only the even ranges are filled:
33
33
  # a + b = N a: positive steps, b: negative steps
34
34
  # a - b = s s: the end point after the N steps
@@ -18,7 +18,7 @@ p r.max
18
18
  p r.get
19
19
  p r.get
20
20
 
21
- r2 = GSL::Random::Rng.alloc
21
+ r2 = GSL::Random::Rng.alloc
22
22
  p r2.uniform
23
23
  p r2.uniform
24
24
 
@@ -8,7 +8,7 @@ f = GSL::Function.alloc { |x, params|
8
8
  }
9
9
  f.set_params(1, 0, -5)
10
10
  expected = Math::sqrt(5.0)
11
- printf("%5s [%9s, %9s] %9s %10s %9s\n",
11
+ printf("%5s [%9s, %9s] %9s %10s %9s\n",
12
12
  "iter", "lower", "upper", "root", "err", "err(est)")
13
13
  solver.set(f, 0.0, 5.0)
14
14
  iter = 0
@@ -19,7 +19,7 @@ f.set_params(1, 0, -5)
19
19
  expected = sqrt(5.0)
20
20
 
21
21
  printf("%5s [%9s, %9s] %9s %10s %9s\n",
22
- "iter", "lower", "upper", "root",
22
+ "iter", "lower", "upper", "root",
23
23
  "err", "err(est)")
24
24
 
25
25
  solver.set(f, 0.0, 5.0)
@@ -7,7 +7,7 @@ f = GSL::Function.alloc { |x, params|
7
7
  }
8
8
  f.set_params(1, 0, -5)
9
9
  expected = Math::sqrt(5.0)
10
- printf("%5s [%9s, %9s] %9s %10s %9s\n", "iter",
10
+ printf("%5s [%9s, %9s] %9s %10s %9s\n", "iter",
11
11
  "lower", "upper", "root", "err", "err(est)")
12
12
  solver.set(f, 0.0, 5.0)
13
13
  iter = 0; status = nil
@@ -2,14 +2,14 @@
2
2
  require("gsl")
3
3
  include Math
4
4
 
5
- f = Proc.new { |x, params|
5
+ f = Proc.new { |x, params|
6
6
  a = params[0]
7
7
  b = params[1]
8
8
  c = params[2]
9
9
  (a*x + b)*x + c
10
10
  }
11
11
 
12
- df = Proc.new { |x, params|
12
+ df = Proc.new { |x, params|
13
13
  a = params[0]
14
14
  b = params[1]
15
15
  2.0*a*x + b
@@ -9,4 +9,3 @@ set ylabel 'Fractional ionization'
9
9
  set x2label 'Temperature of the Universe [K]'
10
10
  set grid
11
11
  plot 'recombination.dat' u 2:3 title '' w l lw 3, 1 title ''
12
-
@@ -1,5 +1,5 @@
1
1
  #!/usr/bin/env ruby
2
- require("gsl")
2
+ require("gsl")
3
3
 
4
4
  p GSL::Sf.hyperg_0F1(2, 3)
5
5
  r = GSL::Sf.hyperg_0F1_e(2, 3)
@@ -19,7 +19,7 @@ File.open("sphbessel.dat", "w") do |file|
19
19
  k0 = bessel_k0_scaled(x)
20
20
  k1 = bessel_k1_scaled(x)
21
21
  k2 = bessel_k2_scaled(x)
22
- file.printf("%e %e %e %e %e %e %e %e %e %e %e %e %e %e\n",
22
+ file.printf("%e %e %e %e %e %e %e %e %e %e %e %e %e %e\n",
23
23
  x, j0, j1, j2, j3, y0, y1, y2, i0, i1, i2, k0, k1, k2)
24
24
  x += 0.1
25
25
  end
@@ -1,5 +1,5 @@
1
1
  #!/usr/bin/env ruby
2
- require("gsl")
2
+ require("gsl")
3
3
 
4
4
  N = 100000
5
5
  k = 5
@@ -1,16 +1,16 @@
1
1
  #!/usr/bin/env ruby
2
2
  require("gsl")
3
3
 
4
- data = GSL::Vector[88.60,73.20,91.40,68.00,75.20,63.00,53.90,69.20,
5
- 50.10,71.50,44.90,59.50,40.20,56.30,38.70,31.00,
4
+ data = GSL::Vector[88.60,73.20,91.40,68.00,75.20,63.00,53.90,69.20,
5
+ 50.10,71.50,44.90,59.50,40.20,56.30,38.70,31.00,
6
6
  39.60,45.30,25.20,22.70]
7
7
  factor = GSL::Vector::Int[1,1,1,1,1,2,2,2,2,2,3,3,3,3,3,4,4,4,4,4]
8
8
 
9
9
  table = GSL::TAMU_ANOVA::Table.oneway(data, factor, 4)
10
10
  table.print
11
11
 
12
- data = GSL::Vector[45.50,45.30,45.40,44.40,44.60,43.90,44.60,44.00,44.20,
13
- 43.90,44.70,44.20,44.00,43.80,44.60,43.10,46.00,45.90,
12
+ data = GSL::Vector[45.50,45.30,45.40,44.40,44.60,43.90,44.60,44.00,44.20,
13
+ 43.90,44.70,44.20,44.00,43.80,44.60,43.10,46.00,45.90,
14
14
  44.80,46.20,45.10,45.50]
15
15
  factor = GSL::Vector::Int[1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,3,3,3,3,3,3]
16
16
 
@@ -23,7 +23,7 @@ a += b
23
23
  a
24
24
 
25
25
  # Subtract b from a
26
- a -= b
26
+ a -= b
27
27
 
28
28
  # Show a
29
29
  a