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
data/ext/{sf_airy.c → gsl_native/sf_airy.c} RENAMED
@@ -10,8 +10,8 @@
10
10
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
11
11
  */
12
12
 
13
- #include "rb_gsl_sf.h"
14
- #include "rb_gsl_array.h"
13
+ #include "include/rb_gsl_sf.h"
14
+ #include "include/rb_gsl_array.h"
15
15
 
16
16
  /* m: precision, d(double), s(single), a(apporox) */
17
17
  static VALUE rb_gsl_sf_airy_Ai(int argc, VALUE *argv, VALUE obj)
data/ext/{sf_bessel.c → gsl_native/sf_bessel.c} RENAMED
@@ -10,7 +10,7 @@
10
10
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
11
11
  */
12
12
 
13
- #include "rb_gsl_sf.h"
13
+ #include "include/rb_gsl_sf.h"
14
14
  EXTERN VALUE cgsl_vector;
15
15
 
16
16
  /* Cylindrical Bessel Functions */
@@ -45,7 +45,7 @@ static VALUE rb_gsl_sf_bessel_Jn_e(VALUE obj, VALUE n, VALUE x)
45
45
  }
46
46
 
47
47
  static VALUE rb_gsl_sf_bessel_Xn_array(VALUE obj, VALUE n0, VALUE n1, VALUE x,
48
- int (*f)(int, int, double, double[]))
48
+ int (*f)(int, int, double, double[]))
49
49
  {
50
50
  int nmin, nmax, n;
51
51
  gsl_vector *v = NULL;
@@ -167,7 +167,7 @@ static VALUE rb_gsl_sf_bessel_In_scaled_e(VALUE obj, VALUE n, VALUE x)
167
167
  }
168
168
 
169
169
  static VALUE rb_gsl_sf_bessel_In_scaled_array(VALUE obj, VALUE n0, VALUE n1,
170
- VALUE x)
170
+ VALUE x)
171
171
  {
172
172
  return rb_gsl_sf_bessel_Xn_array(obj, n0, n1, x, gsl_sf_bessel_In_scaled_array);
173
173
  }
@@ -239,7 +239,7 @@ static VALUE rb_gsl_sf_bessel_Kn_scaled_e(VALUE obj, VALUE n, VALUE x)
239
239
  }
240
240
 
241
241
  static VALUE rb_gsl_sf_bessel_Kn_scaled_array(VALUE obj, VALUE n0, VALUE n1,
242
- VALUE x)
242
+ VALUE x)
243
243
  {
244
244
  return rb_gsl_sf_bessel_Xn_array(obj, n0, n1, x, gsl_sf_bessel_Kn_scaled_array);
245
245
  }
@@ -286,7 +286,7 @@ static VALUE rb_gsl_sf_bessel_jl_e(VALUE obj, VALUE n, VALUE x)
286
286
  }
287
287
 
288
288
  static VALUE rb_gsl_sf_bessel_xl_array(VALUE obj, VALUE n1, VALUE x,
289
- int (*f)(int, double, double[]))
289
+ int (*f)(int, double, double[]))
290
290
  {
291
291
  int nmax, n;
292
292
  // local variable "status" declared and set, but never used
@@ -446,7 +446,7 @@ static VALUE rb_gsl_sf_bessel_kl_scaled_e(VALUE obj, VALUE n, VALUE x)
446
446
  }
447
447
 
448
448
  static VALUE rb_gsl_sf_bessel_kl_scaled_array(VALUE obj, VALUE n1,
449
- VALUE x)
449
+ VALUE x)
450
450
  {
451
451
  return rb_gsl_sf_bessel_xl_array(obj, n1, x, gsl_sf_bessel_kl_scaled_array);
452
452
  }
@@ -497,7 +497,7 @@ static VALUE rb_gsl_sf_bessel_sequence_Jnu_e(int argc, VALUE *argv, VALUE obj)
497
497
  break;
498
498
  default:
499
499
  rb_raise(rb_eTypeError, "wrong argument type %s (String or Fixnum expected)",
500
- rb_class2name(CLASS_OF(m)));
500
+ rb_class2name(CLASS_OF(m)));
501
501
  break;
502
502
  }
503
503
  default:
data/ext/{sf_clausen.c → gsl_native/sf_clausen.c} RENAMED
@@ -9,7 +9,7 @@
9
9
  WITHOUT ANY WARRANTY.
10
10
  */
11
11
 
12
- #include "rb_gsl_sf.h"
12
+ #include "include/rb_gsl_sf.h"
13
13
 
14
14
  static VALUE rb_gsl_sf_clausen(VALUE obj, VALUE x)
15
15
  {
data/ext/{sf_coulomb.c → gsl_native/sf_coulomb.c} RENAMED
@@ -10,7 +10,7 @@
10
10
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
11
11
  */
12
12
 
13
- #include "rb_gsl_sf.h"
13
+ #include "include/rb_gsl_sf.h"
14
14
  EXTERN VALUE cgsl_vector;
15
15
 
16
16
  static VALUE rb_gsl_sf_hydrogenicR_1(VALUE obj, VALUE Z, VALUE r)
@@ -23,15 +23,15 @@ static VALUE rb_gsl_sf_hydrogenicR_1_e(VALUE obj, VALUE Z, VALUE r)
23
23
  return rb_gsl_sf_eval_e_double2(gsl_sf_hydrogenicR_1_e, Z, r);
24
24
  }
25
25
 
26
- static VALUE rb_gsl_sf_hydrogenicR(VALUE obj, VALUE n, VALUE l,
27
- VALUE Z, VALUE r)
26
+ static VALUE rb_gsl_sf_hydrogenicR(VALUE obj, VALUE n, VALUE l,
27
+ VALUE Z, VALUE r)
28
28
  {
29
29
  return rb_float_new(gsl_sf_hydrogenicR(FIX2INT(n), FIX2INT(l),
30
- NUM2DBL(Z), NUM2DBL(r)));
30
+ NUM2DBL(Z), NUM2DBL(r)));
31
31
  }
32
32
 
33
- static VALUE rb_gsl_sf_hydrogenicR_e(VALUE obj, VALUE n, VALUE l,
34
- VALUE Z, VALUE r)
33
+ static VALUE rb_gsl_sf_hydrogenicR_e(VALUE obj, VALUE n, VALUE l,
34
+ VALUE Z, VALUE r)
35
35
  {
36
36
  gsl_sf_result *rslt = NULL;
37
37
  VALUE v;
@@ -41,12 +41,12 @@ static VALUE rb_gsl_sf_hydrogenicR_e(VALUE obj, VALUE n, VALUE l,
41
41
  Need_Float(Z); Need_Float(r);
42
42
  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
43
43
  /*status =*/ gsl_sf_hydrogenicR_e(FIX2INT(n), FIX2INT(l),
44
- NUM2DBL(Z), NUM2DBL(r), rslt);
44
+ NUM2DBL(Z), NUM2DBL(r), rslt);
45
45
  return v;
46
46
  }
47
47
 
48
48
  static VALUE rb_gsl_sf_coulomb_wave_FG_e(VALUE obj, VALUE eta, VALUE x,
49
- VALUE L_F, VALUE k)
49
+ VALUE L_F, VALUE k)
50
50
 
51
51
  {
52
52
  gsl_sf_result *F, *Fp, *G, *Gp;
@@ -60,13 +60,13 @@ static VALUE rb_gsl_sf_coulomb_wave_FG_e(VALUE obj, VALUE eta, VALUE x,
60
60
  vG = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, G);
61
61
  vGp = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, Gp);
62
62
  status = gsl_sf_coulomb_wave_FG_e(NUM2DBL(eta), NUM2DBL(x), NUM2DBL(L_F),
63
- FIX2INT(k), F, Fp, G, Gp, &exp_F, &exp_G);
63
+ FIX2INT(k), F, Fp, G, Gp, &exp_F, &exp_G);
64
64
  return rb_ary_new3(7, vF, vFp, vG, vGp,
65
- rb_float_new(exp_F), rb_float_new(exp_G), INT2FIX(status));
65
+ rb_float_new(exp_F), rb_float_new(exp_G), INT2FIX(status));
66
66
  }
67
67
 
68
68
  static VALUE rb_gsl_sf_coulomb_wave_F_array(VALUE obj, VALUE Lmin, VALUE kmax,
69
- VALUE eta, VALUE x)
69
+ VALUE eta, VALUE x)
70
70
  {
71
71
  double F_exponent;
72
72
  int status;
@@ -76,15 +76,15 @@ static VALUE rb_gsl_sf_coulomb_wave_F_array(VALUE obj, VALUE Lmin, VALUE kmax,
76
76
  Need_Float(Lmin); Need_Float(eta); Need_Float(x);
77
77
  size = FIX2INT(kmax);
78
78
  v = gsl_vector_alloc(size);
79
- status = gsl_sf_coulomb_wave_F_array(NUM2DBL(Lmin), size, NUM2DBL(eta),
80
- NUM2DBL(x), v->data, &F_exponent);
81
-
82
- return rb_ary_new3(3, Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v),
83
- rb_float_new(F_exponent), INT2FIX(status));
79
+ status = gsl_sf_coulomb_wave_F_array(NUM2DBL(Lmin), size, NUM2DBL(eta),
80
+ NUM2DBL(x), v->data, &F_exponent);
81
+
82
+ return rb_ary_new3(3, Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v),
83
+ rb_float_new(F_exponent), INT2FIX(status));
84
84
  }
85
85
 
86
86
  static VALUE rb_gsl_sf_coulomb_wave_FG_array(VALUE obj, VALUE Lmin, VALUE kmax,
87
- VALUE eta, VALUE x)
87
+ VALUE eta, VALUE x)
88
88
  {
89
89
  double F_exponent, G_exponent;
90
90
  int status;
@@ -97,18 +97,18 @@ static VALUE rb_gsl_sf_coulomb_wave_FG_array(VALUE obj, VALUE Lmin, VALUE kmax,
97
97
  vf = gsl_vector_alloc(size);
98
98
  vg = gsl_vector_alloc(size);
99
99
 
100
- status = gsl_sf_coulomb_wave_FG_array(NUM2DBL(Lmin), size, NUM2DBL(eta),
101
- NUM2DBL(x), vf->data, vg->data,
102
- &F_exponent, &G_exponent);
100
+ status = gsl_sf_coulomb_wave_FG_array(NUM2DBL(Lmin), size, NUM2DBL(eta),
101
+ NUM2DBL(x), vf->data, vg->data,
102
+ &F_exponent, &G_exponent);
103
103
  fary = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vf);
104
104
  gary = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vf);
105
105
  return rb_ary_new3(5, fary, gary,
106
- rb_float_new(F_exponent), rb_float_new(G_exponent),
107
- INT2FIX(status));
106
+ rb_float_new(F_exponent), rb_float_new(G_exponent),
107
+ INT2FIX(status));
108
108
  }
109
109
 
110
110
  static VALUE rb_gsl_sf_coulomb_wave_FGp_array(VALUE obj, VALUE Lmin, VALUE kmax,
111
- VALUE eta, VALUE x)
111
+ VALUE eta, VALUE x)
112
112
  {
113
113
  double F_exponent, G_exponent;
114
114
  int status;
@@ -123,21 +123,21 @@ static VALUE rb_gsl_sf_coulomb_wave_FGp_array(VALUE obj, VALUE Lmin, VALUE kmax,
123
123
  vg = gsl_vector_alloc(size);
124
124
  vgp = gsl_vector_alloc(size);
125
125
 
126
- status = gsl_sf_coulomb_wave_FGp_array(NUM2DBL(Lmin), size, NUM2DBL(eta),
127
- NUM2DBL(x), vf->data, vfp->data,
128
- vg->data, vgp->data,
129
- &F_exponent, &G_exponent);
126
+ status = gsl_sf_coulomb_wave_FGp_array(NUM2DBL(Lmin), size, NUM2DBL(eta),
127
+ NUM2DBL(x), vf->data, vfp->data,
128
+ vg->data, vgp->data,
129
+ &F_exponent, &G_exponent);
130
130
  fary = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vf);
131
- fpary =Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vfp);
131
+ fpary = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vfp);
132
132
  gary = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vg);
133
- gpary =Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vgp);
133
+ gpary = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vgp);
134
134
  return rb_ary_new3(7, fary, fpary, gary, gpary,
135
- rb_float_new(F_exponent), rb_float_new(G_exponent),
136
- INT2FIX(status));
135
+ rb_float_new(F_exponent), rb_float_new(G_exponent),
136
+ INT2FIX(status));
137
137
  }
138
138
 
139
139
  static VALUE rb_gsl_sf_coulomb_wave_sphF_array(VALUE obj, VALUE Lmin, VALUE kmax,
140
- VALUE eta, VALUE x)
140
+ VALUE eta, VALUE x)
141
141
  {
142
142
  int status;
143
143
  size_t size;
@@ -147,11 +147,11 @@ static VALUE rb_gsl_sf_coulomb_wave_sphF_array(VALUE obj, VALUE Lmin, VALUE kmax
147
147
  size = FIX2INT(kmax);
148
148
  v = gsl_vector_alloc(size);
149
149
  v2 = gsl_vector_alloc(size);
150
- status = gsl_sf_coulomb_wave_sphF_array(NUM2DBL(Lmin), size, NUM2DBL(eta),
151
- NUM2DBL(x), v->data, v2->data);
152
- return rb_ary_new3(3, Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v),
153
- Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v2),
154
- INT2FIX(status));
150
+ status = gsl_sf_coulomb_wave_sphF_array(NUM2DBL(Lmin), size, NUM2DBL(eta),
151
+ NUM2DBL(x), v->data, v2->data);
152
+ return rb_ary_new3(3, Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v),
153
+ Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v2),
154
+ INT2FIX(status));
155
155
  }
156
156
 
157
157
  static VALUE rb_gsl_sf_coulomb_CL_e(VALUE obj, VALUE L, VALUE eta)
@@ -160,14 +160,14 @@ static VALUE rb_gsl_sf_coulomb_CL_e(VALUE obj, VALUE L, VALUE eta)
160
160
  }
161
161
 
162
162
  static VALUE rb_gsl_sf_coulomb_CL_array(VALUE obj, VALUE Lmin, VALUE kmax,
163
- VALUE eta)
163
+ VALUE eta)
164
164
  {
165
165
  gsl_vector *v = NULL;
166
166
  size_t size;
167
167
  // local variable "status" declared and set, but never used
168
168
  //int status;
169
169
  CHECK_FIXNUM(kmax);
170
- Need_Float(Lmin); Need_Float(eta);
170
+ Need_Float(Lmin); Need_Float(eta);
171
171
  size = FIX2INT(kmax);
172
172
  v = gsl_vector_alloc(size);
173
173
  /*status =*/ gsl_sf_coulomb_CL_array(NUM2DBL(Lmin), size, NUM2DBL(eta), v->data);
@@ -191,7 +191,7 @@ void Init_gsl_sf_coulomb(VALUE module)
191
191
  rb_define_module_function(module, "coulomb_CL_array", rb_gsl_sf_coulomb_CL_array, 3);
192
192
 
193
193
  mgsl_sf_coulomb = rb_define_module_under(module, "Coulomb");
194
-
194
+
195
195
  rb_define_module_function(mgsl_sf_coulomb, "hydrogenicR_1", rb_gsl_sf_hydrogenicR_1, 2);
196
196
  rb_define_module_function(mgsl_sf_coulomb, "hydrogenicR_1_e", rb_gsl_sf_hydrogenicR_1_e, 2);
197
197
  rb_define_module_function(mgsl_sf_coulomb, "hydrogenicR", rb_gsl_sf_hydrogenicR, 4);
data/ext/{sf_coupling.c → gsl_native/sf_coupling.c} RENAMED
@@ -10,18 +10,18 @@
10
10
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
11
11
  */
12
12
 
13
- #include "rb_gsl_sf.h"
13
+ #include "include/rb_gsl_sf.h"
14
14
 
15
- static VALUE rb_gsl_sf_coupling_3j(VALUE obj, VALUE two_ja, VALUE two_jb, VALUE two_jc, VALUE two_ma, VALUE two_mb, VALUE two_mc)
15
+ static VALUE rb_gsl_sf_coupling_3j(VALUE obj, VALUE two_ja, VALUE two_jb, VALUE two_jc, VALUE two_ma, VALUE two_mb, VALUE two_mc)
16
16
  {
17
17
  CHECK_FIXNUM(two_ja); CHECK_FIXNUM(two_jb); CHECK_FIXNUM(two_jc);
18
18
  CHECK_FIXNUM(two_ma); CHECK_FIXNUM(two_mb); CHECK_FIXNUM(two_mc);
19
19
  return rb_float_new(gsl_sf_coupling_3j(FIX2INT(two_ja), FIX2INT(two_jb),
20
- FIX2INT(two_jc), FIX2INT(two_ma),
21
- FIX2INT(two_mb), FIX2INT(two_mc)));
20
+ FIX2INT(two_jc), FIX2INT(two_ma),
21
+ FIX2INT(two_mb), FIX2INT(two_mc)));
22
22
  }
23
23
 
24
- static VALUE rb_gsl_sf_coupling_3j_e(VALUE obj, VALUE two_ja, VALUE two_jb, VALUE two_jc, VALUE two_ma, VALUE two_mb, VALUE two_mc)
24
+ static VALUE rb_gsl_sf_coupling_3j_e(VALUE obj, VALUE two_ja, VALUE two_jb, VALUE two_jc, VALUE two_ma, VALUE two_mb, VALUE two_mc)
25
25
  {
26
26
  gsl_sf_result *rslt = NULL;
27
27
  VALUE v;
@@ -31,23 +31,23 @@ static VALUE rb_gsl_sf_coupling_3j_e(VALUE obj, VALUE two_ja, VALUE two_jb, VALU
31
31
  CHECK_FIXNUM(two_ma); CHECK_FIXNUM(two_mb); CHECK_FIXNUM(two_mc);
32
32
  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
33
33
  /*status =*/ gsl_sf_coupling_3j_e(FIX2INT(two_ja), FIX2INT(two_jb),
34
- FIX2INT(two_jc), FIX2INT(two_ma),
35
- FIX2INT(two_mb), FIX2INT(two_mc),
36
- rslt);
34
+ FIX2INT(two_jc), FIX2INT(two_ma),
35
+ FIX2INT(two_mb), FIX2INT(two_mc),
36
+ rslt);
37
37
 
38
38
  return v;
39
39
  }
40
40
 
41
- static VALUE rb_gsl_sf_coupling_6j(VALUE obj, VALUE two_ja, VALUE two_jb, VALUE two_jc, VALUE two_jd, VALUE two_je, VALUE two_jf)
41
+ static VALUE rb_gsl_sf_coupling_6j(VALUE obj, VALUE two_ja, VALUE two_jb, VALUE two_jc, VALUE two_jd, VALUE two_je, VALUE two_jf)
42
42
  {
43
43
  CHECK_FIXNUM(two_ja); CHECK_FIXNUM(two_jb); CHECK_FIXNUM(two_jc);
44
44
  CHECK_FIXNUM(two_jd); CHECK_FIXNUM(two_je); CHECK_FIXNUM(two_jf);
45
45
  return rb_float_new(gsl_sf_coupling_6j(FIX2INT(two_ja), FIX2INT(two_jb),
46
- FIX2INT(two_jc), FIX2INT(two_jd),
47
- FIX2INT(two_je), FIX2INT(two_jf)));
46
+ FIX2INT(two_jc), FIX2INT(two_jd),
47
+ FIX2INT(two_je), FIX2INT(two_jf)));
48
48
  }
49
49
 
50
- static VALUE rb_gsl_sf_coupling_6j_e(VALUE obj, VALUE two_ja, VALUE two_jb, VALUE two_jc, VALUE two_jd, VALUE two_je, VALUE two_jf)
50
+ static VALUE rb_gsl_sf_coupling_6j_e(VALUE obj, VALUE two_ja, VALUE two_jb, VALUE two_jc, VALUE two_jd, VALUE two_je, VALUE two_jf)
51
51
  {
52
52
  gsl_sf_result *rslt = NULL;
53
53
  VALUE v;
@@ -57,31 +57,31 @@ static VALUE rb_gsl_sf_coupling_6j_e(VALUE obj, VALUE two_ja, VALUE two_jb, VALU
57
57
  CHECK_FIXNUM(two_jd); CHECK_FIXNUM(two_je); CHECK_FIXNUM(two_jf);
58
58
  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
59
59
  /*status =*/ gsl_sf_coupling_6j_e(FIX2INT(two_ja), FIX2INT(two_jb),
60
- FIX2INT(two_jc), FIX2INT(two_jd),
61
- FIX2INT(two_je), FIX2INT(two_jf),
62
- rslt);
60
+ FIX2INT(two_jc), FIX2INT(two_jd),
61
+ FIX2INT(two_je), FIX2INT(two_jf),
62
+ rslt);
63
63
  return v;
64
64
  }
65
65
 
66
66
  static VALUE rb_gsl_sf_coupling_9j(VALUE obj, VALUE two_ja, VALUE two_jb,
67
- VALUE two_jc, VALUE two_jd, VALUE two_je,
68
- VALUE two_jf, VALUE two_jg, VALUE two_jh,
69
- VALUE two_ji)
67
+ VALUE two_jc, VALUE two_jd, VALUE two_je,
68
+ VALUE two_jf, VALUE two_jg, VALUE two_jh,
69
+ VALUE two_ji)
70
70
  {
71
71
  CHECK_FIXNUM(two_ja); CHECK_FIXNUM(two_jb); CHECK_FIXNUM(two_jc);
72
72
  CHECK_FIXNUM(two_jd); CHECK_FIXNUM(two_je); CHECK_FIXNUM(two_jf);
73
73
  CHECK_FIXNUM(two_jg); CHECK_FIXNUM(two_jh); CHECK_FIXNUM(two_ji);
74
74
  return rb_float_new(gsl_sf_coupling_9j(FIX2INT(two_ja), FIX2INT(two_jb),
75
- FIX2INT(two_jc), FIX2INT(two_jd),
76
- FIX2INT(two_je), FIX2INT(two_jf),
77
- FIX2INT(two_jg), FIX2INT(two_jh),
78
- FIX2INT(two_ji)));
75
+ FIX2INT(two_jc), FIX2INT(two_jd),
76
+ FIX2INT(two_je), FIX2INT(two_jf),
77
+ FIX2INT(two_jg), FIX2INT(two_jh),
78
+ FIX2INT(two_ji)));
79
79
  }
80
80
 
81
81
  static VALUE rb_gsl_sf_coupling_9j_e(VALUE obj, VALUE two_ja, VALUE two_jb,
82
- VALUE two_jc, VALUE two_jd, VALUE two_je,
83
- VALUE two_jf, VALUE two_jg, VALUE two_jh,
84
- VALUE two_ji)
82
+ VALUE two_jc, VALUE two_jd, VALUE two_je,
83
+ VALUE two_jf, VALUE two_jg, VALUE two_jh,
84
+ VALUE two_ji)
85
85
  {
86
86
  gsl_sf_result *rslt = NULL;
87
87
  VALUE v;
@@ -92,10 +92,10 @@ static VALUE rb_gsl_sf_coupling_9j_e(VALUE obj, VALUE two_ja, VALUE two_jb,
92
92
  CHECK_FIXNUM(two_jg); CHECK_FIXNUM(two_jh); CHECK_FIXNUM(two_ji);
93
93
  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
94
94
  /*status =*/ gsl_sf_coupling_9j_e(FIX2INT(two_ja), FIX2INT(two_jb),
95
- FIX2INT(two_jc), FIX2INT(two_jd),
96
- FIX2INT(two_je), FIX2INT(two_jf),
97
- FIX2INT(two_jg), FIX2INT(two_jh),
98
- FIX2INT(two_ji), rslt);
95
+ FIX2INT(two_jc), FIX2INT(two_jd),
96
+ FIX2INT(two_je), FIX2INT(two_jf),
97
+ FIX2INT(two_jg), FIX2INT(two_jh),
98
+ FIX2INT(two_ji), rslt);
99
99
  return v;
100
100
  }
101
101
 
@@ -111,7 +111,7 @@ void Init_gsl_sf_coupling(VALUE module)
111
111
  rb_define_module_function(module, "coupling_9j_e", rb_gsl_sf_coupling_9j_e, 9);
112
112
 
113
113
  mgsl_sf_coupling = rb_define_module_under(module, "Coupling");
114
-
114
+
115
115
  rb_define_module_function(mgsl_sf_coupling, "3j", rb_gsl_sf_coupling_3j, 6);
116
116
  rb_define_module_function(mgsl_sf_coupling, "3j_e", rb_gsl_sf_coupling_3j_e, 6);
117
117
  rb_define_module_function(mgsl_sf_coupling, "6j", rb_gsl_sf_coupling_6j, 6);
data/ext/{sf_dawson.c → gsl_native/sf_dawson.c} RENAMED
@@ -10,7 +10,7 @@
10
10
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
11
11
  */
12
12
 
13
- #include "rb_gsl_sf.h"
13
+ #include "include/rb_gsl_sf.h"
14
14
 
15
15
  static VALUE rb_gsl_sf_dawson(VALUE obj, VALUE x)
16
16
  {
data/ext/{sf_debye.c → gsl_native/sf_debye.c} RENAMED
@@ -10,7 +10,7 @@
10
10
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
11
11
  */
12
12
 
13
- #include "rb_gsl_sf.h"
13
+ #include "include/rb_gsl_sf.h"
14
14
 
15
15
  static VALUE rb_gsl_sf_debye_1(VALUE obj, VALUE x)
16
16
  {
@@ -83,21 +83,18 @@ static VALUE rb_gsl_sf_debye_n(int argc, VALUE *argv, VALUE obj)
83
83
  case 4:
84
84
  return rb_gsl_sf_eval1(gsl_sf_debye_4, x);
85
85
  break;
86
- #ifdef GSL_1_8_LATER
87
86
  case 5:
88
87
  return rb_gsl_sf_eval1(gsl_sf_debye_5, x);
89
88
  break;
90
89
  case 6:
91
90
  return rb_gsl_sf_eval1(gsl_sf_debye_6, x);
92
91
  break;
93
- #endif
94
92
  default:
95
93
  rb_raise(rb_eRuntimeError, "n must be 1, 2, 3, or 4");
96
94
  break;
97
95
  }
98
96
  }
99
97
 
100
- #ifdef GSL_1_8_LATER
101
98
  static VALUE rb_gsl_sf_debye_5(VALUE obj, VALUE x)
102
99
  {
103
100
  return rb_gsl_sf_eval1(gsl_sf_debye_5, x);
@@ -117,8 +114,6 @@ static VALUE rb_gsl_sf_debye_6_e(VALUE obj, VALUE x)
117
114
  return rb_gsl_sf_eval_e(gsl_sf_debye_6_e, x);
118
115
  }
119
116
 
120
- #endif
121
-
122
117
  void Init_gsl_sf_debye(VALUE module)
123
118
  {
124
119
  VALUE mgsl_sf_debye;
@@ -130,12 +125,10 @@ void Init_gsl_sf_debye(VALUE module)
130
125
  rb_define_module_function(module, "debye_3_e", rb_gsl_sf_debye_3_e, 1);
131
126
  rb_define_module_function(module, "debye_4", rb_gsl_sf_debye_4, 1);
132
127
  rb_define_module_function(module, "debye_4_e", rb_gsl_sf_debye_4_e, 1);
133
- #ifdef GSL_1_8_LATER
134
128
  rb_define_module_function(module, "debye_5", rb_gsl_sf_debye_5, 1);
135
129
  rb_define_module_function(module, "debye_5_e", rb_gsl_sf_debye_5_e, 1);
136
130
  rb_define_module_function(module, "debye_6", rb_gsl_sf_debye_6, 1);
137
131
  rb_define_module_function(module, "debye_6_e", rb_gsl_sf_debye_6_e, 1);
138
- #endif
139
132
  rb_define_module_function(module, "debye_n", rb_gsl_sf_debye_n, -1);
140
133
 
141
134
  mgsl_sf_debye = rb_define_module_under(module, "Debye");
@@ -147,11 +140,9 @@ void Init_gsl_sf_debye(VALUE module)
147
140
  rb_define_module_function(mgsl_sf_debye, "three_e", rb_gsl_sf_debye_3_e, 1);
148
141
  rb_define_module_function(mgsl_sf_debye, "four", rb_gsl_sf_debye_4, 1);
149
142
  rb_define_module_function(mgsl_sf_debye, "four_e", rb_gsl_sf_debye_4_e, 1);
150
- #ifdef GSL_1_8_LATER
151
143
  rb_define_module_function(mgsl_sf_debye, "five", rb_gsl_sf_debye_5, 1);
152
144
  rb_define_module_function(mgsl_sf_debye, "five_e", rb_gsl_sf_debye_5_e, 1);
153
145
  rb_define_module_function(mgsl_sf_debye, "six", rb_gsl_sf_debye_6, 1);
154
146
  rb_define_module_function(mgsl_sf_debye, "six_e", rb_gsl_sf_debye_6_e, 1);
155
- #endif
156
147
  rb_define_module_function(mgsl_sf_debye, "n", rb_gsl_sf_debye_n, -1);
157
148
  }