rb-gsl 1.16.0.5 → 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 +4 -4
- data/rb-gsl.gemspec +5 -10
- metadata +10 -713
- data/.gitignore +0 -12
- data/.travis.yml +0 -24
- data/AUTHORS +0 -12
- data/COPYING +0 -341
- data/ChangeLog +0 -621
- data/Gemfile +0 -4
- data/README.md +0 -77
- data/Rakefile +0 -20
- data/THANKS +0 -17
- data/examples/alf/alf.gp +0 -15
- data/examples/alf/alf.rb +0 -32
- data/examples/blas/blas.rb +0 -13
- data/examples/blas/dnrm2.rb +0 -16
- data/examples/blas/level1.rb +0 -81
- data/examples/blas/level2.rb +0 -11
- data/examples/blas/level3.rb +0 -12
- data/examples/bspline.rb +0 -57
- data/examples/cdf.rb +0 -16
- data/examples/cheb.rb +0 -21
- data/examples/combination.rb +0 -23
- data/examples/complex/RC-lpf.rb +0 -47
- data/examples/complex/add.rb +0 -36
- data/examples/complex/coerce.rb +0 -14
- data/examples/complex/complex.rb +0 -25
- data/examples/complex/fpmi.rb +0 -70
- data/examples/complex/functions.rb +0 -77
- data/examples/complex/michelson.rb +0 -36
- data/examples/complex/mul.rb +0 -28
- data/examples/complex/oscillator.rb +0 -17
- data/examples/complex/set.rb +0 -37
- data/examples/const/physconst.rb +0 -151
- data/examples/const/travel.rb +0 -45
- data/examples/deriv/demo.rb +0 -13
- data/examples/deriv/deriv.rb +0 -36
- data/examples/deriv/diff.rb +0 -35
- data/examples/dht.rb +0 -42
- data/examples/dirac.rb +0 -56
- data/examples/eigen/eigen.rb +0 -34
- data/examples/eigen/herm.rb +0 -22
- data/examples/eigen/narray.rb +0 -9
- data/examples/eigen/nonsymm.rb +0 -37
- data/examples/eigen/nonsymmv.rb +0 -43
- data/examples/eigen/qhoscillator.gp +0 -35
- data/examples/eigen/qhoscillator.rb +0 -90
- data/examples/eigen/vander.rb +0 -41
- data/examples/fft/fft.rb +0 -17
- data/examples/fft/fft2.rb +0 -17
- data/examples/fft/forward.rb +0 -25
- data/examples/fft/forward2.rb +0 -26
- data/examples/fft/radix2.rb +0 -18
- data/examples/fft/real-halfcomplex.rb +0 -33
- data/examples/fft/real-halfcomplex2.rb +0 -30
- data/examples/fft/realradix2.rb +0 -19
- data/examples/fft/sunspot.dat +0 -256
- data/examples/fft/sunspot.rb +0 -16
- data/examples/fit/expdata.dat +0 -20
- data/examples/fit/expfit.rb +0 -31
- data/examples/fit/gaussfit.rb +0 -29
- data/examples/fit/gaussian_2peaks.rb +0 -34
- data/examples/fit/hillfit.rb +0 -40
- data/examples/fit/lognormal.rb +0 -26
- data/examples/fit/lorentzfit.rb +0 -22
- data/examples/fit/multifit.rb +0 -72
- data/examples/fit/ndlinear.rb +0 -133
- data/examples/fit/nonlinearfit.rb +0 -89
- data/examples/fit/plot.gp +0 -36
- data/examples/fit/polyfit.rb +0 -9
- data/examples/fit/powerfit.rb +0 -21
- data/examples/fit/sigmoidfit.rb +0 -40
- data/examples/fit/sinfit.rb +0 -22
- data/examples/fit/wlinear.rb +0 -46
- data/examples/fresnel.rb +0 -11
- data/examples/function/function.rb +0 -36
- data/examples/function/log.rb +0 -7
- data/examples/function/min.rb +0 -33
- data/examples/function/sin.rb +0 -10
- data/examples/function/synchrotron.rb +0 -18
- data/examples/gallery/butterfly.rb +0 -7
- data/examples/gallery/cayley.rb +0 -12
- data/examples/gallery/cornu.rb +0 -23
- data/examples/gallery/eight.rb +0 -11
- data/examples/gallery/koch.rb +0 -40
- data/examples/gallery/lemniscate.rb +0 -11
- data/examples/gallery/polar.rb +0 -11
- data/examples/gallery/rgplot/cossin.rb +0 -35
- data/examples/gallery/rgplot/rgplot.replaced +0 -0
- data/examples/gallery/rgplot/roesller.rb +0 -55
- data/examples/gallery/roesller.rb +0 -39
- data/examples/gallery/scarabaeus.rb +0 -14
- data/examples/histogram/cauchy.rb +0 -27
- data/examples/histogram/cauchy.sh +0 -2
- data/examples/histogram/exponential.rb +0 -19
- data/examples/histogram/gauss.rb +0 -16
- data/examples/histogram/gsl-histogram.rb +0 -40
- data/examples/histogram/histo2d.rb +0 -31
- data/examples/histogram/histo3d.rb +0 -34
- data/examples/histogram/histogram-pdf.rb +0 -27
- data/examples/histogram/histogram.rb +0 -26
- data/examples/histogram/integral.rb +0 -28
- data/examples/histogram/poisson.rb +0 -27
- data/examples/histogram/power.rb +0 -25
- data/examples/histogram/rebin.rb +0 -17
- data/examples/histogram/smp.dat +0 -5
- data/examples/histogram/xexp.rb +0 -21
- data/examples/integration/ahmed.rb +0 -21
- data/examples/integration/cosmology.rb +0 -75
- data/examples/integration/friedmann.gp +0 -16
- data/examples/integration/friedmann.rb +0 -35
- data/examples/integration/gamma-zeta.rb +0 -35
- data/examples/integration/integration.rb +0 -22
- data/examples/integration/qag.rb +0 -8
- data/examples/integration/qag2.rb +0 -14
- data/examples/integration/qag3.rb +0 -8
- data/examples/integration/qagi.rb +0 -28
- data/examples/integration/qagi2.rb +0 -49
- data/examples/integration/qagiu.rb +0 -29
- data/examples/integration/qagp.rb +0 -20
- data/examples/integration/qags.rb +0 -14
- data/examples/integration/qawc.rb +0 -18
- data/examples/integration/qawf.rb +0 -41
- data/examples/integration/qawo.rb +0 -29
- data/examples/integration/qaws.rb +0 -30
- data/examples/integration/qng.rb +0 -17
- data/examples/interp/demo.gp +0 -20
- data/examples/interp/demo.rb +0 -45
- data/examples/interp/interp.rb +0 -37
- data/examples/interp/points +0 -10
- data/examples/interp/spline.rb +0 -20
- data/examples/jacobi/deriv.rb +0 -40
- data/examples/jacobi/integrate.rb +0 -34
- data/examples/jacobi/interp.rb +0 -43
- data/examples/jacobi/jacobi.rb +0 -11
- data/examples/linalg/HH.rb +0 -15
- data/examples/linalg/HH_narray.rb +0 -13
- data/examples/linalg/LQ_solve.rb +0 -73
- data/examples/linalg/LU.rb +0 -84
- data/examples/linalg/LU2.rb +0 -31
- data/examples/linalg/LU_narray.rb +0 -24
- data/examples/linalg/PTLQ.rb +0 -47
- data/examples/linalg/QR.rb +0 -18
- data/examples/linalg/QRPT.rb +0 -47
- data/examples/linalg/QR_solve.rb +0 -78
- data/examples/linalg/QR_solve_narray.rb +0 -13
- data/examples/linalg/SV.rb +0 -16
- data/examples/linalg/SV_narray.rb +0 -12
- data/examples/linalg/SV_solve.rb +0 -49
- data/examples/linalg/chol.rb +0 -29
- data/examples/linalg/chol_narray.rb +0 -15
- data/examples/linalg/complex.rb +0 -57
- data/examples/linalg/invert_narray.rb +0 -10
- data/examples/math/const.rb +0 -67
- data/examples/math/elementary.rb +0 -35
- data/examples/math/functions.rb +0 -41
- data/examples/math/inf_nan.rb +0 -34
- data/examples/math/minmax.rb +0 -22
- data/examples/math/power.rb +0 -18
- data/examples/math/test.rb +0 -31
- data/examples/matrix/a.dat +0 -0
- data/examples/matrix/add.rb +0 -45
- data/examples/matrix/b.dat +0 -4
- data/examples/matrix/cat.rb +0 -31
- data/examples/matrix/colvectors.rb +0 -24
- data/examples/matrix/complex.rb +0 -41
- data/examples/matrix/det.rb +0 -29
- data/examples/matrix/diagonal.rb +0 -23
- data/examples/matrix/get_all.rb +0 -159
- data/examples/matrix/hilbert.rb +0 -31
- data/examples/matrix/iterator.rb +0 -19
- data/examples/matrix/matrix.rb +0 -57
- data/examples/matrix/minmax.rb +0 -53
- data/examples/matrix/mul.rb +0 -39
- data/examples/matrix/rand.rb +0 -20
- data/examples/matrix/read.rb +0 -29
- data/examples/matrix/rowcol.rb +0 -47
- data/examples/matrix/set.rb +0 -41
- data/examples/matrix/set_all.rb +0 -100
- data/examples/matrix/view.rb +0 -32
- data/examples/matrix/view_all.rb +0 -148
- data/examples/matrix/write.rb +0 -23
- data/examples/min.rb +0 -29
- data/examples/monte/miser.rb +0 -47
- data/examples/monte/monte.rb +0 -47
- data/examples/monte/plain.rb +0 -47
- data/examples/monte/vegas.rb +0 -46
- data/examples/multimin/bundle.rb +0 -66
- data/examples/multimin/cqp.rb +0 -109
- data/examples/multimin/fdfminimizer.rb +0 -40
- data/examples/multimin/fminimizer.rb +0 -41
- data/examples/multiroot/demo.rb +0 -36
- data/examples/multiroot/fdfsolver.rb +0 -50
- data/examples/multiroot/fsolver.rb +0 -33
- data/examples/multiroot/fsolver2.rb +0 -32
- data/examples/multiroot/fsolver3.rb +0 -26
- data/examples/narray/histogram.rb +0 -14
- data/examples/narray/mandel.rb +0 -27
- data/examples/narray/narray.rb +0 -28
- data/examples/narray/narray2.rb +0 -44
- data/examples/narray/sf.rb +0 -26
- data/examples/ntuple/create.rb +0 -17
- data/examples/ntuple/project.rb +0 -31
- data/examples/odeiv/binarysystem.gp +0 -23
- data/examples/odeiv/binarysystem.rb +0 -104
- data/examples/odeiv/demo.gp +0 -24
- data/examples/odeiv/demo.rb +0 -69
- data/examples/odeiv/demo2.gp +0 -26
- data/examples/odeiv/duffing.rb +0 -45
- data/examples/odeiv/frei1.rb +0 -109
- data/examples/odeiv/frei2.rb +0 -76
- data/examples/odeiv/legendre.rb +0 -52
- data/examples/odeiv/odeiv.rb +0 -32
- data/examples/odeiv/odeiv2.rb +0 -45
- data/examples/odeiv/oscillator.rb +0 -42
- data/examples/odeiv/sedov.rb +0 -97
- data/examples/odeiv/whitedwarf.gp +0 -40
- data/examples/odeiv/whitedwarf.rb +0 -158
- data/examples/ool/conmin.rb +0 -100
- data/examples/ool/gencan.rb +0 -99
- data/examples/ool/pgrad.rb +0 -100
- data/examples/ool/spg.rb +0 -100
- data/examples/pdf/bernoulli.rb +0 -5
- data/examples/pdf/beta.rb +0 -7
- data/examples/pdf/binomiral.rb +0 -10
- data/examples/pdf/cauchy.rb +0 -6
- data/examples/pdf/chisq.rb +0 -8
- data/examples/pdf/exponential.rb +0 -7
- data/examples/pdf/exppow.rb +0 -6
- data/examples/pdf/fdist.rb +0 -7
- data/examples/pdf/flat.rb +0 -7
- data/examples/pdf/gamma.rb +0 -8
- data/examples/pdf/gauss-tail.rb +0 -5
- data/examples/pdf/gauss.rb +0 -6
- data/examples/pdf/geometric.rb +0 -5
- data/examples/pdf/gumbel.rb +0 -6
- data/examples/pdf/hypergeometric.rb +0 -11
- data/examples/pdf/landau.rb +0 -5
- data/examples/pdf/laplace.rb +0 -7
- data/examples/pdf/logarithmic.rb +0 -5
- data/examples/pdf/logistic.rb +0 -6
- data/examples/pdf/lognormal.rb +0 -6
- data/examples/pdf/neg-binomiral.rb +0 -10
- data/examples/pdf/pareto.rb +0 -7
- data/examples/pdf/pascal.rb +0 -10
- data/examples/pdf/poisson.rb +0 -5
- data/examples/pdf/rayleigh-tail.rb +0 -6
- data/examples/pdf/rayleigh.rb +0 -6
- data/examples/pdf/tdist.rb +0 -6
- data/examples/pdf/weibull.rb +0 -8
- data/examples/permutation/ex1.rb +0 -22
- data/examples/permutation/permutation.rb +0 -16
- data/examples/poly/bell.rb +0 -6
- data/examples/poly/bessel.rb +0 -6
- data/examples/poly/cheb.rb +0 -6
- data/examples/poly/cheb_II.rb +0 -6
- data/examples/poly/cubic.rb +0 -9
- data/examples/poly/demo.rb +0 -20
- data/examples/poly/eval.rb +0 -28
- data/examples/poly/eval_derivs.rb +0 -14
- data/examples/poly/fit.rb +0 -21
- data/examples/poly/hermite.rb +0 -6
- data/examples/poly/poly.rb +0 -13
- data/examples/poly/quadratic.rb +0 -25
- data/examples/random/diffusion.rb +0 -34
- data/examples/random/gaussian.rb +0 -9
- data/examples/random/generator.rb +0 -27
- data/examples/random/hdsobol.rb +0 -21
- data/examples/random/poisson.rb +0 -9
- data/examples/random/qrng.rb +0 -19
- data/examples/random/randomwalk.rb +0 -37
- data/examples/random/randomwalk2d.rb +0 -19
- data/examples/random/rayleigh.rb +0 -36
- data/examples/random/rng.rb +0 -33
- data/examples/random/rngextra.rb +0 -14
- data/examples/roots/bisection.rb +0 -25
- data/examples/roots/brent.rb +0 -43
- data/examples/roots/demo.rb +0 -30
- data/examples/roots/newton.rb +0 -46
- data/examples/roots/recombination.gp +0 -11
- data/examples/roots/recombination.rb +0 -61
- data/examples/roots/steffenson.rb +0 -48
- data/examples/sf/ShiChi.rb +0 -6
- data/examples/sf/SiCi.rb +0 -6
- data/examples/sf/airy_Ai.rb +0 -8
- data/examples/sf/airy_Bi.rb +0 -8
- data/examples/sf/bessel_IK.rb +0 -12
- data/examples/sf/bessel_JY.rb +0 -13
- data/examples/sf/beta_inc.rb +0 -9
- data/examples/sf/clausen.rb +0 -6
- data/examples/sf/dawson.rb +0 -5
- data/examples/sf/debye.rb +0 -9
- data/examples/sf/dilog.rb +0 -6
- data/examples/sf/ellint.rb +0 -6
- data/examples/sf/expint.rb +0 -8
- data/examples/sf/fermi.rb +0 -10
- data/examples/sf/gamma_inc_P.rb +0 -9
- data/examples/sf/gegenbauer.rb +0 -8
- data/examples/sf/hyperg.rb +0 -7
- data/examples/sf/laguerre.rb +0 -19
- data/examples/sf/lambertW.rb +0 -5
- data/examples/sf/legendre_P.rb +0 -10
- data/examples/sf/lngamma.rb +0 -5
- data/examples/sf/psi.rb +0 -54
- data/examples/sf/sphbessel.gp +0 -27
- data/examples/sf/sphbessel.rb +0 -30
- data/examples/sf/synchrotron.rb +0 -5
- data/examples/sf/transport.rb +0 -10
- data/examples/sf/zetam1.rb +0 -5
- data/examples/siman.rb +0 -44
- data/examples/sort/heapsort.rb +0 -23
- data/examples/sort/heapsort_vector_complex.rb +0 -21
- data/examples/sort/sort.rb +0 -23
- data/examples/sort/sort2.rb +0 -16
- data/examples/stats/mean.rb +0 -17
- data/examples/stats/statistics.rb +0 -18
- data/examples/stats/test.rb +0 -9
- data/examples/sum.rb +0 -34
- data/examples/tamu_anova.rb +0 -18
- data/examples/vector/a.dat +0 -0
- data/examples/vector/add.rb +0 -56
- data/examples/vector/b.dat +0 -4
- data/examples/vector/c.dat +0 -3
- data/examples/vector/collect.rb +0 -26
- data/examples/vector/compare.rb +0 -28
- data/examples/vector/complex.rb +0 -51
- data/examples/vector/complex_get_all.rb +0 -85
- data/examples/vector/complex_set_all.rb +0 -131
- data/examples/vector/complex_view_all.rb +0 -77
- data/examples/vector/connect.rb +0 -22
- data/examples/vector/decimate.rb +0 -38
- data/examples/vector/diff.rb +0 -31
- data/examples/vector/filescan.rb +0 -17
- data/examples/vector/floor.rb +0 -23
- data/examples/vector/get_all.rb +0 -82
- data/examples/vector/gnuplot.rb +0 -38
- data/examples/vector/graph.rb +0 -28
- data/examples/vector/histogram.rb +0 -22
- data/examples/vector/linspace.rb +0 -24
- data/examples/vector/log.rb +0 -17
- data/examples/vector/logic.rb +0 -33
- data/examples/vector/logspace.rb +0 -25
- data/examples/vector/minmax.rb +0 -47
- data/examples/vector/mul.rb +0 -49
- data/examples/vector/narray.rb +0 -46
- data/examples/vector/read.rb +0 -29
- data/examples/vector/set.rb +0 -35
- data/examples/vector/set_all.rb +0 -121
- data/examples/vector/smpv.dat +0 -15
- data/examples/vector/test.rb +0 -43
- data/examples/vector/test_gslblock.rb +0 -58
- data/examples/vector/vector.rb +0 -110
- data/examples/vector/view.rb +0 -35
- data/examples/vector/view_all.rb +0 -73
- data/examples/vector/where.rb +0 -29
- data/examples/vector/write.rb +0 -24
- data/examples/vector/zip.rb +0 -34
- data/examples/wavelet/ecg.dat +0 -256
- data/examples/wavelet/wavelet1.rb +0 -50
- data/ext/gsl_native/alf.c +0 -206
- data/ext/gsl_native/array.c +0 -553
- data/ext/gsl_native/array_complex.c +0 -245
- data/ext/gsl_native/blas.c +0 -28
- data/ext/gsl_native/blas1.c +0 -733
- data/ext/gsl_native/blas2.c +0 -1088
- data/ext/gsl_native/blas3.c +0 -880
- data/ext/gsl_native/block.c +0 -40
- data/ext/gsl_native/block_source.h +0 -885
- data/ext/gsl_native/bspline.c +0 -122
- data/ext/gsl_native/bundle.c +0 -3
- data/ext/gsl_native/cdf.c +0 -740
- data/ext/gsl_native/cheb.c +0 -531
- data/ext/gsl_native/combination.c +0 -275
- data/ext/gsl_native/common.c +0 -293
- data/ext/gsl_native/complex.c +0 -1002
- data/ext/gsl_native/const.c +0 -331
- data/ext/gsl_native/const_additional.c +0 -99
- data/ext/gsl_native/cqp.c +0 -283
- data/ext/gsl_native/deriv.c +0 -187
- data/ext/gsl_native/dht.c +0 -353
- data/ext/gsl_native/diff.c +0 -164
- data/ext/gsl_native/dirac.c +0 -388
- data/ext/gsl_native/eigen.c +0 -2322
- data/ext/gsl_native/error.c +0 -193
- data/ext/gsl_native/extconf.rb +0 -118
- data/ext/gsl_native/fft.c +0 -1095
- data/ext/gsl_native/fit.c +0 -204
- data/ext/gsl_native/fresnel.c +0 -312
- data/ext/gsl_native/function.c +0 -518
- data/ext/gsl_native/geometry.c +0 -139
- data/ext/gsl_native/graph.c +0 -1590
- data/ext/gsl_native/gsl.c +0 -259
- data/ext/gsl_native/gsl_narray.c +0 -794
- data/ext/gsl_native/histogram.c +0 -1964
- data/ext/gsl_native/histogram2d.c +0 -1042
- data/ext/gsl_native/histogram3d.c +0 -884
- data/ext/gsl_native/histogram3d_source.c +0 -749
- data/ext/gsl_native/histogram_find.c +0 -99
- data/ext/gsl_native/histogram_oper.c +0 -150
- data/ext/gsl_native/ieee.c +0 -88
- data/ext/gsl_native/include/rb_gsl.h +0 -136
- data/ext/gsl_native/include/rb_gsl_array.h +0 -214
- data/ext/gsl_native/include/rb_gsl_cheb.h +0 -19
- data/ext/gsl_native/include/rb_gsl_common.h +0 -348
- data/ext/gsl_native/include/rb_gsl_complex.h +0 -25
- data/ext/gsl_native/include/rb_gsl_const.h +0 -23
- data/ext/gsl_native/include/rb_gsl_dirac.h +0 -6
- data/ext/gsl_native/include/rb_gsl_eigen.h +0 -17
- data/ext/gsl_native/include/rb_gsl_fft.h +0 -49
- data/ext/gsl_native/include/rb_gsl_fit.h +0 -23
- data/ext/gsl_native/include/rb_gsl_function.h +0 -23
- data/ext/gsl_native/include/rb_gsl_graph.h +0 -68
- data/ext/gsl_native/include/rb_gsl_histogram.h +0 -63
- data/ext/gsl_native/include/rb_gsl_histogram3d.h +0 -97
- data/ext/gsl_native/include/rb_gsl_integration.h +0 -17
- data/ext/gsl_native/include/rb_gsl_interp.h +0 -41
- data/ext/gsl_native/include/rb_gsl_linalg.h +0 -21
- data/ext/gsl_native/include/rb_gsl_math.h +0 -20
- data/ext/gsl_native/include/rb_gsl_odeiv.h +0 -18
- data/ext/gsl_native/include/rb_gsl_poly.h +0 -67
- data/ext/gsl_native/include/rb_gsl_rational.h +0 -30
- data/ext/gsl_native/include/rb_gsl_rng.h +0 -20
- data/ext/gsl_native/include/rb_gsl_root.h +0 -22
- data/ext/gsl_native/include/rb_gsl_sf.h +0 -110
- data/ext/gsl_native/include/rb_gsl_statistics.h +0 -17
- data/ext/gsl_native/include/rb_gsl_tensor.h +0 -43
- data/ext/gsl_native/include/rb_gsl_with_narray.h +0 -31
- data/ext/gsl_native/include/templates_off.h +0 -87
- data/ext/gsl_native/include/templates_on.h +0 -241
- data/ext/gsl_native/integration.c +0 -1154
- data/ext/gsl_native/interp.c +0 -499
- data/ext/gsl_native/jacobi.c +0 -733
- data/ext/gsl_native/linalg.c +0 -3915
- data/ext/gsl_native/linalg_complex.c +0 -726
- data/ext/gsl_native/math.c +0 -706
- data/ext/gsl_native/matrix.c +0 -36
- data/ext/gsl_native/matrix_complex.c +0 -1733
- data/ext/gsl_native/matrix_double.c +0 -557
- data/ext/gsl_native/matrix_int.c +0 -255
- data/ext/gsl_native/matrix_source.h +0 -2708
- data/ext/gsl_native/min.c +0 -219
- data/ext/gsl_native/monte.c +0 -978
- data/ext/gsl_native/multifit.c +0 -1862
- data/ext/gsl_native/multimin.c +0 -778
- data/ext/gsl_native/multimin_fsdf.c +0 -156
- data/ext/gsl_native/multiroots.c +0 -952
- data/ext/gsl_native/multiset.c +0 -210
- data/ext/gsl_native/ndlinear.c +0 -320
- data/ext/gsl_native/nmf.c +0 -171
- data/ext/gsl_native/nmf_wrap.c +0 -75
- data/ext/gsl_native/ntuple.c +0 -469
- data/ext/gsl_native/odeiv.c +0 -947
- data/ext/gsl_native/ool.c +0 -879
- data/ext/gsl_native/permutation.c +0 -598
- data/ext/gsl_native/poly.c +0 -39
- data/ext/gsl_native/poly2.c +0 -265
- data/ext/gsl_native/poly_source.h +0 -1871
- data/ext/gsl_native/qrng.c +0 -160
- data/ext/gsl_native/randist.c +0 -1848
- data/ext/gsl_native/rational.c +0 -480
- data/ext/gsl_native/rng.c +0 -595
- data/ext/gsl_native/root.c +0 -407
- data/ext/gsl_native/sf.c +0 -1446
- data/ext/gsl_native/sf_airy.c +0 -200
- data/ext/gsl_native/sf_bessel.c +0 -871
- data/ext/gsl_native/sf_clausen.c +0 -28
- data/ext/gsl_native/sf_coulomb.c +0 -206
- data/ext/gsl_native/sf_coupling.c +0 -121
- data/ext/gsl_native/sf_dawson.c +0 -29
- data/ext/gsl_native/sf_debye.c +0 -148
- data/ext/gsl_native/sf_dilog.c +0 -43
- data/ext/gsl_native/sf_elementary.c +0 -46
- data/ext/gsl_native/sf_ellint.c +0 -206
- data/ext/gsl_native/sf_elljac.c +0 -30
- data/ext/gsl_native/sf_erfc.c +0 -89
- data/ext/gsl_native/sf_exp.c +0 -169
- data/ext/gsl_native/sf_expint.c +0 -201
- data/ext/gsl_native/sf_fermi_dirac.c +0 -148
- data/ext/gsl_native/sf_gamma.c +0 -343
- data/ext/gsl_native/sf_gegenbauer.c +0 -97
- data/ext/gsl_native/sf_hyperg.c +0 -203
- data/ext/gsl_native/sf_laguerre.c +0 -113
- data/ext/gsl_native/sf_lambert.c +0 -47
- data/ext/gsl_native/sf_legendre.c +0 -368
- data/ext/gsl_native/sf_log.c +0 -105
- data/ext/gsl_native/sf_mathieu.c +0 -235
- data/ext/gsl_native/sf_power.c +0 -47
- data/ext/gsl_native/sf_psi.c +0 -89
- data/ext/gsl_native/sf_synchrotron.c +0 -48
- data/ext/gsl_native/sf_transport.c +0 -76
- data/ext/gsl_native/sf_trigonometric.c +0 -210
- data/ext/gsl_native/sf_zeta.c +0 -115
- data/ext/gsl_native/signal.c +0 -303
- data/ext/gsl_native/siman.c +0 -713
- data/ext/gsl_native/sort.c +0 -207
- data/ext/gsl_native/spline.c +0 -377
- data/ext/gsl_native/stats.c +0 -787
- data/ext/gsl_native/sum.c +0 -168
- data/ext/gsl_native/tamu_anova.c +0 -56
- data/ext/gsl_native/tensor.c +0 -35
- data/ext/gsl_native/tensor_source.h +0 -1122
- data/ext/gsl_native/vector.c +0 -35
- data/ext/gsl_native/vector_complex.c +0 -2241
- data/ext/gsl_native/vector_double.c +0 -1433
- data/ext/gsl_native/vector_int.c +0 -202
- data/ext/gsl_native/vector_source.h +0 -3321
- data/ext/gsl_native/wavelet.c +0 -923
- data/lib/gsl.rb +0 -8
- data/lib/gsl/gnuplot.rb +0 -41
- data/lib/gsl/oper.rb +0 -43
- data/lib/gsl/version.rb +0 -3
- data/lib/ool.rb +0 -22
- data/lib/ool/conmin.rb +0 -30
- data/lib/rbgsl.rb +0 -1
- data/rdoc/alf.rdoc +0 -77
- data/rdoc/blas.rdoc +0 -269
- data/rdoc/bspline.rdoc +0 -42
- data/rdoc/changes.rdoc +0 -159
- data/rdoc/cheb.rdoc +0 -99
- data/rdoc/cholesky_complex.rdoc +0 -46
- data/rdoc/combi.rdoc +0 -125
- data/rdoc/complex.rdoc +0 -210
- data/rdoc/const.rdoc +0 -546
- data/rdoc/dht.rdoc +0 -122
- data/rdoc/diff.rdoc +0 -133
- data/rdoc/ehandling.rdoc +0 -50
- data/rdoc/eigen.rdoc +0 -401
- data/rdoc/fft.rdoc +0 -535
- data/rdoc/fit.rdoc +0 -284
- data/rdoc/function.rdoc +0 -94
- data/rdoc/graph.rdoc +0 -137
- data/rdoc/hist.rdoc +0 -409
- data/rdoc/hist2d.rdoc +0 -279
- data/rdoc/hist3d.rdoc +0 -112
- data/rdoc/integration.rdoc +0 -398
- data/rdoc/interp.rdoc +0 -231
- data/rdoc/intro.rdoc +0 -27
- data/rdoc/linalg.rdoc +0 -681
- data/rdoc/linalg_complex.rdoc +0 -88
- data/rdoc/math.rdoc +0 -276
- data/rdoc/matrix.rdoc +0 -1093
- data/rdoc/min.rdoc +0 -189
- data/rdoc/monte.rdoc +0 -234
- data/rdoc/multimin.rdoc +0 -312
- data/rdoc/multiroot.rdoc +0 -293
- data/rdoc/narray.rdoc +0 -177
- data/rdoc/ndlinear.rdoc +0 -250
- data/rdoc/nonlinearfit.rdoc +0 -348
- data/rdoc/ntuple.rdoc +0 -88
- data/rdoc/odeiv.rdoc +0 -378
- data/rdoc/perm.rdoc +0 -221
- data/rdoc/poly.rdoc +0 -335
- data/rdoc/qrng.rdoc +0 -90
- data/rdoc/randist.rdoc +0 -233
- data/rdoc/ref.rdoc +0 -93
- data/rdoc/rng.rdoc +0 -203
- data/rdoc/roots.rdoc +0 -305
- data/rdoc/sf.rdoc +0 -1622
- data/rdoc/siman.rdoc +0 -89
- data/rdoc/sort.rdoc +0 -94
- data/rdoc/start.rdoc +0 -16
- data/rdoc/stats.rdoc +0 -219
- data/rdoc/sum.rdoc +0 -65
- data/rdoc/tensor.rdoc +0 -251
- data/rdoc/tut.rdoc +0 -5
- data/rdoc/use.rdoc +0 -177
- data/rdoc/vector.rdoc +0 -1243
- data/rdoc/vector_complex.rdoc +0 -347
- data/rdoc/wavelet.rdoc +0 -218
- data/test/gsl/blas_test.rb +0 -79
- data/test/gsl/bspline_test.rb +0 -63
- data/test/gsl/cdf_test.rb +0 -1512
- data/test/gsl/cheb_test.rb +0 -80
- data/test/gsl/combination_test.rb +0 -100
- data/test/gsl/complex_test.rb +0 -20
- data/test/gsl/const_test.rb +0 -29
- data/test/gsl/deriv_test.rb +0 -62
- data/test/gsl/dht_test.rb +0 -79
- data/test/gsl/diff_test.rb +0 -53
- data/test/gsl/eigen_test.rb +0 -563
- data/test/gsl/err_test.rb +0 -23
- data/test/gsl/fit_test.rb +0 -101
- data/test/gsl/histo_test.rb +0 -14
- data/test/gsl/index_test.rb +0 -61
- data/test/gsl/integration_test.rb +0 -274
- data/test/gsl/interp_test.rb +0 -27
- data/test/gsl/linalg_test.rb +0 -463
- data/test/gsl/matrix_nmf_test.rb +0 -37
- data/test/gsl/matrix_test.rb +0 -98
- data/test/gsl/min_test.rb +0 -89
- data/test/gsl/monte_test.rb +0 -77
- data/test/gsl/multifit_test.rb +0 -753
- data/test/gsl/multimin_test.rb +0 -157
- data/test/gsl/multiroot_test.rb +0 -135
- data/test/gsl/multiset_test.rb +0 -52
- data/test/gsl/odeiv_test.rb +0 -275
- data/test/gsl/oper_test.rb +0 -98
- data/test/gsl/poly_test.rb +0 -338
- data/test/gsl/qrng_test.rb +0 -94
- data/test/gsl/quartic_test.rb +0 -28
- data/test/gsl/randist_test.rb +0 -122
- data/test/gsl/rng_test.rb +0 -303
- data/test/gsl/roots_test.rb +0 -78
- data/test/gsl/sf_test.rb +0 -2079
- data/test/gsl/stats_test.rb +0 -122
- data/test/gsl/sum_test.rb +0 -69
- data/test/gsl/tensor_test.rb +0 -396
- data/test/gsl/vector_test.rb +0 -223
- data/test/gsl/wavelet_test.rb +0 -130
- data/test/gsl_test.rb +0 -321
- data/test/test_helper.rb +0 -42
- data/uncrustify.cfg +0 -1693
data/rdoc/hist.rdoc
DELETED
@@ -1,409 +0,0 @@
|
|
1
|
-
#
|
2
|
-
# = Histograms
|
3
|
-
# 1. {Histogram allocation}[link:rdoc/hist_rdoc.html#label-Histogram+allocation]
|
4
|
-
# 1. {Copying histograms}[link:rdoc/hist_rdoc.html#label-Copying+Histograms]
|
5
|
-
# 1. {Updating and accessing histogram elements}[link:rdoc/hist_rdoc.html#label-Updating+and+accessing+histogram+elements]
|
6
|
-
# 1. {Searching histogram ranges}[link:rdoc/hist_rdoc.html#label-Searching+histogram+ranges]
|
7
|
-
# 1. {Histogram Statistics}[link:rdoc/hist_rdoc.html#label-Histogram+Statistics]
|
8
|
-
# 1. {Histogram Operations}[link:rdoc/hist_rdoc.html#label-Histogram+Operations]
|
9
|
-
# 1. {Reading and writing histograms}[link:rdoc/hist_rdoc.html#label-Reading+and+writing+histograms]
|
10
|
-
# 1. {Extensions}[link:rdoc/hist_rdoc.html#label-Extentions]
|
11
|
-
# 1. {Histogram Operations}[link:rdoc/hist_rdoc.html#label-Histogram+operations]
|
12
|
-
# 1. {Graph interface}[link:rdoc/hist_rdoc.html#label-Graphics]
|
13
|
-
# 1. {Histogram Fittings}[link:rdoc/hist_rdoc.html#label-Fitting]
|
14
|
-
# 1. {The histogram probability distribution}[link:rdoc/hist_rdoc.html#label-The+histogram+probability+distribution]
|
15
|
-
#
|
16
|
-
# == Histogram allocation
|
17
|
-
# ---
|
18
|
-
# * GSL::Histogram.alloc(n)
|
19
|
-
# * GSL::Histogram.alloc(n, [xmin, xmax])
|
20
|
-
# * GSL::Histogram.alloc(n, xmin, xmax)
|
21
|
-
# * GSL::Histogram.alloc(n)
|
22
|
-
# * GSL::Histogram.alloc(array)
|
23
|
-
# * GSL::Histogram.alloc(vector)
|
24
|
-
#
|
25
|
-
# Constructor for a histogram object with <tt>n</tt> bins.
|
26
|
-
#
|
27
|
-
# Examples:
|
28
|
-
#
|
29
|
-
# 1. With an integer:
|
30
|
-
# h = Histogram.alloc(4) <--- Histogram of 4 bins.
|
31
|
-
# The range is not defined yet.
|
32
|
-
#
|
33
|
-
# [ bin[0] )[ bin[1] )[ bin[2] )[ bin[3] )
|
34
|
-
# |---------|---------|---------|---------|
|
35
|
-
# range[0] range[1] range[2] range[3] range[4]
|
36
|
-
#
|
37
|
-
# 1. With an array or a vector:
|
38
|
-
# h = Histogram.alloc([1, 3, 7, 9, 20]) <--- Histogram of 4 bins.
|
39
|
-
# The range is initialized as
|
40
|
-
# range[0] = 1, range[1] = 3, ..., range[4] = 20.
|
41
|
-
#
|
42
|
-
# 1. With size and the range [min, max]:
|
43
|
-
#
|
44
|
-
# >> h = Histogram.alloc(5, [0, 5])
|
45
|
-
# >> h.range
|
46
|
-
# => GSL::Histogram::Range:
|
47
|
-
# [ 0.000e+00 1.000e+00 2.000e+00 3.000e+00 4.000e+00 5.000e+00 ]
|
48
|
-
# >> h.bin
|
49
|
-
# => GSL::Histogram::Bin:
|
50
|
-
# [ 0.000e+00 0.000e+00 0.000e+00 0.000e+00 0.000e+00 ]
|
51
|
-
# >> h.increment(2.5)
|
52
|
-
# >> h.bin
|
53
|
-
# => GSL::Histogram::Bin:
|
54
|
-
# [ 0.000e+00 0.000e+00 1.000e+00 0.000e+00 0.000e+00 ]
|
55
|
-
#
|
56
|
-
# ---
|
57
|
-
# * GSL::Histogram.alloc_uniform(n, min, max)
|
58
|
-
# * GSL::Histogram.alloc_uniform(n, [min, max])
|
59
|
-
# * GSL::Histogram.equal_bins_p(h1, h2)
|
60
|
-
# * GSL::Histogram.equal_bins(h1, h2)
|
61
|
-
#
|
62
|
-
# Return 1 if the all of the individual bin ranges of the two histograms
|
63
|
-
# are identical, and 0 otherwise.
|
64
|
-
# ---
|
65
|
-
# * GSL::Histogram.equal_bins_p?(h1, h2)
|
66
|
-
# * GSL::Histogram.equal_bins?(h1, h2)
|
67
|
-
#
|
68
|
-
# Return <tt>true</tt> if the all of the individual bin ranges of the two histograms
|
69
|
-
# are identical, and <tt>false</tt> otherwise.
|
70
|
-
#
|
71
|
-
# ---
|
72
|
-
# * GSL::Histogram#set_ranges(v)
|
73
|
-
#
|
74
|
-
# This sets the ranges of the existing histogram using a {GSL::Vector}[link:rdoc/vector_rdoc.html] object.
|
75
|
-
# ---
|
76
|
-
# * GSL::Histogram#set_ranges_uniform(xmin, xmax)
|
77
|
-
# * GSL::Histogram#set_ranges_uniform([xmin, xmax])
|
78
|
-
#
|
79
|
-
# This method sets the ranges of the existing histogram <tt>self</tt>
|
80
|
-
# to cover the range <tt>xmin</tt> to <tt>xmax</tt> uniformly.
|
81
|
-
# The values of the histogram bins are reset to zero.
|
82
|
-
# The bin ranges are shown as below,
|
83
|
-
# bin[0] corresponds to xmin <= x < xmin + d
|
84
|
-
# bin[1] corresponds to xmin + d <= x < xmin + 2 d
|
85
|
-
# ......
|
86
|
-
# bin[n-1] corresponds to xmin + (n-1)d <= x < xmax
|
87
|
-
# where d is the bin spacing, d = (xmax-xmin)/n.
|
88
|
-
#
|
89
|
-
# == Copying Histograms
|
90
|
-
# ---
|
91
|
-
# * GSL::Histogram.memcpy(dest, src)
|
92
|
-
#
|
93
|
-
# Copies the histogram <tt>src</tt> into the pre-existing histogram <tt>dest</tt>,
|
94
|
-
# making dest into an exact copy of <tt>src</tt>.
|
95
|
-
# The two histograms must be of the same size.
|
96
|
-
# ---
|
97
|
-
# * GSL::Histogram#clone
|
98
|
-
#
|
99
|
-
# Returns a newly created histogram which is an exact copy of the histogram
|
100
|
-
# <tt>self</tt>.
|
101
|
-
#
|
102
|
-
# == Updating and accessing histogram elements
|
103
|
-
# ---
|
104
|
-
# * GSL::Histogram#increment(x, weight = 1)
|
105
|
-
# * GSL::Histogram#fill(x, weight = 1)
|
106
|
-
# * GSL::Histogram#accumulate(x, weight = 1)
|
107
|
-
#
|
108
|
-
# These methods updates the histogram <tt>self</tt> by adding <tt>weight</tt>
|
109
|
-
# (default = 1) to the bin whose range contains the coordinate <tt>x</tt>.
|
110
|
-
# If <tt>x</tt> is an instance of <tt>GSL::Vector</tt> or <tt>Array</tt>,
|
111
|
-
# all the elements are filled into the histogram.
|
112
|
-
# If <tt>x</tt> is less than (greater than) the lower limit (upper limit)
|
113
|
-
# of the histogram then none of bins are modified.
|
114
|
-
#
|
115
|
-
# ---
|
116
|
-
# * GSL::Histogram#increment2(x, weight = 1)
|
117
|
-
# * GSL::Histogram#fill2(x, weight = 1)
|
118
|
-
# * GSL::Histogram#accumulate2(x, weight = 1)
|
119
|
-
#
|
120
|
-
# These methods updates the histogram <tt>self</tt> by adding <tt>weight</tt>
|
121
|
-
# to the bin whose range contains the coordinate <tt>x</tt>.
|
122
|
-
# If <tt>x</tt> is less than the lower limit, the lowest bin is incremented.
|
123
|
-
# If <tt>x</tt> is greater than the upper limit, the highest bin is incremented.
|
124
|
-
#
|
125
|
-
# ---
|
126
|
-
# * GSL::Histogram#get(i)
|
127
|
-
# * \GSL::Histogram#[i]
|
128
|
-
#
|
129
|
-
# These methods return the contents of the <tt>i</tt>-th bin of the histogram
|
130
|
-
# <tt>self</tt>.
|
131
|
-
#
|
132
|
-
# ---
|
133
|
-
# * GSL::Hiatogram#get_range(i)
|
134
|
-
#
|
135
|
-
# This method finds the upper and lower range limits of the <tt>i</tt>-th bin
|
136
|
-
# of the histogram <tt>self</tt>, and returns an array [<tt>lower, upper</tt>].
|
137
|
-
#
|
138
|
-
# ---
|
139
|
-
# * GSL::Histogram#range
|
140
|
-
#
|
141
|
-
# This returns a <tt>Vector::View</tt> object as a reference to the pointer
|
142
|
-
# <tt>double *range</tt> in the <tt>gsl_histogram</tt> struct.
|
143
|
-
#
|
144
|
-
# ---
|
145
|
-
# * GSL::Histogram#bin
|
146
|
-
#
|
147
|
-
# This returns a <tt>Vector::View</tt> object to access the pointer <tt>double *bin</tt> in the <tt>gsl_histogram</tt> struct.
|
148
|
-
#
|
149
|
-
# ---
|
150
|
-
# * GSL::Histogram#max
|
151
|
-
# * GSL::Histogram#min
|
152
|
-
# * GSL::Histogram#bins
|
153
|
-
#
|
154
|
-
# These methods return the maximum upper and minimum lower range
|
155
|
-
# limits and the number of bins of the histogram <tt>self</tt>.
|
156
|
-
#
|
157
|
-
# ---
|
158
|
-
# * GSL::Histogram#reset
|
159
|
-
#
|
160
|
-
# This method resets all the bins in the histogram <tt>self</tt> to zero.
|
161
|
-
#
|
162
|
-
# == Searching histogram ranges
|
163
|
-
# ---
|
164
|
-
# * GSL::Histogram#find(x)
|
165
|
-
#
|
166
|
-
# This method finds and sets the index i to the bin number which
|
167
|
-
# covers the coordinate <tt>x</tt> in the histogram <tt>self</tt>.
|
168
|
-
#
|
169
|
-
# == Histogram Statistics
|
170
|
-
# ---
|
171
|
-
# * GSL::Histogram#max_val
|
172
|
-
#
|
173
|
-
# This returns the maximum value contained in the histogram bins.
|
174
|
-
#
|
175
|
-
# ---
|
176
|
-
# * GSL::Histogram#max_bin
|
177
|
-
#
|
178
|
-
# This returns the index of the bin containing the maximum value.
|
179
|
-
# In the case where several bins contain the same maximum value the
|
180
|
-
# smallest index is returned.
|
181
|
-
#
|
182
|
-
# ---
|
183
|
-
# * GSL::Histogram#min_val
|
184
|
-
#
|
185
|
-
# This returns the minimum value contained in the histogram bins.
|
186
|
-
#
|
187
|
-
# ---
|
188
|
-
# * GSL::Histogram#min_bin
|
189
|
-
#
|
190
|
-
# This returns the index of the bin containing the minimum value.
|
191
|
-
# In the case where several bins contain the same maximum value
|
192
|
-
# the smallest index is returned.
|
193
|
-
#
|
194
|
-
# ---
|
195
|
-
# * GSL::Histogram#mean
|
196
|
-
#
|
197
|
-
# This returns the mean of the histogrammed variable,
|
198
|
-
# where the histogram is regarded as a probability distribution.
|
199
|
-
# Negative bin values are ignored for the purposes of this calculation.
|
200
|
-
# The accuracy of the result is limited by the bin width.
|
201
|
-
#
|
202
|
-
# ---
|
203
|
-
# * GSL::Histogram#sigma
|
204
|
-
#
|
205
|
-
# This function returns the standard deviation of the histogrammed variable,
|
206
|
-
# where the histogram is regarded as a probability distribution.
|
207
|
-
# Negative bin values are ignored for the purposes of this calculation.
|
208
|
-
# The accuracy of the result is limited by the bin width.
|
209
|
-
#
|
210
|
-
# ---
|
211
|
-
# * GSL::Histogram#sum(istart = 0, iend = n-1)
|
212
|
-
#
|
213
|
-
# The sum of values of the histogram <tt>self</tt> from the <tt>istart</tt>-th bin
|
214
|
-
# to the <tt>iend</tt>-th bin.
|
215
|
-
#
|
216
|
-
#
|
217
|
-
# == Histogram Operations
|
218
|
-
#
|
219
|
-
# ---
|
220
|
-
# * GSL::Histogram#add(h2)
|
221
|
-
# * GSL::Histogram#sub(h2)
|
222
|
-
# * GSL::Histogram#mul(h2)
|
223
|
-
# * GSL::Histogram#div(h2)
|
224
|
-
# * GSL::Histogram#scale(val)
|
225
|
-
# * GSL::Histogram#shift(val)
|
226
|
-
#
|
227
|
-
#
|
228
|
-
#
|
229
|
-
# == Reading and writing histograms
|
230
|
-
# ---
|
231
|
-
# * GSL::Histogram#fwrite(io)
|
232
|
-
# * GSL::Histogram#fwrite(filename)
|
233
|
-
# * GSL::Histogram#fread(io)
|
234
|
-
# * GSL::Histogram#fread(filename)
|
235
|
-
# * GSL::Histogram#fprintf(io, range_format = "%e", bin_format = "%e")
|
236
|
-
# * GSL::Histogram#fprintf(filename, range_format = "%e", bin_format = "%e")
|
237
|
-
# * GSL::Histogram#fscanf(io)
|
238
|
-
# * GSL::Histogram#fscanf(filename)
|
239
|
-
#
|
240
|
-
#
|
241
|
-
# == Extentions
|
242
|
-
# === Histogram operations
|
243
|
-
# ---
|
244
|
-
# * GSL::Histogram#normalize
|
245
|
-
#
|
246
|
-
# This methods scales the contents of the bins of histogram <tt>self</tt> by its
|
247
|
-
# maximum value.
|
248
|
-
#
|
249
|
-
# ---
|
250
|
-
# * GSL::Histogram#rebin(m = 2)
|
251
|
-
#
|
252
|
-
# This method creates a new histogram merging <tt>m</tt> bins in one in the
|
253
|
-
# histogram <tt>self</tt>. This method cannot be used for histograms of
|
254
|
-
# non-uniform bin size. If <tt>m</tt> is not an exact divider of the number
|
255
|
-
# of bins of <tt>self</tt>, the range of the rebinned histogram is extended
|
256
|
-
# not to lose the entries in the last <tt>m-1</tt> (at most) bins.
|
257
|
-
#
|
258
|
-
# Example: a histogram <tt>h</tt> of size 5 with the range [0, 5), binned as
|
259
|
-
#
|
260
|
-
# GSL::Histogram::Range:
|
261
|
-
# [ 0.000e+00 1.000e+00 2.000e+00 3.000e+00 4.000e+00 5.000e+00 ]
|
262
|
-
# GSL::Histogram::Bin:
|
263
|
-
# [ 0.000e+00 3.000e+00 1.000e+00 1.000e+00 3.000e+00 ]
|
264
|
-
#
|
265
|
-
# When a new histogram is created merging two bins into one as
|
266
|
-
# <tt>h2 = h.rebin</tt>, then <tt>h2</tt> looks like
|
267
|
-
#
|
268
|
-
# GSL::Histogram::Range:
|
269
|
-
# [ 0.000e+00 2.000e+00 4.000e+00 6.000e+00 ]
|
270
|
-
# GSL::Histogram::Bin:
|
271
|
-
# [ 3.000e+00 2.000e+00 3.000e+00 ]
|
272
|
-
#
|
273
|
-
# ---
|
274
|
-
# * GSL::Histogram#reverse
|
275
|
-
#
|
276
|
-
# This method create a new histogram reversing the order of the range and the bin of
|
277
|
-
# histogram <tt>self</tt>.
|
278
|
-
#
|
279
|
-
# ---
|
280
|
-
# * GSL::Histogram#integrate(istart = 0, iend = n-1)
|
281
|
-
# * GSL::Histogram#integrate([istart, iend])
|
282
|
-
# * GSL::Histogram#integrate(direction = 1 or -1)
|
283
|
-
#
|
284
|
-
# This method calculates cumulative counts of the histogram <tt>self</tt>
|
285
|
-
# from the <tt>istart</tt>-th bin to the <tt>iend</tt>-th bin (<tt>iend</tt> inclusive),
|
286
|
-
# and returns a <tt>GSL::Histogram::Integral</tt>
|
287
|
-
# object. If <tt>istart <= iend</tt> (or <tt>direction == 1</tt>),
|
288
|
-
# the <tt>i</tt>-th bin value of a
|
289
|
-
# <tt>GSL::Histogram::Integral</tt> object <tt>hi</tt> created from a
|
290
|
-
# <tt>GSL::Histogram</tt> <tt>h</tt> is given by <tt>hi[i] = hi[i-1] + h[i]</tt>.
|
291
|
-
# If <tt>istart > iend</tt> (or <tt>direction == -1</tt>), <tt>hi[i] = hi[i+1] = h[i]</tt>.
|
292
|
-
#
|
293
|
-
# ---
|
294
|
-
# * GSL::Histogram::Integral#differentiate
|
295
|
-
# * GSL::Histogram::Integral#diff
|
296
|
-
#
|
297
|
-
# === Graphics
|
298
|
-
# ---
|
299
|
-
# * GSL::Histogram#graph(options)
|
300
|
-
#
|
301
|
-
# This method uses the GNU plotutils <tt>graph</tt> to draw the histogram <tt>self</tt>.
|
302
|
-
# The options as "-T X -C -l x" etc are given by a String.
|
303
|
-
#
|
304
|
-
# === Fitting
|
305
|
-
# ---
|
306
|
-
# * GSL::Histogram#fit_exponential(binstart = 0, binend = n-1)
|
307
|
-
#
|
308
|
-
# This method fits the histogram <tt>self</tt> to an exponential model
|
309
|
-
# <tt>h[n] = a exp(b x[n])</tt> using the bins of indices from <tt>binstart</tt>
|
310
|
-
# to <tt>binend</tt>. The result is returned as an Array of 6 elements,
|
311
|
-
# <tt>[a, b, erra, errb, sumsq, dof]</tt>, where
|
312
|
-
# * <tt>a</tt>: scale factor
|
313
|
-
# * <tt>b</tt>: exponent
|
314
|
-
# * <tt>erra, errb</tt>: fitting errors
|
315
|
-
# * <tt>sumsq</tt>: fitting chi-squared (not reduced chi-squared)
|
316
|
-
# * <tt>dof</tt>: degree-of-freedom, the number of bins used minus the number of parameters (2)
|
317
|
-
#
|
318
|
-
# ---
|
319
|
-
# * GSL::Histogram#fit_power(binstart = 0, binend = n-1)
|
320
|
-
#
|
321
|
-
# This method fits the histogram <tt>self</tt> to a power-law model
|
322
|
-
# <tt>h[n] = a x[n]^b</tt> using the bins of indices from <tt>binstart</tt>
|
323
|
-
# to <tt>binend</tt>. The result is returned as an Array of 6 elements,
|
324
|
-
# <tt>[a, b, erra, errb, sumsq, dof]</tt>.
|
325
|
-
#
|
326
|
-
# ---
|
327
|
-
# * GSL::Histogram#fit_gaussian(binstart = 0, binend = n-1)
|
328
|
-
#
|
329
|
-
# This method fits the histogram <tt>self</tt> to Gaussian distribution
|
330
|
-
# using the bins of indices from <tt>binstart</tt> to <tt>binend</tt>,
|
331
|
-
# and returns an Array of 8 elements,
|
332
|
-
# <tt>[sigma, mean, height, errsig, errmean, errhei, sumsq, dof]</tt>.
|
333
|
-
#
|
334
|
-
# Example:
|
335
|
-
# #!/usr/bin/env ruby
|
336
|
-
# require("gsl")
|
337
|
-
#
|
338
|
-
# N = 10000
|
339
|
-
# MAX = 8
|
340
|
-
# rng = Rng.alloc
|
341
|
-
#
|
342
|
-
# data = Ran.gaussian(rng, 1.5, N) + 2
|
343
|
-
# h = Histogram.alloc(100, [-MAX, MAX])
|
344
|
-
# h.increment(data)
|
345
|
-
#
|
346
|
-
# sigma, mean, height, = h.fit_gaussian
|
347
|
-
# x = Vector.linspace(-MAX, MAX, 100)
|
348
|
-
# y = height*Ran::gaussian_pdf(x-mean, sigma)
|
349
|
-
# GSL::graph(h, [x, y], "-T X -C -g 3")
|
350
|
-
#
|
351
|
-
# == The histogram probability distribution
|
352
|
-
# The probability distribution function for a histogram consists of a set of bins
|
353
|
-
# which measure the probability of an event falling into a given range of
|
354
|
-
# a continuous variable x. A probability distribution function is defined
|
355
|
-
# by the following class, which actually stores the cumulative probability
|
356
|
-
# distribution function. This is the natural quantity for generating samples
|
357
|
-
# via the inverse transform method, because there is a one-to-one mapping
|
358
|
-
# between the cumulative probability distribution and the range [0,1].
|
359
|
-
# It can be shown that by taking a uniform random number in this
|
360
|
-
# range and finding its corresponding coordinate in the cumulative probability
|
361
|
-
# distribution we obtain samples with the desired probability distribution.
|
362
|
-
#
|
363
|
-
# === Pdf class
|
364
|
-
# ---
|
365
|
-
# * GSL::Histogram::Pdf.alloc(n)
|
366
|
-
# * GSL::Histogram::Pdf.alloc(h)
|
367
|
-
#
|
368
|
-
# Constructors. If a histogram <tt>h</tt> is given,
|
369
|
-
# the probability distribution is initialized with the contents of <tt>h</tt>.
|
370
|
-
#
|
371
|
-
# ---
|
372
|
-
# * GSL::Histogram::Pdf#init(h)
|
373
|
-
#
|
374
|
-
# This initializes the probability distribution <tt>self</tt> with the contents
|
375
|
-
# of the histogram <tt>h</tt>.
|
376
|
-
#
|
377
|
-
# ---
|
378
|
-
# * GSL::Histogram::Pdf#sample(r)
|
379
|
-
#
|
380
|
-
# This method uses <tt>r</tt>, a uniform random number between zero and one,
|
381
|
-
# to compute a single random sample from the probability distribution <tt>self</tt>.
|
382
|
-
# The algorithm used to compute the sample s is given by the following formula,
|
383
|
-
# s = range[i] + delta * (range[i+1] - range[i])
|
384
|
-
# where i is the index which satisfies
|
385
|
-
# <tt>sum[i] <= r < sum[i+1]</tt> and <tt>delta</tt> is <tt>(r - sum[i])/(sum[i+1] - sum[i])</tt>.
|
386
|
-
#
|
387
|
-
# ---
|
388
|
-
# * GSL::Histogram::Pdf#n
|
389
|
-
#
|
390
|
-
# This returns the number of bins of the probability distribution function.
|
391
|
-
#
|
392
|
-
# ---
|
393
|
-
# * GSL::Histogram:Pdf#range
|
394
|
-
#
|
395
|
-
# This returns a <tt>Vector::View</tt> object as a reference to the pointer
|
396
|
-
# <tt>double *range</tt> in the <tt>gsl_histogram_pdf</tt> struct.
|
397
|
-
# ---
|
398
|
-
# * GSL::Histogram:Pdf#sum
|
399
|
-
#
|
400
|
-
# This returns a <tt>Vector::View</tt> object as a reference to the pointer
|
401
|
-
# <tt>double *sum</tt> in the <tt>gsl_histogram_pdf</tt> struct.
|
402
|
-
#
|
403
|
-
# {prev}[link:rdoc/stats_rdoc.html]
|
404
|
-
# {next}[link:rdoc/hist2d_rdoc.html]
|
405
|
-
#
|
406
|
-
# {Reference index}[link:rdoc/ref_rdoc.html]
|
407
|
-
# {top}[link:index.html]
|
408
|
-
#
|
409
|
-
#
|
data/rdoc/hist2d.rdoc
DELETED
@@ -1,279 +0,0 @@
|
|
1
|
-
#
|
2
|
-
# = GSL::Histogram2d class
|
3
|
-
#
|
4
|
-
# == Class methods
|
5
|
-
# ---
|
6
|
-
# * GSL::Histogram2d.alloc(nx, ny)
|
7
|
-
# * GSL::Histogram2d.alloc(xrange, yrange)
|
8
|
-
# * GSL::Histogram2d.alloc(nx, xmin, xmax, ny, ymin, ymax)
|
9
|
-
# * GSL::Histogram2d.alloc(nx, [xmin, xmax], ny, [ymin, ymax])
|
10
|
-
#
|
11
|
-
# Constructors.
|
12
|
-
#
|
13
|
-
# 1. With sizes
|
14
|
-
# >> h2 = GSL::Histogram2d.alloc(2, 3) # size 6
|
15
|
-
# >> h2.xrange
|
16
|
-
# => GSL::Histogram::Range:
|
17
|
-
# [ 0.000e+00 1.000e+00 2.000e+00 ]
|
18
|
-
# >> h2.yrange
|
19
|
-
# => GSL::Histogram::Range:
|
20
|
-
# [ 0.000e+00 1.000e+00 2.000e+00 3.000e+00 ]
|
21
|
-
# >> h2.bin
|
22
|
-
# => GSL::Histogram::Bin:
|
23
|
-
# [ 0.000e+00 0.000e+00 0.000e+00 0.000e+00 0.000e+00 0.000e+00 ]
|
24
|
-
# >> h2.xrange.class.superclass
|
25
|
-
# => GSL::Vector::View::ReadOnly
|
26
|
-
# >> h2.bin.class.superclass
|
27
|
-
# => GSL::Vector::View::ReadOnly
|
28
|
-
#
|
29
|
-
# 1. With ranges
|
30
|
-
# >> h2 = Histogram2d.alloc([3, 6, 7], [2, 3, 8])
|
31
|
-
# >> h2.xrange
|
32
|
-
# => GSL::Histogram::Range:
|
33
|
-
# [ 3.000e+00 6.000e+00 7.000e+00 ]
|
34
|
-
# >> h2.yrange
|
35
|
-
# => GSL::Histogram::Range:
|
36
|
-
# [ 2.000e+00 3.000e+00 8.000e+00 ]
|
37
|
-
# >> h2.bin
|
38
|
-
# => GSL::Histogram::Bin:
|
39
|
-
# [ 0.000e+00 0.000e+00 0.000e+00 0.000e+00 ]
|
40
|
-
#
|
41
|
-
# 1. With sizes and ranges
|
42
|
-
# >> h2 = Histogram2d.alloc(4, [0, 4], 3, [1, 5])
|
43
|
-
# >> h2.xrange
|
44
|
-
# => GSL::Histogram::Range:
|
45
|
-
# [ 0.000e+00 1.000e+00 2.000e+00 3.000e+00 4.000e+00 ]
|
46
|
-
# >> h2.yrange
|
47
|
-
# => GSL::Histogram::Range:
|
48
|
-
# [ 1.000e+00 2.333e+00 3.667e+00 5.000e+00 ]
|
49
|
-
#
|
50
|
-
# ---
|
51
|
-
# * GSL::Histogram2d.alloc_uniform(nx, xmin, xmax, ny, ymin, ymax)
|
52
|
-
# * GSL::Histogram2d.equal_bins_p(h1, h2)
|
53
|
-
#
|
54
|
-
#
|
55
|
-
# == Methods
|
56
|
-
# ---
|
57
|
-
# * GSL::Histogram2d#set_ranges(vx, vy)
|
58
|
-
#
|
59
|
-
# This sets the ranges of the existing histogram using
|
60
|
-
# two {GSL::Vector}[link:rdoc/vector_rdoc.html] objects or arrays.
|
61
|
-
# ---
|
62
|
-
# * GSL::Histogram2d#set_ranges_uniform(xmin, xmax, ymin, ymax)
|
63
|
-
# * GSL::Histogram2d#set_ranges_uniform([xmin, xmax], [ymin, ymax])
|
64
|
-
#
|
65
|
-
# This sets the ranges of the existing histogram <tt>self</tt> to cover
|
66
|
-
# the ranges <tt>xmin</tt> to <tt>xmax</tt> and <tt>ymin</tt> to <tt>ymax</tt> uniformly.
|
67
|
-
# The values of the histogram bins are reset to zero.
|
68
|
-
#
|
69
|
-
# ---
|
70
|
-
# * GSL::Histogram2d#clone
|
71
|
-
# * GSL::Histogram2d#duplicate
|
72
|
-
#
|
73
|
-
# Returns a newly created histogram which is an exact copy of the
|
74
|
-
# histogram <tt>self</tt>.
|
75
|
-
#
|
76
|
-
# ---
|
77
|
-
# * GSL::Histogram2d#increment(x, y, weight = 1)
|
78
|
-
# * GSL::Histogram2d#accumulate(x, y, weight = 1)
|
79
|
-
# * GSL::Histogram2d#fill(x, y, weight = 1)
|
80
|
-
#
|
81
|
-
# These method update the histogram <tt>self</tt> by adding <tt>weight</tt>
|
82
|
-
# to the bin whose <tt>x</tt> and <tt>y</tt> ranges contain the coordinates (x,y).
|
83
|
-
# If (x,y) lies outside the limits of the histogram then none of the
|
84
|
-
# bins are modified.
|
85
|
-
#
|
86
|
-
# ---
|
87
|
-
# * GSL::Histogram2d#increment2(x, y, weight = 1)
|
88
|
-
# * GSL::Histogram2d#accumulate2(x, y, weight = 1)
|
89
|
-
# * GSL::Histogram2d#fill2(x, y, weight = 1)
|
90
|
-
#
|
91
|
-
# These are similar to <tt>GSL::Histogram#increment/accumulate/fill</tt>,
|
92
|
-
# but when (x,y) lies outside the limits of the histogram then the edge
|
93
|
-
# bin is incremented.
|
94
|
-
#
|
95
|
-
# ---
|
96
|
-
# * GSL::Histogram2d#get(i[, j])
|
97
|
-
# * GSL::Histogram2d#get([i, j])
|
98
|
-
# * GSL::Histogram2d#[i, j]
|
99
|
-
#
|
100
|
-
# Return the contents of the <tt>(i,j)</tt>-th bin of the histogram <tt>self</tt>.
|
101
|
-
#
|
102
|
-
# ---
|
103
|
-
# * GSL::Histogram2d#get_xrange(i)
|
104
|
-
# * GSL::Histogram2d#get_yrange(j)
|
105
|
-
#
|
106
|
-
# These methods find the upper and lower range limits of the <tt>i</tt>-th
|
107
|
-
# and <tt>j</tt>-th bins in the x and y directions of the histogram <tt>self</tt>.
|
108
|
-
#
|
109
|
-
# Ex:
|
110
|
-
# >> h2 = Histogram2d.alloc(2, [0, 2], 3, [1, 4])
|
111
|
-
# >> h2.get_xrange(1)
|
112
|
-
# => [1.0, 2.0]
|
113
|
-
# >> h2.get_yrange(1)
|
114
|
-
# => [2.0, 3.0]
|
115
|
-
#
|
116
|
-
# ---
|
117
|
-
# * GSL::Histogram2d#xrange
|
118
|
-
# * GSL::Histogram2d#yrange
|
119
|
-
#
|
120
|
-
# These methods returns the range of histogram as <tt>Vector::View</tt> objects.
|
121
|
-
#
|
122
|
-
# ---
|
123
|
-
# * GSL::Histogram2d#xmax
|
124
|
-
# * GSL::Histogram2d#xmin
|
125
|
-
# * GSL::Histogram2d#ymax
|
126
|
-
# * GSL::Histogram2d#ymin
|
127
|
-
# * GSL::Histogram2d#nx
|
128
|
-
# * GSL::Histogram2d#ny
|
129
|
-
#
|
130
|
-
# These methodss return the maximum upper and minimum lower range limits
|
131
|
-
# and the number of bins for the x and y directions of the histogram <tt>self</tt>.
|
132
|
-
#
|
133
|
-
# ---
|
134
|
-
# * GSL::Histogram2d#reset
|
135
|
-
#
|
136
|
-
# Resets all the bins of the histogram <tt>self</tt> to zero.
|
137
|
-
#
|
138
|
-
# ---
|
139
|
-
# * GSL::Histogram2d#find(x, y)
|
140
|
-
#
|
141
|
-
# Finds the indices i and j to the to the bin which covers the
|
142
|
-
# coordinates (<tt>x,y</tt>).
|
143
|
-
#
|
144
|
-
# ---
|
145
|
-
# * GSL::Histogram2d#max_val
|
146
|
-
#
|
147
|
-
# Returns the maximum value contained in the histogram bins.
|
148
|
-
#
|
149
|
-
# ---
|
150
|
-
# * GSL::Histogram2d#max_bin
|
151
|
-
#
|
152
|
-
# Returns the indices (i,j) of the bin containing the maximum value in
|
153
|
-
# the histogram <tt>self</tt>.
|
154
|
-
#
|
155
|
-
# ---
|
156
|
-
# * GSL::Histogram2d#min_val
|
157
|
-
#
|
158
|
-
# Returns the minimum value contained in the histogram bins.
|
159
|
-
# ---
|
160
|
-
# * GSL::Histogram2d#min_bin
|
161
|
-
#
|
162
|
-
# Returns the indices (i,j) of the bin containing the minimum value
|
163
|
-
# in the histogram <tt>self</tt>.
|
164
|
-
#
|
165
|
-
# ---
|
166
|
-
# * GSL::Histogram2d#xmean
|
167
|
-
#
|
168
|
-
# Returns the mean of the histogrammed x variable, where the histogram
|
169
|
-
# is regarded as a probability distribution.
|
170
|
-
# Negative bin values are ignored for the purposes of this calculation.
|
171
|
-
#
|
172
|
-
# ---
|
173
|
-
# * GSL::Histogram2d#xsigma
|
174
|
-
#
|
175
|
-
# Returns the standard deviation of the histogrammed x variable,
|
176
|
-
# where the histogram is regarded as a probability distribution.
|
177
|
-
# Negative bin values are ignored for the purposes of this calculation.
|
178
|
-
#
|
179
|
-
# ---
|
180
|
-
# * GSL::Histogram2d#ymean
|
181
|
-
#
|
182
|
-
# Returns the mean of the histogrammed y variable, where the histogram
|
183
|
-
# is regarded as a probability distribution. Negative bin values are
|
184
|
-
# ignored for the purposes of this calculation.
|
185
|
-
#
|
186
|
-
# ---
|
187
|
-
# * GSL::Histogram2d#ysigma
|
188
|
-
#
|
189
|
-
# Returns the standard deviation of the histogrammed y variable,
|
190
|
-
# where the histogram is regarded as a probability distribution.
|
191
|
-
# Negative bin values are ignored for the purposes of this calculation.
|
192
|
-
#
|
193
|
-
# ---
|
194
|
-
# * GSL::Histogram2d#cov
|
195
|
-
#
|
196
|
-
# Returns the covariance of the histogrammed x and y variables,
|
197
|
-
# where the histogram is regarded as a probability distribution.
|
198
|
-
# Negative bin values are ignored for the purposes of this calculation.
|
199
|
-
#
|
200
|
-
# ---
|
201
|
-
# * GSL::Histogram2d#sum
|
202
|
-
# * GSL::Histogram2d#integral
|
203
|
-
#
|
204
|
-
# Return the sum of all bin values. Negative bin values are included in the sum.
|
205
|
-
# ---
|
206
|
-
# * GSL::Histogram2d#add
|
207
|
-
# * GSL::Histogram2d#sub
|
208
|
-
# * GSL::Histogram2d#mul
|
209
|
-
# * GSL::Histogram2d#div
|
210
|
-
# * GSL::Histogram2d#scale(val)
|
211
|
-
# * GSL::Histogram2d#shift(val)
|
212
|
-
# * GSL::Histogram2d#fwrite(io)
|
213
|
-
# * GSL::Histogram2d#fwrite(filename)
|
214
|
-
# * GSL::Histogram2d#fread(io)
|
215
|
-
# * GSL::Histogram2d#fread(filename)
|
216
|
-
# * GSL::Histogram2d#fprintf(io, range_format, bin_format)
|
217
|
-
# * GSL::Histogram2d#fprintf(filename, range_format, bin_format)
|
218
|
-
# * GSL::Histogram2d#fscanf(io)
|
219
|
-
# * GSL::Histogram2d#fscanf(filename)
|
220
|
-
#
|
221
|
-
#
|
222
|
-
# ---
|
223
|
-
# * GSL::Histogram2d#xproject(jstart = 0, jend = ny-1)
|
224
|
-
#
|
225
|
-
# Creates a <tt>GSL::Histogram</tt> object by projecting the 2D histogram
|
226
|
-
# <tt>self</tt> onto the x-axis over the y-range from <tt>jstart</tt> to <tt>jend</tt>.
|
227
|
-
#
|
228
|
-
# ---
|
229
|
-
# * GSL::Histogram2d#yproject(istart = 0, iend = nx-1)
|
230
|
-
#
|
231
|
-
# Creates a <tt>GSL::Histogram</tt> object by projecting the 2D histogram
|
232
|
-
# <tt>self</tt> onto the y-axis over the x-range from <tt>istart</tt> to <tt>iend</tt>.
|
233
|
-
#
|
234
|
-
# == Pdf class
|
235
|
-
# ---
|
236
|
-
# * GSL::Histogram2d::Pdf.alloc(nx, ny)
|
237
|
-
#
|
238
|
-
# Constructors
|
239
|
-
#
|
240
|
-
# ---
|
241
|
-
# * GSL::Histogram2d::Pdf#init(h)
|
242
|
-
#
|
243
|
-
# Initializer with a GSL::Histogram2d object <tt>h</tt>.
|
244
|
-
#
|
245
|
-
# ---
|
246
|
-
# * GSL::Histogram2d::Pdf#sample(r1, r2)
|
247
|
-
#
|
248
|
-
#
|
249
|
-
# == Example
|
250
|
-
# #!/usr/bin/env ruby
|
251
|
-
# require("gsl")
|
252
|
-
#
|
253
|
-
# N = 10000
|
254
|
-
# BINS = 100
|
255
|
-
#
|
256
|
-
# rng = Rng.new(2)
|
257
|
-
# h2 = Histogram2d.alloc(BINS, [-5, 5], BINS, [-8, 8])
|
258
|
-
#
|
259
|
-
# sig1 = 1.0
|
260
|
-
# sig2 = 2.0
|
261
|
-
#
|
262
|
-
# for i in 0...N do
|
263
|
-
# r1 = rng.gaussian(sig1)
|
264
|
-
# r2 = rng.gaussian(sig2)
|
265
|
-
# h2.increment(r1, r2)
|
266
|
-
# end
|
267
|
-
#
|
268
|
-
# hx = h2.xproject
|
269
|
-
# hy = h2.yproject
|
270
|
-
# printf("%f %f %f %f\n", h2.xsigma, h2.ysigma, hx.sigma, hy.sigma)
|
271
|
-
# GSL::graph(hx, hy, "-T X -C -g 3")
|
272
|
-
#
|
273
|
-
#
|
274
|
-
# {prev}[link:rdoc/hist_rdoc.html]
|
275
|
-
# {next}[link:rdoc/hist3d_rdoc.html]
|
276
|
-
#
|
277
|
-
# {Reference index}[link:rdoc/ref_rdoc.html]
|
278
|
-
# {top}[link:index.html]
|
279
|
-
#
|