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
@@ -1,23 +1,23 @@
|
|
1
1
|
/* histogram/find.c
|
2
|
-
*
|
2
|
+
*
|
3
3
|
* Copyright (C) 1996, 1997, 1998, 1999, 2000 Brian Gough
|
4
|
-
*
|
4
|
+
*
|
5
5
|
* This program is free software; you can redistribute it and/or modify
|
6
6
|
* it under the terms of the GNU General Public License as published by
|
7
7
|
* the Free Software Foundation; either version 2 of the License, or (at
|
8
8
|
* your option) any later version.
|
9
|
-
*
|
9
|
+
*
|
10
10
|
* This program is distributed in the hope that it will be useful, but
|
11
11
|
* WITHOUT ANY WARRANTY; without even the implied warranty of
|
12
12
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
13
13
|
* General Public License for more details.
|
14
|
-
*
|
14
|
+
*
|
15
15
|
* You should have received a copy of the GNU General Public License
|
16
16
|
* along with this program; if not, write to the Free Software
|
17
17
|
* Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
|
18
18
|
*/
|
19
19
|
|
20
|
-
#include "rb_gsl_histogram.h"
|
20
|
+
#include "include/rb_gsl_histogram.h"
|
21
21
|
/* determines whether to optimize for linear ranges */
|
22
22
|
#define LINEAR_OPT 1
|
23
23
|
|
@@ -26,7 +26,6 @@ int mygsl_find (const size_t n, const double range[], const double x, size_t * i
|
|
26
26
|
size_t i_linear, lower, upper, mid;
|
27
27
|
if (x < range[0]) return -1;
|
28
28
|
if (x >= range[n]) return +1;
|
29
|
-
|
30
29
|
/* optimize for linear case */
|
31
30
|
|
32
31
|
#ifdef LINEAR_OPT
|
@@ -36,47 +35,46 @@ int mygsl_find (const size_t n, const double range[], const double x, size_t * i
|
|
36
35
|
}
|
37
36
|
|
38
37
|
if (x >= range[i_linear] && x < range[i_linear + 1])
|
39
|
-
|
40
|
-
|
41
|
-
|
42
|
-
|
38
|
+
{
|
39
|
+
*i = i_linear;
|
40
|
+
return 0;
|
41
|
+
}
|
43
42
|
#endif
|
44
43
|
|
45
44
|
/* perform binary search */
|
46
45
|
|
47
|
-
upper = n
|
48
|
-
lower = 0
|
46
|
+
upper = n;
|
47
|
+
lower = 0;
|
49
48
|
|
50
49
|
while (upper - lower > 1)
|
50
|
+
{
|
51
|
+
mid = (upper + lower) / 2;
|
52
|
+
|
53
|
+
if (x >= range[mid])
|
51
54
|
{
|
52
|
-
|
53
|
-
|
54
|
-
if (x >= range[mid])
|
55
|
-
{
|
56
|
-
lower = mid ;
|
57
|
-
}
|
58
|
-
else
|
59
|
-
{
|
60
|
-
upper = mid ;
|
61
|
-
}
|
55
|
+
lower = mid;
|
62
56
|
}
|
57
|
+
else
|
58
|
+
{
|
59
|
+
upper = mid;
|
60
|
+
}
|
61
|
+
}
|
63
62
|
|
64
|
-
*i = lower
|
63
|
+
*i = lower;
|
65
64
|
|
66
65
|
/* sanity check the result */
|
67
66
|
|
68
67
|
if (x < range[lower] || x >= range[lower + 1])
|
69
|
-
|
70
|
-
|
71
|
-
|
72
|
-
|
68
|
+
{
|
69
|
+
GSL_ERROR ("x not found in range", GSL_ESANITY);
|
70
|
+
}
|
73
71
|
return 0;
|
74
72
|
}
|
75
73
|
|
76
74
|
int mygsl_find2d (const size_t nx, const double xrange[],
|
77
|
-
|
78
|
-
|
79
|
-
|
75
|
+
const size_t ny, const double yrange[],
|
76
|
+
const double x, const double y,
|
77
|
+
size_t * i, size_t * j)
|
80
78
|
{
|
81
79
|
int status = mygsl_find (nx, xrange, x, i);
|
82
80
|
if (status) return status;
|
@@ -86,10 +84,10 @@ int mygsl_find2d (const size_t nx, const double xrange[],
|
|
86
84
|
}
|
87
85
|
|
88
86
|
int mygsl_find3d (const size_t nx, const double xrange[],
|
89
|
-
|
90
|
-
|
91
|
-
|
92
|
-
|
87
|
+
const size_t ny, const double yrange[],
|
88
|
+
const size_t nz, const double zrange[],
|
89
|
+
const double x, const double y, const double z,
|
90
|
+
size_t * i, size_t * j, size_t *k)
|
93
91
|
{
|
94
92
|
int status = mygsl_find (nx, xrange, x, i);
|
95
93
|
if (status) return status;
|
@@ -18,11 +18,11 @@
|
|
18
18
|
*/
|
19
19
|
/***************************************************************
|
20
20
|
*
|
21
|
-
* File gsl_histogram_oper.c:
|
22
|
-
* Routine to make operation on histograms.
|
21
|
+
* File gsl_histogram_oper.c:
|
22
|
+
* Routine to make operation on histograms.
|
23
23
|
* Need GSL library and header.
|
24
24
|
* Contains the routines:
|
25
|
-
* gsl_histogram_same_binning check if two histograms have the same binning
|
25
|
+
* gsl_histogram_same_binning check if two histograms have the same binning
|
26
26
|
* gsl_histogram_add add two histograms
|
27
27
|
* gsl_histogram_sub subctract two histograms
|
28
28
|
* gsl_histogram_mult multiply two histograms
|
@@ -38,7 +38,7 @@
|
|
38
38
|
#include <gsl/gsl_math.h>
|
39
39
|
#include <gsl/gsl_histogram.h>
|
40
40
|
|
41
|
-
/*
|
41
|
+
/*
|
42
42
|
* gsl_histogram_same_binning:
|
43
43
|
* control if two histograms have the
|
44
44
|
* same binning
|
@@ -48,112 +48,103 @@ int
|
|
48
48
|
mygsl_histogram_equal_bins_p (const gsl_histogram * h1, const gsl_histogram * h2)
|
49
49
|
{
|
50
50
|
if (h1->n != h2->n)
|
51
|
-
|
52
|
-
|
53
|
-
|
54
|
-
|
51
|
+
{
|
52
|
+
return 0;
|
53
|
+
}
|
55
54
|
{
|
56
55
|
size_t i;
|
57
56
|
/* init ranges */
|
58
57
|
|
59
58
|
for (i = 0; i <= h1->n; i++)
|
59
|
+
{
|
60
|
+
if (gsl_fcmp(h1->range[i],h2->range[i], 1e-12))
|
60
61
|
{
|
61
|
-
|
62
|
-
{
|
63
|
-
return 0;
|
64
|
-
}
|
62
|
+
return 0;
|
65
63
|
}
|
64
|
+
}
|
66
65
|
}
|
67
66
|
|
68
67
|
return 1;
|
69
68
|
}
|
70
69
|
|
71
|
-
/*
|
70
|
+
/*
|
72
71
|
* gsl_histogram_add:
|
73
72
|
* add two histograms
|
74
73
|
*/
|
75
|
-
int
|
74
|
+
int
|
76
75
|
mygsl_histogram_add (gsl_histogram * h1, const gsl_histogram * h2)
|
77
76
|
{
|
78
77
|
size_t i;
|
79
78
|
|
80
79
|
if (!mygsl_histogram_equal_bins_p (h1, h2))
|
81
|
-
|
82
|
-
|
83
|
-
|
84
|
-
|
80
|
+
{
|
81
|
+
GSL_ERROR ("histograms have different binning", GSL_EINVAL);
|
82
|
+
}
|
85
83
|
for (i = 0; i < h1->n; i++)
|
86
|
-
|
87
|
-
|
88
|
-
|
89
|
-
|
84
|
+
{
|
85
|
+
h1->bin[i] += h2->bin[i];
|
86
|
+
}
|
90
87
|
return GSL_SUCCESS;
|
91
88
|
}
|
92
89
|
|
93
|
-
/*
|
90
|
+
/*
|
94
91
|
* gsl_histogram_sub:
|
95
92
|
* subtract two histograms
|
96
93
|
*/
|
97
94
|
|
98
|
-
int
|
95
|
+
int
|
99
96
|
mygsl_histogram_sub (gsl_histogram * h1, const gsl_histogram * h2)
|
100
97
|
{
|
101
98
|
size_t i;
|
102
99
|
|
103
100
|
if (!mygsl_histogram_equal_bins_p (h1, h2))
|
104
|
-
|
105
|
-
|
106
|
-
|
107
|
-
|
101
|
+
{
|
102
|
+
GSL_ERROR ("histograms have different binning", GSL_EINVAL);
|
103
|
+
}
|
108
104
|
for (i = 0; i < h1->n; i++)
|
109
|
-
|
110
|
-
|
111
|
-
|
112
|
-
|
105
|
+
{
|
106
|
+
h1->bin[i] -= h2->bin[i];
|
107
|
+
}
|
113
108
|
return GSL_SUCCESS;
|
114
109
|
|
115
110
|
}
|
116
111
|
|
117
|
-
/*
|
112
|
+
/*
|
118
113
|
* gsl_histogram_mult:
|
119
114
|
* multiply two histograms
|
120
115
|
*/
|
121
116
|
|
122
|
-
int
|
117
|
+
int
|
123
118
|
mygsl_histogram_mul (gsl_histogram * h1, const gsl_histogram * h2)
|
124
119
|
{
|
125
120
|
size_t i;
|
126
121
|
|
127
122
|
if (!mygsl_histogram_equal_bins_p (h1, h2))
|
128
|
-
|
129
|
-
|
130
|
-
|
131
|
-
|
123
|
+
{
|
124
|
+
GSL_ERROR ("histograms have different binning", GSL_EINVAL);
|
125
|
+
}
|
132
126
|
for (i = 0; i < h1->n; i++)
|
133
|
-
|
134
|
-
|
135
|
-
|
136
|
-
|
127
|
+
{
|
128
|
+
h1->bin[i] *= h2->bin[i];
|
129
|
+
}
|
137
130
|
return GSL_SUCCESS;
|
138
131
|
}
|
139
|
-
/*
|
132
|
+
/*
|
140
133
|
* gsl_histogram_div:
|
141
134
|
* divide two histograms
|
142
135
|
*/
|
143
|
-
int
|
136
|
+
int
|
144
137
|
mygsl_histogram_div (gsl_histogram * h1, const gsl_histogram * h2)
|
145
138
|
{
|
146
139
|
size_t i;
|
147
140
|
|
148
141
|
if (!mygsl_histogram_equal_bins_p (h1, h2))
|
149
|
-
|
150
|
-
|
151
|
-
|
152
|
-
|
142
|
+
{
|
143
|
+
GSL_ERROR ("histograms have different binning", GSL_EINVAL);
|
144
|
+
}
|
153
145
|
for (i = 0; i < h1->n; i++)
|
154
|
-
|
155
|
-
|
156
|
-
|
157
|
-
|
146
|
+
{
|
147
|
+
h1->bin[i] /= h2->bin[i];
|
148
|
+
}
|
158
149
|
return GSL_SUCCESS;
|
159
150
|
}
|
@@ -9,7 +9,7 @@
|
|
9
9
|
WITHOUT ANY WARRANTY.
|
10
10
|
*/
|
11
11
|
|
12
|
-
#include "rb_gsl.h"
|
12
|
+
#include "include/rb_gsl.h"
|
13
13
|
|
14
14
|
static VALUE rb_gsl_ieee_env_setup(VALUE obj)
|
15
15
|
{
|
@@ -19,15 +19,11 @@ static VALUE rb_gsl_ieee_env_setup(VALUE obj)
|
|
19
19
|
|
20
20
|
static VALUE rb_gsl_ieee_fprintf_double(int argc, VALUE *argv, VALUE obj)
|
21
21
|
{
|
22
|
-
#ifdef RUBY_1_9_LATER
|
23
22
|
rb_io_t *fptr = NULL;
|
24
|
-
double ftmp;
|
25
|
-
#else
|
26
|
-
OpenFile *fptr = NULL;
|
27
|
-
#endif
|
28
23
|
FILE *fp = NULL;
|
29
24
|
int flag = 0;
|
30
25
|
VALUE vtmp;
|
26
|
+
double ftmp;
|
31
27
|
|
32
28
|
switch (argc) {
|
33
29
|
case 2:
|
@@ -39,15 +35,11 @@ static VALUE rb_gsl_ieee_fprintf_double(int argc, VALUE *argv, VALUE obj)
|
|
39
35
|
case T_FILE:
|
40
36
|
GetOpenFile(argv[0], fptr);
|
41
37
|
rb_io_check_writable(fptr);
|
42
|
-
#ifdef RUBY_1_9_LATER
|
43
38
|
fp = rb_io_stdio_file(fptr);
|
44
|
-
#else
|
45
|
-
fp = GetWriteFile(fptr);
|
46
|
-
#endif
|
47
39
|
break;
|
48
40
|
default:
|
49
41
|
rb_raise(rb_eTypeError, "wrong type argument %s (IO or String expected)",
|
50
|
-
|
42
|
+
rb_class2name(CLASS_OF(argv[0])));
|
51
43
|
}
|
52
44
|
vtmp = argv[1];
|
53
45
|
break;
|
@@ -60,13 +52,9 @@ static VALUE rb_gsl_ieee_fprintf_double(int argc, VALUE *argv, VALUE obj)
|
|
60
52
|
}
|
61
53
|
if (TYPE(vtmp) != T_FLOAT)
|
62
54
|
rb_raise(rb_eTypeError, "wrong argument type %s (Float expected)",
|
63
|
-
|
64
|
-
#ifdef RUBY_1_9_LATER
|
55
|
+
rb_class2name(CLASS_OF(vtmp)));
|
65
56
|
ftmp = RFLOAT_VALUE(vtmp);
|
66
57
|
gsl_ieee_fprintf_double(fp, &ftmp);
|
67
|
-
#else
|
68
|
-
gsl_ieee_fprintf_double(fp, &(RFLOAT(vtmp)->value));
|
69
|
-
#endif
|
70
58
|
if (fp == stdout) fprintf(stdout, "\n");
|
71
59
|
if (flag == 1) fclose(fp);
|
72
60
|
return obj;
|
@@ -86,15 +74,15 @@ void Init_gsl_ieee(VALUE module)
|
|
86
74
|
mgsl_ieee = rb_define_module_under(module, "IEEE");
|
87
75
|
|
88
76
|
rb_define_singleton_method(mgsl_ieee, "env_setup",
|
89
|
-
|
77
|
+
rb_gsl_ieee_env_setup, 0);
|
90
78
|
rb_define_module_function(module, "ieee_env_setup", rb_gsl_ieee_env_setup, 0);
|
91
79
|
rb_define_singleton_method(mgsl_ieee, "fprintf_double",
|
92
|
-
|
80
|
+
rb_gsl_ieee_fprintf_double, -1);
|
93
81
|
rb_define_singleton_method(mgsl_ieee, "fprintf",
|
94
|
-
|
82
|
+
rb_gsl_ieee_fprintf_double, -1);
|
95
83
|
rb_define_singleton_method(mgsl_ieee, "printf",
|
96
|
-
|
84
|
+
rb_gsl_ieee_printf_double, -1);
|
97
85
|
rb_define_singleton_method(mgsl_ieee, "printf_double",
|
98
|
-
|
86
|
+
rb_gsl_ieee_printf_double, -1);
|
99
87
|
|
100
88
|
}
|
@@ -12,40 +12,25 @@
|
|
12
12
|
#ifndef ___RB_GSL_H___
|
13
13
|
#define ___RB_GSL_H___
|
14
14
|
|
15
|
-
#include <stdio.h>
|
16
15
|
#include <string.h>
|
17
|
-
#include <ctype.h>
|
18
16
|
#include <stdlib.h>
|
19
17
|
#include <gsl/gsl_mode.h>
|
20
|
-
#include "rb_gsl_common.h"
|
21
|
-
#include "rb_gsl_math.h"
|
22
|
-
#include "rb_gsl_complex.h"
|
23
|
-
#include "rb_gsl_array.h"
|
24
|
-
#include "rb_gsl_function.h"
|
25
18
|
#include "rb_gsl_poly.h"
|
26
|
-
#include "rb_gsl_sf.h"
|
27
19
|
#include "rb_gsl_linalg.h"
|
20
|
+
#include "rb_gsl_math.h"
|
21
|
+
#include "rb_gsl_root.h"
|
22
|
+
#include "rb_gsl_histogram.h"
|
23
|
+
#include "rb_gsl_sf.h"
|
28
24
|
#include "rb_gsl_eigen.h"
|
29
25
|
#include "rb_gsl_fft.h"
|
30
26
|
#include "rb_gsl_integration.h"
|
31
27
|
#include "rb_gsl_rng.h"
|
32
28
|
#include "rb_gsl_statistics.h"
|
33
|
-
#include "rb_gsl_histogram.h"
|
34
29
|
#include "rb_gsl_odeiv.h"
|
35
30
|
#include "rb_gsl_interp.h"
|
36
31
|
#include "rb_gsl_cheb.h"
|
37
|
-
#include "rb_gsl_root.h"
|
38
32
|
#include "rb_gsl_fit.h"
|
39
33
|
#include "rb_gsl_const.h"
|
40
|
-
#include "rb_gsl_config.h"
|
41
|
-
|
42
|
-
|
43
|
-
#include "ruby.h"
|
44
|
-
#ifdef RUBY_1_9_LATER
|
45
|
-
#include "ruby/io.h"
|
46
|
-
#else
|
47
|
-
#include "rubyio.h"
|
48
|
-
#endif
|
49
34
|
|
50
35
|
void Init_gsl_error(VALUE module);
|
51
36
|
void Init_gsl_math(VALUE module);
|
@@ -82,9 +67,7 @@ void Init_gsl_odeiv(VALUE module);
|
|
82
67
|
void Init_gsl_interp(VALUE module);
|
83
68
|
void Init_gsl_spline(VALUE module);
|
84
69
|
void Init_gsl_diff(VALUE module);
|
85
|
-
#ifdef GSL_1_4_9_LATER
|
86
70
|
void Init_gsl_deriv(VALUE module);
|
87
|
-
#endif
|
88
71
|
|
89
72
|
void Init_gsl_cheb(VALUE module);
|
90
73
|
void Init_gsl_sum(VALUE module);
|
@@ -138,9 +121,7 @@ void Init_cqp(VALUE module);
|
|
138
121
|
|
139
122
|
void Init_fresnel(VALUE module);
|
140
123
|
|
141
|
-
#ifdef GSL_1_9_LATER
|
142
124
|
void Init_bspline(VALUE module);
|
143
|
-
#endif
|
144
125
|
|
145
126
|
#ifdef HAVE_ALF_ALF_H
|
146
127
|
#include "alf/alf.h"
|
@@ -148,11 +129,8 @@ void Init_bspline(VALUE module);
|
|
148
129
|
void Init_alf(VALUE module);
|
149
130
|
void Init_geometry(VALUE module);
|
150
131
|
|
151
|
-
#ifdef GSL_1_14_LATER
|
152
132
|
#include <gsl/gsl_multiset.h>
|
153
133
|
extern VALUE cMultiset;
|
154
134
|
void Init_multiset(VALUE module);
|
155
135
|
|
156
136
|
#endif
|
157
|
-
|
158
|
-
#endif
|
@@ -13,15 +13,10 @@
|
|
13
13
|
#define ___RB_GSL_ARRAY_H___
|
14
14
|
|
15
15
|
#include <math.h>
|
16
|
-
#include <gsl/gsl_errno.h>
|
17
|
-
#include <gsl/gsl_complex.h>
|
18
16
|
#include <gsl/gsl_block.h>
|
19
|
-
#include <gsl/gsl_vector.h>
|
20
17
|
#include <gsl/gsl_vector_complex.h>
|
21
18
|
#include <gsl/gsl_vector_complex_double.h>
|
22
|
-
#include <gsl/gsl_matrix.h>
|
23
19
|
#include <gsl/gsl_matrix_complex_double.h>
|
24
|
-
#include <gsl/gsl_permutation.h>
|
25
20
|
#include <gsl/gsl_combination.h>
|
26
21
|
#include <gsl/gsl_sort.h>
|
27
22
|
#include <gsl/gsl_sort_vector.h>
|
@@ -29,15 +24,6 @@
|
|
29
24
|
#include <gsl/gsl_eigen.h>
|
30
25
|
#include <gsl/gsl_odeiv.h>
|
31
26
|
#include <gsl/gsl_blas.h>
|
32
|
-
#include "rb_gsl_config.h"
|
33
|
-
|
34
|
-
#include "ruby.h"
|
35
|
-
#ifdef RUBY_1_9_LATER
|
36
|
-
#include "ruby/io.h"
|
37
|
-
#else
|
38
|
-
#include "rubyio.h"
|
39
|
-
#endif
|
40
|
-
|
41
27
|
#include "rb_gsl_common.h"
|
42
28
|
|
43
29
|
typedef gsl_permutation gsl_index;
|
@@ -123,14 +109,14 @@ double* get_vector_ptr(VALUE ary, size_t *stride, size_t *n);
|
|
123
109
|
|
124
110
|
gsl_matrix_complex* matrix_to_complex(const gsl_matrix *m);
|
125
111
|
|
126
|
-
void gsl_matrix_complex_mul(gsl_matrix_complex *mnew, const gsl_matrix_complex *m,
|
127
|
-
|
112
|
+
void gsl_matrix_complex_mul(gsl_matrix_complex *mnew, const gsl_matrix_complex *m,
|
113
|
+
const gsl_matrix_complex *mb);
|
128
114
|
void gsl_matrix_mul(gsl_matrix *mnew, gsl_matrix *m, gsl_matrix *b);
|
129
|
-
void gsl_matrix_complex_mul_vector(gsl_vector_complex *vnew,
|
130
|
-
|
131
|
-
|
132
|
-
void gsl_matrix_mul_vector(gsl_vector *vnew,
|
133
|
-
|
115
|
+
void gsl_matrix_complex_mul_vector(gsl_vector_complex *vnew,
|
116
|
+
const gsl_matrix_complex *m,
|
117
|
+
const gsl_vector_complex *v);
|
118
|
+
void gsl_matrix_mul_vector(gsl_vector *vnew,
|
119
|
+
const gsl_matrix *m, const gsl_vector *v);
|
134
120
|
gsl_vector_complex* vector_to_complex(const gsl_vector *v);
|
135
121
|
|
136
122
|
gsl_vector* make_vector_clone(const gsl_vector *v);
|
@@ -164,8 +150,8 @@ gsl_matrix_int_view* rb_gsl_matrix_int_view_alloc();
|
|
164
150
|
void rb_gsl_matrix_int_view_free(gsl_matrix_int_view *v);
|
165
151
|
VALUE rb_gsl_matrix_to_i(VALUE obj);
|
166
152
|
VALUE rb_gsl_matrix_int_to_f(VALUE obj);
|
167
|
-
void gsl_matrix_int_mul_vector(gsl_vector_int *vnew,
|
168
|
-
|
153
|
+
void gsl_matrix_int_mul_vector(gsl_vector_int *vnew,
|
154
|
+
const gsl_matrix_int *m, const gsl_vector_int *v);
|
169
155
|
VALUE rb_gsl_vector_to_i(VALUE obj);
|
170
156
|
VALUE make_rarray_from_cvector_int(const gsl_vector_int *v);
|
171
157
|
VALUE rb_gsl_vector_int_to_f(VALUE obj);
|
@@ -174,34 +160,24 @@ VALUE rb_gsl_vector_print(VALUE obj);
|
|
174
160
|
void gsl_vector_print(const gsl_vector *v, VALUE klass);
|
175
161
|
int rbgsl_vector_equal(const gsl_vector *v1, const gsl_vector *v2, double eps);
|
176
162
|
|
177
|
-
#ifndef GSL_1_2_LATER
|
178
|
-
int gsl_matrix_complex_add(gsl_matrix_complex * a, const gsl_matrix_complex * b);
|
179
|
-
int gsl_matrix_complex_sub(gsl_matrix_complex * a, const gsl_matrix_complex * b);
|
180
|
-
int gsl_matrix_complex_mul_elements(gsl_matrix_complex * a, const gsl_matrix_complex * b);
|
181
|
-
int gsl_matrix_complex_div_elements(gsl_matrix_complex * a, const gsl_matrix_complex * b);
|
182
|
-
int gsl_matrix_complex_scale(gsl_matrix_complex * a, const gsl_complex x);
|
183
|
-
int gsl_matrix_complex_add_constant(gsl_matrix_complex * a, const gsl_complex x);
|
184
|
-
int gsl_matrix_complex_add_diagonal(gsl_matrix_complex * a, const gsl_complex x);
|
185
|
-
#endif
|
186
|
-
|
187
163
|
void Init_gsl_vector_init(VALUE module);
|
188
164
|
void Init_gsl_vector_int_init(VALUE module);
|
189
165
|
void Init_gsl_matrix_init(VALUE module);
|
190
166
|
void Init_gsl_matrix_int_init(VALUE module);
|
191
167
|
|
192
|
-
gsl_matrix* gsl_matrix_alloc_from_array_sizes(VALUE ary,
|
193
|
-
|
168
|
+
gsl_matrix* gsl_matrix_alloc_from_array_sizes(VALUE ary,
|
169
|
+
VALUE nn1, VALUE nn2);
|
194
170
|
gsl_matrix* gsl_matrix_alloc_from_arrays(int argc, VALUE *argv);
|
195
171
|
|
196
172
|
gsl_matrix* gsl_matrix_alloc_from_vector_sizes(VALUE ary,
|
197
|
-
|
173
|
+
VALUE nn1, VALUE nn2);
|
198
174
|
gsl_matrix* gsl_matrix_alloc_from_vectors(int argc, VALUE *argv);
|
199
|
-
gsl_matrix_int* gsl_matrix_int_alloc_from_array_sizes(VALUE ary,
|
200
|
-
|
175
|
+
gsl_matrix_int* gsl_matrix_int_alloc_from_array_sizes(VALUE ary,
|
176
|
+
VALUE nn1, VALUE nn2);
|
201
177
|
gsl_matrix_int* gsl_matrix_int_alloc_from_arrays(int argc, VALUE *argv);
|
202
178
|
|
203
179
|
gsl_matrix_int* gsl_matrix_int_alloc_from_vector_sizes(VALUE ary,
|
204
|
-
|
180
|
+
VALUE nn1, VALUE nn2);
|
205
181
|
gsl_matrix_int* gsl_matrix_int_alloc_from_vectors(int argc, VALUE *argv);
|
206
182
|
|
207
183
|
VALUE rb_gsl_matrix_do_something(VALUE obj, void (*f)(gsl_matrix *));
|