gsl 1.15.3 → 1.16.0.6

data/examples/multiroot/fsolver3.rb

@@ -9,7 +9,7 @@ func = GSL::MultiRoot::Function.alloc(2) { |x, f|
   f[0] = -2.0*x0*x0 + 3.0*x0*x1 + 4.0*sin(x1) - 6.0
   f[1] = 3.0*x0*x0 - 2.0*x0*x1*x1 + 3.0*cos(x0) + 4.0
 }
-                            
+
 p func.solve([1.0, 2.0].to_gv, 1000, 1e-7, "hybrids")
 p func.solve([1.0, 2.0].to_gv, 1000, "broyden")
 p func.solve([1.0, 2.0], "hybrid")

data/examples/odeiv/binarysystem.rb

@@ -2,13 +2,13 @@
 #      19/Apr/2004         by Yoshiki Tsunesada
 #
 #   This is an example to calculate the orbital evolution of
-# a double neutron star (binary) system. General relativity predicts 
-# that the binary orbital decays by radiating gravitational waves, 
-# and the two stars will coalesce in time scale of 100-1000 Mega-years. 
-#   The values used here are of the binary system J0730-3039 discovered 
-# in 2003 (Burgay et al., Nature 2003). The result shows that the two 
-# neutron stars will merge after about 85 Mega-years. From the age of 
-# the system 100 Mega-year, the lifetime of the system is estimated 
+# a double neutron star (binary) system. General relativity predicts
+# that the binary orbital decays by radiating gravitational waves,
+# and the two stars will coalesce in time scale of 100-1000 Mega-years.
+#   The values used here are of the binary system J0730-3039 discovered
+# in 2003 (Burgay et al., Nature 2003). The result shows that the two
+# neutron stars will merge after about 85 Mega-years. From the age of
+# the system 100 Mega-year, the lifetime of the system is estimated
 # about 185 Mega-years.
 #
 # References:
@@ -16,7 +16,7 @@
 #   2. Shapiro & Teukolsky, "Black holes, white dwarfs and neutron stars"
 #        John Wiley and Sans (1983)
 #
- 
+
 require("gsl")
 include Math
 
@@ -33,7 +33,7 @@ GMsolarC3 = 4.925490947e-6
 MegaYear = 3600*24*365*1e6
 
 # Time evolution of the binary orbital period and the eccentricity
-# due to gravitational radiation. 
+# due to gravitational radiation.
 # The calculation is based on general relativity (See e.g. Ref.2).
 #     y[0]: orbital period (pb)
 #     y[1]: eccentricity (e)
@@ -45,7 +45,7 @@ deriv = Proc.new { |t, y, dydt, binary|
   e = y[1]             # eccentricity
   m1 = binary.m1       # neutron star masses
   m2 = binary.m2
-  totalM = m1 + m2     # total mass                           
+  totalM = m1 + m2     # total mass
   mu = m1*m2/totalM    # reduced mass
   mm = mu*GSL::pow(totalM, 2.0/3.0)
   f_e = GSL::pow(1.0 - e*e, -3.5)*(1.0 + (73.0/24.0 + 37.0/96.0*e*e)*e*e);
@@ -85,8 +85,8 @@ tend = 2500*MegaYear
 h = 1.0*MegaYear
 
 begin
-  file = File.open("binarysystem.dat", "w") 
-  while t < tend 
+  file = File.open("binarysystem.dat", "w")
+  while t < tend
     t, h, status = solver.apply(t, tend, h, y)
     break if status != GSL::SUCCESS
     break if GSL::isnan?(y[0])

data/examples/odeiv/demo.rb

@@ -1,5 +1,5 @@
 #!/usr/bin/env ruby
-# Solve 
+# Solve
 #   dydt = -2y     ---> y(t) = exp(-2t)
 
 require("gsl")
@@ -40,11 +40,11 @@ def odeiv_compare_algorithm(solver, steptype, t0, tend, h0, y0, outfile)
 end
 
 if GSL::VERSION >= "1.5.90"
-  ALGORITHMS = ["rk2", "rk4", "rkf45", "rkck", "rk8pd", "rk2imp", "rk4imp", 
+  ALGORITHMS = ["rk2", "rk4", "rkf45", "rkck", "rk8pd", "rk2imp", "rk4imp",
                 "bsimp", "gear1", "gear2", "rk2simp"]
   gpfile = "demo2.gp"
 else
-  ALGORITHMS = ["rk2", "rk4", "rkf45", "rkck", "rk8pd", "rk2imp", "rk4imp", 
+  ALGORITHMS = ["rk2", "rk4", "rkf45", "rkck", "rk8pd", "rk2imp", "rk4imp",
                 "bsimp", "gear1", "gear2"]
   gpfile = "demo.gp"
 end

data/examples/odeiv/frei1.rb

@@ -1,13 +1,13 @@
 #!/usr/bin/env ruby
 #  Solve Schroedinger equation
 #
-#  This example is taken from frei1.cpp 
+#  This example is taken from frei1.cpp
 #  in "Numerische Physik" p201-204 (Springer),
 #  which simulates the time evolution of a probability density.
 #
 #  Name: frei1.cpp
-#  Zweck: Simuliert ein quantenmechanisches freies Teilchen 
-#  Gleichung: Schroedingergleichung ohne Potential verwendete 
+#  Zweck: Simuliert ein quantenmechanisches freies Teilchen
+#  Gleichung: Schroedingergleichung ohne Potential verwendete
 #  Bibiliothek: GSL
 #
 #  Reference:
@@ -20,7 +20,7 @@ require("gsl")
 #NMAX = 8192
 NMAX = 256
 
-# The wave equation: 
+# The wave equation:
 # calculate time derivative of the wave function.
 # The second spatial derivative is approximated by
 #   d2_psi/dx2 ~ (psi[n+1] - 2*psi[n] + pxi[n-1])/(dx*dx)
@@ -68,7 +68,7 @@ sum = 0.0
 for n in 0...NMAX do
   x = (n-NMAX/2) * $dx
   psi[n] = Math::exp(-GSL::pow_2(x/alpha)/2)
-  sum += GSL::pow_2(psi[n]) 
+  sum += GSL::pow_2(psi[n])
 end
 sum = 1.0/Math::sqrt(sum)
 
@@ -89,12 +89,12 @@ IO.popen("graph -T X -C -g 3", "w") do |io|
   c = GSL::Odeiv::Control.y_new(atol, rtol)
   evolve = GSL::Odeiv::Evolve.alloc(2*NMAX)
   sys = GSL::Odeiv::System.alloc(f, 2*NMAX)
-  
+
   t = 0.0
   for n in 1..n_out do
     t1 = n*dt
     STDOUT.printf("t = %2.1f (%2d/%2d)\n", t1-dt, n, n_out)
-    while t < t1 
+    while t < t1
       t, h, status = evolve.apply(c, step, sys, t, t1, h, psi)
       break if status != GSL::SUCCESS
     end

data/examples/odeiv/frei2.rb

@@ -1,5 +1,5 @@
 #!/usr/bin/env ruby
-# This example is taken from frei2.cpp 
+# This example is taken from frei2.cpp
 # in "Numerische Physik" p205-206 (Springer).
 #
 # Reference:
@@ -45,10 +45,10 @@ IO.popen("graph -T X -C -g 3", "w") do |io|
     io.printf("%e %e\n", x, psi[n1].abs)
   end
   io.printf("\n")
-  
+
   psi_p.radix2_forward!
-  
-  
+
+
   t = 0.0
   for n in 1..n_out do
     t1 = n*dt

data/examples/odeiv/oscillator.rb

@@ -5,7 +5,7 @@
 #   b: resist
 #   f: external force
 
-require("gsl") 
+require("gsl")
 
 dim = 2
 

data/examples/odeiv/sedov.rb

@@ -18,16 +18,16 @@ include Math
 Sedov = Proc.new { |logx, y, dydlogx, sh|
   a = -5.0*(6.0*y[2] - 15.0*y[2]*y[0]*sh + 2.0*y[0]*y[1] - 7.0*y[0]*y[0]*y[1] + 5.0*GSL::pow_3(y[0])*y[1])
   b = -25.0*y[2]*sh + 4.0*y[1] - 20.0*y[0]*y[1] + 25.0*y[0]*y[0]*y[1]
-  dydlogx[0] = a/b             
+  dydlogx[0] = a/b
 
   bb = -5.0*(6.0*y[2] - 15.0*y[2]*y[0]*sh + 2.0*y[0]*y[1] - 7.0*y[0]*y[0]*y[1])
 
   a = -5.0*(-30.0*y[2]*y[1] + 2.0*y[0]*y[1]*y[1] - 25.0*GSL::pow_2(y[0]*y[1]) + 50.0*GSL::pow_3(y[0])*GSL::pow_2(y[1]))
   bb = (-2.0 + 5.0*y[0])*(-25.0*y[2]*sh + 4.0*y[1] - 20.0*y[0]*y[1] + 25.0*y[0]*y[0]*y[1])
-  dydlogx[1] = a/bb       
+  dydlogx[1] = a/bb
 
   a = -5.0*y[2]*(-10.0*y[2]*sh + 4.0*y[1] - 14.0*y[0]*y[1] + 10.0*y[0]*y[0]*y[1] - y[0]*sh*y[1] + 10.0*y[0]*y[0]*sh*y[1])
-  dydlogx[2] = a/b             
+  dydlogx[2] = a/b
 }
 
 DIM = 3

data/examples/odeiv/whitedwarf.rb

@@ -15,8 +15,8 @@ module Degenerate
   hbar = PLANCKS_CONSTANT_HBAR    # Planck's constant
   me = MASS_ELECTRON              # Electron mass
   mn = MASS_NEUTRON               # Neutron mass
-  mu = UNIFIED_ATOMIC_MASS        
-  c = SPEED_OF_LIGHT            
+  mu = UNIFIED_ATOMIC_MASS
+  c = SPEED_OF_LIGHT
   ELambda = hbar/me/c             # Compton length of electron
   NLambda = hbar/mn/c
   MeC2 = me*c*c                   # Electron rest mass energy
@@ -29,24 +29,24 @@ module Degenerate
     tmp = sqrt(1.0 + x*x)
     (x*tmp*(2.0*x*x/3.0 - 1.0) + log(x + tmp))/8/PI/PI
   end
-  
+
   def chi(x)
     tmp = sqrt(1.0 + x*x)
     (x*tmp*(1.0 + 2*x*x) - log(x + tmp))/8/PI/PI
   end
-  
+
   def xe_rho_mue(rho, mue)
     Factor_xe*GSL::pow(rho/mue, 1.0/3.0)
   end
- 
+
   def xn_rho(rho)
     Factor_xn*GSL::pow(rho, 1.0/3.0)
   end
- 
+
 end
 
 # Polytrope gas sphere
-module Polytrope 
+module Polytrope
 
 # Lane-Emden equation
 #   n: polytrope index
@@ -59,7 +59,7 @@ module Polytrope
     dim = 2
     y = GSL::Vector[1.0, 0.0]
     dydx = GSL::Vector.alloc(dim)
-    
+
     solver = Solver.alloc(Step::RKF45, [1e-6, 0], EmdenEq, dim)
     solver.set_params(n)
     solver.reset
@@ -71,7 +71,7 @@ module Polytrope
     x = 0.0001
     xend = 10.0
     h = 1e-6
-    
+
     file = File.open("polytrope.dat", "w")
     i = 0
     while x < xend
@@ -94,10 +94,10 @@ module Polytrope
     vdy2 = vdy.subvector(0, i)
     spline = GSL::Spline.alloc(GSL::Interp::AKIMA, i)
     spline.init(vy2.reverse, vx2.reverse)
-    
+
 # Find star surface:
 # Star sufrace is defined as the zero point of density structure function
-    x1 = spline.eval(0.0)      
+    x1 = spline.eval(0.0)
     spline.init(vx2, vdy2)
     yx2 = spline.eval(x1).abs
     return [x1, yx2*x1*x1]

data/examples/permutation/ex1.rb

@@ -11,12 +11,12 @@ pp.init
 pp.fprintf(STDOUT, " %u")
 printf("\n")
 
-puts(" random permutation:");  
+puts(" random permutation:");
 r.shuffle(pp)
 pp.fprintf(STDOUT, " %u")
 printf("\n");
 
-puts("inverse permutation:");  
+puts("inverse permutation:");
 q = pp.inverse
 q.fprintf(STDOUT, " %u")
 printf ("\n");

data/examples/permutation/permutation.rb

@@ -1,5 +1,5 @@
 #!/usr/bin/env ruby
-require("gsl") 
+require("gsl")
 
 p = GSL::Permutation::alloc(10)
 p.init

data/examples/poly/demo.rb

@@ -6,7 +6,7 @@ include GSL
 poly = Poly[1.5, -1.25, -3.75, 0, 1]
 # Solve the equation p(x) == 0
 root = poly.solve    # Vector::Complex
-# Extract only the real parts 
+# Extract only the real parts
 # (imaginary parts are zero for this case)
 re = root.real       # Vector::View
 

data/examples/random/diffusion.rb

@@ -19,7 +19,7 @@ h[2] = GSL::Histogram.alloc(61, -30, 30)
 
 i = 0
 for n in [6, 50, 100] do
-  M.times do 
+  M.times do
     s = 0
     n.times do
       ds = rng.get%2 == 0 ? 1 : -1

data/examples/random/generator.rb

@@ -7,8 +7,8 @@ names = ["default", "mt19937", "mt19937_1999", "mt19937_1998", "ranlxs0", "ranlx
 "gfsr4", "rand", "random_bsd", "random8_bsd", "random32_bsd", "random64_bsd",
 "random128_bsd", "random256_bsd", "random_libc5", "random_glibc2", "rand48", "ran0",
 "ran1", "ran2", "ran3", "ranf", "ranmar", "r250", "tt800", "vax", "transputer",
-"randu", "minstd", "uni", "uni32", "slatec", "zuf", "borosh13", "coveyou", 
-"fishman18", "fishman20", "fishman2x", "knuthran2", "knuthran", "lecuyer21", 
+"randu", "minstd", "uni", "uni32", "slatec", "zuf", "borosh13", "coveyou",
+"fishman18", "fishman20", "fishman2x", "knuthran2", "knuthran", "lecuyer21",
 "waterman14"]
 
 names.each do |name|

data/examples/random/randomwalk.rb

@@ -17,7 +17,7 @@ sigma = Math::sqrt(N).to_i
 
 h = GSL::Histogram.alloc(8*sigma+1, [-4*sigma-0.5, 4*sigma+0.5])
 
-M.times do 
+M.times do
   s = 0
   N.times do
     ds = rng.get%2 == 0 ? 1 : -1
@@ -27,8 +27,8 @@ M.times do
 end
 
 x = GSL::Vector.linspace(-40, 40, 80)
-y = GSL::Ran::gaussian_pdf(x, sigma)*M*2  
-# Factor 2 is not important, but necessary 
+y = GSL::Ran::gaussian_pdf(x, sigma)*M*2
+# Factor 2 is not important, but necessary
 # because only the even ranges are filled:
 #   a + b = N     a: positive steps, b: negative steps
 #   a - b = s     s: the end point after the N steps

data/examples/random/rng.rb

@@ -18,7 +18,7 @@ p r.max
 p r.get
 p r.get
 
-r2 = GSL::Random::Rng.alloc 
+r2 = GSL::Random::Rng.alloc
 p r2.uniform
 p r2.uniform
 

data/examples/roots/bisection.rb

@@ -8,7 +8,7 @@ f = GSL::Function.alloc { |x, params|
 }
 f.set_params(1, 0, -5)
 expected = Math::sqrt(5.0)
-printf("%5s [%9s, %9s] %9s %10s %9s\n", 
+printf("%5s [%9s, %9s] %9s %10s %9s\n",
        "iter", "lower", "upper", "root",  "err", "err(est)")
 solver.set(f, 0.0, 5.0)
 iter = 0

data/examples/roots/brent.rb

@@ -19,7 +19,7 @@ f.set_params(1, 0, -5)
 expected = sqrt(5.0)
 
 printf("%5s [%9s, %9s] %9s %10s %9s\n",
-          "iter", "lower", "upper", "root", 
+          "iter", "lower", "upper", "root",
           "err", "err(est)")
 
 solver.set(f, 0.0, 5.0)

data/examples/roots/demo.rb

@@ -7,7 +7,7 @@ f = GSL::Function.alloc { |x, params|
 }
 f.set_params(1, 0, -5)
 expected = Math::sqrt(5.0)
-printf("%5s [%9s, %9s] %9s %10s %9s\n", "iter", 
+printf("%5s [%9s, %9s] %9s %10s %9s\n", "iter",
        "lower", "upper", "root", "err", "err(est)")
 solver.set(f, 0.0, 5.0)
 iter = 0; status = nil

data/examples/roots/newton.rb

@@ -2,14 +2,14 @@
 require("gsl")
 include Math
 
-f = Proc.new { |x, params| 
+f = Proc.new { |x, params|
   a = params[0]
   b = params[1]
   c = params[2]
   (a*x + b)*x + c
 }
 
-df = Proc.new { |x, params| 
+df = Proc.new { |x, params|
   a = params[0]
   b = params[1]
   2.0*a*x + b

data/examples/roots/recombination.gp

@@ -9,4 +9,3 @@ set ylabel 'Fractional ionization'
 set x2label 'Temperature of the Universe [K]'
 set grid
 plot 'recombination.dat' u 2:3 title '' w l lw 3, 1 title ''
- 

data/examples/sf/hyperg.rb

@@ -1,5 +1,5 @@
 #!/usr/bin/env ruby
-require("gsl") 
+require("gsl")
 
 p GSL::Sf.hyperg_0F1(2, 3)
 r = GSL::Sf.hyperg_0F1_e(2, 3)

data/examples/sf/sphbessel.rb

@@ -19,7 +19,7 @@ File.open("sphbessel.dat", "w") do |file|
     k0 = bessel_k0_scaled(x)
     k1 = bessel_k1_scaled(x)
     k2 = bessel_k2_scaled(x)
-    file.printf("%e %e %e %e %e %e %e %e %e %e %e %e %e %e\n", 
+    file.printf("%e %e %e %e %e %e %e %e %e %e %e %e %e %e\n",
                 x, j0, j1, j2, j3, y0, y1, y2, i0, i1, i2, k0, k1, k2)
     x += 0.1
   end

data/examples/sort/sort.rb

@@ -1,5 +1,5 @@
 #!/usr/bin/env ruby
-require("gsl") 
+require("gsl")
 
 N = 100000
 k = 5

data/examples/tamu_anova.rb

@@ -1,16 +1,16 @@
 #!/usr/bin/env ruby
 require("gsl")
 
-data = GSL::Vector[88.60,73.20,91.40,68.00,75.20,63.00,53.90,69.20, 
-50.10,71.50,44.90,59.50,40.20,56.30,38.70,31.00, 
+data = GSL::Vector[88.60,73.20,91.40,68.00,75.20,63.00,53.90,69.20,
+50.10,71.50,44.90,59.50,40.20,56.30,38.70,31.00,
 39.60,45.30,25.20,22.70]
 factor = GSL::Vector::Int[1,1,1,1,1,2,2,2,2,2,3,3,3,3,3,4,4,4,4,4]
 
 table = GSL::TAMU_ANOVA::Table.oneway(data, factor, 4)
 table.print
 
-data = GSL::Vector[45.50,45.30,45.40,44.40,44.60,43.90,44.60,44.00,44.20, 
-43.90,44.70,44.20,44.00,43.80,44.60,43.10,46.00,45.90, 
+data = GSL::Vector[45.50,45.30,45.40,44.40,44.60,43.90,44.60,44.00,44.20,
+43.90,44.70,44.20,44.00,43.80,44.60,43.10,46.00,45.90,
 44.80,46.20,45.10,45.50]
 factor = GSL::Vector::Int[1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,3,3,3,3,3,3]
 

data/examples/vector/add.rb

@@ -23,7 +23,7 @@ a += b
 a
 
 # Subtract b from a
-a -= b 
+a -= b
 
 # Show a
 a