rb-gsl 1.16.0.5 → 1.16.0.6
Sign up to get free protection for your applications and to get access to all the features.
- 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
|
-
#
|