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
@@ -13,9 +13,7 @@
13
13
  #ifndef ___RB_GSL_CHEB_H___
14
14
  #define ___RB_GSL_CHEB_H___
15
15
 
16
- #include "rb_gsl.h"
17
16
  #include "rb_gsl_function.h"
18
- #include <gsl/gsl_math.h>
19
17
  #include <gsl/gsl_chebyshev.h>
20
18
 
21
19
  #endif
@@ -12,26 +12,23 @@
12
12
  #ifndef ___RB_GSL_COMMON_H___
13
13
  #define ___RB_GSL_COMMON_H___
14
14
 
15
- #include "rb_gsl_config.h"
16
- #include "ruby.h"
15
+ #include <ruby.h>
16
+ #include <ruby/io.h>
17
+
17
18
  #include <ctype.h>
18
19
  #include <gsl/gsl_errno.h>
19
20
  #include <gsl/gsl_version.h>
20
21
  #include <gsl/gsl_math.h>
21
22
  #include <gsl/gsl_complex.h>
22
- #include <gsl/gsl_vector.h>
23
23
  #include <gsl/gsl_matrix.h>
24
24
  #include <gsl/gsl_permutation.h>
25
25
  #include <gsl/gsl_ieee_utils.h>
26
- #ifdef HAVE_NARRAY_H
27
- #include "narray.h"
28
26
  #include "rb_gsl_with_narray.h"
29
- #endif
30
27
 
31
28
  EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
32
29
 
33
30
  #ifndef CHECK_FIXNUM
34
- #define CHECK_FIXNUM(x) if(!FIXNUM_P(x))rb_raise(rb_eTypeError,"Fixnum expected");
31
+ #define CHECK_FIXNUM(x) if(!FIXNUM_P(x)) rb_raise(rb_eTypeError,"Fixnum expected");
35
32
  #endif
36
33
 
37
34
  #ifndef Need_Float
@@ -39,9 +36,9 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
39
36
  #endif
40
37
 
41
38
  #ifndef Need_Float2
42
- #define Need_Float2(x,y) do {\
43
- Need_Float(x);\
44
- Need_Float(y);} while (0)
39
+ #define Need_Float2(x,y) do { \
40
+ Need_Float(x); \
41
+ Need_Float(y); } while (0)
45
42
  #endif
46
43
 
47
44
  #ifndef COMPLEX_P
@@ -49,12 +46,12 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
49
46
  #endif
50
47
 
51
48
  #ifndef CHECK_RNG
52
- #define CHECK_RNG(x) if(!rb_obj_is_kind_of(x,cgsl_rng))\
49
+ #define CHECK_RNG(x) if(!rb_obj_is_kind_of(x,cgsl_rng)) \
53
50
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Rng expected)");
54
51
  #endif
55
52
 
56
53
  #ifndef CHECK_COMPLEX
57
- #define CHECK_COMPLEX(x) if(!rb_obj_is_kind_of(x,cgsl_complex))\
54
+ #define CHECK_COMPLEX(x) if(!rb_obj_is_kind_of(x,cgsl_complex)) \
58
55
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Complex expected)");
59
56
  #endif
60
57
 
@@ -63,12 +60,12 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
63
60
  #endif
64
61
 
65
62
  #ifndef CHECK_POLY
66
- #define CHECK_POLY(x) if(!rb_obj_is_kind_of(x,cgsl_poly))\
63
+ #define CHECK_POLY(x) if(!rb_obj_is_kind_of(x,cgsl_poly)) \
67
64
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Poly expected)");
68
65
  #endif
69
66
 
70
67
  #ifndef CHECK_POLY_INT
71
- #define CHECK_POLY_INT(x) if(!rb_obj_is_kind_of(x,cgsl_poly_int))\
68
+ #define CHECK_POLY_INT(x) if(!rb_obj_is_kind_of(x,cgsl_poly_int)) \
72
69
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Poly::Int expected)");
73
70
  #endif
74
71
 
@@ -77,7 +74,7 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
77
74
  #endif
78
75
 
79
76
  #ifndef CHECK_RATIONAL
80
- #define CHECK_RATIONAL(x) if(!rb_obj_is_kind_of(x,cgsl_rational))\
77
+ #define CHECK_RATIONAL(x) if(!rb_obj_is_kind_of(x,cgsl_rational)) \
81
78
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Rational expected)");
82
79
  #endif
83
80
 
@@ -87,7 +84,7 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
87
84
  #endif
88
85
 
89
86
  #ifndef CHECK_BLOCK
90
- #define CHECK_BLOCK(x) if(!rb_obj_is_kind_of(x,cgsl_block))\
87
+ #define CHECK_BLOCK(x) if(!rb_obj_is_kind_of(x,cgsl_block)) \
91
88
  rb_raise(rb_eTypeError, "wrong argument type %s (GSL::Block expected)", rb_class2name(CLASS_OF(x)));
92
89
  #endif
93
90
  #ifndef BLOCK_INT_P
@@ -95,7 +92,7 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
95
92
  #endif
96
93
 
97
94
  #ifndef CHECK_BLOCK_INT
98
- #define CHECK_BLOCK_INT(x) if(!rb_obj_is_kind_of(x,cgsl_block_int))\
95
+ #define CHECK_BLOCK_INT(x) if(!rb_obj_is_kind_of(x,cgsl_block_int)) \
99
96
  rb_raise(rb_eTypeError, "wrong argument type %s (GSL::Block::Int expected)", rb_class2name(CLASS_OF(x)));
100
97
  #endif
101
98
  #ifndef BLOCK_UCHAR_P
@@ -103,7 +100,7 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
103
100
  #endif
104
101
 
105
102
  #ifndef CHECK_BLOCK_UCHAR
106
- #define CHECK_BLOCK_UCHAR(x) if(!rb_obj_is_kind_of(x,cgsl_block_uchar))\
103
+ #define CHECK_BLOCK_UCHAR(x) if(!rb_obj_is_kind_of(x,cgsl_block_uchar)) \
107
104
  rb_raise(rb_eTypeError, "wrong argument type %s (GSL::Block::Byte expected)", rb_class2name(CLASS_OF(x)));
108
105
  #endif
109
106
  /*****/
@@ -125,27 +122,27 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
125
122
  #endif
126
123
 
127
124
  #ifndef VECTOR_ROW_COL
128
- #define VECTOR_ROW_COL(x) ((rb_obj_is_kind_of(x,cgsl_vector_col)||rb_obj_is_kind_of(x,cgsl_vector_int_col))?cgsl_vector_col:cgsl_vector)
125
+ #define VECTOR_ROW_COL(x) ((rb_obj_is_kind_of(x,cgsl_vector_col)||rb_obj_is_kind_of(x,cgsl_vector_int_col)) ? cgsl_vector_col : cgsl_vector)
129
126
  #endif
130
127
 
131
128
  #ifndef CHECK_VECTOR
132
- #define CHECK_VECTOR(x) if(!rb_obj_is_kind_of(x,cgsl_vector))\
129
+ #define CHECK_VECTOR(x) if(!rb_obj_is_kind_of(x,cgsl_vector)) \
133
130
  rb_raise(rb_eTypeError, "wrong argument type %s (GSL::Vector expected)", rb_class2name(CLASS_OF(x)));
134
131
  #endif
135
132
 
136
133
  #ifdef HAVE_NARRAY_H
137
- #define Data_Get_Vector(obj,sval) do {\
138
- if (NA_IsNArray(obj)) {\
139
- /* Convert obj to GSL::Vector::View */\
140
- obj = rb_gsl_na_to_gsl_vector_view_method(obj);\
141
- }\
142
- CHECK_VECTOR(obj);\
143
- Data_Get_Struct(obj,gsl_vector,sval);\
134
+ #define Data_Get_Vector(obj,sval) do { \
135
+ if (NA_IsNArray(obj)) { \
136
+ /* Convert obj to GSL::Vector::View */ \
137
+ obj = rb_gsl_na_to_gsl_vector_view_method(obj); \
138
+ } \
139
+ CHECK_VECTOR(obj); \
140
+ Data_Get_Struct(obj,gsl_vector,sval); \
144
141
  } while (0)
145
142
  #else
146
- #define Data_Get_Vector(obj,sval) do {\
147
- CHECK_VECTOR(obj);\
148
- Data_Get_Struct(obj,gsl_vector,sval);\
143
+ #define Data_Get_Vector(obj,sval) do { \
144
+ CHECK_VECTOR(obj); \
145
+ Data_Get_Struct(obj,gsl_vector,sval); \
149
146
  } while (0)
150
147
  #endif
151
148
 
@@ -167,11 +164,11 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
167
164
  #endif
168
165
 
169
166
  #ifndef VECTOR_INT_ROW_COL
170
- #define VECTOR_INT_ROW_COL(x) (VECTOR_INT_ROW_P(x)?cgsl_vector_int:cgsl_vector_int_col)
167
+ #define VECTOR_INT_ROW_COL(x) (VECTOR_INT_ROW_P(x) ? cgsl_vector_int : cgsl_vector_int_col)
171
168
  #endif
172
169
 
173
170
  #ifndef CHECK_VECTOR_INT
174
- #define CHECK_VECTOR_INT(x) if(!rb_obj_is_kind_of(x,cgsl_vector_int))\
171
+ #define CHECK_VECTOR_INT(x) if(!rb_obj_is_kind_of(x,cgsl_vector_int)) \
175
172
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Vector::Int expected)");
176
173
  #endif
177
174
 
@@ -189,11 +186,11 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
189
186
  #endif
190
187
 
191
188
  #ifndef VECTOR_COMPLEX_ROW_COL
192
- #define VECTOR_COMPLEX_ROW_COL(x) (VECTOR_COMPLEX_ROW_P(x)?cgsl_vector_complex:cgsl_vector_complex_col)
189
+ #define VECTOR_COMPLEX_ROW_COL(x) (VECTOR_COMPLEX_ROW_P(x) ? cgsl_vector_complex : cgsl_vector_complex_col)
193
190
  #endif
194
191
 
195
192
  #ifndef CHECK_VECTOR_COMPLEX
196
- #define CHECK_VECTOR_COMPLEX(x) if(!rb_obj_is_kind_of(x,cgsl_vector_complex))\
193
+ #define CHECK_VECTOR_COMPLEX(x) if(!rb_obj_is_kind_of(x,cgsl_vector_complex)) \
197
194
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Vector::Complex expected)");
198
195
  #endif
199
196
 
@@ -202,13 +199,13 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
202
199
  #endif
203
200
 
204
201
  #ifndef CHECK_MATRIX
205
- #define CHECK_MATRIX(x) if(!rb_obj_is_kind_of(x,cgsl_matrix))\
202
+ #define CHECK_MATRIX(x) if(!rb_obj_is_kind_of(x,cgsl_matrix)) \
206
203
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Matrix expected)");
207
204
  #endif
208
205
 
209
- #define Data_Get_Matrix(obj,sval) do {\
210
- CHECK_MATRIX(obj);\
211
- Data_Get_Struct(obj,gsl_matrix,sval);\
206
+ #define Data_Get_Matrix(obj,sval) do { \
207
+ CHECK_MATRIX(obj); \
208
+ Data_Get_Struct(obj,gsl_matrix,sval); \
212
209
  } while (0)
213
210
 
214
211
 
@@ -217,7 +214,7 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
217
214
  #endif
218
215
 
219
216
  #ifndef CHECK_MATRIX_INT
220
- #define CHECK_MATRIX_INT(x) if(!rb_obj_is_kind_of(x,cgsl_matrix_int))\
217
+ #define CHECK_MATRIX_INT(x) if(!rb_obj_is_kind_of(x,cgsl_matrix_int)) \
221
218
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Matrix::Int expected)");
222
219
  #endif
223
220
 
@@ -226,7 +223,7 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
226
223
  #endif
227
224
 
228
225
  #ifndef CHECK_MATRIX_COMPLEX
229
- #define CHECK_MATRIX_COMPLEX(x) if(!rb_obj_is_kind_of(x,cgsl_matrix_complex))\
226
+ #define CHECK_MATRIX_COMPLEX(x) if(!rb_obj_is_kind_of(x,cgsl_matrix_complex)) \
230
227
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Matrix::Complex expected)");
231
228
  #endif
232
229
 
@@ -235,7 +232,7 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
235
232
  #endif
236
233
 
237
234
  #ifndef CHECK_TENSOR
238
- #define CHECK_TENSOR(x) if(CLASS_OF(x)!=cgsl_tensor)\
235
+ #define CHECK_TENSOR(x) if(CLASS_OF(x)!=cgsl_tensor) \
239
236
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Tensor expected)");
240
237
  #endif
241
238
 
@@ -244,7 +241,7 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
244
241
  #endif
245
242
 
246
243
  #ifndef CHECK_TENSOR_INT
247
- #define CHECK_TENSOR_INT(x) if(CLASS_OF(x)!=cgsl_tensor_int)\
244
+ #define CHECK_TENSOR_INT(x) if(CLASS_OF(x)!=cgsl_tensor_int) \
248
245
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Tensor::Int expected)");
249
246
  #endif
250
247
 
@@ -253,7 +250,7 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
253
250
  #endif
254
251
 
255
252
  #ifndef CHECK_PERMUTATION
256
- #define CHECK_PERMUTATION(x) if(!rb_obj_is_kind_of(x,cgsl_permutation))\
253
+ #define CHECK_PERMUTATION(x) if(!rb_obj_is_kind_of(x,cgsl_permutation)) \
257
254
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Permutation expected)");
258
255
  #endif
259
256
 
@@ -262,7 +259,7 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
262
259
  #endif
263
260
 
264
261
  #ifndef CHECK_PROC
265
- #define CHECK_PROC(x) if(!rb_obj_is_kind_of(x,rb_cProc))\
262
+ #define CHECK_PROC(x) if(!rb_obj_is_kind_of(x,rb_cProc)) \
266
263
  rb_raise(rb_eTypeError, "wrong argument type (Proc expected)");
267
264
  #endif
268
265
 
@@ -271,7 +268,7 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
271
268
  #endif
272
269
 
273
270
  #ifndef CHECK_FUNCTION
274
- #define CHECK_FUNCTION(x) if(!rb_obj_is_kind_of(x,cgsl_function))\
271
+ #define CHECK_FUNCTION(x) if(!rb_obj_is_kind_of(x,cgsl_function)) \
275
272
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Function expected)");
276
273
  #endif
277
274
 
@@ -280,7 +277,7 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
280
277
  #endif
281
278
 
282
279
  #ifndef CHECK_FUNCTION_FDF
283
- #define CHECK_FUNCTION_FDF(x) if(!rb_obj_is_kind_of(x,cgsl_function_fdf))\
280
+ #define CHECK_FUNCTION_FDF(x) if(!rb_obj_is_kind_of(x,cgsl_function_fdf)) \
284
281
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Function_fdf expected)");
285
282
  #endif
286
283
 
@@ -289,12 +286,24 @@ EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
289
286
  #endif
290
287
 
291
288
  #ifndef CHECK_HISTOGRAM
292
- #define CHECK_HISTOGRAM(x) if(!rb_obj_is_kind_of(x,cgsl_histogram))\
289
+ #define CHECK_HISTOGRAM(x) if(!rb_obj_is_kind_of(x,cgsl_histogram)) \
293
290
  rb_raise(rb_eTypeError, "wrong argument type (GSL::Histogram expected)");
294
291
  #endif
295
292
 
293
+ #ifndef RBGSL_SET_CLASS
294
+ #ifdef RB_OBJ_WRITE
295
+ #define RBGSL_SET_CLASS0(obj0, cls) RB_OBJ_WRITE(obj0, &(RBASIC_CLASS(obj0)), cls)
296
+ #else
297
+ #define RBGSL_SET_CLASS0(obj0, cls) RBASIC(obj0)->klass = cls
298
+ #endif
299
+ #define RBGSL_SET_CLASS(obj, cls) do { \
300
+ VALUE _obj_ = (obj); \
301
+ RBGSL_SET_CLASS0(_obj_, cls); \
302
+ } while (0)
303
+ #endif
304
+
296
305
  void rb_gsl_error_handler(const char *reason, const char *file,
297
- int line, int gsl_errno);
306
+ int line, int gsl_errno);
298
307
 
299
308
  FILE* rb_gsl_open_writefile(VALUE io, int *flag);
300
309
  FILE* rb_gsl_open_readfile(VALUE io, int *flag);
@@ -304,17 +313,12 @@ VALUE rb_gsl_obj_read_only(int argc, VALUE *argv, VALUE obj);
304
313
  int str_tail_grep(const char *s0, const char *s1);
305
314
  int str_head_grep(const char *s0, const char *s1);
306
315
 
307
- #ifdef RUBY_1_8_LATER
308
- #define RB_GSL_MAKE_PROC rb_block_proc()
309
- #define STR2CHARPTR StringValuePtr
310
- #else
311
- #define RB_GSL_MAKE_PROC rb_f_lambda()
312
- #define STR2CHARPTR STR2CSTR
316
+ #ifndef STR2CSTR
317
+ #define STR2CSTR StringValuePtr
313
318
  #endif
314
319
 
315
- // Added 2010/Sep/29
316
- #ifdef RUBY_1_9_2_LATER
317
- #define STR2CSTR StringValuePtr
320
+ #ifndef RCLASS_SUPER
321
+ #define RCLASS_SUPER(cls) RCLASS(cls)->super
318
322
  #endif
319
323
 
320
324
  void make_graphcommand(char *command, VALUE hash);
@@ -328,10 +332,6 @@ gsl_vector_int* mygsl_vector_int_down(gsl_vector_int *p);
328
332
  void mygsl_vector_int_up2(gsl_vector_int *pnew, gsl_vector_int *p);
329
333
  gsl_vector_int* mygsl_vector_int_up(gsl_vector_int *p);
330
334
 
331
- #ifndef GSL_1_3_LATER
332
- int gsl_fcmp(const double x1, const double x2, const double epsilon);
333
- #endif
334
-
335
335
  size_t count_columns(const char *str);
336
336
  char* str_scan_double(const char *str, double *val);
337
337
  char* str_scan_int(const char *str, int *val);
@@ -14,7 +14,7 @@
14
14
  #define ___RB_GSL_COMPLEX_H___
15
15
 
16
16
  #include <math.h>
17
- #include "ruby.h"
17
+ #include <ruby.h>
18
18
  #include <gsl/gsl_complex.h>
19
19
  #include <gsl/gsl_complex_math.h>
20
20
 
@@ -13,15 +13,9 @@
13
13
  #ifndef ___RB_GSL_CONST_H___
14
14
  #define ___RB_GSL_CONST_H___
15
15
 
16
- #include "ruby.h"
17
16
  #include "rb_gsl.h"
18
- #ifdef GSL_CONST_OLD
19
- #include <gsl/gsl_const_mks.h>
20
- #include <gsl/gsl_const_cgs.h>
21
- #else
22
17
  #include <gsl/gsl_const_mksa.h>
23
18
  #include <gsl/gsl_const_cgsm.h>
24
- #endif
25
19
  #include <gsl/gsl_const_num.h>
26
20
 
27
21
  EXTERN VALUE mgsl_const_mks, mgsl_const_cgs;
@@ -0,0 +1,6 @@
1
+ #ifndef ___RB_GSL_DIRAC_H___
2
+ #define ___RB_GSL_DIRAC_H___
3
+
4
+ #include "rb_gsl.h"
5
+
6
+ #endif
@@ -12,6 +12,6 @@
12
12
  #ifndef ___RB_GSL_EIGEN_H___
13
13
  #define ___RB_GSL_EIGEN_H___
14
14
 
15
- #include "gsl/gsl_eigen.h"
15
+ #include <gsl/gsl_eigen.h>
16
16
 
17
17
  #endif
@@ -13,25 +13,12 @@
13
13
  #ifndef ___RB_GSL_FFT_H___
14
14
  #define ___RB_GSL_FFT_H___
15
15
 
16
-
17
- #include <stdio.h>
18
- #include <math.h>
19
- #include "ruby.h"
20
- #include <gsl/gsl_vector.h>
21
- #include <gsl/gsl_errno.h>
22
16
  #include <gsl/gsl_fft.h>
23
17
  #include <gsl/gsl_fft_complex.h>
24
18
  #include <gsl/gsl_fft_real.h>
25
19
  #include <gsl/gsl_fft_halfcomplex.h>
26
- #include "rb_gsl_array.h"
27
20
  #include "rb_gsl.h"
28
21
 
29
- #include "rb_gsl_complex.h"
30
- #include "rb_gsl_array.h"
31
- #ifdef HAVE_NARRAY_H
32
- #include "narray.h"
33
- #endif
34
-
35
22
  typedef struct
36
23
  {
37
24
  size_t n;
@@ -13,8 +13,6 @@
13
13
  #ifndef ___RB_GSL_FIT_H___
14
14
  #define ___RB_GSL_FIT_H___
15
15
 
16
- #include "ruby.h"
17
- #include <gsl/gsl_blas.h>
18
16
  #include <gsl/gsl_fit.h>
19
17
  #include <gsl/gsl_multifit.h>
20
18
  #include <gsl/gsl_multifit_nlin.h>
@@ -13,10 +13,6 @@
13
13
  #ifndef ___RB_GSL_FUNCTION_H___
14
14
  #define ___RB_GSL_FUNCTION_H___
15
15
 
16
- #include "ruby.h"
17
- #include <gsl/gsl_vector.h>
18
- #include <gsl/gsl_integration.h>
19
- #include <gsl/gsl_math.h>
20
16
  #include "rb_gsl.h"
21
17
 
22
18
  EXTERN VALUE cgsl_function;
@@ -14,8 +14,6 @@
14
14
  #define ___RB_GSL_GRAPH_H___
15
15
 
16
16
  #include "rb_gsl.h"
17
- #include "rb_gsl_array.h"
18
- #include "rb_gsl_common.h"
19
17
 
20
18
  typedef struct __rb_gsl_plot {
21
19
  VALUE xdata, ydata;
@@ -60,8 +58,8 @@ typedef struct __rb_gsl_plot {
60
58
  VALUE reposition;
61
59
  VALUE blankout;
62
60
  VALUE O;
63
- VALUE X, Y;
64
-
61
+ VALUE X, Y;
62
+
65
63
  } gsl_graph;
66
64
 
67
65
  gsl_graph* gsl_graph_new();
@@ -29,21 +29,21 @@ typedef struct {
29
29
  } mygsl_histogram2d_view;
30
30
 
31
31
  #ifndef HISTOGRAM2D_P
32
- #define HISTOGRAM2D_P(x) (rb_obj_is_kind_of(x,cgsl_histogram2d)?1:0)
32
+ #define HISTOGRAM2D_P(x) (rb_obj_is_kind_of(x,cgsl_histogram2d) ? 1 : 0)
33
33
  #endif
34
34
 
35
35
  #ifndef CHECK_HISTOGRAM2D
36
- #define CHECK_HISTOGRAM2D(x) if(!rb_obj_is_kind_of(x,cgsl_histogram2d))\
36
+ #define CHECK_HISTOGRAM2D(x) if(!rb_obj_is_kind_of(x,cgsl_histogram2d)) \
37
37
  rb_raise(rb_eTypeError, "wrong type (Histogram2d expected)");
38
38
  #endif
39
39
 
40
40
 
41
41
  #ifndef HISTOGRAM3D_P
42
- #define HISTOGRAM3D_P(x) (rb_obj_is_kind_of(x,cgsl_histogram3d)?1:0)
42
+ #define HISTOGRAM3D_P(x) (rb_obj_is_kind_of(x,cgsl_histogram3d) ? 1 : 0)
43
43
  #endif
44
44
 
45
45
  #ifndef CHECK_HISTOGRAM3D
46
- #define CHECK_HISTOGRAM3D(x) if(!rb_obj_is_kind_of(x,cgsl_histogram3d))\
46
+ #define CHECK_HISTOGRAM3D(x) if(!rb_obj_is_kind_of(x,cgsl_histogram3d)) \
47
47
  rb_raise(rb_eTypeError, "wrong type (Histogram3d expected)");
48
48
  #endif
49
49
 
@@ -51,13 +51,13 @@ typedef struct {
51
51
 
52
52
  int
53
53
  mygsl_histogram_equal_bins_p (const gsl_histogram * h1, const gsl_histogram * h2);
54
- int
54
+ int
55
55
  mygsl_histogram_add (gsl_histogram * h1, const gsl_histogram * h2);
56
- int
56
+ int
57
57
  mygsl_histogram_sub (gsl_histogram * h1, const gsl_histogram * h2);
58
- int
58
+ int
59
59
  mygsl_histogram_mul (gsl_histogram * h1, const gsl_histogram * h2);
60
- int
60
+ int
61
61
  mygsl_histogram_div (gsl_histogram * h1, const gsl_histogram * h2);
62
62
 
63
63
  #endif