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
@@ -1,23 +1,23 @@
1
1
  /* histogram/find.c
2
- *
2
+ *
3
3
  * Copyright (C) 1996, 1997, 1998, 1999, 2000 Brian Gough
4
- *
4
+ *
5
5
  * This program is free software; you can redistribute it and/or modify
6
6
  * it under the terms of the GNU General Public License as published by
7
7
  * the Free Software Foundation; either version 2 of the License, or (at
8
8
  * your option) any later version.
9
- *
9
+ *
10
10
  * This program is distributed in the hope that it will be useful, but
11
11
  * WITHOUT ANY WARRANTY; without even the implied warranty of
12
12
  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13
13
  * General Public License for more details.
14
- *
14
+ *
15
15
  * You should have received a copy of the GNU General Public License
16
16
  * along with this program; if not, write to the Free Software
17
17
  * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
18
18
  */
19
19
 
20
- #include "rb_gsl_histogram.h"
20
+ #include "include/rb_gsl_histogram.h"
21
21
  /* determines whether to optimize for linear ranges */
22
22
  #define LINEAR_OPT 1
23
23
 
@@ -26,7 +26,6 @@ int mygsl_find (const size_t n, const double range[], const double x, size_t * i
26
26
  size_t i_linear, lower, upper, mid;
27
27
  if (x < range[0]) return -1;
28
28
  if (x >= range[n]) return +1;
29
-
30
29
  /* optimize for linear case */
31
30
 
32
31
  #ifdef LINEAR_OPT
@@ -36,47 +35,46 @@ int mygsl_find (const size_t n, const double range[], const double x, size_t * i
36
35
  }
37
36
 
38
37
  if (x >= range[i_linear] && x < range[i_linear + 1])
39
- {
40
- *i = i_linear;
41
- return 0;
42
- }
38
+ {
39
+ *i = i_linear;
40
+ return 0;
41
+ }
43
42
  #endif
44
43
 
45
44
  /* perform binary search */
46
45
 
47
- upper = n ;
48
- lower = 0 ;
46
+ upper = n;
47
+ lower = 0;
49
48
 
50
49
  while (upper - lower > 1)
50
+ {
51
+ mid = (upper + lower) / 2;
52
+
53
+ if (x >= range[mid])
51
54
  {
52
- mid = (upper + lower) / 2 ;
53
-
54
- if (x >= range[mid])
55
- {
56
- lower = mid ;
57
- }
58
- else
59
- {
60
- upper = mid ;
61
- }
55
+ lower = mid;
62
56
  }
57
+ else
58
+ {
59
+ upper = mid;
60
+ }
61
+ }
63
62
 
64
- *i = lower ;
63
+ *i = lower;
65
64
 
66
65
  /* sanity check the result */
67
66
 
68
67
  if (x < range[lower] || x >= range[lower + 1])
69
- {
70
- GSL_ERROR ("x not found in range", GSL_ESANITY);
71
- }
72
-
68
+ {
69
+ GSL_ERROR ("x not found in range", GSL_ESANITY);
70
+ }
73
71
  return 0;
74
72
  }
75
73
 
76
74
  int mygsl_find2d (const size_t nx, const double xrange[],
77
- const size_t ny, const double yrange[],
78
- const double x, const double y,
79
- size_t * i, size_t * j)
75
+ const size_t ny, const double yrange[],
76
+ const double x, const double y,
77
+ size_t * i, size_t * j)
80
78
  {
81
79
  int status = mygsl_find (nx, xrange, x, i);
82
80
  if (status) return status;
@@ -86,10 +84,10 @@ int mygsl_find2d (const size_t nx, const double xrange[],
86
84
  }
87
85
 
88
86
  int mygsl_find3d (const size_t nx, const double xrange[],
89
- const size_t ny, const double yrange[],
90
- const size_t nz, const double zrange[],
91
- const double x, const double y, const double z,
92
- size_t * i, size_t * j, size_t *k)
87
+ const size_t ny, const double yrange[],
88
+ const size_t nz, const double zrange[],
89
+ const double x, const double y, const double z,
90
+ size_t * i, size_t * j, size_t *k)
93
91
  {
94
92
  int status = mygsl_find (nx, xrange, x, i);
95
93
  if (status) return status;
@@ -18,11 +18,11 @@
18
18
  */
19
19
  /***************************************************************
20
20
  *
21
- * File gsl_histogram_oper.c:
22
- * Routine to make operation on histograms.
21
+ * File gsl_histogram_oper.c:
22
+ * Routine to make operation on histograms.
23
23
  * Need GSL library and header.
24
24
  * Contains the routines:
25
- * gsl_histogram_same_binning check if two histograms have the same binning
25
+ * gsl_histogram_same_binning check if two histograms have the same binning
26
26
  * gsl_histogram_add add two histograms
27
27
  * gsl_histogram_sub subctract two histograms
28
28
  * gsl_histogram_mult multiply two histograms
@@ -38,7 +38,7 @@
38
38
  #include <gsl/gsl_math.h>
39
39
  #include <gsl/gsl_histogram.h>
40
40
 
41
- /*
41
+ /*
42
42
  * gsl_histogram_same_binning:
43
43
  * control if two histograms have the
44
44
  * same binning
@@ -48,112 +48,103 @@ int
48
48
  mygsl_histogram_equal_bins_p (const gsl_histogram * h1, const gsl_histogram * h2)
49
49
  {
50
50
  if (h1->n != h2->n)
51
- {
52
- return 0;
53
- }
54
-
51
+ {
52
+ return 0;
53
+ }
55
54
  {
56
55
  size_t i;
57
56
  /* init ranges */
58
57
 
59
58
  for (i = 0; i <= h1->n; i++)
59
+ {
60
+ if (gsl_fcmp(h1->range[i],h2->range[i], 1e-12))
60
61
  {
61
- if (gsl_fcmp(h1->range[i],h2->range[i], 1e-12))
62
- {
63
- return 0;
64
- }
62
+ return 0;
65
63
  }
64
+ }
66
65
  }
67
66
 
68
67
  return 1;
69
68
  }
70
69
 
71
- /*
70
+ /*
72
71
  * gsl_histogram_add:
73
72
  * add two histograms
74
73
  */
75
- int
74
+ int
76
75
  mygsl_histogram_add (gsl_histogram * h1, const gsl_histogram * h2)
77
76
  {
78
77
  size_t i;
79
78
 
80
79
  if (!mygsl_histogram_equal_bins_p (h1, h2))
81
- {
82
- GSL_ERROR ("histograms have different binning", GSL_EINVAL);
83
- }
84
-
80
+ {
81
+ GSL_ERROR ("histograms have different binning", GSL_EINVAL);
82
+ }
85
83
  for (i = 0; i < h1->n; i++)
86
- {
87
- h1->bin[i] += h2->bin[i];
88
- }
89
-
84
+ {
85
+ h1->bin[i] += h2->bin[i];
86
+ }
90
87
  return GSL_SUCCESS;
91
88
  }
92
89
 
93
- /*
90
+ /*
94
91
  * gsl_histogram_sub:
95
92
  * subtract two histograms
96
93
  */
97
94
 
98
- int
95
+ int
99
96
  mygsl_histogram_sub (gsl_histogram * h1, const gsl_histogram * h2)
100
97
  {
101
98
  size_t i;
102
99
 
103
100
  if (!mygsl_histogram_equal_bins_p (h1, h2))
104
- {
105
- GSL_ERROR ("histograms have different binning", GSL_EINVAL);
106
- }
107
-
101
+ {
102
+ GSL_ERROR ("histograms have different binning", GSL_EINVAL);
103
+ }
108
104
  for (i = 0; i < h1->n; i++)
109
- {
110
- h1->bin[i] -= h2->bin[i];
111
- }
112
-
105
+ {
106
+ h1->bin[i] -= h2->bin[i];
107
+ }
113
108
  return GSL_SUCCESS;
114
109
 
115
110
  }
116
111
 
117
- /*
112
+ /*
118
113
  * gsl_histogram_mult:
119
114
  * multiply two histograms
120
115
  */
121
116
 
122
- int
117
+ int
123
118
  mygsl_histogram_mul (gsl_histogram * h1, const gsl_histogram * h2)
124
119
  {
125
120
  size_t i;
126
121
 
127
122
  if (!mygsl_histogram_equal_bins_p (h1, h2))
128
- {
129
- GSL_ERROR ("histograms have different binning", GSL_EINVAL);
130
- }
131
-
123
+ {
124
+ GSL_ERROR ("histograms have different binning", GSL_EINVAL);
125
+ }
132
126
  for (i = 0; i < h1->n; i++)
133
- {
134
- h1->bin[i] *= h2->bin[i];
135
- }
136
-
127
+ {
128
+ h1->bin[i] *= h2->bin[i];
129
+ }
137
130
  return GSL_SUCCESS;
138
131
  }
139
- /*
132
+ /*
140
133
  * gsl_histogram_div:
141
134
  * divide two histograms
142
135
  */
143
- int
136
+ int
144
137
  mygsl_histogram_div (gsl_histogram * h1, const gsl_histogram * h2)
145
138
  {
146
139
  size_t i;
147
140
 
148
141
  if (!mygsl_histogram_equal_bins_p (h1, h2))
149
- {
150
- GSL_ERROR ("histograms have different binning", GSL_EINVAL);
151
- }
152
-
142
+ {
143
+ GSL_ERROR ("histograms have different binning", GSL_EINVAL);
144
+ }
153
145
  for (i = 0; i < h1->n; i++)
154
- {
155
- h1->bin[i] /= h2->bin[i];
156
- }
157
-
146
+ {
147
+ h1->bin[i] /= h2->bin[i];
148
+ }
158
149
  return GSL_SUCCESS;
159
150
  }
@@ -9,7 +9,7 @@
9
9
  WITHOUT ANY WARRANTY.
10
10
  */
11
11
 
12
- #include "rb_gsl.h"
12
+ #include "include/rb_gsl.h"
13
13
 
14
14
  static VALUE rb_gsl_ieee_env_setup(VALUE obj)
15
15
  {
@@ -19,15 +19,11 @@ static VALUE rb_gsl_ieee_env_setup(VALUE obj)
19
19
 
20
20
  static VALUE rb_gsl_ieee_fprintf_double(int argc, VALUE *argv, VALUE obj)
21
21
  {
22
- #ifdef RUBY_1_9_LATER
23
22
  rb_io_t *fptr = NULL;
24
- double ftmp;
25
- #else
26
- OpenFile *fptr = NULL;
27
- #endif
28
23
  FILE *fp = NULL;
29
24
  int flag = 0;
30
25
  VALUE vtmp;
26
+ double ftmp;
31
27
 
32
28
  switch (argc) {
33
29
  case 2:
@@ -39,15 +35,11 @@ static VALUE rb_gsl_ieee_fprintf_double(int argc, VALUE *argv, VALUE obj)
39
35
  case T_FILE:
40
36
  GetOpenFile(argv[0], fptr);
41
37
  rb_io_check_writable(fptr);
42
- #ifdef RUBY_1_9_LATER
43
38
  fp = rb_io_stdio_file(fptr);
44
- #else
45
- fp = GetWriteFile(fptr);
46
- #endif
47
39
  break;
48
40
  default:
49
41
  rb_raise(rb_eTypeError, "wrong type argument %s (IO or String expected)",
50
- rb_class2name(CLASS_OF(argv[0])));
42
+ rb_class2name(CLASS_OF(argv[0])));
51
43
  }
52
44
  vtmp = argv[1];
53
45
  break;
@@ -60,13 +52,9 @@ static VALUE rb_gsl_ieee_fprintf_double(int argc, VALUE *argv, VALUE obj)
60
52
  }
61
53
  if (TYPE(vtmp) != T_FLOAT)
62
54
  rb_raise(rb_eTypeError, "wrong argument type %s (Float expected)",
63
- rb_class2name(CLASS_OF(vtmp)));
64
- #ifdef RUBY_1_9_LATER
55
+ rb_class2name(CLASS_OF(vtmp)));
65
56
  ftmp = RFLOAT_VALUE(vtmp);
66
57
  gsl_ieee_fprintf_double(fp, &ftmp);
67
- #else
68
- gsl_ieee_fprintf_double(fp, &(RFLOAT(vtmp)->value));
69
- #endif
70
58
  if (fp == stdout) fprintf(stdout, "\n");
71
59
  if (flag == 1) fclose(fp);
72
60
  return obj;
@@ -86,15 +74,15 @@ void Init_gsl_ieee(VALUE module)
86
74
  mgsl_ieee = rb_define_module_under(module, "IEEE");
87
75
 
88
76
  rb_define_singleton_method(mgsl_ieee, "env_setup",
89
- rb_gsl_ieee_env_setup, 0);
77
+ rb_gsl_ieee_env_setup, 0);
90
78
  rb_define_module_function(module, "ieee_env_setup", rb_gsl_ieee_env_setup, 0);
91
79
  rb_define_singleton_method(mgsl_ieee, "fprintf_double",
92
- rb_gsl_ieee_fprintf_double, -1);
80
+ rb_gsl_ieee_fprintf_double, -1);
93
81
  rb_define_singleton_method(mgsl_ieee, "fprintf",
94
- rb_gsl_ieee_fprintf_double, -1);
82
+ rb_gsl_ieee_fprintf_double, -1);
95
83
  rb_define_singleton_method(mgsl_ieee, "printf",
96
- rb_gsl_ieee_printf_double, -1);
84
+ rb_gsl_ieee_printf_double, -1);
97
85
  rb_define_singleton_method(mgsl_ieee, "printf_double",
98
- rb_gsl_ieee_printf_double, -1);
86
+ rb_gsl_ieee_printf_double, -1);
99
87
 
100
88
  }
@@ -12,40 +12,25 @@
12
12
  #ifndef ___RB_GSL_H___
13
13
  #define ___RB_GSL_H___
14
14
 
15
- #include <stdio.h>
16
15
  #include <string.h>
17
- #include <ctype.h>
18
16
  #include <stdlib.h>
19
17
  #include <gsl/gsl_mode.h>
20
- #include "rb_gsl_common.h"
21
- #include "rb_gsl_math.h"
22
- #include "rb_gsl_complex.h"
23
- #include "rb_gsl_array.h"
24
- #include "rb_gsl_function.h"
25
18
  #include "rb_gsl_poly.h"
26
- #include "rb_gsl_sf.h"
27
19
  #include "rb_gsl_linalg.h"
20
+ #include "rb_gsl_math.h"
21
+ #include "rb_gsl_root.h"
22
+ #include "rb_gsl_histogram.h"
23
+ #include "rb_gsl_sf.h"
28
24
  #include "rb_gsl_eigen.h"
29
25
  #include "rb_gsl_fft.h"
30
26
  #include "rb_gsl_integration.h"
31
27
  #include "rb_gsl_rng.h"
32
28
  #include "rb_gsl_statistics.h"
33
- #include "rb_gsl_histogram.h"
34
29
  #include "rb_gsl_odeiv.h"
35
30
  #include "rb_gsl_interp.h"
36
31
  #include "rb_gsl_cheb.h"
37
- #include "rb_gsl_root.h"
38
32
  #include "rb_gsl_fit.h"
39
33
  #include "rb_gsl_const.h"
40
- #include "rb_gsl_config.h"
41
-
42
-
43
- #include "ruby.h"
44
- #ifdef RUBY_1_9_LATER
45
- #include "ruby/io.h"
46
- #else
47
- #include "rubyio.h"
48
- #endif
49
34
 
50
35
  void Init_gsl_error(VALUE module);
51
36
  void Init_gsl_math(VALUE module);
@@ -82,9 +67,7 @@ void Init_gsl_odeiv(VALUE module);
82
67
  void Init_gsl_interp(VALUE module);
83
68
  void Init_gsl_spline(VALUE module);
84
69
  void Init_gsl_diff(VALUE module);
85
- #ifdef GSL_1_4_9_LATER
86
70
  void Init_gsl_deriv(VALUE module);
87
- #endif
88
71
 
89
72
  void Init_gsl_cheb(VALUE module);
90
73
  void Init_gsl_sum(VALUE module);
@@ -138,9 +121,7 @@ void Init_cqp(VALUE module);
138
121
 
139
122
  void Init_fresnel(VALUE module);
140
123
 
141
- #ifdef GSL_1_9_LATER
142
124
  void Init_bspline(VALUE module);
143
- #endif
144
125
 
145
126
  #ifdef HAVE_ALF_ALF_H
146
127
  #include "alf/alf.h"
@@ -148,11 +129,8 @@ void Init_bspline(VALUE module);
148
129
  void Init_alf(VALUE module);
149
130
  void Init_geometry(VALUE module);
150
131
 
151
- #ifdef GSL_1_14_LATER
152
132
  #include <gsl/gsl_multiset.h>
153
133
  extern VALUE cMultiset;
154
134
  void Init_multiset(VALUE module);
155
135
 
156
136
  #endif
157
-
158
- #endif
@@ -13,15 +13,10 @@
13
13
  #define ___RB_GSL_ARRAY_H___
14
14
 
15
15
  #include <math.h>
16
- #include <gsl/gsl_errno.h>
17
- #include <gsl/gsl_complex.h>
18
16
  #include <gsl/gsl_block.h>
19
- #include <gsl/gsl_vector.h>
20
17
  #include <gsl/gsl_vector_complex.h>
21
18
  #include <gsl/gsl_vector_complex_double.h>
22
- #include <gsl/gsl_matrix.h>
23
19
  #include <gsl/gsl_matrix_complex_double.h>
24
- #include <gsl/gsl_permutation.h>
25
20
  #include <gsl/gsl_combination.h>
26
21
  #include <gsl/gsl_sort.h>
27
22
  #include <gsl/gsl_sort_vector.h>
@@ -29,15 +24,6 @@
29
24
  #include <gsl/gsl_eigen.h>
30
25
  #include <gsl/gsl_odeiv.h>
31
26
  #include <gsl/gsl_blas.h>
32
- #include "rb_gsl_config.h"
33
-
34
- #include "ruby.h"
35
- #ifdef RUBY_1_9_LATER
36
- #include "ruby/io.h"
37
- #else
38
- #include "rubyio.h"
39
- #endif
40
-
41
27
  #include "rb_gsl_common.h"
42
28
 
43
29
  typedef gsl_permutation gsl_index;
@@ -123,14 +109,14 @@ double* get_vector_ptr(VALUE ary, size_t *stride, size_t *n);
123
109
 
124
110
  gsl_matrix_complex* matrix_to_complex(const gsl_matrix *m);
125
111
 
126
- void gsl_matrix_complex_mul(gsl_matrix_complex *mnew, const gsl_matrix_complex *m,
127
- const gsl_matrix_complex *mb);
112
+ void gsl_matrix_complex_mul(gsl_matrix_complex *mnew, const gsl_matrix_complex *m,
113
+ const gsl_matrix_complex *mb);
128
114
  void gsl_matrix_mul(gsl_matrix *mnew, gsl_matrix *m, gsl_matrix *b);
129
- void gsl_matrix_complex_mul_vector(gsl_vector_complex *vnew,
130
- const gsl_matrix_complex *m,
131
- const gsl_vector_complex *v);
132
- void gsl_matrix_mul_vector(gsl_vector *vnew,
133
- const gsl_matrix *m, const gsl_vector *v);
115
+ void gsl_matrix_complex_mul_vector(gsl_vector_complex *vnew,
116
+ const gsl_matrix_complex *m,
117
+ const gsl_vector_complex *v);
118
+ void gsl_matrix_mul_vector(gsl_vector *vnew,
119
+ const gsl_matrix *m, const gsl_vector *v);
134
120
  gsl_vector_complex* vector_to_complex(const gsl_vector *v);
135
121
 
136
122
  gsl_vector* make_vector_clone(const gsl_vector *v);
@@ -164,8 +150,8 @@ gsl_matrix_int_view* rb_gsl_matrix_int_view_alloc();
164
150
  void rb_gsl_matrix_int_view_free(gsl_matrix_int_view *v);
165
151
  VALUE rb_gsl_matrix_to_i(VALUE obj);
166
152
  VALUE rb_gsl_matrix_int_to_f(VALUE obj);
167
- void gsl_matrix_int_mul_vector(gsl_vector_int *vnew,
168
- const gsl_matrix_int *m, const gsl_vector_int *v);
153
+ void gsl_matrix_int_mul_vector(gsl_vector_int *vnew,
154
+ const gsl_matrix_int *m, const gsl_vector_int *v);
169
155
  VALUE rb_gsl_vector_to_i(VALUE obj);
170
156
  VALUE make_rarray_from_cvector_int(const gsl_vector_int *v);
171
157
  VALUE rb_gsl_vector_int_to_f(VALUE obj);
@@ -174,34 +160,24 @@ VALUE rb_gsl_vector_print(VALUE obj);
174
160
  void gsl_vector_print(const gsl_vector *v, VALUE klass);
175
161
  int rbgsl_vector_equal(const gsl_vector *v1, const gsl_vector *v2, double eps);
176
162
 
177
- #ifndef GSL_1_2_LATER
178
- int gsl_matrix_complex_add(gsl_matrix_complex * a, const gsl_matrix_complex * b);
179
- int gsl_matrix_complex_sub(gsl_matrix_complex * a, const gsl_matrix_complex * b);
180
- int gsl_matrix_complex_mul_elements(gsl_matrix_complex * a, const gsl_matrix_complex * b);
181
- int gsl_matrix_complex_div_elements(gsl_matrix_complex * a, const gsl_matrix_complex * b);
182
- int gsl_matrix_complex_scale(gsl_matrix_complex * a, const gsl_complex x);
183
- int gsl_matrix_complex_add_constant(gsl_matrix_complex * a, const gsl_complex x);
184
- int gsl_matrix_complex_add_diagonal(gsl_matrix_complex * a, const gsl_complex x);
185
- #endif
186
-
187
163
  void Init_gsl_vector_init(VALUE module);
188
164
  void Init_gsl_vector_int_init(VALUE module);
189
165
  void Init_gsl_matrix_init(VALUE module);
190
166
  void Init_gsl_matrix_int_init(VALUE module);
191
167
 
192
- gsl_matrix* gsl_matrix_alloc_from_array_sizes(VALUE ary,
193
- VALUE nn1, VALUE nn2);
168
+ gsl_matrix* gsl_matrix_alloc_from_array_sizes(VALUE ary,
169
+ VALUE nn1, VALUE nn2);
194
170
  gsl_matrix* gsl_matrix_alloc_from_arrays(int argc, VALUE *argv);
195
171
 
196
172
  gsl_matrix* gsl_matrix_alloc_from_vector_sizes(VALUE ary,
197
- VALUE nn1, VALUE nn2);
173
+ VALUE nn1, VALUE nn2);
198
174
  gsl_matrix* gsl_matrix_alloc_from_vectors(int argc, VALUE *argv);
199
- gsl_matrix_int* gsl_matrix_int_alloc_from_array_sizes(VALUE ary,
200
- VALUE nn1, VALUE nn2);
175
+ gsl_matrix_int* gsl_matrix_int_alloc_from_array_sizes(VALUE ary,
176
+ VALUE nn1, VALUE nn2);
201
177
  gsl_matrix_int* gsl_matrix_int_alloc_from_arrays(int argc, VALUE *argv);
202
178
 
203
179
  gsl_matrix_int* gsl_matrix_int_alloc_from_vector_sizes(VALUE ary,
204
- VALUE nn1, VALUE nn2);
180
+ VALUE nn1, VALUE nn2);
205
181
  gsl_matrix_int* gsl_matrix_int_alloc_from_vectors(int argc, VALUE *argv);
206
182
 
207
183
  VALUE rb_gsl_matrix_do_something(VALUE obj, void (*f)(gsl_matrix *));