fibonacci_rng 0.0.3 → 0.0.4

checksums.yaml

@@ -1,15 +1,15 @@
 ---
 !binary "U0hBMQ==":
   metadata.gz: !binary |-
-    NTI1ZTdjMjg4MTMyN2Q0YjBiNjg4NzViYThlZGExOGM4ZmFkYjc4MQ==
+    NTJlYjNlYzYxMGMxMzAyOTk2MmEyYWM2OWI5N2RlN2JhMWQ5ZDQ3OA==
   data.tar.gz: !binary |-
-    MGM3NTk5MDk0MDM5ZWE1MGZiYmY3NTM3NjAxZDRiZGExMWFlZjI1OA==
+    OGQxYWNhZjViOTA5MjYzNDI0OTNkZTgxMTllNmYyOWRkNmQ5MDUwMA==
 SHA512:
   metadata.gz: !binary |-
-    NTdmMTkzNzQ5OGU2YThkNzVlN2NmYjJkYjU5YjBjZTE2NTg2M2FmOTllOWI5
-    MjdlMzZlY2UwODMwODkxYjYyNWNlZTZhNTkwODYzNWI3ZDA0OTM2ZDgyNWVm
-    Y2VkYzA0YzQ3MzFlZTg2NDRjZmIyNWE4NmFlZmU0NmFiMzFiYWY=
+    ZTlhNWYzYzRjOTQ0OTljMjE5ZGM4NjA1NTZlNDkxNTc4ZDNjYmE3NGJiM2Yy
+    YzFjZTA2YzM1ZWJhZDdiMWNlZWIxY2YyYWU1ODNiNjY4ZWVmNWZhOTQyNmE2
+    NjI1N2RjMzZlNWY1NThkMWZjZTJhMTRjMDgyN2QwMjhhYmQxZmY=
   data.tar.gz: !binary |-
-    YWQwOGZkZWMzNzZjNWY0MGVhNDE1Mzk2ZmZhNzM5MDlkZGNmYjY0MmY2ODNi
-    NDhhNDNjNDQ1NjMxYzRiMTVkZGM5ZmU0YWRmYTk1YmQzNzkwYjMyZDNmODlm
-    MTFiOWJkYjE4YzJmOTFkMGU1OGZkM2Y5YTljMDZhMDMyMGY1NDc=
+    ZmU1MGRmYzE4OTk2MTUyOGUzNTExNjMwNzdjYWFkM2ZkYTIxNmZlZmQxOGJm
+    MTkxY2NjYzIyNzYxNzBjN2M5NzMwMmM5YjY5M2FmMDVhNjZjNGU5NDMzOTVm
+    MDM4ZWIwNThhNzJhN2NhNGViYzQ1MjY4NzUwMThmYmM1NjdlZGM=

data/lib/fibonacci_rng/version.rb

@@ -2,5 +2,5 @@
 
 #The class of Fibonacci inspired random number generators.
 class FibonacciRng
-  VERSION = "0.0.3"
+  VERSION = "0.0.4"
 end

data/lib/fibonacci_rng.rb

@@ -13,10 +13,10 @@ class FibonacciRng
     Time.now.to_s + @tickle.succ!
   end
 
-  CHOP = 0x3FFFFFFF
+  CHOP = 0x1FFFFFFF
   BYTE = 0xFF
   WORD = 0xFFFF
-  BASE = 1073741824.0
+  BASE = (CHOP+1).to_f
 
   #An accessor for the depth!
   attr_reader :depth
@@ -38,7 +38,12 @@ class FibonacciRng
 
   #Roll a dice.
   def dice(sides)
-    do_spin
+    limit = ((CHOP+1) / sides) * sides
+
+    begin
+      do_spin
+    end until @buffer[0] < limit
+
     @buffer[0] % sides
   end
 

data/prng.rb

@@ -0,0 +1,79 @@
+# coding: utf-8
+
+#The Pseudo Random Number Generation Test Suite.
+module PrngTestSuite
+
+  #Display usage info and exit
+  def self.usage(msg=nil)
+    puts
+    puts "PrngTestSuite Version 0.0.2"
+
+    if msg
+      puts
+      puts msg
+    end
+
+    puts
+    puts "Usage Info:"
+    puts
+    puts "$ ruby prng.rb <locn> <gen> <test> <max> <count>"
+    puts
+    puts "Where:"
+    puts "  <locn>  fib code location, either 'gem' or 'local'"
+    puts "  <gen>   select a generator, either 'fib' or 'int'"
+    puts "  <test>  select a test, either 'chisq' or 'auto'"
+    puts "  <max>   the number of bins to test    (2..65535)"
+    puts "  <count> the number of samples to test (1..1000000)"
+    puts
+
+    exit
+  end
+
+end
+
+if ARGV[0] == 'gem'
+  require 'fibonacci_rng'
+elsif ARGV[0] == 'local'
+  require_relative 'lib/fibonacci_rng'
+elsif ARGV[0] == 'help'
+  PrngTestSuite.usage
+else
+  PrngTestSuite.usage "Error: missing or invalid locn parameter."
+end
+
+require_relative 'tools/internal_rng'
+require_relative 'tools/chi_squared'
+require_relative 'tools/auto_corr'
+
+module PrngTestSuite
+
+  def self.main
+    if ARGV[1] == 'fib'
+      gen = FibonacciRng.new
+    elsif ARGV[1] == 'int'
+      gen = InternalRng.new
+    else
+      PrngTestSuite.usage "Error: missing or invalid gen parameter."
+    end
+
+    if ARGV[2] == 'chisq'
+      tester = ChiSquaredTest.new
+    elsif ARGV[2] == 'auto'
+      tester = AutoCorrelator.new
+    else
+      PrngTestSuite.usage "Error: missing or invalid test parameter."
+    end
+
+    unless (2..65535) === (max = ARGV[3].to_i)
+      PrngTestSuite.usage "Error: missing or invalid max parameter."
+    end
+
+    unless (1..10000000) === (count = ARGV[4].to_i)
+      PrngTestSuite.usage "Error: missing or invalid count parameter."
+    end
+
+    tester.run(gen, max, count)
+  end
+end
+
+PrngTestSuite.main

data/tests/fibinacci_rng_tests.rb

@@ -44,6 +44,16 @@ class CoreTester < MiniTest::Unit::TestCase
 
   end
 
+  def test_that_it_creates_floats
+    prng = FibonacciRng.new
+
+    100.times do
+      value = prng.float
+      assert((value >= 0.0) && (value < 1.0))
+    end
+
+  end
+
   def test_that_it_makes_unique_sequnces
     prnga = FibonacciRng.new
     prngb = FibonacciRng.new

data/tools/auto_corr.rb

@@ -0,0 +1,29 @@
+# coding: utf-8
+
+class AutoCorrelator
+
+  #Run a auto correlation test on the random number generator.
+  def run(gen, _width, count)
+    data   = Array.new(2*count) { gen.float - 0.5}
+    result = Array.new(count)
+
+    (0..count).each do |lag|
+      result[lag] = correlate(data, count, lag)
+    end
+
+    puts result
+
+  end
+
+  #Perform a single correlation
+  def correlate(data, count, lag)
+    result = 0.0
+
+    (0...count).each do |idx|
+      result += data[idx] * data[idx+lag]
+    end
+
+    result
+  end
+
+end

data/tools/chi_squared.rb

@@ -0,0 +1,35 @@
+# coding: utf-8
+
+#A Chi Squared tester
+class ChiSquaredTest
+
+  #Run a chi squared test on the random number generator.
+  def run(gen, bin_size, count)
+    max  = bin_size
+    bins = Array.new(bin_size, 0)
+
+    start_time = Time.now
+    puts "Starting test"
+
+    count.times do
+      bins[gen.dice(max)] += 1
+    end
+
+    puts "Test completed (#{Time.now - start_time} elapsed)"
+    puts "Raw Bins = #{bins.inspect}"
+
+    normalized = bins.collect {|value| value.to_f/count }
+
+    #puts "Nrm Bins = #{normalized.inspect}"
+
+    expected = 1/max.to_f
+    #puts "Expected = #{expected}"
+
+    err_sqr = normalized.collect{|value| diff = expected - value; diff*diff}
+    #puts "Esq Bins = #{err_sqr.inspect}"
+
+    chi_sq = err_sqr.inject {|sum, value| sum + value}
+    puts "Chi Squared = #{chi_sq}"
+  end
+
+end

data/tools/internal_rng.rb

@@ -0,0 +1,16 @@
+# coding: utf-8
+
+#A wrapper for the internal random number generator
+class InternalRng
+
+  #Roll the dice!
+  def dice(sides)
+    rand(sides)
+  end
+
+  #Get a float!
+  def float
+    rand(0)
+  end
+
+end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: fibonacci_rng
 version: !ruby/object:Gem::Version
-  version: 0.0.3
+  version: 0.0.4
 platform: ruby
 authors:
 - Peter Camilleri
@@ -82,8 +82,12 @@ files:
 - fibonacci_rng.gemspec
 - lib/fibonacci_rng.rb
 - lib/fibonacci_rng/version.rb
+- prng.rb
 - reek.txt
 - tests/fibinacci_rng_tests.rb
+- tools/auto_corr.rb
+- tools/chi_squared.rb
+- tools/internal_rng.rb
 homepage: ''
 licenses:
 - MIT