Contents

#!/usr/bin/env ruby
require("gsl")
include GSL
include Monte
include Math

proc_f = Proc.new { |k, dim, params|
  pi = Math::PI
  a = 1.0/(pi*pi*pi)
  a/(1.0 - cos(k[0])*cos(k[1])*cos(k[2]))
}

def display_results(title, result, error)
  exact = 1.3932039296856768591842462603255

  diff = result - exact
  printf("%s ==================\n", title);
  printf("result = % .6f\n", result);
  printf("sigma  = % .6f\n", error);
  printf("exact  = % .6f\n", exact);
  printf("error  = % .6f = %.1g sigma\n", diff, diff.abs/error)
end

dim = 3
xl = Vector.alloc([0, 0, 0])
xu = Vector.alloc([PI, PI, PI])
G = Monte::Function.alloc(proc_f, dim)
calls = 10000
r = GSL::Rng.alloc(Rng::DEFAULT)

plain = Monte::Plain.alloc(dim)

result, error = Monte::Plain::integrate(G, xl, xu, dim, calls, r, plain)
#result, error = Monte::Plain::integrate(G, xl, xu, calls, r, plain)
#result, error = Monte::Plain::integrate(G, xl, xu, calls, plain)
#result, error = Monte::Plain::integrate(G, xl, xu, dim, calls, r, "plain")

#result, error = G.integrate(xl, xu, dim, calls, r, plain)
#result, error = G.integrate(xl, xu, calls, r, plain)
#result, error = G.integrate(xl, xu, calls, plain)

#result, error = plain.integrate(G, xl, xu, dim, calls, r)
#result, error = plain.integrate(G, xl, xu, calls, r)
#result, error = plain.integrate(G, xl, xu, dim, calls)
#result, error = plain.integrate(G, xl, xu, calls)

display_results("plain", result, error)

Version data entries

22 entries across 22 versions & 4 rubygems

Version	Path
gsl-2.1.0.3	examples/monte/plain.rb
gsl-2.1.0.2	examples/monte/plain.rb
gsl-2.1.0.1	examples/monte/plain.rb
gsl-2.1.0	examples/monte/plain.rb
gsl-1.16.0.6	examples/monte/plain.rb
rb-gsl-1.16.0.5	examples/monte/plain.rb
rb-gsl-1.16.0.4	examples/monte/plain.rb
rb-gsl-1.16.0.3	examples/monte/plain.rb
rb-gsl-1.16.0.3.rc1	examples/monte/plain.rb
gsl-nmatrix-1.15.3.2	examples/monte/plain.rb
gsl-nmatrix-1.15.3.1	examples/monte/plain.rb
rb-gsl-1.16.0.2	examples/monte/plain.rb
rb-gsl-1.16.0.1	examples/monte/plain.rb
rb-gsl-1.16.0	examples/monte/plain.rb
rb-gsl-1.15.3.2	examples/monte/plain.rb
rb-gsl-1.15.3.1	examples/monte/plain.rb
blackwinter-gsl-1.15.3.2	examples/monte/plain.rb
gsl-nmatrix-1.15.3.0	examples/monte/plain.rb
gsl-1.15.3	examples/monte/plain.rb
gsl-1.14.7	examples/monte/plain.rb