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