Sha256: d4894957094238ef60c05d0e783492f10b406d0f5edd314e04a5b586d5e2201b

Contents?: true

Size: 1.4 KB

Versions: 22

Compression:

Stored size: 1.4 KB

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)

miser = Monte::Miser.alloc(dim)

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

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

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

display_results("miser", result, error)

Version data entries

22 entries across 22 versions & 4 rubygems

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