RubyGems - gsl - Versions diffs - 1.15.3 → 1.16.0.6 - Mend

gsl 1.15.3 → 1.16.0.6

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.

Files changed (446) hide show

checksums.yaml +7 -0
data/.gitignore +12 -0
data/.travis.yml +24 -0
data/AUTHORS +10 -2
data/COPYING +341 -339
data/ChangeLog +612 -554
data/Gemfile +4 -0
data/README.md +77 -0
data/Rakefile +14 -93
data/THANKS +1 -1
data/examples/blas/blas.rb +1 -1
data/examples/bspline.rb +3 -3
data/examples/complex/functions.rb +4 -4
data/examples/complex/mul.rb +1 -1
data/examples/const/physconst.rb +4 -4
data/examples/const/travel.rb +4 -4
data/examples/eigen/nonsymmv.rb +1 -1
data/examples/eigen/qhoscillator.rb +3 -3
data/examples/fft/radix2.rb +1 -1
data/examples/fft/real-halfcomplex.rb +3 -3
data/examples/fft/real-halfcomplex2.rb +3 -3
data/examples/fit/expfit.rb +1 -1
data/examples/fit/multifit.rb +1 -1
data/examples/fit/ndlinear.rb +44 -44
data/examples/fit/nonlinearfit.rb +1 -1
data/examples/fit/wlinear.rb +3 -3
data/examples/function/function.rb +1 -1
data/examples/function/min.rb +1 -1
data/examples/function/synchrotron.rb +2 -2
data/examples/gallery/koch.rb +1 -1
data/examples/histogram/cauchy.rb +2 -2
data/examples/histogram/exponential.rb +1 -1
data/examples/histogram/histo3d.rb +1 -1
data/examples/histogram/histogram-pdf.rb +2 -2
data/examples/histogram/xexp.rb +1 -1
data/examples/integration/ahmed.rb +2 -2
data/examples/integration/cosmology.rb +7 -7
data/examples/integration/friedmann.rb +4 -4
data/examples/integration/qng.rb +1 -1
data/examples/interp/demo.rb +2 -2
data/examples/linalg/LQ_solve.rb +1 -1
data/examples/linalg/LU.rb +1 -1
data/examples/linalg/LU2.rb +1 -1
data/examples/linalg/LU_narray.rb +1 -1
data/examples/linalg/PTLQ.rb +1 -1
data/examples/linalg/QRPT.rb +1 -1
data/examples/linalg/chol.rb +1 -1
data/examples/linalg/chol_narray.rb +1 -1
data/examples/linalg/complex.rb +1 -1
data/examples/math/elementary.rb +1 -1
data/examples/math/functions.rb +1 -1
data/examples/math/inf_nan.rb +1 -1
data/examples/math/minmax.rb +1 -1
data/examples/math/power.rb +1 -1
data/examples/math/test.rb +1 -1
data/examples/min.rb +1 -1
data/examples/multimin/bundle.rb +1 -1
data/examples/multimin/cqp.rb +17 -17
data/examples/multiroot/fsolver3.rb +1 -1
data/examples/odeiv/binarysystem.rb +12 -12
data/examples/odeiv/demo.rb +3 -3
data/examples/odeiv/frei1.rb +7 -7
data/examples/odeiv/frei2.rb +4 -4
data/examples/odeiv/oscillator.rb +1 -1
data/examples/odeiv/sedov.rb +3 -3
data/examples/odeiv/whitedwarf.rb +11 -11
data/examples/permutation/ex1.rb +2 -2
data/examples/permutation/permutation.rb +1 -1
data/examples/poly/demo.rb +1 -1
data/examples/random/diffusion.rb +1 -1
data/examples/random/generator.rb +2 -2
data/examples/random/randomwalk.rb +3 -3
data/examples/random/rng.rb +1 -1
data/examples/roots/bisection.rb +1 -1
data/examples/roots/brent.rb +1 -1
data/examples/roots/demo.rb +1 -1
data/examples/roots/newton.rb +2 -2
data/examples/roots/recombination.gp +0 -1
data/examples/sf/hyperg.rb +1 -1
data/examples/sf/sphbessel.rb +1 -1
data/examples/sort/sort.rb +1 -1
data/examples/tamu_anova.rb +4 -4
data/examples/vector/add.rb +1 -1
data/examples/vector/decimate.rb +1 -1
data/examples/vector/gnuplot.rb +8 -8
data/examples/vector/vector.rb +2 -2
data/examples/wavelet/wavelet1.rb +1 -1
data/ext/{alf.c → gsl_native/alf.c} +10 -10
data/ext/{array.c → gsl_native/array.c} +70 -159
data/ext/{array_complex.c → gsl_native/array_complex.c} +63 -66
data/ext/{blas.c → gsl_native/blas.c} +2 -3
data/ext/{blas1.c → gsl_native/blas1.c} +35 -36
data/ext/{blas2.c → gsl_native/blas2.c} +57 -62
data/ext/{blas3.c → gsl_native/blas3.c} +57 -58
data/ext/{block.c → gsl_native/block.c} +14 -18
data/ext/{block_source.c → gsl_native/block_source.h} +110 -112
data/ext/gsl_native/bspline.c +122 -0
data/ext/{bundle.c → gsl_native/bundle.c} +0 -0
data/ext/{cdf.c → gsl_native/cdf.c} +79 -93
data/ext/{cheb.c → gsl_native/cheb.c} +78 -89
data/ext/{combination.c → gsl_native/combination.c} +11 -19
data/ext/{common.c → gsl_native/common.c} +9 -41
data/ext/{complex.c → gsl_native/complex.c} +116 -118
data/ext/gsl_native/const.c +331 -0
data/ext/{const_additional.c → gsl_native/const_additional.c} +13 -34
data/ext/gsl_native/cqp.c +283 -0
data/ext/{deriv.c → gsl_native/deriv.c} +25 -33
data/ext/{dht.c → gsl_native/dht.c} +23 -31
data/ext/{diff.c → gsl_native/diff.c} +26 -28
data/ext/{dirac.c → gsl_native/dirac.c} +45 -46
data/ext/{eigen.c → gsl_native/eigen.c} +1044 -1095
data/ext/{error.c → gsl_native/error.c} +18 -18
data/ext/gsl_native/extconf.rb +118 -0
data/ext/{fft.c → gsl_native/fft.c} +197 -204
data/ext/{fit.c → gsl_native/fit.c} +17 -18
data/ext/gsl_native/fresnel.c +312 -0
data/ext/{function.c → gsl_native/function.c} +37 -43
data/ext/{geometry.c → gsl_native/geometry.c} +16 -16
data/ext/{graph.c → gsl_native/graph.c} +39 -89
data/ext/{gsl.c → gsl_native/gsl.c} +12 -33
data/ext/{gsl_narray.c → gsl_native/gsl_narray.c} +20 -30
data/ext/{histogram.c → gsl_native/histogram.c} +133 -160
data/ext/{histogram2d.c → gsl_native/histogram2d.c} +78 -104
data/ext/{histogram3d.c → gsl_native/histogram3d.c} +76 -76
data/ext/{histogram3d_source.c → gsl_native/histogram3d_source.c} +196 -197
data/ext/{histogram_find.c → gsl_native/histogram_find.c} +32 -34
data/ext/{histogram_oper.c → gsl_native/histogram_oper.c} +43 -52
data/ext/{ieee.c → gsl_native/ieee.c} +9 -21
data/{include → ext/gsl_native/include}/rb_gsl.h +4 -26
data/{include → ext/gsl_native/include}/rb_gsl_array.h +15 -39
data/{include → ext/gsl_native/include}/rb_gsl_cheb.h +0 -2
data/{include → ext/gsl_native/include}/rb_gsl_common.h +61 -61
data/{include → ext/gsl_native/include}/rb_gsl_complex.h +1 -1
data/{include → ext/gsl_native/include}/rb_gsl_const.h +0 -6
data/ext/gsl_native/include/rb_gsl_dirac.h +6 -0
data/{include → ext/gsl_native/include}/rb_gsl_eigen.h +1 -1
data/{include → ext/gsl_native/include}/rb_gsl_fft.h +0 -13
data/{include → ext/gsl_native/include}/rb_gsl_fit.h +0 -2
data/{include → ext/gsl_native/include}/rb_gsl_function.h +0 -4
data/{include → ext/gsl_native/include}/rb_gsl_graph.h +2 -4
data/{include → ext/gsl_native/include}/rb_gsl_histogram.h +8 -8
data/{include → ext/gsl_native/include}/rb_gsl_histogram3d.h +50 -50
data/{include → ext/gsl_native/include}/rb_gsl_integration.h +1 -1
data/{include → ext/gsl_native/include}/rb_gsl_interp.h +0 -5
data/{include → ext/gsl_native/include}/rb_gsl_linalg.h +2 -6
data/{include → ext/gsl_native/include}/rb_gsl_math.h +0 -6
data/{include → ext/gsl_native/include}/rb_gsl_odeiv.h +0 -3
data/{include → ext/gsl_native/include}/rb_gsl_poly.h +3 -7
data/{include → ext/gsl_native/include}/rb_gsl_rational.h +1 -8
data/{include → ext/gsl_native/include}/rb_gsl_rng.h +0 -1
data/{include → ext/gsl_native/include}/rb_gsl_root.h +1 -1
data/{include → ext/gsl_native/include}/rb_gsl_sf.h +39 -48
data/{include → ext/gsl_native/include}/rb_gsl_statistics.h +1 -1
data/{include → ext/gsl_native/include}/rb_gsl_tensor.h +0 -2
data/{include → ext/gsl_native/include}/rb_gsl_with_narray.h +3 -1
data/{include → ext/gsl_native/include}/templates_off.h +0 -0
data/{include → ext/gsl_native/include}/templates_on.h +1 -1
data/ext/{integration.c → gsl_native/integration.c} +164 -189
data/ext/{interp.c → gsl_native/interp.c} +25 -38
data/ext/gsl_native/jacobi.c +733 -0
data/ext/{linalg.c → gsl_native/linalg.c} +462 -589
data/ext/{linalg_complex.c → gsl_native/linalg_complex.c} +93 -106
data/ext/{math.c → gsl_native/math.c} +48 -67
data/ext/{matrix.c → gsl_native/matrix.c} +13 -16
data/ext/{matrix_complex.c → gsl_native/matrix_complex.c} +119 -123
data/ext/{matrix_double.c → gsl_native/matrix_double.c} +79 -82
data/ext/{matrix_int.c → gsl_native/matrix_int.c} +53 -54
data/ext/{matrix_source.c → gsl_native/matrix_source.h} +292 -318
data/ext/{min.c → gsl_native/min.c} +45 -76
data/ext/{monte.c → gsl_native/monte.c} +50 -64
data/ext/{multifit.c → gsl_native/multifit.c} +142 -151
data/ext/{multimin.c → gsl_native/multimin.c} +64 -92
data/ext/{multimin_fsdf.c → gsl_native/multimin_fsdf.c} +16 -16
data/ext/{multiroots.c → gsl_native/multiroots.c} +73 -76
data/ext/{multiset.c → gsl_native/multiset.c} +4 -8
data/ext/{ndlinear.c → gsl_native/ndlinear.c} +320 -321
data/ext/{nmf.c → gsl_native/nmf.c} +11 -11
data/ext/{nmf_wrap.c → gsl_native/nmf_wrap.c} +1 -1
data/ext/{ntuple.c → gsl_native/ntuple.c} +23 -23
data/ext/{odeiv.c → gsl_native/odeiv.c} +101 -116
data/ext/gsl_native/ool.c +879 -0
data/ext/{permutation.c → gsl_native/permutation.c} +39 -37
data/ext/{poly.c → gsl_native/poly.c} +10 -13
data/ext/{poly2.c → gsl_native/poly2.c} +16 -16
data/ext/{poly_source.c → gsl_native/poly_source.h} +249 -293
data/ext/{qrng.c → gsl_native/qrng.c} +9 -20
data/ext/{randist.c → gsl_native/randist.c} +222 -247
data/ext/{rational.c → gsl_native/rational.c} +12 -12
data/ext/{rng.c → gsl_native/rng.c} +30 -47
data/ext/{root.c → gsl_native/root.c} +47 -48
data/ext/{sf.c → gsl_native/sf.c} +196 -244
data/ext/{sf_airy.c → gsl_native/sf_airy.c} +2 -2
data/ext/{sf_bessel.c → gsl_native/sf_bessel.c} +7 -7
data/ext/{sf_clausen.c → gsl_native/sf_clausen.c} +1 -1
data/ext/{sf_coulomb.c → gsl_native/sf_coulomb.c} +40 -40
data/ext/{sf_coupling.c → gsl_native/sf_coupling.c} +30 -30
data/ext/{sf_dawson.c → gsl_native/sf_dawson.c} +1 -1
data/ext/{sf_debye.c → gsl_native/sf_debye.c} +1 -10
data/ext/{sf_dilog.c → gsl_native/sf_dilog.c} +1 -1
data/ext/{sf_elementary.c → gsl_native/sf_elementary.c} +3 -3
data/ext/{sf_ellint.c → gsl_native/sf_ellint.c} +43 -43
data/ext/{sf_elljac.c → gsl_native/sf_elljac.c} +3 -3
data/ext/{sf_erfc.c → gsl_native/sf_erfc.c} +1 -5
data/ext/{sf_exp.c → gsl_native/sf_exp.c} +3 -3
data/ext/{sf_expint.c → gsl_native/sf_expint.c} +2 -12
data/ext/{sf_fermi_dirac.c → gsl_native/sf_fermi_dirac.c} +1 -1
data/ext/{sf_gamma.c → gsl_native/sf_gamma.c} +2 -6
data/ext/{sf_gegenbauer.c → gsl_native/sf_gegenbauer.c} +1 -1
data/ext/{sf_hyperg.c → gsl_native/sf_hyperg.c} +1 -1
data/ext/{sf_laguerre.c → gsl_native/sf_laguerre.c} +4 -4
data/ext/{sf_lambert.c → gsl_native/sf_lambert.c} +1 -1
data/ext/{sf_legendre.c → gsl_native/sf_legendre.c} +1 -1
data/ext/{sf_log.c → gsl_native/sf_log.c} +4 -4
data/ext/gsl_native/sf_mathieu.c +235 -0
data/ext/{sf_power.c → gsl_native/sf_power.c} +1 -1
data/ext/{sf_psi.c → gsl_native/sf_psi.c} +3 -12
data/ext/{sf_synchrotron.c → gsl_native/sf_synchrotron.c} +1 -1
data/ext/{sf_transport.c → gsl_native/sf_transport.c} +1 -1
data/ext/{sf_trigonometric.c → gsl_native/sf_trigonometric.c} +4 -4
data/ext/{sf_zeta.c → gsl_native/sf_zeta.c} +1 -5
data/ext/{signal.c → gsl_native/signal.c} +63 -68
data/ext/{siman.c → gsl_native/siman.c} +45 -49
data/ext/{sort.c → gsl_native/sort.c} +6 -7
data/ext/{spline.c → gsl_native/spline.c} +28 -46
data/ext/{stats.c → gsl_native/stats.c} +105 -118
data/ext/{sum.c → gsl_native/sum.c} +34 -34
data/ext/{tamu_anova.c → gsl_native/tamu_anova.c} +1 -1
data/ext/{tensor.c → gsl_native/tensor.c} +8 -11
data/ext/{tensor_source.c → gsl_native/tensor_source.h} +147 -148
data/ext/{vector.c → gsl_native/vector.c} +11 -14
data/ext/{vector_complex.c → gsl_native/vector_complex.c} +179 -184
data/ext/{vector_double.c → gsl_native/vector_double.c} +178 -183
data/ext/{vector_int.c → gsl_native/vector_int.c} +27 -29
data/ext/{vector_source.c → gsl_native/vector_source.h} +428 -443
data/ext/{wavelet.c → gsl_native/wavelet.c} +224 -246
data/gsl.gemspec +29 -0
data/lib/gsl.rb +8 -3
data/lib/gsl/gnuplot.rb +3 -3
data/lib/gsl/oper.rb +35 -60
data/lib/gsl/version.rb +3 -0
data/lib/rbgsl.rb +1 -3
data/rdoc/alf.rdoc +5 -5
data/rdoc/blas.rdoc +9 -9
data/rdoc/bspline.rdoc +17 -17
data/rdoc/changes.rdoc +4 -9
data/rdoc/cheb.rdoc +25 -25
data/rdoc/cholesky_complex.rdoc +21 -21
data/rdoc/combi.rdoc +37 -37
data/rdoc/complex.rdoc +22 -22
data/rdoc/const.rdoc +47 -47
data/rdoc/dht.rdoc +49 -49
data/rdoc/diff.rdoc +42 -42
data/rdoc/ehandling.rdoc +6 -6
data/rdoc/eigen.rdoc +153 -153
data/rdoc/fft.rdoc +146 -146
data/rdoc/fit.rdoc +109 -109
data/rdoc/function.rdoc +11 -11
data/rdoc/graph.rdoc +17 -17
data/rdoc/hist.rdoc +103 -103
data/rdoc/hist2d.rdoc +42 -42
data/rdoc/hist3d.rdoc +9 -9
data/rdoc/integration.rdoc +110 -110
data/rdoc/interp.rdoc +71 -71
data/rdoc/intro.rdoc +8 -8
data/rdoc/linalg.rdoc +188 -188
data/rdoc/linalg_complex.rdoc +1 -1
data/rdoc/math.rdoc +58 -58
data/rdoc/matrix.rdoc +275 -275
data/rdoc/min.rdoc +57 -57
data/rdoc/monte.rdoc +22 -22
data/rdoc/multimin.rdoc +95 -95
data/rdoc/multiroot.rdoc +80 -80
data/rdoc/narray.rdoc +32 -32
data/rdoc/ndlinear.rdoc +54 -54
data/rdoc/nonlinearfit.rdoc +100 -100
data/rdoc/ntuple.rdoc +31 -31
data/rdoc/odeiv.rdoc +88 -88
data/rdoc/perm.rdoc +90 -90
data/rdoc/poly.rdoc +66 -66
data/rdoc/qrng.rdoc +21 -21
data/rdoc/randist.rdoc +82 -82
data/rdoc/ref.rdoc +57 -57
data/rdoc/rng.rdoc +85 -85
data/rdoc/roots.rdoc +57 -57
data/rdoc/sf.rdoc +428 -428
data/rdoc/siman.rdoc +19 -19
data/rdoc/sort.rdoc +30 -30
data/rdoc/start.rdoc +8 -8
data/rdoc/stats.rdoc +52 -52
data/rdoc/sum.rdoc +12 -12
data/rdoc/tensor.rdoc +31 -31
data/rdoc/tut.rdoc +1 -1
data/rdoc/use.rdoc +39 -39
data/rdoc/vector.rdoc +188 -188
data/rdoc/vector_complex.rdoc +24 -24
data/rdoc/wavelet.rdoc +46 -46
data/test/gsl/blas_test.rb +79 -0
data/test/gsl/bspline_test.rb +63 -0
data/test/gsl/cdf_test.rb +1512 -0
data/test/gsl/cheb_test.rb +80 -0
data/test/gsl/combination_test.rb +100 -0
data/test/gsl/complex_test.rb +20 -0
data/test/gsl/const_test.rb +29 -0
data/test/gsl/deriv_test.rb +62 -0
data/test/gsl/dht_test.rb +79 -0
data/test/gsl/diff_test.rb +53 -0
data/test/gsl/eigen_test.rb +563 -0
data/test/gsl/err_test.rb +23 -0
data/test/gsl/fit_test.rb +101 -0
data/test/gsl/histo_test.rb +14 -0
data/test/gsl/index_test.rb +61 -0
data/test/gsl/integration_test.rb +274 -0
data/test/gsl/interp_test.rb +27 -0
data/test/gsl/linalg_test.rb +463 -0
data/test/gsl/matrix_nmf_test.rb +37 -0
data/test/gsl/matrix_test.rb +98 -0
data/test/gsl/min_test.rb +89 -0
data/test/gsl/monte_test.rb +77 -0
data/test/gsl/multifit_test.rb +753 -0
data/test/gsl/multimin_test.rb +157 -0
data/test/gsl/multiroot_test.rb +135 -0
data/test/gsl/multiset_test.rb +52 -0
data/test/gsl/odeiv_test.rb +275 -0
data/test/gsl/oper_test.rb +98 -0
data/test/gsl/poly_test.rb +338 -0
data/test/gsl/qrng_test.rb +94 -0
data/test/gsl/quartic_test.rb +28 -0
data/test/gsl/randist_test.rb +122 -0
data/test/gsl/rng_test.rb +303 -0
data/test/gsl/roots_test.rb +78 -0
data/test/gsl/sf_test.rb +2079 -0
data/test/gsl/stats_test.rb +122 -0
data/test/gsl/sum_test.rb +69 -0
data/test/gsl/tensor_test.rb +396 -0
data/test/gsl/vector_test.rb +223 -0
data/test/gsl/wavelet_test.rb +130 -0
data/test/gsl_test.rb +321 -0
data/test/test_helper.rb +42 -0
data/uncrustify.cfg +1693 -0
metadata +337 -378
data/README +0 -32
data/VERSION +0 -1
data/ext/bspline.c +0 -130
data/ext/const.c +0 -673
data/ext/cqp.c +0 -283
data/ext/extconf.rb +0 -295
data/ext/fcmp.c +0 -66
data/ext/fresnel.c +0 -312
data/ext/jacobi.c +0 -739
data/ext/ool.c +0 -879
data/ext/oper_complex_source.c +0 -253
data/ext/sf_mathieu.c +0 -238
data/include/rb_gsl_config.h +0 -62
data/include/rb_gsl_dirac.h +0 -13
data/rdoc/index.rdoc +0 -62
data/rdoc/rngextra.rdoc +0 -11
data/rdoc/screenshot.rdoc +0 -40
data/setup.rb +0 -1585
data/tests/blas/amax.rb +0 -14
data/tests/blas/asum.rb +0 -16
data/tests/blas/axpy.rb +0 -25
data/tests/blas/copy.rb +0 -23
data/tests/blas/dot.rb +0 -23
data/tests/bspline.rb +0 -53
data/tests/cdf.rb +0 -1388
data/tests/cheb.rb +0 -112
data/tests/combination.rb +0 -123
data/tests/complex.rb +0 -17
data/tests/const.rb +0 -24
data/tests/deriv.rb +0 -85
data/tests/dht/dht1.rb +0 -17
data/tests/dht/dht2.rb +0 -23
data/tests/dht/dht3.rb +0 -23
data/tests/dht/dht4.rb +0 -23
data/tests/diff.rb +0 -78
data/tests/eigen/eigen.rb +0 -220
data/tests/eigen/gen.rb +0 -105
data/tests/eigen/genherm.rb +0 -66
data/tests/eigen/gensymm.rb +0 -68
data/tests/eigen/nonsymm.rb +0 -53
data/tests/eigen/nonsymmv.rb +0 -53
data/tests/eigen/symm-herm.rb +0 -74
data/tests/err.rb +0 -58
data/tests/fit.rb +0 -124
data/tests/gsl_test.rb +0 -118
data/tests/gsl_test2.rb +0 -110
data/tests/histo.rb +0 -12
data/tests/integration/integration1.rb +0 -72
data/tests/integration/integration2.rb +0 -71
data/tests/integration/integration3.rb +0 -71
data/tests/integration/integration4.rb +0 -71
data/tests/interp.rb +0 -45
data/tests/linalg/HH.rb +0 -64
data/tests/linalg/LU.rb +0 -47
data/tests/linalg/QR.rb +0 -77
data/tests/linalg/SV.rb +0 -24
data/tests/linalg/TDN.rb +0 -116
data/tests/linalg/TDS.rb +0 -122
data/tests/linalg/bidiag.rb +0 -73
data/tests/linalg/cholesky.rb +0 -20
data/tests/linalg/linalg.rb +0 -158
data/tests/matrix/matrix_complex_test.rb +0 -36
data/tests/matrix/matrix_nmf_test.rb +0 -39
data/tests/matrix/matrix_test.rb +0 -48
data/tests/min.rb +0 -99
data/tests/monte/miser.rb +0 -31
data/tests/monte/vegas.rb +0 -45
data/tests/multifit/test_2dgauss.rb +0 -112
data/tests/multifit/test_brown.rb +0 -90
data/tests/multifit/test_enso.rb +0 -246
data/tests/multifit/test_filip.rb +0 -155
data/tests/multifit/test_gauss.rb +0 -97
data/tests/multifit/test_longley.rb +0 -110
data/tests/multifit/test_multifit.rb +0 -52
data/tests/multimin.rb +0 -139
data/tests/multiroot.rb +0 -131
data/tests/multiset.rb +0 -52
data/tests/narray/blas_dnrm2.rb +0 -20
data/tests/odeiv.rb +0 -353
data/tests/poly/poly.rb +0 -290
data/tests/poly/special.rb +0 -65
data/tests/qrng.rb +0 -131
data/tests/quartic.rb +0 -29
data/tests/randist.rb +0 -134
data/tests/rng.rb +0 -305
data/tests/roots.rb +0 -76
data/tests/run-test.sh +0 -17
data/tests/sf/gsl_test_sf.rb +0 -249
data/tests/sf/test_airy.rb +0 -83
data/tests/sf/test_bessel.rb +0 -306
data/tests/sf/test_coulomb.rb +0 -17
data/tests/sf/test_dilog.rb +0 -25
data/tests/sf/test_gamma.rb +0 -209
data/tests/sf/test_hyperg.rb +0 -356
data/tests/sf/test_legendre.rb +0 -227
data/tests/sf/test_mathieu.rb +0 -59
data/tests/sf/test_mode.rb +0 -19
data/tests/sf/test_sf.rb +0 -839
data/tests/stats.rb +0 -174
data/tests/stats_mt.rb +0 -16
data/tests/sum.rb +0 -98
data/tests/sys.rb +0 -323
data/tests/tensor.rb +0 -419
data/tests/vector/vector_complex_test.rb +0 -101
data/tests/vector/vector_test.rb +0 -141
data/tests/wavelet.rb +0 -142

data/ext/{rational.c → gsl_native/rational.c} RENAMED

@@ -8,9 +8,9 @@
   This library is distributed in the hope that it will be useful, but
   WITHOUT ANY WARRANTY.
 */
-#include "rb_gsl_config.h"
-#include "rb_gsl_rational.h"
-#include "rb_gsl_array.h"
+#include "include/rb_gsl_rational.h"
+#include "include/rb_gsl_array.h"
 VALUE cgsl_rational;
 static gsl_rational* gsl_rational_div_poly(const gsl_rational *r1, const gsl_poly *p);
@@ -199,7 +199,7 @@ static VALUE rb_gsl_poly_make_rational(VALUE obj, VALUE other)
     case T_ARRAY:
       p2 = gsl_vector_alloc(RARRAY_LEN(other));
       for (i = 0; i < p2->size; i++)
-	gsl_vector_set(p2, i, NUM2DBL(rb_ary_entry(other, i)));
+        gsl_vector_set(p2, i, NUM2DBL(rb_ary_entry(other, i)));
       rnew = gsl_rational_new(p, p2);
       gsl_vector_free(p2);
       break;
@@ -210,8 +210,8 @@ static VALUE rb_gsl_poly_make_rational(VALUE obj, VALUE other)
       return Data_Wrap_Struct(cgsl_poly, 0, gsl_vector_free, p2);
       break;
     default:
-      rb_raise(rb_eTypeError, "wrong argument type %s",
-	       rb_class2name(CLASS_OF(other)));
+      rb_raise(rb_eTypeError, "wrong argument type %s",
+               rb_class2name(CLASS_OF(other)));
       break;
     }
   }
@@ -280,7 +280,7 @@ static VALUE rb_gsl_rational_uminus(VALUE obj)
   if (RATIONAL_P(obj)) {
     Data_Get_Struct(obj, gsl_rational, r);
     rnew = gsl_rational_new(r->pnum, r->pden);
-    for (i = 0; i < rnew->pnum->size; i++)
+    for (i = 0; i < rnew->pnum->size; i++)
       gsl_vector_set(rnew->pnum, i, -gsl_vector_get(r->pnum, i));
     return Data_Wrap_Struct(cgsl_rational, gsl_rational_mark, gsl_rational_free, rnew);
   } else {
@@ -342,7 +342,7 @@ static VALUE rb_gsl_rational_div(VALUE obj, VALUE other)
     case T_ARRAY:
       p = gsl_vector_alloc(RARRAY_LEN(other));
       for (i = 0; i < p->size; i++)
-	gsl_vector_set(p, i, NUM2DBL(rb_ary_entry(other, i)));
+        gsl_vector_set(p, i, NUM2DBL(rb_ary_entry(other, i)));
       rnew = gsl_rational_div_poly(r, p);
       gsl_vector_free(p);
       break;
@@ -352,8 +352,8 @@ static VALUE rb_gsl_rational_div(VALUE obj, VALUE other)
       gsl_vector_scale(rnew->pnum, 1.0/NUM2DBL(other));
       break;
     default:
-      rb_raise(rb_eTypeError, "wrong argument type %s",
-	       rb_class2name(CLASS_OF(other)));
+      rb_raise(rb_eTypeError, "wrong argument type %s",
+               rb_class2name(CLASS_OF(other)));
       break;
     }
   }
@@ -407,8 +407,8 @@ static VALUE rb_gsl_rational_coerce(VALUE obj, VALUE other)
   ptmp = gsl_vector_alloc(1);
   gsl_vector_set(ptmp, 0, 1.0);
   r = gsl_rational_new2(p, ptmp);
-  return rb_ary_new3(2,
-		     Data_Wrap_Struct(cgsl_rational, gsl_rational_mark, gsl_rational_free, r), obj);
+  return rb_ary_new3(2,
+                     Data_Wrap_Struct(cgsl_rational, gsl_rational_mark, gsl_rational_free, r), obj);
 }
 static VALUE rb_gsl_rational_zero(VALUE obj)

data/ext/{rng.c → gsl_native/rng.c} RENAMED

@@ -9,13 +9,12 @@
   WITHOUT ANY WARRANTY.
 */
-/*
+/*
    Document-class: <i>GSL::Rng</i>
    Random number generator
 */
-#include "rb_gsl_config.h"
-#include "rb_gsl_rng.h"
+#include "include/rb_gsl_rng.h"
 #ifdef HAVE_RNGEXTRA_RNGEXTRA_H
 #include "rngextra/rngextra.h"
 #endif
@@ -25,22 +24,22 @@ VALUE cgsl_rng;
 enum rb_gsl_rng_generator {
   GSL_RNG_DEFAULT,
-  GSL_RNG_MT19937, GSL_RNG_MT19937_1999, GSL_RNG_MT19937_1998,
-  GSL_RNG_RANLXS0, GSL_RNG_RANLXS1, GSL_RNG_RANLXS2,
-  GSL_RNG_RANLXD1, GSL_RNG_RANLXD2,
+  GSL_RNG_MT19937, GSL_RNG_MT19937_1999, GSL_RNG_MT19937_1998,
+  GSL_RNG_RANLXS0, GSL_RNG_RANLXS1, GSL_RNG_RANLXS2,
+  GSL_RNG_RANLXD1, GSL_RNG_RANLXD2,
   GSL_RNG_RANLUX, GSL_RNG_RANLUX389,
-  GSL_RNG_CMRG, GSL_RNG_MRG,
+  GSL_RNG_CMRG, GSL_RNG_MRG,
   GSL_RNG_TAUS, GSL_RNG_TAUS2, GSL_RNG_TAUS113,  GSL_RNG_GFSR4,
-  GSL_RNG_RAND,
+  GSL_RNG_RAND,
   GSL_RNG_RANDOM_BSD, GSL_RNG_RANDOM_GLIBC2,
-  GSL_RNG_RANDOM8_GLIBC2, GSL_RNG_RANDOM32_GLIBC2, GSL_RNG_RANDOM64_GLIBC2,
+  GSL_RNG_RANDOM8_GLIBC2, GSL_RNG_RANDOM32_GLIBC2, GSL_RNG_RANDOM64_GLIBC2,
   GSL_RNG_RANDOM128_GLIBC2, GSL_RNG_RANDOM256_GLIBC2,
-  GSL_RNG_RANDOM8_BSD, GSL_RNG_RANDOM32_BSD, GSL_RNG_RANDOM64_BSD,
+  GSL_RNG_RANDOM8_BSD, GSL_RNG_RANDOM32_BSD, GSL_RNG_RANDOM64_BSD,
   GSL_RNG_RANDOM128_BSD, GSL_RNG_RANDOM256_BSD,
-  GSL_RNG_RANDOM_LIBC5, GSL_RNG_RANDOM8_LIBC5, GSL_RNG_RANDOM32_LIBC5,
-  GSL_RNG_RANDOM64_LIBC5, GSL_RNG_RANDOM128_LIBC5, GSL_RNG_RANDOM256_LIBC5,
-  GSL_RNG_RAND48,
-  GSL_RNG_RAN0, GSL_RNG_RAN1, GSL_RNG_RAN2, GSL_RNG_RAN3,
+  GSL_RNG_RANDOM_LIBC5, GSL_RNG_RANDOM8_LIBC5, GSL_RNG_RANDOM32_LIBC5,
+  GSL_RNG_RANDOM64_LIBC5, GSL_RNG_RANDOM128_LIBC5, GSL_RNG_RANDOM256_LIBC5,
+  GSL_RNG_RAND48,
+  GSL_RNG_RAN0, GSL_RNG_RAN1, GSL_RNG_RAN2, GSL_RNG_RAN3,
   GSL_RNG_RANF, GSL_RNG_RANMAR, GSL_RNG_R250, GSL_RNG_TT800, GSL_RNG_VAX,
   GSL_RNG_TRANSPUTER, GSL_RNG_RANDU, GSL_RNG_MINSTD,
   GSL_RNG_UNI, GSL_RNG_UNI32, GSL_RNG_SLATEC, GSL_RNG_ZUF,
@@ -71,16 +70,16 @@ static VALUE rb_gsl_rng_alloc(int argc, VALUE *argv, VALUE klass)
     T = gsl_rng_default;
     seed = gsl_rng_default_seed;
   } else {
-    T = get_gsl_rng_type(argv[0]);
+    T = get_gsl_rng_type(argv[0]);
     if (argc == 1) {
       seed = gsl_rng_default_seed;
     } else if (argc == 2) {
       itype = TYPE(argv[1]);
       if (itype == T_FIXNUM || itype == T_BIGNUM) {
-	seed = FIX2INT(argv[1]);
+        seed = FIX2INT(argv[1]);
       } else {
-	rb_raise(rb_eArgError,
-		 "bad argument 2, seed must be an integer.");
+        rb_raise(rb_eArgError,
+                 "bad argument 2, seed must be an integer.");
       }
     } else {
       rb_raise(rb_eArgError, "too many arguments (%d for 0 or 1)", argc);
@@ -111,7 +110,6 @@ static const gsl_rng_type* get_gsl_rng_type_name(char *name)
 {
   if (str_tail_grep(name, "default") == 0) return gsl_rng_default;
   else if (str_tail_grep(name, "mt19937") == 0) return gsl_rng_mt19937;
-#ifdef GSL_1_1_LATER
   else if (str_tail_grep(name, "borosh13") == 0) return gsl_rng_borosh13;
   else if (str_tail_grep(name, "coveyou") == 0) return gsl_rng_coveyou;
   else if (str_tail_grep(name, "fishman18") == 0) return gsl_rng_fishman18;
@@ -121,15 +119,12 @@ static const gsl_rng_type* get_gsl_rng_type_name(char *name)
   else if (str_tail_grep(name, "waterman14") == 0) return gsl_rng_waterman14;
   else if (str_tail_grep(name, "knuthran") == 0) return gsl_rng_knuthran;
   else if (str_tail_grep(name, "knuthran2") == 0) return gsl_rng_knuthran2;
-#endif
-#ifdef GSL_1_2_LATER
   else if (str_tail_grep(name, "mt19937_1999") == 0) return gsl_rng_mt19937_1999;
   else if (str_tail_grep(name, "mt19937-1999") == 0) return gsl_rng_mt19937_1999;
-  else if (str_tail_grep(name, "mt19937_1998") == 0) return gsl_rng_mt19937_1998;
-  else if (str_tail_grep(name, "mt19937-1998") == 0) return gsl_rng_mt19937_1998;
+  else if (str_tail_grep(name, "mt19937_1998") == 0) return gsl_rng_mt19937_1998;
+  else if (str_tail_grep(name, "mt19937-1998") == 0) return gsl_rng_mt19937_1998;
   else if (str_tail_grep(name, "taus113") == 0) return gsl_rng_taus113;
   else if (str_tail_grep(name, "taus2") == 0) return gsl_rng_taus2;
-#endif
   else if (str_tail_grep(name, "mt19937") == 0) return gsl_rng_mt19937;
   else if (str_tail_grep(name, "ranlxs0") == 0) return gsl_rng_ranlxs0;
   else if (str_tail_grep(name, "ranlxs1") == 0) return gsl_rng_ranlxs1;
@@ -148,7 +143,7 @@ static const gsl_rng_type* get_gsl_rng_type_name(char *name)
   else if (str_tail_grep(name, "random32_libc5") == 0) return gsl_rng_random32_libc5;
   else if (str_tail_grep(name, "random64_libc5") == 0) return gsl_rng_random64_libc5;
   else if (str_tail_grep(name, "random128_libc5") == 0) return gsl_rng_random128_libc5;
-  else if (str_tail_grep(name, "random256_libc5") == 0) return gsl_rng_random256_libc5;
+  else if (str_tail_grep(name, "random256_libc5") == 0) return gsl_rng_random256_libc5;
   else if (str_tail_grep(name, "random-libc5") == 0) return gsl_rng_random_libc5;
   else if (str_tail_grep(name, "random8-libc5") == 0) return gsl_rng_random8_libc5;
   else if (str_tail_grep(name, "random32-libc5") == 0) return gsl_rng_random32_libc5;
@@ -202,18 +197,16 @@ static const gsl_rng_type* get_gsl_rng_type_name(char *name)
   else if (str_tail_grep(name, "rngextra-rng1") == 0) return rngextra_rng1;
   else if (str_tail_grep(name, "rngextra-rng2") == 0) return rngextra_rng2;
 #else
-  else if (str_tail_grep(name, "rngextra_rng1")*str_tail_grep(name, "rngextra_rng2") == 0)
+  else if (str_tail_grep(name, "rngextra_rng1")*str_tail_grep(name, "rngextra_rng2") == 0)
     rb_raise(rb_eNotImpError, "Install the rngextra package found at <http://www.network-theory.co.uk/download/rngextra/>.");
-  else if (str_tail_grep(name, "rngextra_rng2")*str_tail_grep(name, "rngextra_rng2") == 0)
+  else if (str_tail_grep(name, "rngextra_rng2")*str_tail_grep(name, "rngextra_rng2") == 0)
     rb_raise(rb_eNotImpError, "Install the rngextra package found at <http://www.network-theory.co.uk/download/rngextra/>.");
-  else if (str_tail_grep(name, "rngextra-rng1")*str_tail_grep(name, "rngextra_rng2") == 0)
+  else if (str_tail_grep(name, "rngextra-rng1")*str_tail_grep(name, "rngextra_rng2") == 0)
     rb_raise(rb_eNotImpError, "Install the rngextra package found at <http://www.network-theory.co.uk/download/rngextra/>.");
-  else if (str_tail_grep(name, "rngextra-rng2")*str_tail_grep(name, "rngextra_rng2") == 0)
+  else if (str_tail_grep(name, "rngextra-rng2")*str_tail_grep(name, "rngextra_rng2") == 0)
     rb_raise(rb_eNotImpError, "Install the rngextra package found at <http://www.network-theory.co.uk/download/rngextra/>.");
 #endif
-#ifdef GSL_1_9_LATER
   else if (str_tail_grep(name, "knuthran2002") == 0) return gsl_rng_knuthran2002;
-#endif
   else
     rb_raise(rb_eArgError, "unknown generator type \"%s\"", name);
 }
@@ -225,12 +218,10 @@ static const gsl_rng_type* get_gsl_rng_type_int(int itype)
   switch (itype) {
   case GSL_RNG_DEFAULT: T = gsl_rng_default; break;
   case GSL_RNG_MT19937: T = gsl_rng_mt19937; break; /* default */
-#ifdef GSL_1_2_LATER
   case GSL_RNG_MT19937_1999: T = gsl_rng_mt19937_1999; break;
   case GSL_RNG_MT19937_1998: T = gsl_rng_mt19937_1998; break;
   case GSL_RNG_TAUS113: T = gsl_rng_taus113; break;
   case GSL_RNG_TAUS2: T = gsl_rng_taus2; break;
-#endif
   case GSL_RNG_RANLXS0: T = gsl_rng_ranlxs0; break;
   case GSL_RNG_RANLXS1: T = gsl_rng_ranlxs1; break;
   case GSL_RNG_RANLXS2: T = gsl_rng_ranlxs2; break;
@@ -278,7 +269,6 @@ static const gsl_rng_type* get_gsl_rng_type_int(int itype)
   case GSL_RNG_UNI32: T = gsl_rng_uni32; break;
   case GSL_RNG_SLATEC: T = gsl_rng_slatec; break;
   case GSL_RNG_ZUF: T = gsl_rng_zuf; break;
-#ifdef GSL_1_1_LATER
   case GSL_RNG_BOROSH13: T = gsl_rng_borosh13; break;
   case GSL_RNG_COVEYOU: T = gsl_rng_coveyou; break;
   case GSL_RNG_FISHMAN18: T = gsl_rng_fishman18; break;
@@ -288,7 +278,6 @@ static const gsl_rng_type* get_gsl_rng_type_int(int itype)
   case GSL_RNG_KNUTHRAN2: T = gsl_rng_knuthran2; break;
   case GSL_RNG_LECUYER21: T = gsl_rng_lecuyer21; break;
   case GSL_RNG_WATERMAN14: T = gsl_rng_waterman14; break;
-#endif
 #ifdef HAVE_RNGEXTRA_RNGEXTRA_H
   case GSL_RNGEXTRA_RNG1: T = rngextra_rng1; break;
   case GSL_RNGEXTRA_RNG2: T = rngextra_rng2; break;
@@ -298,14 +287,12 @@ static const gsl_rng_type* get_gsl_rng_type_int(int itype)
     rb_raise(rb_eNotImpError, "Install the rngextra package found at <http://www.network-theory.co.uk/download/rngextra/>.");
     break;
 #endif
-#ifdef GSL_1_9_LATER
   case GSL_RNG_KNUTHRAN2002: T = gsl_rng_knuthran2002; break;
-#endif
   default:
     rb_raise(rb_eTypeError, "wrong generator type");
   }
-  return  T;
+  return T;
 }
 static void rb_gsl_rng_define_const_type(VALUE module)
@@ -405,10 +392,10 @@ static VALUE rb_gsl_rng_set(VALUE obj, VALUE s)
 /*
   Document-method: <i>GSL::Rng#get</i>
-    Returns a random integer from the generator.
-    The minimum and maximum values depend on the algorithm used,
-    but all integers in the range [min,max] are equally likely.
-    The values of min and max can determined using the auxiliary
+    Returns a random integer from the generator.
+    The minimum and maximum values depend on the algorithm used,
+    but all integers in the range [min,max] are equally likely.
+    The values of min and max can determined using the auxiliary
     methodss GSL::Rng#max and GSL::Rng#min.
 */
 static VALUE rb_gsl_rng_get(int argc, VALUE *argv, VALUE obj)
@@ -534,7 +521,6 @@ static VALUE rb_gsl_rng_print_state(VALUE obj)
   return obj;
 }
-#ifdef GSL_1_4_LATER
 static VALUE rb_gsl_rng_fwrite(VALUE obj, VALUE io)
 {
   gsl_rng *h = NULL;
@@ -558,7 +544,6 @@ static VALUE rb_gsl_rng_fread(VALUE obj, VALUE io)
   if (flag == 1) fclose(f);
   return INT2FIX(status);
 }
-#endif
 static VALUE rb_gsl_rng_memcpy(VALUE obj, VALUE dst, VALUE org)
 {
@@ -575,7 +560,7 @@ void Init_gsl_rng(VALUE module)
   cgsl_rng = rb_define_class_under(module, "Rng", cGSL_Object);
   rb_gsl_rng_define_const_type(module);
   rb_define_singleton_method(cgsl_rng, "alloc", rb_gsl_rng_alloc, -1);
   rb_define_singleton_method(cgsl_rng, "default_seed", rb_gsl_rng_default_seed, 0);
@@ -604,9 +589,7 @@ void Init_gsl_rng(VALUE module)
   rb_define_alias(cgsl_rng, "duplicate", "clone");
   rb_define_method(cgsl_rng, "print_state", rb_gsl_rng_print_state, 0);
-#ifdef GSL_1_4_LATER
   rb_define_method(cgsl_rng, "fwrite", rb_gsl_rng_fwrite, 1);
   rb_define_method(cgsl_rng, "fread", rb_gsl_rng_fread, 1);
-#endif
   rb_define_singleton_method(cgsl_rng, "memcpy", rb_gsl_rng_memcpy, 2);
 }

data/ext/{root.c → gsl_native/root.c} RENAMED

@@ -2,18 +2,17 @@
   root.c
   Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
   (C) Copyright 2004 by Yoshiki Tsunesada
   Ruby/GSL is free software: you can redistribute it and/or modify it
   under the terms of the GNU General Public License.
   This library is distributed in the hope that it will be useful, but
   WITHOUT ANY WARRANTY.
 */
-#include "rb_gsl_config.h"
-#include "rb_gsl_array.h"
-#include "rb_gsl_common.h"
-#include "rb_gsl_function.h"
-#include "rb_gsl_root.h"
+#include "include/rb_gsl_array.h"
+#include "include/rb_gsl_common.h"
+#include "include/rb_gsl_function.h"
+#include "include/rb_gsl_root.h"
 EXTERN VALUE cgsl_function_fdf;
@@ -43,8 +42,8 @@ static VALUE rb_gsl_fsolver_new(VALUE klass, VALUE t)
     } else if (!str_tail_grep(name, "brent")) {
       T = gsl_root_fsolver_brent;
     } else {
-      rb_raise(rb_eTypeError,
-	       "type must be \"bisection\" or \"falsepos\", or \"brent\".");
+      rb_raise(rb_eTypeError,
+               "type must be \"bisection\" or \"falsepos\", or \"brent\".");
     }
     break;
   case T_FIXNUM:
@@ -65,7 +64,7 @@ static VALUE rb_gsl_fsolver_new(VALUE klass, VALUE t)
     break;
   default:
     rb_raise(rb_eTypeError, "wrong argument type %s (String or Fixnum expected)",
-	     rb_class2name(CLASS_OF(t)));
+             rb_class2name(CLASS_OF(t)));
     break;
   }
   s = gsl_root_fsolver_alloc(T);
@@ -125,28 +124,28 @@ static VALUE rb_gsl_fsolver_name(VALUE obj)
 static VALUE rb_gsl_fsolver_test_interval(VALUE obj, VALUE eabs, VALUE erel)
 {
   gsl_root_fsolver *s = NULL;
-  Need_Float(eabs); Need_Float(erel);
+  Need_Float(eabs); Need_Float(erel);
   Data_Get_Struct(obj, gsl_root_fsolver, s);
   return INT2FIX(gsl_root_test_interval(s->x_lower, s->x_upper,
-					NUM2DBL(eabs), NUM2DBL(erel)));
+                                        NUM2DBL(eabs), NUM2DBL(erel)));
 }
 static VALUE rb_gsl_root_test_interval(VALUE obj, VALUE xl, VALUE xu, VALUE eabs,
-				       VALUE erel)
+                                       VALUE erel)
 {
   Need_Float(xl); Need_Float(xu);
-  Need_Float(eabs); Need_Float(erel);
+  Need_Float(eabs); Need_Float(erel);
   return INT2FIX(gsl_root_test_interval(NUM2DBL(xl), NUM2DBL(xu),
-					NUM2DBL(eabs), NUM2DBL(erel)));
+                                        NUM2DBL(eabs), NUM2DBL(erel)));
 }
 static VALUE rb_gsl_root_test_delta(VALUE obj, VALUE xl, VALUE xu, VALUE eabs,
-				       VALUE erel)
+                                    VALUE erel)
 {
   Need_Float(xl); Need_Float(xu);
-  Need_Float(eabs); Need_Float(erel);
+  Need_Float(eabs); Need_Float(erel);
   return INT2FIX(gsl_root_test_delta(NUM2DBL(xl), NUM2DBL(xu),
-				     NUM2DBL(eabs), NUM2DBL(erel)));
+                                     NUM2DBL(eabs), NUM2DBL(erel)));
 }
 static VALUE rb_gsl_root_test_residual(VALUE obj, VALUE xl,VALUE eabs)
@@ -155,7 +154,7 @@ static VALUE rb_gsl_root_test_residual(VALUE obj, VALUE xl,VALUE eabs)
   return INT2FIX(gsl_root_test_residual(NUM2DBL(xl),  NUM2DBL(eabs)));
 }
-static VALUE rb_gsl_fsolver_solve(int argc, VALUE *argv, VALUE *obj)
+static VALUE rb_gsl_fsolver_solve(int argc, VALUE *argv, VALUE obj)
 {
   gsl_root_fsolver *s = NULL;
   gsl_function *F = NULL;
@@ -166,15 +165,15 @@ static VALUE rb_gsl_fsolver_solve(int argc, VALUE *argv, VALUE *obj)
     Check_Type(argv[2], T_ARRAY);
     epsabs = NUM2DBL(rb_ary_entry(argv[2], 0));
     epsrel = NUM2DBL(rb_ary_entry(argv[2], 1));
-    /* no break */
+  /* no break */
   case 2:
     Check_Type(argv[1], T_ARRAY);
     xl = NUM2DBL(rb_ary_entry(argv[1], 0));
     xh = NUM2DBL(rb_ary_entry(argv[1], 1));
     break;
   default:
-    rb_raise(rb_eArgError,
-	     "Usage: solve(f = Function, range = Array, eps = Array)");
+    rb_raise(rb_eArgError,
+             "Usage: solve(f = Function, range = Array, eps = Array)");
     break;
   }
   CHECK_FUNCTION(argv[0]);
@@ -229,7 +228,7 @@ static VALUE rb_gsl_fdfsolver_new(VALUE klass, VALUE t)
     break;
   default:
     rb_raise(rb_eTypeError, "wrong argument type %s (String or Fixnum expected)",
-	     rb_class2name(CLASS_OF(t)));
+             rb_class2name(CLASS_OF(t)));
     break;
   }
   s = gsl_root_fdfsolver_alloc(T);
@@ -270,7 +269,7 @@ static VALUE rb_gsl_fdfsolver_name(VALUE obj)
   return rb_str_new2(gsl_root_fdfsolver_name(s));
 }
-static VALUE rb_gsl_fdfsolver_solve(int argc, VALUE *argv, VALUE *obj)
+static VALUE rb_gsl_fdfsolver_solve(int argc, VALUE *argv, VALUE obj)
 {
   gsl_root_fdfsolver *s = NULL;
   double x = 0.0, x0, epsabs = 0.0, epsrel = 1e-6;
@@ -281,7 +280,7 @@ static VALUE rb_gsl_fdfsolver_solve(int argc, VALUE *argv, VALUE *obj)
     Check_Type(argv[2], T_ARRAY);
     epsabs = NUM2DBL(rb_ary_entry(argv[2], 0));
     epsrel = NUM2DBL(rb_ary_entry(argv[2], 1));
-    /* no break */
+  /* no break */
   case 2:
     Need_Float(argv[1]);
     x0 = NUM2DBL(argv[1]);
@@ -293,7 +292,7 @@ static VALUE rb_gsl_fdfsolver_solve(int argc, VALUE *argv, VALUE *obj)
   CHECK_FUNCTION_FDF(argv[0]);
   Data_Get_Struct(argv[0], gsl_function_fdf, F);
   Data_Get_Struct(obj, gsl_root_fdfsolver, s);
-  gsl_root_fdfsolver_set(s, F, x0);
+  gsl_root_fdfsolver_set(s, F, x0);
   do {
     iter++;
     status = gsl_root_fdfsolver_iterate (s);
@@ -362,32 +361,32 @@ void Init_gsl_root(VALUE module)
   mgsl_root = rb_define_module_under(module, "Root");
   cgsl_fsolver = rb_define_class_under(mgsl_root, "FSolver", cGSL_Object);
-  rb_define_singleton_method(cgsl_fsolver, "alloc", rb_gsl_fsolver_new, 1);
-  rb_define_method(cgsl_fsolver, "set", rb_gsl_fsolver_set, 3);
-  rb_define_method(cgsl_fsolver, "iterate", rb_gsl_fsolver_iterate, 0);
-  rb_define_method(cgsl_fsolver, "root", rb_gsl_fsolver_root, 0);
-  rb_define_method(cgsl_fsolver, "name", rb_gsl_fsolver_name, 0);
-  rb_define_method(cgsl_fsolver, "x_lower", rb_gsl_fsolver_x_lower, 0);
-  rb_define_method(cgsl_fsolver, "x_upper", rb_gsl_fsolver_x_upper, 0);
-  rb_define_method(cgsl_fsolver, "test_interval", rb_gsl_fsolver_test_interval, 2);
-  rb_define_method(cgsl_fsolver, "solve", rb_gsl_fsolver_solve, -1);
-  rb_define_singleton_method(mgsl_root, "test_interval",
-			    rb_gsl_root_test_interval, 4);
-  rb_define_singleton_method(mgsl_root, "test_delta",
-			    rb_gsl_root_test_delta, 4);
-  rb_define_singleton_method(mgsl_root, "test_residual",
-			    rb_gsl_root_test_residual, 2);
+  rb_define_singleton_method(cgsl_fsolver, "alloc", rb_gsl_fsolver_new, 1);
+  rb_define_method(cgsl_fsolver, "set", rb_gsl_fsolver_set, 3);
+  rb_define_method(cgsl_fsolver, "iterate", rb_gsl_fsolver_iterate, 0);
+  rb_define_method(cgsl_fsolver, "root", rb_gsl_fsolver_root, 0);
+  rb_define_method(cgsl_fsolver, "name", rb_gsl_fsolver_name, 0);
+  rb_define_method(cgsl_fsolver, "x_lower", rb_gsl_fsolver_x_lower, 0);
+  rb_define_method(cgsl_fsolver, "x_upper", rb_gsl_fsolver_x_upper, 0);
+  rb_define_method(cgsl_fsolver, "test_interval", rb_gsl_fsolver_test_interval, 2);
+  rb_define_method(cgsl_fsolver, "solve", rb_gsl_fsolver_solve, -1);
+  rb_define_singleton_method(mgsl_root, "test_interval",
+                             rb_gsl_root_test_interval, 4);
+  rb_define_singleton_method(mgsl_root, "test_delta",
+                             rb_gsl_root_test_delta, 4);
+  rb_define_singleton_method(mgsl_root, "test_residual",
+                             rb_gsl_root_test_residual, 2);
   cgsl_fdfsolver = rb_define_class_under(mgsl_root, "FdfSolver", cGSL_Object);
-  rb_define_singleton_method(cgsl_fdfsolver, "alloc", rb_gsl_fdfsolver_new, 1);
+  rb_define_singleton_method(cgsl_fdfsolver, "alloc", rb_gsl_fdfsolver_new, 1);
-  rb_define_method(cgsl_fdfsolver, "set", rb_gsl_fdfsolver_set, 2);
-  rb_define_method(cgsl_fdfsolver, "iterate", rb_gsl_fdfsolver_iterate, 0);
-  rb_define_method(cgsl_fdfsolver, "root", rb_gsl_fdfsolver_root, 0);
-  rb_define_method(cgsl_fdfsolver, "name", rb_gsl_fdfsolver_name, 0);
-  rb_define_method(cgsl_fdfsolver, "solve", rb_gsl_fdfsolver_solve, -1);
+  rb_define_method(cgsl_fdfsolver, "set", rb_gsl_fdfsolver_set, 2);
+  rb_define_method(cgsl_fdfsolver, "iterate", rb_gsl_fdfsolver_iterate, 0);
+  rb_define_method(cgsl_fdfsolver, "root", rb_gsl_fdfsolver_root, 0);
+  rb_define_method(cgsl_fdfsolver, "name", rb_gsl_fdfsolver_name, 0);
+  rb_define_method(cgsl_fdfsolver, "solve", rb_gsl_fdfsolver_solve, -1);
   rb_define_method(cgsl_function, "fsolve", rb_gsl_function_rootfinder, -1);
   rb_define_alias(cgsl_function, "solve", "fsolve");