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
@@ -0,0 +1,122 @@
1
+ #include "include/rb_gsl.h"
2
+ #include "gsl/gsl_bspline.h"
3
+
4
+ static VALUE cBSWS;
5
+
6
+ static VALUE rb_gsl_bspline_alloc(VALUE klass, VALUE k, VALUE n)
7
+ {
8
+ gsl_bspline_workspace *w;
9
+ w = gsl_bspline_alloc(FIX2INT(k), FIX2INT(n));
10
+ return Data_Wrap_Struct(klass, 0, gsl_bspline_free, w);
11
+ }
12
+
13
+ static VALUE rb_gsl_bspline_ncoeffs(VALUE obj)
14
+ {
15
+ gsl_bspline_workspace *w;
16
+ Data_Get_Struct(obj, gsl_bspline_workspace, w);
17
+ return INT2FIX((int)gsl_bspline_ncoeffs(w));
18
+ }
19
+ static VALUE rb_gsl_bspline_order(VALUE obj)
20
+ {
21
+ gsl_bspline_workspace *w;
22
+ Data_Get_Struct(obj, gsl_bspline_workspace, w);
23
+ return INT2FIX((int)gsl_bspline_order(w));
24
+ }
25
+ static VALUE rb_gsl_bspline_nbreak(VALUE obj)
26
+ {
27
+ gsl_bspline_workspace *w;
28
+ Data_Get_Struct(obj, gsl_bspline_workspace, w);
29
+ return INT2FIX((int)gsl_bspline_nbreak(w));
30
+ }
31
+ static VALUE rb_gsl_bspline_breakpoint(VALUE obj, VALUE i)
32
+ {
33
+ gsl_bspline_workspace *w;
34
+ Data_Get_Struct(obj, gsl_bspline_workspace, w);
35
+ return rb_float_new(gsl_bspline_breakpoint(FIX2INT(i), w));
36
+ }
37
+ static VALUE rb_gsl_bspline_knots(VALUE obj, VALUE b)
38
+ {
39
+ gsl_bspline_workspace *w;
40
+ gsl_vector *bpts;
41
+ CHECK_VECTOR(b);
42
+ Data_Get_Struct(obj, gsl_bspline_workspace, w);
43
+ Data_Get_Struct(b, gsl_vector, bpts);
44
+ gsl_bspline_knots(bpts, w);
45
+ return Data_Wrap_Struct(cgsl_vector_view_ro, 0, NULL, w->knots);
46
+ }
47
+ static VALUE rb_gsl_bspline_knots_uniform(int argc, VALUE *argv, VALUE obj)
48
+ {
49
+ gsl_bspline_workspace *w;
50
+ int argc2;
51
+ switch (TYPE(obj)) {
52
+ case T_MODULE:
53
+ case T_CLASS:
54
+ case T_OBJECT:
55
+ if (!rb_obj_is_kind_of(argv[argc-1], cBSWS)) {
56
+ rb_raise(rb_eTypeError, "Wrong argument type %s (GSL::BSpline expected)",
57
+ rb_class2name(CLASS_OF(argv[argc-1])));
58
+ }
59
+ Data_Get_Struct(argv[argc-1], gsl_bspline_workspace, w);
60
+ argc2 = argc-1;
61
+ break;
62
+ default:
63
+ Data_Get_Struct(obj, gsl_bspline_workspace, w);
64
+ argc2 = argc;
65
+ }
66
+ if (argc2 != 2) rb_raise(rb_eArgError, "Wrong number of arguments.");
67
+ gsl_bspline_knots_uniform(NUM2DBL(argv[0]), NUM2DBL(argv[1]), w);
68
+ return Data_Wrap_Struct(cgsl_vector_view_ro, 0, NULL, w->knots);
69
+ }
70
+ static VALUE rb_gsl_bspline_eval(int argc, VALUE *argv, VALUE obj)
71
+ {
72
+ gsl_bspline_workspace *w;
73
+ double x;
74
+ gsl_vector *B;
75
+ VALUE vB;
76
+
77
+ Data_Get_Struct(obj, gsl_bspline_workspace, w);
78
+
79
+ switch (argc) {
80
+ case 2:
81
+ CHECK_VECTOR(argv[1]);
82
+ Data_Get_Struct(argv[1], gsl_vector, B);
83
+ vB = argv[1];
84
+ x = NUM2DBL(argv[0]);
85
+ break;
86
+ case 1:
87
+ x = NUM2DBL(argv[0]);
88
+ B = gsl_vector_alloc(w->nbreak+w->k-2);
89
+ vB = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, B);
90
+ break;
91
+ default:
92
+ rb_raise(rb_eArgError, "Wrong number of arguments (%d for 1 or 2)", argc);
93
+ }
94
+
95
+ gsl_bspline_eval(x, B, w);
96
+
97
+ return vB;
98
+ }
99
+
100
+ static VALUE rb_gsl_bspline_greville_abscissa(VALUE obj, VALUE i)
101
+ {
102
+ gsl_bspline_workspace *w;
103
+ Data_Get_Struct(obj, gsl_bspline_workspace, w);
104
+ return rb_float_new(gsl_bspline_greville_abscissa(i, w));
105
+ }
106
+
107
+ void Init_bspline(VALUE module)
108
+ {
109
+ cBSWS = rb_define_class_under(module, "BSpline", cGSL_Object);
110
+
111
+ rb_define_singleton_method(cBSWS, "alloc", rb_gsl_bspline_alloc, 2);
112
+
113
+ rb_define_method(cBSWS, "ncoeffs", rb_gsl_bspline_ncoeffs, 0);
114
+ rb_define_method(cBSWS, "order", rb_gsl_bspline_order, 0);
115
+ rb_define_method(cBSWS, "nbreak", rb_gsl_bspline_nbreak, 0);
116
+ rb_define_method(cBSWS, "breakpoint", rb_gsl_bspline_breakpoint, 1);
117
+ rb_define_method(cBSWS, "knots", rb_gsl_bspline_knots, 1);
118
+ rb_define_method(cBSWS, "knots_uniform", rb_gsl_bspline_knots_uniform, -1);
119
+ rb_define_singleton_method(cBSWS, "knots_uniform", rb_gsl_bspline_knots_uniform, -1);
120
+ rb_define_method(cBSWS, "eval", rb_gsl_bspline_eval, -1);
121
+ rb_define_method(cBSWS, "greville_abscissa", rb_gsl_bspline_greville_abscissa, 1);
122
+ }
File without changes
@@ -9,21 +9,18 @@
9
9
  WITHOUT ANY WARRANTY.
10
10
  */
11
11
 
12
- #include "rb_gsl_config.h"
13
-
14
- #ifdef GSL_1_4_LATER
15
- #include "rb_gsl_array.h"
16
- #include "rb_gsl_common.h"
17
- #include "rb_gsl_rng.h"
12
+ #include "include/rb_gsl_array.h"
13
+ #include "include/rb_gsl_common.h"
14
+ #include "include/rb_gsl_rng.h"
18
15
  #include <gsl/gsl_randist.h>
19
16
  #include <gsl/gsl_cdf.h>
20
17
 
21
18
  VALUE rb_gsl_eval_pdf_cdf(VALUE xx, double (*f)(double));
22
19
  VALUE rb_gsl_eval_pdf_cdf2(VALUE xx, VALUE aa, double (*f)(double, double));
23
- VALUE rb_gsl_eval_pdf_cdf3(VALUE xx, VALUE aa, VALUE bb,
24
- double (*f)(double, double, double));
25
- VALUE rb_gsl_eval_pdf_cdf2_uint(VALUE xx, VALUE aa,
26
- double (*f)(unsigned int, double));
20
+ VALUE rb_gsl_eval_pdf_cdf3(VALUE xx, VALUE aa, VALUE bb,
21
+ double (*f)(double, double, double));
22
+ VALUE rb_gsl_eval_pdf_cdf2_uint(VALUE xx, VALUE aa,
23
+ double (*f)(unsigned int, double));
27
24
 
28
25
 
29
26
  static VALUE rb_gsl_cdf_gaussian_P(int argc, VALUE *argv, VALUE obj)
@@ -72,316 +69,313 @@ static VALUE rb_gsl_cdf_gaussian_Qinv(int argc, VALUE *argv, VALUE obj)
72
69
  static VALUE rb_gsl_cdf_exponential_P(VALUE obj, VALUE x, VALUE mu)
73
70
  {
74
71
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_exponential_P);
75
- }
72
+ }
76
73
 
77
74
  static VALUE rb_gsl_cdf_exponential_Q(VALUE obj, VALUE x, VALUE mu)
78
75
  {
79
76
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_exponential_Q);
80
- }
77
+ }
81
78
 
82
79
  static VALUE rb_gsl_cdf_exponential_Pinv(VALUE obj, VALUE x, VALUE mu)
83
80
  {
84
81
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_exponential_Pinv);
85
- }
82
+ }
86
83
 
87
84
  static VALUE rb_gsl_cdf_exponential_Qinv(VALUE obj, VALUE x, VALUE mu)
88
85
  {
89
86
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_exponential_Qinv);
90
- }
87
+ }
91
88
 
92
89
  static VALUE rb_gsl_cdf_laplace_P(VALUE obj, VALUE x, VALUE mu)
93
90
  {
94
91
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_laplace_P);
95
- }
92
+ }
96
93
 
97
94
  static VALUE rb_gsl_cdf_laplace_Q(VALUE obj, VALUE x, VALUE mu)
98
95
  {
99
96
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_laplace_Q);
100
- }
97
+ }
101
98
 
102
99
  static VALUE rb_gsl_cdf_laplace_Pinv(VALUE obj, VALUE x, VALUE mu)
103
100
  {
104
101
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_laplace_Pinv);
105
- }
102
+ }
106
103
 
107
104
  static VALUE rb_gsl_cdf_laplace_Qinv(VALUE obj, VALUE x, VALUE mu)
108
105
  {
109
106
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_laplace_Qinv);
110
- }
107
+ }
111
108
 
112
109
  static VALUE rb_gsl_cdf_cauchy_P(VALUE obj, VALUE x, VALUE mu)
113
110
  {
114
111
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_cauchy_P);
115
- }
112
+ }
116
113
 
117
114
  static VALUE rb_gsl_cdf_cauchy_Q(VALUE obj, VALUE x, VALUE mu)
118
115
  {
119
116
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_cauchy_Q);
120
- }
117
+ }
121
118
 
122
119
  static VALUE rb_gsl_cdf_cauchy_Pinv(VALUE obj, VALUE x, VALUE mu)
123
120
  {
124
121
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_cauchy_Pinv);
125
- }
122
+ }
126
123
 
127
124
  static VALUE rb_gsl_cdf_cauchy_Qinv(VALUE obj, VALUE x, VALUE mu)
128
125
  {
129
126
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_cauchy_Qinv);
130
- }
127
+ }
131
128
 
132
129
  static VALUE rb_gsl_cdf_rayleigh_P(VALUE obj, VALUE x, VALUE mu)
133
130
  {
134
131
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_rayleigh_P);
135
- }
132
+ }
136
133
 
137
134
  static VALUE rb_gsl_cdf_rayleigh_Q(VALUE obj, VALUE x, VALUE mu)
138
135
  {
139
136
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_rayleigh_Q);
140
- }
137
+ }
141
138
 
142
139
  static VALUE rb_gsl_cdf_rayleigh_Pinv(VALUE obj, VALUE x, VALUE mu)
143
140
  {
144
141
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_rayleigh_Pinv);
145
- }
142
+ }
146
143
 
147
144
  static VALUE rb_gsl_cdf_rayleigh_Qinv(VALUE obj, VALUE x, VALUE mu)
148
145
  {
149
146
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_rayleigh_Qinv);
150
- }
147
+ }
151
148
 
152
149
  static VALUE rb_gsl_cdf_gamma_P(VALUE obj, VALUE x, VALUE a, VALUE b)
153
150
  {
154
151
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gamma_P);
155
- }
152
+ }
156
153
 
157
154
  static VALUE rb_gsl_cdf_gamma_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
158
155
  {
159
156
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gamma_Q);
160
- }
157
+ }
161
158
 
162
159
  static VALUE rb_gsl_cdf_gamma_Pinv(VALUE obj, VALUE x, VALUE a, VALUE b)
163
160
  {
164
161
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gamma_Pinv);
165
- }
162
+ }
166
163
 
167
164
  static VALUE rb_gsl_cdf_gamma_Qinv(VALUE obj, VALUE x, VALUE a, VALUE b)
168
165
  {
169
166
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gamma_Qinv);
170
- }
167
+ }
171
168
 
172
169
  static VALUE rb_gsl_cdf_flat_P(VALUE obj, VALUE x, VALUE a, VALUE b)
173
170
  {
174
171
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_flat_P);
175
- }
172
+ }
176
173
 
177
174
  static VALUE rb_gsl_cdf_flat_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
178
175
  {
179
176
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_flat_Q);
180
- }
177
+ }
181
178
 
182
179
  static VALUE rb_gsl_cdf_flat_Pinv(VALUE obj, VALUE x, VALUE a, VALUE b)
183
180
  {
184
181
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_flat_Pinv);
185
- }
182
+ }
186
183
 
187
184
  static VALUE rb_gsl_cdf_flat_Qinv(VALUE obj, VALUE x, VALUE a, VALUE b)
188
185
  {
189
186
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_flat_Qinv);
190
- }
187
+ }
191
188
 
192
189
  static VALUE rb_gsl_cdf_lognormal_P(VALUE obj, VALUE x, VALUE a, VALUE b)
193
190
  {
194
191
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_lognormal_P);
195
- }
192
+ }
196
193
 
197
194
  static VALUE rb_gsl_cdf_lognormal_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
198
195
  {
199
196
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_lognormal_Q);
200
- }
197
+ }
201
198
 
202
199
  static VALUE rb_gsl_cdf_lognormal_Pinv(VALUE obj, VALUE x, VALUE a, VALUE b)
203
200
  {
204
201
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_lognormal_Pinv);
205
- }
202
+ }
206
203
 
207
204
  static VALUE rb_gsl_cdf_lognormal_Qinv(VALUE obj, VALUE x, VALUE a, VALUE b)
208
205
  {
209
206
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_lognormal_Qinv);
210
- }
207
+ }
211
208
 
212
- #ifdef GSL_1_6_LATER
213
209
  static VALUE rb_gsl_cdf_exppow_P(VALUE obj, VALUE x, VALUE a, VALUE b)
214
210
  {
215
211
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_exppow_P);
216
- }
212
+ }
217
213
 
218
214
  static VALUE rb_gsl_cdf_exppow_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
219
215
  {
220
216
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_exppow_Q);
221
- }
222
- #endif
217
+ }
223
218
 
224
219
  static VALUE rb_gsl_cdf_chisq_P(VALUE obj, VALUE x, VALUE mu)
225
220
  {
226
221
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_chisq_P);
227
- }
222
+ }
228
223
 
229
224
  static VALUE rb_gsl_cdf_chisq_Q(VALUE obj, VALUE x, VALUE mu)
230
225
  {
231
226
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_chisq_Q);
232
- }
227
+ }
233
228
 
234
229
  static VALUE rb_gsl_cdf_chisq_Pinv(VALUE obj, VALUE x, VALUE mu)
235
230
  {
236
231
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_chisq_Pinv);
237
- }
232
+ }
238
233
 
239
234
  static VALUE rb_gsl_cdf_chisq_Qinv(VALUE obj, VALUE x, VALUE mu)
240
235
  {
241
236
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_chisq_Qinv);
242
- }
237
+ }
243
238
 
244
239
  static VALUE rb_gsl_cdf_fdist_P(VALUE obj, VALUE x, VALUE nu1, VALUE nu2)
245
240
  {
246
241
  return rb_gsl_eval_pdf_cdf3(x, nu1, nu2, gsl_cdf_fdist_P);
247
- }
242
+ }
248
243
 
249
244
  static VALUE rb_gsl_cdf_fdist_Q(VALUE obj, VALUE x, VALUE nu1, VALUE nu2)
250
245
  {
251
246
  return rb_gsl_eval_pdf_cdf3(x, nu1, nu2, gsl_cdf_fdist_Q);
252
- }
247
+ }
253
248
 
254
249
  static VALUE rb_gsl_cdf_tdist_P(VALUE obj, VALUE x, VALUE mu)
255
250
  {
256
251
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_tdist_P);
257
- }
252
+ }
258
253
 
259
254
  static VALUE rb_gsl_cdf_tdist_Q(VALUE obj, VALUE x, VALUE mu)
260
255
  {
261
256
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_tdist_Q);
262
- }
257
+ }
263
258
 
264
259
  static VALUE rb_gsl_cdf_tdist_Pinv(VALUE obj, VALUE x, VALUE mu)
265
260
  {
266
261
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_tdist_Pinv);
267
- }
262
+ }
268
263
 
269
264
  static VALUE rb_gsl_cdf_tdist_Qinv(VALUE obj, VALUE x, VALUE mu)
270
265
  {
271
266
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_tdist_Qinv);
272
- }
267
+ }
273
268
 
274
269
  static VALUE rb_gsl_cdf_beta_P(VALUE obj, VALUE x, VALUE nu1, VALUE nu2)
275
270
  {
276
271
  return rb_gsl_eval_pdf_cdf3(x, nu1, nu2, gsl_cdf_beta_P);
277
- }
272
+ }
278
273
 
279
274
  static VALUE rb_gsl_cdf_beta_Q(VALUE obj, VALUE x, VALUE nu1, VALUE nu2)
280
275
  {
281
276
  return rb_gsl_eval_pdf_cdf3(x, nu1, nu2, gsl_cdf_beta_Q);
282
- }
277
+ }
283
278
 
284
279
  static VALUE rb_gsl_cdf_logistic_P(VALUE obj, VALUE x, VALUE mu)
285
280
  {
286
281
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_logistic_P);
287
- }
282
+ }
288
283
 
289
284
  static VALUE rb_gsl_cdf_logistic_Q(VALUE obj, VALUE x, VALUE mu)
290
285
  {
291
286
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_logistic_Q);
292
- }
287
+ }
293
288
 
294
289
  static VALUE rb_gsl_cdf_logistic_Pinv(VALUE obj, VALUE x, VALUE mu)
295
290
  {
296
291
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_logistic_Pinv);
297
- }
292
+ }
298
293
 
299
294
  static VALUE rb_gsl_cdf_logistic_Qinv(VALUE obj, VALUE x, VALUE mu)
300
295
  {
301
296
  return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_logistic_Qinv);
302
- }
297
+ }
303
298
 
304
299
  static VALUE rb_gsl_cdf_pareto_P(VALUE obj, VALUE x, VALUE a, VALUE b)
305
300
  {
306
301
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_pareto_P);
307
- }
302
+ }
308
303
 
309
304
  static VALUE rb_gsl_cdf_pareto_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
310
305
  {
311
306
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_pareto_Q);
312
- }
307
+ }
313
308
 
314
309
  static VALUE rb_gsl_cdf_pareto_Pinv(VALUE obj, VALUE x, VALUE a, VALUE b)
315
310
  {
316
311
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_pareto_Pinv);
317
- }
312
+ }
318
313
 
319
314
  static VALUE rb_gsl_cdf_pareto_Qinv(VALUE obj, VALUE x, VALUE a, VALUE b)
320
315
  {
321
316
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_pareto_Qinv);
322
- }
317
+ }
323
318
 
324
319
  static VALUE rb_gsl_cdf_weibull_P(VALUE obj, VALUE x, VALUE a, VALUE b)
325
320
  {
326
321
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_weibull_P);
327
- }
322
+ }
328
323
 
329
324
  static VALUE rb_gsl_cdf_weibull_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
330
325
  {
331
326
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_weibull_Q);
332
- }
327
+ }
333
328
 
334
329
  static VALUE rb_gsl_cdf_weibull_Pinv(VALUE obj, VALUE x, VALUE a, VALUE b)
335
330
  {
336
331
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_weibull_Pinv);
337
- }
332
+ }
338
333
 
339
334
  static VALUE rb_gsl_cdf_weibull_Qinv(VALUE obj, VALUE x, VALUE a, VALUE b)
340
335
  {
341
336
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_weibull_Qinv);
342
- }
337
+ }
343
338
 
344
339
  static VALUE rb_gsl_cdf_gumbel1_P(VALUE obj, VALUE x, VALUE a, VALUE b)
345
340
  {
346
341
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel1_P);
347
- }
342
+ }
348
343
 
349
344
  static VALUE rb_gsl_cdf_gumbel1_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
350
345
  {
351
346
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel1_Q);
352
- }
347
+ }
353
348
 
354
349
  static VALUE rb_gsl_cdf_gumbel1_Pinv(VALUE obj, VALUE x, VALUE a, VALUE b)
355
350
  {
356
351
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel1_Pinv);
357
- }
352
+ }
358
353
 
359
354
  static VALUE rb_gsl_cdf_gumbel1_Qinv(VALUE obj, VALUE x, VALUE a, VALUE b)
360
355
  {
361
356
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel1_Qinv);
362
- }
357
+ }
363
358
 
364
359
  static VALUE rb_gsl_cdf_gumbel2_P(VALUE obj, VALUE x, VALUE a, VALUE b)
365
360
  {
366
361
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel2_P);
367
- }
362
+ }
368
363
 
369
364
  static VALUE rb_gsl_cdf_gumbel2_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
370
365
  {
371
366
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel2_Q);
372
- }
367
+ }
373
368
 
374
369
  static VALUE rb_gsl_cdf_gumbel2_Pinv(VALUE obj, VALUE x, VALUE a, VALUE b)
375
370
  {
376
371
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel2_Pinv);
377
- }
372
+ }
378
373
 
379
374
  static VALUE rb_gsl_cdf_gumbel2_Qinv(VALUE obj, VALUE x, VALUE a, VALUE b)
380
375
  {
381
376
  return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel2_Qinv);
382
- }
377
+ }
383
378
 
384
- #ifdef GSL_1_8_LATER
385
379
  static VALUE rb_gsl_cdf_binomial_P(VALUE obj, VALUE kk, VALUE pp, VALUE nn)
386
380
  {
387
381
  unsigned int k, n;
@@ -518,7 +512,6 @@ static VALUE rb_gsl_cdf_fdist_Qinv(VALUE obj, VALUE qq, VALUE aa, VALUE bb)
518
512
  b = NUM2DBL(bb);
519
513
  return rb_float_new(gsl_cdf_fdist_Qinv(Q, a, b));
520
514
  }
521
- #endif
522
515
 
523
516
  void Init_gsl_cdf(VALUE module)
524
517
  {
@@ -537,7 +530,7 @@ void Init_gsl_cdf(VALUE module)
537
530
  rb_define_module_function(module, "cdf_ugaussian_Qinv", rb_gsl_cdf_gaussian_Qinv, -1);
538
531
 
539
532
  rb_define_module_function(module, "cdf_exponential_P", rb_gsl_cdf_exponential_P, 2);
540
- rb_define_module_function(module, "cdf_exponential_Q", rb_gsl_cdf_exponential_Q, 2);
533
+ rb_define_module_function(module, "cdf_exponential_Q", rb_gsl_cdf_exponential_Q, 2);
541
534
  rb_define_module_function(module, "cdf_exponential_Pinv", rb_gsl_cdf_exponential_Pinv, 2);
542
535
  rb_define_module_function(module, "cdf_exponential_Qinv", rb_gsl_cdf_exponential_Qinv, 2);
543
536
 
@@ -575,7 +568,7 @@ void Init_gsl_cdf(VALUE module)
575
568
  rb_define_module_function(module, "cdf_chisq_Q", rb_gsl_cdf_chisq_Q, 2);
576
569
  rb_define_module_function(module, "cdf_chisq_Pinv", rb_gsl_cdf_chisq_Pinv, 2);
577
570
  rb_define_module_function(module, "cdf_chisq_Qinv", rb_gsl_cdf_chisq_Qinv, 2);
578
-
571
+
579
572
  rb_define_module_function(module, "cdf_fdist_P", rb_gsl_cdf_fdist_P, 3);
580
573
  rb_define_module_function(module, "cdf_fdist_Q", rb_gsl_cdf_fdist_Q, 3);
581
574
 
@@ -583,7 +576,7 @@ void Init_gsl_cdf(VALUE module)
583
576
  rb_define_module_function(module, "cdf_tdist_Q", rb_gsl_cdf_tdist_Q, 2);
584
577
  rb_define_module_function(module, "cdf_tdist_Pinv", rb_gsl_cdf_tdist_Pinv, 2);
585
578
  rb_define_module_function(module, "cdf_tdist_Qinv", rb_gsl_cdf_tdist_Qinv, 2);
586
-
579
+
587
580
  rb_define_module_function(module, "cdf_beta_P", rb_gsl_cdf_beta_P, 3);
588
581
  rb_define_module_function(module, "cdf_beta_Q", rb_gsl_cdf_beta_Q, 3);
589
582
 
@@ -591,7 +584,7 @@ void Init_gsl_cdf(VALUE module)
591
584
  rb_define_module_function(module, "cdf_logistic_Q", rb_gsl_cdf_logistic_Q, 2);
592
585
  rb_define_module_function(module, "cdf_logistic_Pinv", rb_gsl_cdf_logistic_Pinv, 2);
593
586
  rb_define_module_function(module, "cdf_logistic_Qinv", rb_gsl_cdf_logistic_Qinv, 2);
594
-
587
+
595
588
  rb_define_module_function(module, "cdf_pareto_P", rb_gsl_cdf_pareto_P, 3);
596
589
  rb_define_module_function(module, "cdf_pareto_Q", rb_gsl_cdf_pareto_Q, 3);
597
590
  rb_define_module_function(module, "cdf_pareto_Pinv", rb_gsl_cdf_pareto_Pinv, 3);
@@ -624,12 +617,12 @@ void Init_gsl_cdf(VALUE module)
624
617
  rb_define_module_function(mgsl_cdf, "ugaussian_Qinv", rb_gsl_cdf_gaussian_Qinv, -1);
625
618
 
626
619
  rb_define_module_function(mgsl_cdf, "exponential_P", rb_gsl_cdf_exponential_P, 2);
627
- rb_define_module_function(mgsl_cdf, "exponential_Q", rb_gsl_cdf_exponential_Q, 2);
620
+ rb_define_module_function(mgsl_cdf, "exponential_Q", rb_gsl_cdf_exponential_Q, 2);
628
621
  rb_define_module_function(mgsl_cdf, "exponential_Pinv", rb_gsl_cdf_exponential_Pinv, 2);
629
622
  rb_define_module_function(mgsl_cdf, "exponential_Qinv", rb_gsl_cdf_exponential_Qinv, 2);
630
623
 
631
- rb_define_module_function(mgsl_cdf, "laplace_P", rb_gsl_cdf_laplace_P, 2);
632
- rb_define_module_function(mgsl_cdf, "laplace_Q", rb_gsl_cdf_laplace_Q, 2);
624
+ rb_define_module_function(mgsl_cdf, "laplace_P", rb_gsl_cdf_laplace_P, 2);
625
+ rb_define_module_function(mgsl_cdf, "laplace_Q", rb_gsl_cdf_laplace_Q, 2);
633
626
  rb_define_module_function(mgsl_cdf, "laplace_Pinv", rb_gsl_cdf_laplace_Pinv, 2);
634
627
  rb_define_module_function(mgsl_cdf, "laplace_Qinv", rb_gsl_cdf_laplace_Qinv, 2);
635
628
 
@@ -662,7 +655,7 @@ void Init_gsl_cdf(VALUE module)
662
655
  rb_define_module_function(mgsl_cdf, "chisq_Q", rb_gsl_cdf_chisq_Q, 2);
663
656
  rb_define_module_function(mgsl_cdf, "chisq_Pinv", rb_gsl_cdf_chisq_Pinv, 2);
664
657
  rb_define_module_function(mgsl_cdf, "chisq_Qinv", rb_gsl_cdf_chisq_Qinv, 2);
665
-
658
+
666
659
  rb_define_module_function(mgsl_cdf, "fdist_P", rb_gsl_cdf_fdist_P, 3);
667
660
  rb_define_module_function(mgsl_cdf, "fdist_Q", rb_gsl_cdf_fdist_Q, 3);
668
661
 
@@ -670,7 +663,7 @@ void Init_gsl_cdf(VALUE module)
670
663
  rb_define_module_function(mgsl_cdf, "tdist_Q", rb_gsl_cdf_tdist_Q, 2);
671
664
  rb_define_module_function(mgsl_cdf, "tdist_Pinv", rb_gsl_cdf_tdist_Pinv, 2);
672
665
  rb_define_module_function(mgsl_cdf, "tdist_Qinv", rb_gsl_cdf_tdist_Qinv, 2);
673
-
666
+
674
667
  rb_define_module_function(mgsl_cdf, "beta_P", rb_gsl_cdf_beta_P, 3);
675
668
  rb_define_module_function(mgsl_cdf, "beta_Q", rb_gsl_cdf_beta_Q, 3);
676
669
 
@@ -678,7 +671,7 @@ void Init_gsl_cdf(VALUE module)
678
671
  rb_define_module_function(mgsl_cdf, "logistic_Q", rb_gsl_cdf_logistic_Q, 2);
679
672
  rb_define_module_function(mgsl_cdf, "logistic_Pinv", rb_gsl_cdf_logistic_Pinv, 2);
680
673
  rb_define_module_function(mgsl_cdf, "logistic_Qinv", rb_gsl_cdf_logistic_Qinv, 2);
681
-
674
+
682
675
  rb_define_module_function(mgsl_cdf, "pareto_P", rb_gsl_cdf_pareto_P, 3);
683
676
  rb_define_module_function(mgsl_cdf, "pareto_Q", rb_gsl_cdf_pareto_Q, 3);
684
677
  rb_define_module_function(mgsl_cdf, "pareto_Pinv", rb_gsl_cdf_pareto_Pinv, 3);
@@ -699,16 +692,12 @@ void Init_gsl_cdf(VALUE module)
699
692
  rb_define_module_function(mgsl_cdf, "gumbel2_Pinv", rb_gsl_cdf_gumbel2_Pinv, 3);
700
693
  rb_define_module_function(mgsl_cdf, "gumbel2_Qinv", rb_gsl_cdf_gumbel2_Qinv, 3);
701
694
 
702
-
703
- #ifdef GSL_1_6_LATER
704
695
  rb_define_module_function(module, "cdf_exppow_P", rb_gsl_cdf_exppow_P, 3);
705
696
  rb_define_module_function(module, "cdf_exppow_Q", rb_gsl_cdf_exppow_Q, 3);
706
697
 
707
698
  rb_define_module_function(mgsl_cdf, "exppow_P", rb_gsl_cdf_exppow_P, 3);
708
699
  rb_define_module_function(mgsl_cdf, "exppow_Q", rb_gsl_cdf_exppow_Q, 3);
709
- #endif
710
700
 
711
- #ifdef GSL_1_8_LATER
712
701
  rb_define_module_function(module, "cdf_binomial_P", rb_gsl_cdf_binomial_P, 3);
713
702
  rb_define_module_function(module, "cdf_binomial_Q", rb_gsl_cdf_binomial_Q, 3);
714
703
  rb_define_module_function(mgsl_cdf, "binomial_P", rb_gsl_cdf_binomial_P, 3);
@@ -748,7 +737,4 @@ void Init_gsl_cdf(VALUE module)
748
737
  rb_define_module_function(mgsl_cdf, "fdist_Pinv", rb_gsl_cdf_fdist_Pinv, 3);
749
738
  rb_define_module_function(module, "cdf_fdist_Qinv", rb_gsl_cdf_fdist_Qinv, 3);
750
739
  rb_define_module_function(mgsl_cdf, "fdist_Qinv", rb_gsl_cdf_fdist_Qinv, 3);
751
- #endif
752
-
753
740
  }
754
- #endif