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.
- checksums.yaml +7 -0
- data/.gitignore +12 -0
- data/.travis.yml +24 -0
- data/AUTHORS +10 -2
- data/COPYING +341 -339
- data/ChangeLog +612 -554
- data/Gemfile +4 -0
- data/README.md +77 -0
- data/Rakefile +14 -93
- data/THANKS +1 -1
- data/examples/blas/blas.rb +1 -1
- data/examples/bspline.rb +3 -3
- data/examples/complex/functions.rb +4 -4
- data/examples/complex/mul.rb +1 -1
- data/examples/const/physconst.rb +4 -4
- data/examples/const/travel.rb +4 -4
- data/examples/eigen/nonsymmv.rb +1 -1
- data/examples/eigen/qhoscillator.rb +3 -3
- data/examples/fft/radix2.rb +1 -1
- data/examples/fft/real-halfcomplex.rb +3 -3
- data/examples/fft/real-halfcomplex2.rb +3 -3
- data/examples/fit/expfit.rb +1 -1
- data/examples/fit/multifit.rb +1 -1
- data/examples/fit/ndlinear.rb +44 -44
- data/examples/fit/nonlinearfit.rb +1 -1
- data/examples/fit/wlinear.rb +3 -3
- data/examples/function/function.rb +1 -1
- data/examples/function/min.rb +1 -1
- data/examples/function/synchrotron.rb +2 -2
- data/examples/gallery/koch.rb +1 -1
- data/examples/histogram/cauchy.rb +2 -2
- data/examples/histogram/exponential.rb +1 -1
- data/examples/histogram/histo3d.rb +1 -1
- data/examples/histogram/histogram-pdf.rb +2 -2
- data/examples/histogram/xexp.rb +1 -1
- data/examples/integration/ahmed.rb +2 -2
- data/examples/integration/cosmology.rb +7 -7
- data/examples/integration/friedmann.rb +4 -4
- data/examples/integration/qng.rb +1 -1
- data/examples/interp/demo.rb +2 -2
- data/examples/linalg/LQ_solve.rb +1 -1
- data/examples/linalg/LU.rb +1 -1
- data/examples/linalg/LU2.rb +1 -1
- data/examples/linalg/LU_narray.rb +1 -1
- data/examples/linalg/PTLQ.rb +1 -1
- data/examples/linalg/QRPT.rb +1 -1
- data/examples/linalg/chol.rb +1 -1
- data/examples/linalg/chol_narray.rb +1 -1
- data/examples/linalg/complex.rb +1 -1
- data/examples/math/elementary.rb +1 -1
- data/examples/math/functions.rb +1 -1
- data/examples/math/inf_nan.rb +1 -1
- data/examples/math/minmax.rb +1 -1
- data/examples/math/power.rb +1 -1
- data/examples/math/test.rb +1 -1
- data/examples/min.rb +1 -1
- data/examples/multimin/bundle.rb +1 -1
- data/examples/multimin/cqp.rb +17 -17
- data/examples/multiroot/fsolver3.rb +1 -1
- data/examples/odeiv/binarysystem.rb +12 -12
- data/examples/odeiv/demo.rb +3 -3
- data/examples/odeiv/frei1.rb +7 -7
- data/examples/odeiv/frei2.rb +4 -4
- data/examples/odeiv/oscillator.rb +1 -1
- data/examples/odeiv/sedov.rb +3 -3
- data/examples/odeiv/whitedwarf.rb +11 -11
- data/examples/permutation/ex1.rb +2 -2
- data/examples/permutation/permutation.rb +1 -1
- data/examples/poly/demo.rb +1 -1
- data/examples/random/diffusion.rb +1 -1
- data/examples/random/generator.rb +2 -2
- data/examples/random/randomwalk.rb +3 -3
- data/examples/random/rng.rb +1 -1
- data/examples/roots/bisection.rb +1 -1
- data/examples/roots/brent.rb +1 -1
- data/examples/roots/demo.rb +1 -1
- data/examples/roots/newton.rb +2 -2
- data/examples/roots/recombination.gp +0 -1
- data/examples/sf/hyperg.rb +1 -1
- data/examples/sf/sphbessel.rb +1 -1
- data/examples/sort/sort.rb +1 -1
- data/examples/tamu_anova.rb +4 -4
- data/examples/vector/add.rb +1 -1
- data/examples/vector/decimate.rb +1 -1
- data/examples/vector/gnuplot.rb +8 -8
- data/examples/vector/vector.rb +2 -2
- data/examples/wavelet/wavelet1.rb +1 -1
- data/ext/{alf.c → gsl_native/alf.c} +10 -10
- data/ext/{array.c → gsl_native/array.c} +70 -159
- data/ext/{array_complex.c → gsl_native/array_complex.c} +63 -66
- data/ext/{blas.c → gsl_native/blas.c} +2 -3
- data/ext/{blas1.c → gsl_native/blas1.c} +35 -36
- data/ext/{blas2.c → gsl_native/blas2.c} +57 -62
- data/ext/{blas3.c → gsl_native/blas3.c} +57 -58
- data/ext/{block.c → gsl_native/block.c} +14 -18
- data/ext/{block_source.c → gsl_native/block_source.h} +110 -112
- data/ext/gsl_native/bspline.c +122 -0
- data/ext/{bundle.c → gsl_native/bundle.c} +0 -0
- data/ext/{cdf.c → gsl_native/cdf.c} +79 -93
- data/ext/{cheb.c → gsl_native/cheb.c} +78 -89
- data/ext/{combination.c → gsl_native/combination.c} +11 -19
- data/ext/{common.c → gsl_native/common.c} +9 -41
- data/ext/{complex.c → gsl_native/complex.c} +116 -118
- data/ext/gsl_native/const.c +331 -0
- data/ext/{const_additional.c → gsl_native/const_additional.c} +13 -34
- data/ext/gsl_native/cqp.c +283 -0
- data/ext/{deriv.c → gsl_native/deriv.c} +25 -33
- data/ext/{dht.c → gsl_native/dht.c} +23 -31
- data/ext/{diff.c → gsl_native/diff.c} +26 -28
- data/ext/{dirac.c → gsl_native/dirac.c} +45 -46
- data/ext/{eigen.c → gsl_native/eigen.c} +1044 -1095
- data/ext/{error.c → gsl_native/error.c} +18 -18
- data/ext/gsl_native/extconf.rb +118 -0
- data/ext/{fft.c → gsl_native/fft.c} +197 -204
- data/ext/{fit.c → gsl_native/fit.c} +17 -18
- data/ext/gsl_native/fresnel.c +312 -0
- data/ext/{function.c → gsl_native/function.c} +37 -43
- data/ext/{geometry.c → gsl_native/geometry.c} +16 -16
- data/ext/{graph.c → gsl_native/graph.c} +39 -89
- data/ext/{gsl.c → gsl_native/gsl.c} +12 -33
- data/ext/{gsl_narray.c → gsl_native/gsl_narray.c} +20 -30
- data/ext/{histogram.c → gsl_native/histogram.c} +133 -160
- data/ext/{histogram2d.c → gsl_native/histogram2d.c} +78 -104
- data/ext/{histogram3d.c → gsl_native/histogram3d.c} +76 -76
- data/ext/{histogram3d_source.c → gsl_native/histogram3d_source.c} +196 -197
- data/ext/{histogram_find.c → gsl_native/histogram_find.c} +32 -34
- data/ext/{histogram_oper.c → gsl_native/histogram_oper.c} +43 -52
- data/ext/{ieee.c → gsl_native/ieee.c} +9 -21
- data/{include → ext/gsl_native/include}/rb_gsl.h +4 -26
- data/{include → ext/gsl_native/include}/rb_gsl_array.h +15 -39
- data/{include → ext/gsl_native/include}/rb_gsl_cheb.h +0 -2
- data/{include → ext/gsl_native/include}/rb_gsl_common.h +61 -61
- data/{include → ext/gsl_native/include}/rb_gsl_complex.h +1 -1
- data/{include → ext/gsl_native/include}/rb_gsl_const.h +0 -6
- data/ext/gsl_native/include/rb_gsl_dirac.h +6 -0
- data/{include → ext/gsl_native/include}/rb_gsl_eigen.h +1 -1
- data/{include → ext/gsl_native/include}/rb_gsl_fft.h +0 -13
- data/{include → ext/gsl_native/include}/rb_gsl_fit.h +0 -2
- data/{include → ext/gsl_native/include}/rb_gsl_function.h +0 -4
- data/{include → ext/gsl_native/include}/rb_gsl_graph.h +2 -4
- data/{include → ext/gsl_native/include}/rb_gsl_histogram.h +8 -8
- data/{include → ext/gsl_native/include}/rb_gsl_histogram3d.h +50 -50
- data/{include → ext/gsl_native/include}/rb_gsl_integration.h +1 -1
- data/{include → ext/gsl_native/include}/rb_gsl_interp.h +0 -5
- data/{include → ext/gsl_native/include}/rb_gsl_linalg.h +2 -6
- data/{include → ext/gsl_native/include}/rb_gsl_math.h +0 -6
- data/{include → ext/gsl_native/include}/rb_gsl_odeiv.h +0 -3
- data/{include → ext/gsl_native/include}/rb_gsl_poly.h +3 -7
- data/{include → ext/gsl_native/include}/rb_gsl_rational.h +1 -8
- data/{include → ext/gsl_native/include}/rb_gsl_rng.h +0 -1
- data/{include → ext/gsl_native/include}/rb_gsl_root.h +1 -1
- data/{include → ext/gsl_native/include}/rb_gsl_sf.h +39 -48
- data/{include → ext/gsl_native/include}/rb_gsl_statistics.h +1 -1
- data/{include → ext/gsl_native/include}/rb_gsl_tensor.h +0 -2
- data/{include → ext/gsl_native/include}/rb_gsl_with_narray.h +3 -1
- data/{include → ext/gsl_native/include}/templates_off.h +0 -0
- data/{include → ext/gsl_native/include}/templates_on.h +1 -1
- data/ext/{integration.c → gsl_native/integration.c} +164 -189
- data/ext/{interp.c → gsl_native/interp.c} +25 -38
- data/ext/gsl_native/jacobi.c +733 -0
- data/ext/{linalg.c → gsl_native/linalg.c} +462 -589
- data/ext/{linalg_complex.c → gsl_native/linalg_complex.c} +93 -106
- data/ext/{math.c → gsl_native/math.c} +48 -67
- data/ext/{matrix.c → gsl_native/matrix.c} +13 -16
- data/ext/{matrix_complex.c → gsl_native/matrix_complex.c} +119 -123
- data/ext/{matrix_double.c → gsl_native/matrix_double.c} +79 -82
- data/ext/{matrix_int.c → gsl_native/matrix_int.c} +53 -54
- data/ext/{matrix_source.c → gsl_native/matrix_source.h} +292 -318
- data/ext/{min.c → gsl_native/min.c} +45 -76
- data/ext/{monte.c → gsl_native/monte.c} +50 -64
- data/ext/{multifit.c → gsl_native/multifit.c} +142 -151
- data/ext/{multimin.c → gsl_native/multimin.c} +64 -92
- data/ext/{multimin_fsdf.c → gsl_native/multimin_fsdf.c} +16 -16
- data/ext/{multiroots.c → gsl_native/multiroots.c} +73 -76
- data/ext/{multiset.c → gsl_native/multiset.c} +4 -8
- data/ext/{ndlinear.c → gsl_native/ndlinear.c} +320 -321
- data/ext/{nmf.c → gsl_native/nmf.c} +11 -11
- data/ext/{nmf_wrap.c → gsl_native/nmf_wrap.c} +1 -1
- data/ext/{ntuple.c → gsl_native/ntuple.c} +23 -23
- data/ext/{odeiv.c → gsl_native/odeiv.c} +101 -116
- data/ext/gsl_native/ool.c +879 -0
- data/ext/{permutation.c → gsl_native/permutation.c} +39 -37
- data/ext/{poly.c → gsl_native/poly.c} +10 -13
- data/ext/{poly2.c → gsl_native/poly2.c} +16 -16
- data/ext/{poly_source.c → gsl_native/poly_source.h} +249 -293
- data/ext/{qrng.c → gsl_native/qrng.c} +9 -20
- data/ext/{randist.c → gsl_native/randist.c} +222 -247
- data/ext/{rational.c → gsl_native/rational.c} +12 -12
- data/ext/{rng.c → gsl_native/rng.c} +30 -47
- data/ext/{root.c → gsl_native/root.c} +47 -48
- data/ext/{sf.c → gsl_native/sf.c} +196 -244
- data/ext/{sf_airy.c → gsl_native/sf_airy.c} +2 -2
- data/ext/{sf_bessel.c → gsl_native/sf_bessel.c} +7 -7
- data/ext/{sf_clausen.c → gsl_native/sf_clausen.c} +1 -1
- data/ext/{sf_coulomb.c → gsl_native/sf_coulomb.c} +40 -40
- data/ext/{sf_coupling.c → gsl_native/sf_coupling.c} +30 -30
- data/ext/{sf_dawson.c → gsl_native/sf_dawson.c} +1 -1
- data/ext/{sf_debye.c → gsl_native/sf_debye.c} +1 -10
- data/ext/{sf_dilog.c → gsl_native/sf_dilog.c} +1 -1
- data/ext/{sf_elementary.c → gsl_native/sf_elementary.c} +3 -3
- data/ext/{sf_ellint.c → gsl_native/sf_ellint.c} +43 -43
- data/ext/{sf_elljac.c → gsl_native/sf_elljac.c} +3 -3
- data/ext/{sf_erfc.c → gsl_native/sf_erfc.c} +1 -5
- data/ext/{sf_exp.c → gsl_native/sf_exp.c} +3 -3
- data/ext/{sf_expint.c → gsl_native/sf_expint.c} +2 -12
- data/ext/{sf_fermi_dirac.c → gsl_native/sf_fermi_dirac.c} +1 -1
- data/ext/{sf_gamma.c → gsl_native/sf_gamma.c} +2 -6
- data/ext/{sf_gegenbauer.c → gsl_native/sf_gegenbauer.c} +1 -1
- data/ext/{sf_hyperg.c → gsl_native/sf_hyperg.c} +1 -1
- data/ext/{sf_laguerre.c → gsl_native/sf_laguerre.c} +4 -4
- data/ext/{sf_lambert.c → gsl_native/sf_lambert.c} +1 -1
- data/ext/{sf_legendre.c → gsl_native/sf_legendre.c} +1 -1
- data/ext/{sf_log.c → gsl_native/sf_log.c} +4 -4
- data/ext/gsl_native/sf_mathieu.c +235 -0
- data/ext/{sf_power.c → gsl_native/sf_power.c} +1 -1
- data/ext/{sf_psi.c → gsl_native/sf_psi.c} +3 -12
- data/ext/{sf_synchrotron.c → gsl_native/sf_synchrotron.c} +1 -1
- data/ext/{sf_transport.c → gsl_native/sf_transport.c} +1 -1
- data/ext/{sf_trigonometric.c → gsl_native/sf_trigonometric.c} +4 -4
- data/ext/{sf_zeta.c → gsl_native/sf_zeta.c} +1 -5
- data/ext/{signal.c → gsl_native/signal.c} +63 -68
- data/ext/{siman.c → gsl_native/siman.c} +45 -49
- data/ext/{sort.c → gsl_native/sort.c} +6 -7
- data/ext/{spline.c → gsl_native/spline.c} +28 -46
- data/ext/{stats.c → gsl_native/stats.c} +105 -118
- data/ext/{sum.c → gsl_native/sum.c} +34 -34
- data/ext/{tamu_anova.c → gsl_native/tamu_anova.c} +1 -1
- data/ext/{tensor.c → gsl_native/tensor.c} +8 -11
- data/ext/{tensor_source.c → gsl_native/tensor_source.h} +147 -148
- data/ext/{vector.c → gsl_native/vector.c} +11 -14
- data/ext/{vector_complex.c → gsl_native/vector_complex.c} +179 -184
- data/ext/{vector_double.c → gsl_native/vector_double.c} +178 -183
- data/ext/{vector_int.c → gsl_native/vector_int.c} +27 -29
- data/ext/{vector_source.c → gsl_native/vector_source.h} +428 -443
- data/ext/{wavelet.c → gsl_native/wavelet.c} +224 -246
- data/gsl.gemspec +29 -0
- data/lib/gsl.rb +8 -3
- data/lib/gsl/gnuplot.rb +3 -3
- data/lib/gsl/oper.rb +35 -60
- data/lib/gsl/version.rb +3 -0
- data/lib/rbgsl.rb +1 -3
- data/rdoc/alf.rdoc +5 -5
- data/rdoc/blas.rdoc +9 -9
- data/rdoc/bspline.rdoc +17 -17
- data/rdoc/changes.rdoc +4 -9
- data/rdoc/cheb.rdoc +25 -25
- data/rdoc/cholesky_complex.rdoc +21 -21
- data/rdoc/combi.rdoc +37 -37
- data/rdoc/complex.rdoc +22 -22
- data/rdoc/const.rdoc +47 -47
- data/rdoc/dht.rdoc +49 -49
- data/rdoc/diff.rdoc +42 -42
- data/rdoc/ehandling.rdoc +6 -6
- data/rdoc/eigen.rdoc +153 -153
- data/rdoc/fft.rdoc +146 -146
- data/rdoc/fit.rdoc +109 -109
- data/rdoc/function.rdoc +11 -11
- data/rdoc/graph.rdoc +17 -17
- data/rdoc/hist.rdoc +103 -103
- data/rdoc/hist2d.rdoc +42 -42
- data/rdoc/hist3d.rdoc +9 -9
- data/rdoc/integration.rdoc +110 -110
- data/rdoc/interp.rdoc +71 -71
- data/rdoc/intro.rdoc +8 -8
- data/rdoc/linalg.rdoc +188 -188
- data/rdoc/linalg_complex.rdoc +1 -1
- data/rdoc/math.rdoc +58 -58
- data/rdoc/matrix.rdoc +275 -275
- data/rdoc/min.rdoc +57 -57
- data/rdoc/monte.rdoc +22 -22
- data/rdoc/multimin.rdoc +95 -95
- data/rdoc/multiroot.rdoc +80 -80
- data/rdoc/narray.rdoc +32 -32
- data/rdoc/ndlinear.rdoc +54 -54
- data/rdoc/nonlinearfit.rdoc +100 -100
- data/rdoc/ntuple.rdoc +31 -31
- data/rdoc/odeiv.rdoc +88 -88
- data/rdoc/perm.rdoc +90 -90
- data/rdoc/poly.rdoc +66 -66
- data/rdoc/qrng.rdoc +21 -21
- data/rdoc/randist.rdoc +82 -82
- data/rdoc/ref.rdoc +57 -57
- data/rdoc/rng.rdoc +85 -85
- data/rdoc/roots.rdoc +57 -57
- data/rdoc/sf.rdoc +428 -428
- data/rdoc/siman.rdoc +19 -19
- data/rdoc/sort.rdoc +30 -30
- data/rdoc/start.rdoc +8 -8
- data/rdoc/stats.rdoc +52 -52
- data/rdoc/sum.rdoc +12 -12
- data/rdoc/tensor.rdoc +31 -31
- data/rdoc/tut.rdoc +1 -1
- data/rdoc/use.rdoc +39 -39
- data/rdoc/vector.rdoc +188 -188
- data/rdoc/vector_complex.rdoc +24 -24
- data/rdoc/wavelet.rdoc +46 -46
- data/test/gsl/blas_test.rb +79 -0
- data/test/gsl/bspline_test.rb +63 -0
- data/test/gsl/cdf_test.rb +1512 -0
- data/test/gsl/cheb_test.rb +80 -0
- data/test/gsl/combination_test.rb +100 -0
- data/test/gsl/complex_test.rb +20 -0
- data/test/gsl/const_test.rb +29 -0
- data/test/gsl/deriv_test.rb +62 -0
- data/test/gsl/dht_test.rb +79 -0
- data/test/gsl/diff_test.rb +53 -0
- data/test/gsl/eigen_test.rb +563 -0
- data/test/gsl/err_test.rb +23 -0
- data/test/gsl/fit_test.rb +101 -0
- data/test/gsl/histo_test.rb +14 -0
- data/test/gsl/index_test.rb +61 -0
- data/test/gsl/integration_test.rb +274 -0
- data/test/gsl/interp_test.rb +27 -0
- data/test/gsl/linalg_test.rb +463 -0
- data/test/gsl/matrix_nmf_test.rb +37 -0
- data/test/gsl/matrix_test.rb +98 -0
- data/test/gsl/min_test.rb +89 -0
- data/test/gsl/monte_test.rb +77 -0
- data/test/gsl/multifit_test.rb +753 -0
- data/test/gsl/multimin_test.rb +157 -0
- data/test/gsl/multiroot_test.rb +135 -0
- data/test/gsl/multiset_test.rb +52 -0
- data/test/gsl/odeiv_test.rb +275 -0
- data/test/gsl/oper_test.rb +98 -0
- data/test/gsl/poly_test.rb +338 -0
- data/test/gsl/qrng_test.rb +94 -0
- data/test/gsl/quartic_test.rb +28 -0
- data/test/gsl/randist_test.rb +122 -0
- data/test/gsl/rng_test.rb +303 -0
- data/test/gsl/roots_test.rb +78 -0
- data/test/gsl/sf_test.rb +2079 -0
- data/test/gsl/stats_test.rb +122 -0
- data/test/gsl/sum_test.rb +69 -0
- data/test/gsl/tensor_test.rb +396 -0
- data/test/gsl/vector_test.rb +223 -0
- data/test/gsl/wavelet_test.rb +130 -0
- data/test/gsl_test.rb +321 -0
- data/test/test_helper.rb +42 -0
- data/uncrustify.cfg +1693 -0
- metadata +337 -378
- data/README +0 -32
- data/VERSION +0 -1
- data/ext/bspline.c +0 -130
- data/ext/const.c +0 -673
- data/ext/cqp.c +0 -283
- data/ext/extconf.rb +0 -295
- data/ext/fcmp.c +0 -66
- data/ext/fresnel.c +0 -312
- data/ext/jacobi.c +0 -739
- data/ext/ool.c +0 -879
- data/ext/oper_complex_source.c +0 -253
- data/ext/sf_mathieu.c +0 -238
- data/include/rb_gsl_config.h +0 -62
- data/include/rb_gsl_dirac.h +0 -13
- data/rdoc/index.rdoc +0 -62
- data/rdoc/rngextra.rdoc +0 -11
- data/rdoc/screenshot.rdoc +0 -40
- data/setup.rb +0 -1585
- data/tests/blas/amax.rb +0 -14
- data/tests/blas/asum.rb +0 -16
- data/tests/blas/axpy.rb +0 -25
- data/tests/blas/copy.rb +0 -23
- data/tests/blas/dot.rb +0 -23
- data/tests/bspline.rb +0 -53
- data/tests/cdf.rb +0 -1388
- data/tests/cheb.rb +0 -112
- data/tests/combination.rb +0 -123
- data/tests/complex.rb +0 -17
- data/tests/const.rb +0 -24
- data/tests/deriv.rb +0 -85
- data/tests/dht/dht1.rb +0 -17
- data/tests/dht/dht2.rb +0 -23
- data/tests/dht/dht3.rb +0 -23
- data/tests/dht/dht4.rb +0 -23
- data/tests/diff.rb +0 -78
- data/tests/eigen/eigen.rb +0 -220
- data/tests/eigen/gen.rb +0 -105
- data/tests/eigen/genherm.rb +0 -66
- data/tests/eigen/gensymm.rb +0 -68
- data/tests/eigen/nonsymm.rb +0 -53
- data/tests/eigen/nonsymmv.rb +0 -53
- data/tests/eigen/symm-herm.rb +0 -74
- data/tests/err.rb +0 -58
- data/tests/fit.rb +0 -124
- data/tests/gsl_test.rb +0 -118
- data/tests/gsl_test2.rb +0 -110
- data/tests/histo.rb +0 -12
- data/tests/integration/integration1.rb +0 -72
- data/tests/integration/integration2.rb +0 -71
- data/tests/integration/integration3.rb +0 -71
- data/tests/integration/integration4.rb +0 -71
- data/tests/interp.rb +0 -45
- data/tests/linalg/HH.rb +0 -64
- data/tests/linalg/LU.rb +0 -47
- data/tests/linalg/QR.rb +0 -77
- data/tests/linalg/SV.rb +0 -24
- data/tests/linalg/TDN.rb +0 -116
- data/tests/linalg/TDS.rb +0 -122
- data/tests/linalg/bidiag.rb +0 -73
- data/tests/linalg/cholesky.rb +0 -20
- data/tests/linalg/linalg.rb +0 -158
- data/tests/matrix/matrix_complex_test.rb +0 -36
- data/tests/matrix/matrix_nmf_test.rb +0 -39
- data/tests/matrix/matrix_test.rb +0 -48
- data/tests/min.rb +0 -99
- data/tests/monte/miser.rb +0 -31
- data/tests/monte/vegas.rb +0 -45
- data/tests/multifit/test_2dgauss.rb +0 -112
- data/tests/multifit/test_brown.rb +0 -90
- data/tests/multifit/test_enso.rb +0 -246
- data/tests/multifit/test_filip.rb +0 -155
- data/tests/multifit/test_gauss.rb +0 -97
- data/tests/multifit/test_longley.rb +0 -110
- data/tests/multifit/test_multifit.rb +0 -52
- data/tests/multimin.rb +0 -139
- data/tests/multiroot.rb +0 -131
- data/tests/multiset.rb +0 -52
- data/tests/narray/blas_dnrm2.rb +0 -20
- data/tests/odeiv.rb +0 -353
- data/tests/poly/poly.rb +0 -290
- data/tests/poly/special.rb +0 -65
- data/tests/qrng.rb +0 -131
- data/tests/quartic.rb +0 -29
- data/tests/randist.rb +0 -134
- data/tests/rng.rb +0 -305
- data/tests/roots.rb +0 -76
- data/tests/run-test.sh +0 -17
- data/tests/sf/gsl_test_sf.rb +0 -249
- data/tests/sf/test_airy.rb +0 -83
- data/tests/sf/test_bessel.rb +0 -306
- data/tests/sf/test_coulomb.rb +0 -17
- data/tests/sf/test_dilog.rb +0 -25
- data/tests/sf/test_gamma.rb +0 -209
- data/tests/sf/test_hyperg.rb +0 -356
- data/tests/sf/test_legendre.rb +0 -227
- data/tests/sf/test_mathieu.rb +0 -59
- data/tests/sf/test_mode.rb +0 -19
- data/tests/sf/test_sf.rb +0 -839
- data/tests/stats.rb +0 -174
- data/tests/stats_mt.rb +0 -16
- data/tests/sum.rb +0 -98
- data/tests/sys.rb +0 -323
- data/tests/tensor.rb +0 -419
- data/tests/vector/vector_complex_test.rb +0 -101
- data/tests/vector/vector_test.rb +0 -141
- 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 "
|
13
|
-
|
14
|
-
#
|
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
|
-
|
25
|
-
VALUE rb_gsl_eval_pdf_cdf2_uint(VALUE xx, VALUE aa,
|
26
|
-
|
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
|
-
|
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
|
-
|
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
|
-
|
632
|
-
|
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
|