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/qrng.rdoc
DELETED
@@ -1,90 +0,0 @@
|
|
1
|
-
#
|
2
|
-
# = Quasi-Random Sequences
|
3
|
-
#
|
4
|
-
# This chapter describes the quasi-random sequence generator <tt>GSL::QRng</tt>
|
5
|
-
# of arbitrary dimensions. A quasi-random sequence progressively covers
|
6
|
-
# a d-dimensional space with a set of points that are uniformly distributed.
|
7
|
-
# Quasi-random sequences are also known as low-discrepancy sequences.
|
8
|
-
# The quasi-random sequence generators use an interface that is similar
|
9
|
-
# to the interface for random number generators.
|
10
|
-
#
|
11
|
-
# Contents:
|
12
|
-
# 1. {Quasi-random number generator initialization}[link:rdoc/qrng_rdoc.html#label-Quasi-random+number+generator+initialization]
|
13
|
-
# 1. {Sampling from a quasi-random number generator}[link:rdoc/qrng_rdoc.html#label-Sampling+from+a+quasi-random+number+generator]
|
14
|
-
# 1. {Auxiliary quasi-random number generator functions}[link:rdoc/qrng_rdoc.html#label-Auxiliary+quasi-random+number+generator+functions]
|
15
|
-
# 1. {Saving and resorting quasi-random number generator state}[link:rdoc/qrng_rdoc.html#label-Saving+and+resorting+quasi-random+number+generator+state]
|
16
|
-
# 1. {Quasi-random number generator algorithms}[link:rdoc/qrng_rdoc.html#label-Quasi-random+number+generator+algorithms]
|
17
|
-
#
|
18
|
-
# == Quasi-random number generator initialization
|
19
|
-
#
|
20
|
-
# ---
|
21
|
-
# * GSL::QRng.alloc(T, d)
|
22
|
-
#
|
23
|
-
# This returns a GSL::QRng object, a quasi-random sequence generator of type <tt>T</tt> and dimension <tt>d</tt>.
|
24
|
-
#
|
25
|
-
# ---
|
26
|
-
# * GSL::QRng::init
|
27
|
-
#
|
28
|
-
# This reinitializes the generator to its starting point.
|
29
|
-
#
|
30
|
-
# == Sampling from a quasi-random number generator
|
31
|
-
#
|
32
|
-
# ---
|
33
|
-
# * GSL::QRng::get(x)
|
34
|
-
#
|
35
|
-
# This calculate the next point <tt>x</tt> from the sequence generator. Here <tt>x</tt> is an instance of the {GSL::Vector}[link:rdoc/vector_rdoc.html] class. The space available for <tt>x</tt> must match the dimension of the generator. The point <tt>x</tt> will lie in the range 0 < x_i < 1 for each x_i.
|
36
|
-
#
|
37
|
-
# This is used as
|
38
|
-
# q = QRng.alloc(QRng::SOBOL, dim)
|
39
|
-
# v = Vector.alloc(dim)
|
40
|
-
# for i in 0..1024 do
|
41
|
-
# q.get(v)
|
42
|
-
# printf("%.5f %.5f\n", v[0], v[1])
|
43
|
-
# end
|
44
|
-
#
|
45
|
-
# == Auxiliary quasi-random number generator functions
|
46
|
-
#
|
47
|
-
# ---
|
48
|
-
# * GSL::QRng::name
|
49
|
-
#
|
50
|
-
# Returns the name of the generator <tt>self</tt>.
|
51
|
-
#
|
52
|
-
# ---
|
53
|
-
# * GSL::QRng::size
|
54
|
-
#
|
55
|
-
#
|
56
|
-
# == Saving and resorting quasi-random number generator state
|
57
|
-
# ---
|
58
|
-
# * GSL::QRng::clone
|
59
|
-
# * GSL::QRng::duplicate
|
60
|
-
#
|
61
|
-
# Return a newly created generator which is an exact copy of the generator <tt>self</tt>.
|
62
|
-
# == Quasi-random number generator algorithms
|
63
|
-
# In creating a generator by the method <tt>GSL::QRng.alloc(T, d)</tt>,
|
64
|
-
# the algorithm type <tt>T</tt> is given by a String or a Fixnum constant.
|
65
|
-
# The following quasi-random sequence algorithms are available,
|
66
|
-
#
|
67
|
-
# * "<tt>niederreiter_2</tt>" (String)
|
68
|
-
# * <tt>GSL::QRng::NIEDERREITER_2</tt> (Fixnum)
|
69
|
-
#
|
70
|
-
# The generator of this type uses the algorithm described in Bratley, Fox, Niederreiter, ACM Trans. Model. Comp. Sim. 2, 195 (1992). It is valid up to 12 dimensions.
|
71
|
-
#
|
72
|
-
# * "<tt>sobol</tt>" (String)
|
73
|
-
# * <tt>GSL::QRng::SOBOL</tt> (Fixnum)
|
74
|
-
#
|
75
|
-
# This generator uses the Sobol sequence described in Antonov, Saleev, USSR Comput. Maths. Math. Phys. 19, 252 (1980). It is valid up to 40 dimensions.
|
76
|
-
#
|
77
|
-
# * "<tt>halton</tt>" (String)
|
78
|
-
# * <tt>GSL::QRng::HALTON</tt> (Fixnum)
|
79
|
-
# * "<tt>reversehalton</tt>" (String)
|
80
|
-
# * <tt>GSL::QRng::REVERSEHALTON</tt> (Fixnum)
|
81
|
-
#
|
82
|
-
# (GSL-1.11 or later) These generators use the Halton and reverse Halton sequences described in J.H. Halton, Numerische Mathematik 2, 84-90 (1960) and B. Vandewoestyne and R. Cools Computational and Applied Mathematics 189, 1&2, 341-361 (2006). They are valid up to 1229 dimensions.
|
83
|
-
#
|
84
|
-
# {prev}[link:rdoc/rng_rdoc.html]
|
85
|
-
# {next}[link:rdoc/randist_rdoc.html]
|
86
|
-
#
|
87
|
-
# {Reference index}[link:rdoc/ref_rdoc.html]
|
88
|
-
# {top}[link:index.html]
|
89
|
-
#
|
90
|
-
#
|
data/rdoc/randist.rdoc
DELETED
@@ -1,233 +0,0 @@
|
|
1
|
-
#
|
2
|
-
# = Random Number Distributions
|
3
|
-
# This chapter describes functions for generating random variates and computing
|
4
|
-
# their probability distributions. Samples from the distributions described
|
5
|
-
# in this chapter can be obtained using any of the random number generators
|
6
|
-
# in the library as an underlying source of randomness.
|
7
|
-
#
|
8
|
-
# In the simplest cases a non-uniform distribution can be obtained analytically
|
9
|
-
# from the uniform distribution of a random number generator by applying an
|
10
|
-
# appropriate transformation. This method uses one call to the random number
|
11
|
-
# generator. More complicated distributions are created by the
|
12
|
-
# acceptance-rejection method, which compares the desired distribution against
|
13
|
-
# a distribution which is similar and known analytically. This usually requires
|
14
|
-
# several samples from the generator.
|
15
|
-
#
|
16
|
-
# The library also provides cumulative distribution functions and inverse
|
17
|
-
# cumulative distribution functions, sometimes referred to as quantile
|
18
|
-
# functions. The cumulative distribution functions and their inverses are
|
19
|
-
# computed separately for the upper and lower tails of the distribution,
|
20
|
-
# allowing full accuracy to be retained for small results.
|
21
|
-
#
|
22
|
-
# Contents:
|
23
|
-
# 1. {Introduction}[link:rdoc/randist_rdoc.html#label-Introduction]
|
24
|
-
# 1. {The Gaussian Distribution}[link:rdoc/randist_rdoc.html#label-The+Gaussian+Distribution]
|
25
|
-
# 1. {The Gaussian Tail Distribution}[link:rdoc/randist_rdoc.html#label-The+Gaussian+Tail+Distribution]
|
26
|
-
# ...
|
27
|
-
# and more, see {the GSL reference}[https://gnu.org/software/gsl/manual/]
|
28
|
-
# 1. {Shuffling and Sampling}[link:rdoc/randist_rdoc.html#label-Shuffling+and+Sampling]
|
29
|
-
#
|
30
|
-
# == Introduction
|
31
|
-
# Continuous random number distributions are defined by a probability density
|
32
|
-
# function, p(x), such that the probability of x occurring in the
|
33
|
-
# infinitesimal range x to x+dx is p dx.
|
34
|
-
#
|
35
|
-
# The cumulative distribution function for the lower tail P(x) is defined by the
|
36
|
-
# integral, and gives the probability of a variate taking a value less than x.
|
37
|
-
#
|
38
|
-
# The cumulative distribution function for the upper tail Q(x) is defined by the
|
39
|
-
# integral, and gives the probability of a variate taking a value greater than
|
40
|
-
# x.
|
41
|
-
#
|
42
|
-
# The upper and lower cumulative distribution functions are related
|
43
|
-
# by P(x) + Q(x) = 1 and satisfy 0 <= P(x) <= 1, 0 <= Q(x) <= 1.
|
44
|
-
#
|
45
|
-
# The inverse cumulative distributions, x=P^{-1}(P) and x=Q^{-1}(Q) give the
|
46
|
-
# values of x which correspond to a specific value of P or Q. They can be used
|
47
|
-
# to find confidence limits from probability values.
|
48
|
-
#
|
49
|
-
# For discrete distributions the probability of sampling the integer value k is
|
50
|
-
# given by p(k), where \sum_k p(k) = 1. The cumulative distribution for the
|
51
|
-
# lower tail P(k) of a discrete distribution is defined as, where the sum is
|
52
|
-
# over the allowed range of the distribution less than or equal to k.
|
53
|
-
#
|
54
|
-
# The cumulative distribution for the upper tail of a discrete distribution Q(k)
|
55
|
-
# is defined as giving the sum of probabilities for all values greater than k.
|
56
|
-
# These two definitions satisfy the identity P(k)+Q(k)=1.
|
57
|
-
#
|
58
|
-
# If the range of the distribution is 1 to n inclusive then P(n)=1, Q(n)=0
|
59
|
-
# while P(1) = p(1), Q(1)=1-p(1).
|
60
|
-
#
|
61
|
-
#
|
62
|
-
#
|
63
|
-
# == The Gaussian Distribution
|
64
|
-
# ---
|
65
|
-
# * GSL::Rng#gaussian(sigma = 1)
|
66
|
-
# * GSL::Ran::gaussian(rng, sigma = 1)
|
67
|
-
# * GSL::Rng#ugaussian
|
68
|
-
# * GSL::Ran::ugaussian
|
69
|
-
#
|
70
|
-
# These return a Gaussian random variate, with mean zero and standard
|
71
|
-
# deviation <tt>sigma</tt>.
|
72
|
-
#
|
73
|
-
# ---
|
74
|
-
# * GSL::Ran::gaussian_pdf(x, sigma = 1)
|
75
|
-
#
|
76
|
-
# Computes the probability density p(x) at <tt>x</tt> for a Gaussian distribution
|
77
|
-
# with standard deviation <tt>sigma</tt>.
|
78
|
-
#
|
79
|
-
# ---
|
80
|
-
# * GSL::Rng#gaussian_ratio_method(sigma = 1)
|
81
|
-
# * GSL::Ran::gaussian_ratio_method(rng, sigma = 1)
|
82
|
-
#
|
83
|
-
# Use Kinderman-Monahan ratio method.
|
84
|
-
#
|
85
|
-
# ---
|
86
|
-
# * GSL::Cdf::gaussian_P(x, sigma = 1)
|
87
|
-
# * GSL::Cdf::gaussian_Q(x, sigma = 1)
|
88
|
-
# * GSL::Cdf::gaussian_Pinv(P, sigma = 1)
|
89
|
-
# * GSL::Cdf::gaussian_Qinv(Q, sigma = 1)
|
90
|
-
# * GSL::Cdf::ugaussian_P(x)
|
91
|
-
# * GSL::Cdf::ugaussian_Q(x)
|
92
|
-
# * GSL::Cdf::ugaussian_Pinv(P)
|
93
|
-
# * GSL::Cdf::ugaussian_Qinv(Q)
|
94
|
-
#
|
95
|
-
# These methods compute the cumulative distribution functions P(x), Q(x)
|
96
|
-
# and their inverses for the Gaussian distribution with standard
|
97
|
-
# deviation <tt>sigma</tt>.
|
98
|
-
#
|
99
|
-
# == The Gaussian Tail Distribution
|
100
|
-
# ---
|
101
|
-
# * GSL::Rng#gaussian_tail(a, sigma = 1)
|
102
|
-
# * GSL::Ran#gaussian_tail(rng, a, sigma = 1)
|
103
|
-
# * GSL::Rng#ugaussian_tail(a)
|
104
|
-
# * GSL::Ran#ugaussian_tail(rng)
|
105
|
-
#
|
106
|
-
# These methods provide random variates from the upper tail of a Gaussian
|
107
|
-
# distribution with standard deviation <tt>sigma</tt>.
|
108
|
-
# The values returned are larger than the lower limit <tt>a</tt>, which must be positive.
|
109
|
-
#
|
110
|
-
# ---
|
111
|
-
# * GSL::Ran::gaussian_tail_pdf(x, a, sigma = 1)
|
112
|
-
# * GSL::Ran::ugaussian_tail_pdf(x, a)
|
113
|
-
#
|
114
|
-
# These methods compute the probability density p(x) at <tt>x</tt> for a Gaussian
|
115
|
-
# tail distribution with standard deviation <tt>sigma</tt>
|
116
|
-
# and lower limit <tt>a</tt>.
|
117
|
-
#
|
118
|
-
# == The Bivariate Gaussian Distribution
|
119
|
-
# ---
|
120
|
-
# * GSL::Rng#bivariate_gaussian(sigma_x, sigma_y, rho)
|
121
|
-
# * GSL::Ran::bivariate_gaussian(rng, sigma_x, sigma_y, rho)
|
122
|
-
#
|
123
|
-
# These methods generate a pair of correlated gaussian variates,
|
124
|
-
# with mean zero, correlation coefficient <tt>rho</tt> and standard deviations
|
125
|
-
# <tt>sigma_x</tt> and <tt>sigma_y</tt> in the x and y directions.
|
126
|
-
#
|
127
|
-
# ---
|
128
|
-
# * GSL::Ran::bivariate_gaussian_pdf(x, y, sigma_x, sigma_y, rho)
|
129
|
-
#
|
130
|
-
# This method computes the probability density p(x,y) at <tt>(x,y)</tt>
|
131
|
-
# for a bivariate gaussian distribution with standard deviations
|
132
|
-
# <tt>sigma_x, sigma_y</tt> and correlation coefficient <tt>rho</tt>.
|
133
|
-
#
|
134
|
-
# == The Exponential Distribution
|
135
|
-
# ---
|
136
|
-
# * GSL::Rng#exponential(mu)
|
137
|
-
# * GSL::Ran::exponential(rng, mu)
|
138
|
-
#
|
139
|
-
# These methods return a random variate from the exponential
|
140
|
-
# distribution with mean <tt>mu</tt>.
|
141
|
-
#
|
142
|
-
# ---
|
143
|
-
# * GSL::Ran::exponential_pdf(x, mu)
|
144
|
-
#
|
145
|
-
# This method computes the probability density p(x) at <tt>x</tt>
|
146
|
-
# for an exponential distribution with mean <tt>mu</tt>.
|
147
|
-
#
|
148
|
-
# ---
|
149
|
-
# * GSL::Cdf::exponential_P(x, mu)
|
150
|
-
# * GSL::Cdf::exponential_Q(x, mu)
|
151
|
-
# * GSL::Cdf::exponential_Pinv(P, mu)
|
152
|
-
# * GSL::Cdf::exponential_Qinv(Q, mu)
|
153
|
-
#
|
154
|
-
# These methods compute the cumulative distribution functions P(x), Q(x)
|
155
|
-
# and their inverses for the exponential distribution with mean <tt>mu</tt>.
|
156
|
-
#
|
157
|
-
# == The Laplace Distribution
|
158
|
-
# ---
|
159
|
-
# * GSL::Rng#laplace(a)
|
160
|
-
# * GSL::Ran::laplace(rng, a)
|
161
|
-
#
|
162
|
-
# These methods return a random variate from the Laplace distribution
|
163
|
-
# with width <tt>a</tt>.
|
164
|
-
#
|
165
|
-
# ---
|
166
|
-
# * GSL::Ran::laplace_pdf(x, a)
|
167
|
-
#
|
168
|
-
# This method computes the probability density p(x) at <tt>x</tt>
|
169
|
-
# for a Laplace distribution with width <tt>a</tt>.
|
170
|
-
#
|
171
|
-
# ---
|
172
|
-
# * GSL::Cdf::laplace_P(x, a)
|
173
|
-
# * GSL::Cdf::laplace_Q(x, a)
|
174
|
-
# * GSL::Cdf::laplace_Pinv(P, a)
|
175
|
-
# * GSL::Cdf::laplace_Qinv(Q, a)
|
176
|
-
#
|
177
|
-
# These methods compute the cumulative distribution functions P(x), Q(x)
|
178
|
-
# and their inverses for the Laplace distribution with width <tt>a</tt>.
|
179
|
-
#
|
180
|
-
# ---
|
181
|
-
# * GSL::Rng#exppow(a, b)
|
182
|
-
# * GSL::Rng#cauchy(a)
|
183
|
-
# * GSL::Rng#rayleigh(sigma)
|
184
|
-
# * GSL::Rng#rayleigh_tail(a, sigma)
|
185
|
-
# * GSL::Rng#landau()
|
186
|
-
# * GSL::Rng#levy(c, alpha)
|
187
|
-
# * GSL::Rng#levy_skew(c, alpha, beta)
|
188
|
-
# * GSL::Rng#gamma(a, b)
|
189
|
-
# * GSL::Rng#flat(a, b)
|
190
|
-
# * GSL::Rng#lognormal(zeta, sigma)
|
191
|
-
# * GSL::Rng#chisq(nu)
|
192
|
-
# * GSL::Rng#fdist(nu1, nu2)
|
193
|
-
# * GSL::Rng#tdist(nu)
|
194
|
-
# * GSL::Rng#beta(a, b)
|
195
|
-
# * GSL::Rng#logistic(a)
|
196
|
-
# * GSL::Rng#pareto(a, b)
|
197
|
-
#
|
198
|
-
#
|
199
|
-
# ...
|
200
|
-
#
|
201
|
-
# and more, see {the GSL reference}[https://gnu.org/software/gsl/manual/gsl-ref_19.html#SEC286].
|
202
|
-
#
|
203
|
-
# == Shuffling and Sampling
|
204
|
-
# ---
|
205
|
-
# * GSL::Rng#shuffle(v, n)
|
206
|
-
#
|
207
|
-
# This randomly shuffles the order of <tt>n</tt> objects,
|
208
|
-
# stored in the {GSL::Vector}[link:rdoc/vector_rdoc.html] object <tt>v</tt>.
|
209
|
-
# ---
|
210
|
-
# * GSL::Rng#choose(v, k)
|
211
|
-
#
|
212
|
-
# This returns a {GSL::Vector}[link:rdoc/vector_rdoc.html] object with <tt>k</tt> objects
|
213
|
-
# taken randomly from the {GSL::Vector}[link:rdoc/vector_rdoc.html] object <tt>v</tt>.
|
214
|
-
#
|
215
|
-
# The objects are sampled without replacement, thus each object can only
|
216
|
-
# appear once in the returned vector. It is required that <tt>k</tt> be less
|
217
|
-
# than or equal to the length of the vector <tt>v</tt>.
|
218
|
-
#
|
219
|
-
# ---
|
220
|
-
# * GSL::Rng#sample(v, k)
|
221
|
-
#
|
222
|
-
# This method is like the method <tt>choose</tt> but samples <tt>k</tt> items
|
223
|
-
# from the original vector <tt>v</tt> with replacement, so the same object
|
224
|
-
# can appear more than once in the output sequence. There is no requirement
|
225
|
-
# that <tt>k</tt> be less than the length of <tt>v</tt>.
|
226
|
-
#
|
227
|
-
# {prev}[link:rdoc/qrng_rdoc.html]
|
228
|
-
# {next}[link:rdoc/stats_rdoc.html]
|
229
|
-
#
|
230
|
-
# {Reference index}[link:rdoc/ref_rdoc.html]
|
231
|
-
# {top}[link:index.html]
|
232
|
-
#
|
233
|
-
#
|
data/rdoc/ref.rdoc
DELETED
@@ -1,93 +0,0 @@
|
|
1
|
-
#
|
2
|
-
# = {Ruby/GSL Reference}[link:rdoc/ref_rdoc.html]
|
3
|
-
# (See also {Gnu Scientific Library -- Reference Manual}[https://gnu.org/software/gsl/manual/html_node/])
|
4
|
-
#
|
5
|
-
# == Front Matter
|
6
|
-
# This document describes the modules, classes and the methods of Ruby/GSL. This
|
7
|
-
# includes cut-and-paste from the
|
8
|
-
# {GNU Scientific Library -- Reference Manual}[https://gnu.org/software/gsl/manual/html_node/],
|
9
|
-
# and documents of the extention libraries.
|
10
|
-
#
|
11
|
-
# === Copyright of the GSL Reference
|
12
|
-
#
|
13
|
-
# Copyright 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006,
|
14
|
-
# 2007, 2008 The GSL Team.
|
15
|
-
#
|
16
|
-
# Permission is granted to copy, distribute and/or modify this document under the
|
17
|
-
# terms of the GNU Free Documentation License, Version 1.3 or any later version
|
18
|
-
# published by the Free Software Foundation; with the Invariant Sections being
|
19
|
-
# "GNU General Public License" and "Free Software Needs Free Documentation", the
|
20
|
-
# Front-Cover text being "A GNU Manual", and with the Back-Cover Text being (a)
|
21
|
-
# (see below). A copy of the license is included in the section entitled "GNU
|
22
|
-
# Free Documentation License". (a) The Back-Cover Text is: "You have the freedom
|
23
|
-
# to copy and modify this GNU Manual."
|
24
|
-
#
|
25
|
-
# === Copyright of this reference
|
26
|
-
#
|
27
|
-
# 2009,2010 Yoshiki Tsunesada, David MacMahon
|
28
|
-
#
|
29
|
-
# Permission is granted to copy, distribute and/or modify this document under
|
30
|
-
# the terms of the GNU Free Documentation License.
|
31
|
-
#
|
32
|
-
# == Ruby/GSL Reference
|
33
|
-
#
|
34
|
-
# 1. {Introduction}[link:rdoc/intro_rdoc.html]
|
35
|
-
# 1. {Using Ruby/GSL}[link:rdoc/use_rdoc.html]
|
36
|
-
# 1. {Error Handling}[link:rdoc/ehandling_rdoc.html]
|
37
|
-
# 1. {Mathematical Functions}[link:rdoc/math_rdoc.html]
|
38
|
-
# 1. {Complex Numbers}[link:rdoc/complex_rdoc.html]
|
39
|
-
# 1. {Polynomials}[link:rdoc/poly_rdoc.html]
|
40
|
-
# 1. {Special Functions}[link:rdoc/sf_rdoc.html]
|
41
|
-
# 1. {Vectors}[link:rdoc/vector_rdoc.html] and {Matrices}[link:rdoc/matrix_rdoc.html]
|
42
|
-
# 1. {Permutations}[link:rdoc/perm_rdoc.html]
|
43
|
-
# 1. {Combinations}[link:rdoc/combi_rdoc.html]
|
44
|
-
--
|
45
|
-
# 1. {Multiset}[link:rdoc/sort_rdoc.html] (GSL-1.14)
|
46
|
-
++
|
47
|
-
# 1. {Sorting}[link:rdoc/sort_rdoc.html]
|
48
|
-
# 1. {BLAS Support}[link:rdoc/blas_rdoc.html]
|
49
|
-
# 1. {Linear Algebra}[link:rdoc/linalg_rdoc.html]
|
50
|
-
# 1. {Eigen Systems}[link:rdoc/eigen_rdoc.html]
|
51
|
-
# 1. {Fast Fourier Transform}[link:rdoc/fft_rdoc.html]
|
52
|
-
# 1. {Numerical Integration}[link:rdoc/integration_rdoc.html]
|
53
|
-
# 1. {Random Number Generation}[link:rdoc/rng_rdoc.html]
|
54
|
-
# 1. {Quasi-Random Sequences}[link:rdoc/qrng_rdoc.html]
|
55
|
-
# 1. {Random Number Distributions}[link:rdoc/randist_rdoc.html]
|
56
|
-
# 1. {Statistics}[link:rdoc/stats_rdoc.html]
|
57
|
-
# 1. {1d-Histograms}[link:rdoc/hist_rdoc.html], {2d-Histograms}[link:rdoc/hist2d_rdoc.html] and {3d-Histograms}[link:rdoc/hist3d_rdoc.html]
|
58
|
-
# 1. {N-tuples}[link:rdoc/ntuple_rdoc.html]
|
59
|
-
# 1. {Monte-Carlo Integration}[link:rdoc/monte_rdoc.html]
|
60
|
-
# 1. {Simulated Annealing}[link:rdoc/siman_rdoc.html]
|
61
|
-
# 1. {Ordinary Differential Equations}[link:rdoc/odeiv_rdoc.html]
|
62
|
-
# 1. {Interpolation}[link:rdoc/interp_rdoc.html]
|
63
|
-
# 1. {Numerical Differentiation}[link:rdoc/diff_rdoc.html]
|
64
|
-
# 1. {Chebyshev Approximations}[link:rdoc/cheb_rdoc.html]
|
65
|
-
# 1. {Series Acceleration}[link:rdoc/sum_rdoc.html]
|
66
|
-
# 1. {Wavelet Transforms}[link:rdoc/wavelet_rdoc.html] (GSL-1.6 feature)
|
67
|
-
# 1. {Discrete Hankel Transforms}[link:rdoc/dht_rdoc.html]
|
68
|
-
# 1. {One dimensional Root-Finding}[link:rdoc/roots_rdoc.html]
|
69
|
-
# 1. {One dimensional Minimization}[link:rdoc/min_rdoc.html]
|
70
|
-
# 1. {Multidimensional Root-Finding}[link:rdoc/multiroot_rdoc.html]
|
71
|
-
# 1. {Multidimensional Minimization}[link:rdoc/multimin_rdoc.html]
|
72
|
-
# 1. {Least-Squares Fitting}[link:rdoc/fit_rdoc.html]
|
73
|
-
# 1. {Nonlinear Least-Squares Fitting}[link:rdoc/nonlinearfit_rdoc.html]
|
74
|
-
# 1. {Basis Splines}[link:rdoc/bspline_rdoc.html]
|
75
|
-
# 1. {Physical Constants}[link:rdoc/const_rdoc.html]
|
76
|
-
# 1. {Graphics}[link:rdoc/graph_rdoc.html]
|
77
|
-
# 1. Supported GSL add-on packages
|
78
|
-
# 1. {rngextra}[link:rdoc/rng_rdoc.html#label-Random+number+generator+initialization]
|
79
|
-
# 1. {Tensor manipulations}[link:rdoc/tensor_rdoc.html]
|
80
|
-
# 1. OOL: Open Optimization library (see examples/ool/*.rb)
|
81
|
-
# 1. CQP and Bundle (see examples/multimin/cqp.rb, bundle.rb)
|
82
|
-
# 1. quartic
|
83
|
-
# 1. jacobi (see examples/jacobi/*.rb)
|
84
|
-
# 1. {NDLINEAR: multi-linear, multi-parameter least squares fitting}[link:rdoc/ndlinear_rdoc.html]
|
85
|
-
# 1. {ALF: associated Legendre polynomials}[link:rdoc/alf_rdoc.html]
|
86
|
-
# 1. {NArray compatibilities}[link:rdoc/narray_rdoc.html]
|
87
|
-
# 1. {Changes since Ruby/GSL 1.10.3}[link:rdoc/changes_rdoc.html]
|
88
|
-
# 1. {GNU Free Documentation Licence}[https://gnu.org/software/gsl/manual/html_node/GNU-Free-Documentation-License.html]
|
89
|
-
#
|
90
|
-
# {next}[link:rdoc/intro_rdoc.html]
|
91
|
-
#
|
92
|
-
# {top}[link:index.html]
|
93
|
-
#
|
data/rdoc/rng.rdoc
DELETED
@@ -1,203 +0,0 @@
|
|
1
|
-
#
|
2
|
-
# = Random Number Generation
|
3
|
-
# The library provides a large collection of random number generators which
|
4
|
-
# can be accessed through a uniform interface. Environment variables allow you
|
5
|
-
# to select different generators and seeds at runtime, so that you can easily
|
6
|
-
# switch between generators without needing to recompile your program.
|
7
|
-
# Each instance of a generator keeps track of its own state, allowing the
|
8
|
-
# generators to be used in multi-threaded programs. Additional functions are
|
9
|
-
# available for transforming uniform random numbers into samples from
|
10
|
-
# continuous or discrete probability distributions such as the Gaussian,
|
11
|
-
# log-normal or Poisson distributions.
|
12
|
-
#
|
13
|
-
#
|
14
|
-
#
|
15
|
-
# Contents:
|
16
|
-
# 1. {General comments on random numbers}[link:rdoc/rng_rdoc.html#label-General+comments+on+random+numbers]
|
17
|
-
# 1. {The Random Number Generator Interface: GSL::Rng class}[link:rdoc/rng_rdoc.html#label-The+Random+Number+Generator+Interface]
|
18
|
-
# 1. {Random number generator initialization}[link:rdoc/rng_rdoc.html#label-Random+number+generator+initialization]
|
19
|
-
# 1. {Sampling from a random number generator}[link:rdoc/rng_rdoc.html#label-Sampling+from+a+random+number+generator]
|
20
|
-
# 1. {Auxiliary random number generator functions}[link:rdoc/rng_rdoc.html#label-Auxiliary+random+number+generator+functions]
|
21
|
-
# 1. {Random number environment variables}[link:rdoc/rng_rdoc.html#label-Random+number+environment+variables]
|
22
|
-
#
|
23
|
-
# == General comments on random numbers
|
24
|
-
# In 1988, Park and Miller wrote a paper entitled "Random number generators:
|
25
|
-
# good ones are hard to find." [Commun. ACM, 31, 1192-1201]. Fortunately, some
|
26
|
-
# excellent random number generators are available, though poor ones are still
|
27
|
-
# in common use. You may be happy with the system-supplied random number
|
28
|
-
# generator on your computer, but you should be aware that as computers get
|
29
|
-
# faster, requirements on random number generators increase. Nowadays, a
|
30
|
-
# simulation that calls a random number generator millions of times can often
|
31
|
-
# finish before you can make it down the hall to the coffee machine and back.
|
32
|
-
#
|
33
|
-
# A very nice review of random number generators was written by Pierre L'Ecuyer,
|
34
|
-
# as Chapter 4 of the book: Handbook on Simulation, Jerry Banks, ed.
|
35
|
-
# (Wiley, 1997). The chapter is available in postscript from L'Ecuyer's
|
36
|
-
# ftp site (see references). Knuth's volume on Seminumerical Algorithms
|
37
|
-
# (originally published in 1968) devotes 170 pages to random number generators,
|
38
|
-
# and has recently been updated in its 3rd edition (1997). It is brilliant,
|
39
|
-
# a classic. If you don't own it, you should stop reading right now, run to the
|
40
|
-
# nearest bookstore, and buy it.
|
41
|
-
#
|
42
|
-
# A good random number generator will satisfy both theoretical and statistical
|
43
|
-
# properties. Theoretical properties are often hard to obtain (they require real
|
44
|
-
# math!), but one prefers a random number generator with a long period,
|
45
|
-
# low serial correlation, and a tendency not to "fall mainly on the planes."
|
46
|
-
# Statistical tests are performed with numerical simulations. Generally,
|
47
|
-
# a random number generator is used to estimate some quantity for which the
|
48
|
-
# theory of probability provides an exact answer. Comparison to this exact
|
49
|
-
# answer provides a measure of "randomness".
|
50
|
-
#
|
51
|
-
# == The Random Number Generator Interface
|
52
|
-
# It is important to remember that a random number generator is not a "real"
|
53
|
-
# function like sine or cosine. Unlike real functions, successive calls to a
|
54
|
-
# random number generator yield different return values. Of course that is just
|
55
|
-
# what you want for a random number generator, but to achieve this effect,
|
56
|
-
# the generator must keep track of some kind of "state" variable. Sometimes this
|
57
|
-
# state is just an integer (sometimes just the value of the previously generated
|
58
|
-
# random number), but often it is more complicated than that and may involve a
|
59
|
-
# whole array of numbers, possibly with some indices thrown in. To use the
|
60
|
-
# random number generators, you do not need to know the details of what
|
61
|
-
# comprises the state, and besides that varies from algorithm to algorithm.
|
62
|
-
#
|
63
|
-
# The random number generator library uses <tt>GSL::Rng</tt> class for the interface.
|
64
|
-
# == Random number generator initialization
|
65
|
-
#
|
66
|
-
# ---
|
67
|
-
# * GSL::Rng.alloc(rng_type[, seed])
|
68
|
-
#
|
69
|
-
# This method returns a GSL::Rng object of a random number generator of type
|
70
|
-
# <tt>rng_type</tt> with a seed <tt>seed</tt>. These two arguments can be omitted,
|
71
|
-
# and the generator 'gsl_rng_mt19937' and a seed 0 are used as defaults.
|
72
|
-
# The GSL library provides a number of random number generator types,
|
73
|
-
# and one can choose with a constant <tt>GSL::RNG_xxx</tt> or a String, as
|
74
|
-
#
|
75
|
-
# * <tt>GSL::Rng::MT19937</tt> or <tt>"gsl_rng_mt19937"</tt> or <tt>"mt19937"</tt>
|
76
|
-
# * <tt>GSL::Rng::RANLXS0</tt> or <tt>"gsl_rng_ranlsx0"</tt> or <tt>"ranlxs0"</tt>
|
77
|
-
# * <tt>GSL::Rng::ZUF</tt> or <tt>"gsl_rng_zuf"</tt> or <tt>"zuf"</tt>
|
78
|
-
# * ...
|
79
|
-
#
|
80
|
-
# See the {GSL reference manual}[https://gnu.org/software/gsl/manual/html_node/Random-number-generator-algorithms.html#Random-number-generator-algorithms] for the complete list.
|
81
|
-
# The following demonstrates how to use this class,
|
82
|
-
#
|
83
|
-
# require 'gsl'
|
84
|
-
#
|
85
|
-
# r = Rng.alloc(Rng::TAUS, 1)
|
86
|
-
# r2 = Rng.alloc("ran0", 2)
|
87
|
-
# p r.name <- "taus"
|
88
|
-
# p r.get <- get an integer
|
89
|
-
# p r2.uniform <- get a float of [0, 1)
|
90
|
-
#
|
91
|
-
# A generator of the type <tt>gsl_rng_taus</tt> is created with seed 1,
|
92
|
-
# and <tt>gsl_rng_ran0</tt> with seed 2. The
|
93
|
-
# method <tt>get</tt> returns a random integer.
|
94
|
-
# The method <tt>uniform</tt> returns
|
95
|
-
# a floating number uniformly distributed in the range [0, 1).
|
96
|
-
#
|
97
|
-
# If the package {rngextra}[http://www.network-theory.co.uk/download/rngextra/] is installed, additional
|
98
|
-
# two generators are available,
|
99
|
-
# * <tt>GSL::Rng::RNGEXTRA_RNG1</tt>, <tt>"rngextra_rng1"</tt>
|
100
|
-
# * <tt>GSL::Rng::RNGEXTRA_RNG2</tt>, <tt>"rngextra_rng2"</tt>
|
101
|
-
#
|
102
|
-
# They are created as
|
103
|
-
#
|
104
|
-
# r1 = Rng.alloc(Rng::RNGEXTRA_RNG1)
|
105
|
-
# p r1.name <- "rng1"
|
106
|
-
#
|
107
|
-
# r2 = Rng.alloc("rngextra_rng2")
|
108
|
-
# p r2.name <- "rng2"
|
109
|
-
#
|
110
|
-
# ---
|
111
|
-
# * GSL::Rng.default_seed
|
112
|
-
#
|
113
|
-
# Returns the default seed
|
114
|
-
# ---
|
115
|
-
# * GSL::Rng.set_default_seed(seed)
|
116
|
-
# * GSL::Rng.default_seed=(seed)
|
117
|
-
#
|
118
|
-
# Override the default seed by <tt>seed</tt>.
|
119
|
-
#
|
120
|
-
# ---
|
121
|
-
# * GSL::Rng.types_setup()
|
122
|
-
#
|
123
|
-
# Returns an array of all the available generators' names.
|
124
|
-
#
|
125
|
-
# ---
|
126
|
-
# * GSK::Rng.memcpy(dest, src)
|
127
|
-
#
|
128
|
-
# Copies the random number generator <tt>src</tt>) into the pre-existing generator
|
129
|
-
# <tt>dest</tt>, making dest into an exact copy of <tt>src</tt>.
|
130
|
-
# The two generators must be of the same type.
|
131
|
-
#
|
132
|
-
# ---
|
133
|
-
# * GSL::Rng#set(s)
|
134
|
-
#
|
135
|
-
# This method initializes the random number generator with a given seed <tt>s</tt>.
|
136
|
-
#
|
137
|
-
# == Sampling from a random number generator
|
138
|
-
#
|
139
|
-
# ---
|
140
|
-
# * GSL::Rng#get
|
141
|
-
#
|
142
|
-
# This returns a random integer from the reciever generator.
|
143
|
-
#
|
144
|
-
# ---
|
145
|
-
# * GSL::Rng#uniform
|
146
|
-
#
|
147
|
-
# This method returns a double precision floating point number uniformly
|
148
|
-
# distributed in the range [0,1).
|
149
|
-
#
|
150
|
-
# ---
|
151
|
-
# * GSL::Rng#uniform_pos
|
152
|
-
#
|
153
|
-
# This returns a positive double precision floating point number uniformly
|
154
|
-
# distributed in the range (0,1), excluding both 0.0 and 1.0.
|
155
|
-
#
|
156
|
-
# ---
|
157
|
-
# * GSL::Rng#uniform_int(n)
|
158
|
-
#
|
159
|
-
# This method returns a random integer from 0 to n-1 inclusive.
|
160
|
-
#
|
161
|
-
# == Auxiliary random number generator functions
|
162
|
-
#
|
163
|
-
# ---
|
164
|
-
# * GSL::Rng#name
|
165
|
-
#
|
166
|
-
# This method returns a Sting object of the name of the generator.
|
167
|
-
#
|
168
|
-
# ---
|
169
|
-
# * GSL::Rng#max
|
170
|
-
# * GSL::Rng#min
|
171
|
-
#
|
172
|
-
# These method return the largest/smallest value that the method
|
173
|
-
# <tt>get</tt> can return.
|
174
|
-
#
|
175
|
-
# ---
|
176
|
-
# * GSL::Rng#clone
|
177
|
-
# * GSL::Rng#duplicate
|
178
|
-
#
|
179
|
-
# Return a newly created generator which is an exact copy of the generator <tt>self</tt>.
|
180
|
-
#
|
181
|
-
# == Random number environment variables
|
182
|
-
# The library allows you to choose a default generator and seed from the
|
183
|
-
# environment variables <tt>GSL_RNG_TYPE</tt> and <tt>GSL_RNG_SEED</tt>
|
184
|
-
# and the method <tt>GSL::Rng::env_setup</tt>.
|
185
|
-
#
|
186
|
-
# ---
|
187
|
-
# * GSL::Rng.env_setup()
|
188
|
-
#
|
189
|
-
# Reads the environment variables <tt>GSL_RNG_TYPE</tt> and
|
190
|
-
# <tt>GSL_RNG_SEED</tt> and uses their values to set the corresponding
|
191
|
-
# library variables.
|
192
|
-
#
|
193
|
-
# If you don't specify a generator for <tt>GSL_RNG_TYPE</tt>
|
194
|
-
# then "mt19937" is used as the default.
|
195
|
-
# The initial value of the default seed is zero.
|
196
|
-
#
|
197
|
-
#
|
198
|
-
# {prev}[link:rdoc/integration_rdoc.html]
|
199
|
-
# {next}[link:rdoc/qrng_rdoc.html]
|
200
|
-
#
|
201
|
-
# {Reference index}[link:rdoc/ref_rdoc.html]
|
202
|
-
# {top}[link:index.html]
|
203
|
-
#
|