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
@@ -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)
|
@@ -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
|
-
|
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
|
-
|
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
|
-
|
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
|
-
|
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
|
-
|
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
|
-
|
500
|
+
rb_class2name(CLASS_OF(m)));
|
501
501
|
break;
|
502
502
|
}
|
503
503
|
default:
|
@@ -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
|
-
|
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
|
-
|
30
|
+
NUM2DBL(Z), NUM2DBL(r)));
|
31
31
|
}
|
32
32
|
|
33
|
-
static VALUE rb_gsl_sf_hydrogenicR_e(VALUE obj, VALUE n, VALUE l,
|
34
|
-
|
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
|
-
|
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
|
-
|
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
|
-
|
63
|
+
FIX2INT(k), F, Fp, G, Gp, &exp_F, &exp_G);
|
64
64
|
return rb_ary_new3(7, vF, vFp, vG, vGp,
|
65
|
-
|
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
|
-
|
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
|
-
|
81
|
-
|
82
|
-
return rb_ary_new3(3, Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v),
|
83
|
-
|
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
|
-
|
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
|
-
|
102
|
-
|
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
|
-
|
107
|
-
|
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
|
-
|
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
|
-
|
128
|
-
|
129
|
-
|
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
|
-
|
136
|
-
|
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
|
-
|
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
|
-
|
152
|
-
return rb_ary_new3(3, Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v),
|
153
|
-
|
154
|
-
|
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
|
-
|
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);
|
@@ -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
|
-
|
21
|
-
|
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
|
-
|
35
|
-
|
36
|
-
|
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
|
-
|
47
|
-
|
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
|
-
|
61
|
-
|
62
|
-
|
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
|
-
|
68
|
-
|
69
|
-
|
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
|
-
|
76
|
-
|
77
|
-
|
78
|
-
|
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
|
-
|
83
|
-
|
84
|
-
|
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
|
-
|
96
|
-
|
97
|
-
|
98
|
-
|
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);
|
@@ -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
|
}
|