simple-random 1.0.2 → 1.0.4

checksums.yaml

@@ -1,7 +1,7 @@
 ---
-SHA1:
-  metadata.gz: 65923da1d24e7888037400e04e8247ce5a65e223
-  data.tar.gz: c1340f384268823745fdaefe3f0a3b8a068539d7
+SHA256:
+  metadata.gz: 2b758006bb39a31380c6505ff4e9ecd52c67ef9b62ad952d63adeb7a2d49411c
+  data.tar.gz: 1ddd254e9a7cfce4b5e7d76dd19f58658501f537f4d46f2c025300d28a0aa471
 SHA512:
-  metadata.gz: 99f369a1a00047572b031eae49326f458a2afa8944c8ccfeb03f02872ea1bedfed70939631ee84ce72aee9eabdcdf062e4911ba71f41d0f8dcaf7e8f8da7d271
-  data.tar.gz: 9750c6caa2986f42a77333c87f9c6d436c49c9fc5282f77d61f8b9cb1eb6a6fe874bf097b72a623c0101fdb295ceb0093ff1c348a6c0dd01f69ec26f7e2fd428
+  metadata.gz: 89da26472354aca342ab707908b2c519389ea92d38b72157c18fcaa71f2911312e46a37f5aaa821a4562a2171ae7c220bcdc4935de8afd2a4676279ff0624ebd
+  data.tar.gz: fcd901b6f50d9925d820f4ced57687716363ea12be33bac361f8734c54b27a6005d8077009814fa7cafd9edc88fa919eb84d41611dd20e6daf7ba43ae1f51c29

data/.travis.yml

@@ -1,8 +1,22 @@
+# single test suite, non-parallel build.
+
+env:
+  global:
+    - CC_TEST_REPORTER_ID=fc1330064eb82caf0f39ecd509e3e1c5b4fc5ca8b02fb9fa63d6e8a92da5ad18
 language: ruby
 rvm:
-  - 2.2
-  - 2.1
-  - 1.9
-  - jruby-19mode
-  - rbx-2
-  - jruby-head
+  - 2.6
+  - 2.5
+  - 2.4
+  - 2.3.6
+  - 2.2.0
+before_install:
+  - gem install bundler  --version 1.17.3
+before_script:
+  - curl -L https://codeclimate.com/downloads/test-reporter/test-reporter-latest-linux-amd64 > ./cc-test-reporter
+  - chmod +x ./cc-test-reporter
+  - ./cc-test-reporter before-build
+script:
+  - bundle exec rake test
+after_script:
+  - ./cc-test-reporter after-build --exit-code $TRAVIS_TEST_RESULT

data/Gemfile

@@ -1,15 +1,11 @@
-source "http://rubygems.org"
-# Add dependencies required to use your gem here.
-# Example:
-#   gem "activesupport", ">= 2.3.5"
+source 'https://rubygems.org'
 
-# Add dependencies to develop your gem here.
-# Include everything needed to run rake, tests, features, etc.
 group :development do
-  gem "minitest", ">= 0"
-  gem "shoulda", ">= 0"
-  gem "rdoc", "~> 3.12"
-  gem "bundler", "~> 1.0"
-  gem "jeweler", "~> 2.0.1"
-  gem "simplecov", ">= 0"
+  gem 'awesome_print', '>= 1.6'
+  gem 'minitest'
+  gem 'shoulda'
+  gem 'jeweler', '~> 2.0.1'
+  gem 'simplecov'
+  gem 'nokogiri', '>= 1.8.5'
+  gem 'codeclimate-test-reporter', require: nil
 end

data/LICENSE

@@ -127,4 +127,4 @@ This License represents the complete agreement concerning subject matter hereof.
 
 10. RESPONSIBILITY FOR CLAIMS.
 
-As between Initial Developer and the Contributors, each party is responsible for claims and damages arising, directly or indirectly, out of its utilization of rights under this License and You agree to work with Initial Developer and Contributors to distribute such responsibility on an equitable basis. Nothing herein is intended or shall be deemed to constitute any admission of liability.
+As between Initial Developer and the Contributors, each party is responsible for claims and damages arising, directly or indirectly, out of its utilization of rights under this License and You agree to work with Initial Developer and Contributors to distribute such responsibility on an equitable basis. Nothing herein is intended or shall be deemed to constitute any admission of liability.

data/README.md

@@ -1,4 +1,4 @@
-# simple-random ![Build status](https://travis-ci.org/ealdent/simple-random.svg?branch=master)
+# simple-random ![Build status](https://travis-ci.com/ealdent/simple-random.svg?branch=master)
 
 Generate random numbers sampled from the following distributions:
 
@@ -34,16 +34,16 @@ Add `gem 'simple-random', '~> 1.0.0'` to your Gemfile and run `bundle install`.
 Some of the methods available:
 
 ``` ruby
-    > @sr = SimpleRandom.new # Initialize a SimpleRandom instance
+    > r = SimpleRandom.new                 # Initialize a SimpleRandom instance
      => #<SimpleRandom:0x007f9e3ad58010 @m_w=521288629, @m_z=362436069>
-    > @sr.set_seed # By default the same random seed is used, so we change it
-    > @sr.uniform(0, 5) # Produce a uniform random sample from the open interval (lower, upper).
+    > r.set_seed                           # By default the same random seed is used, so we change it
+    > r.uniform(0, 5)                      # Produce a uniform random sample from the open interval (lower, upper).
      => 0.6353204359766096
-    > @sr.normal(1000, 200) # Sample normal distribution with given mean and standard deviation
+    > r.normal(1000, 200)                  # Sample normal distribution with given mean and standard deviation
      => 862.5447157384566
-    > @sr.exponential(2) # Get exponential random sample with specified mean
+    > r.exponential(2)                     # Get exponential random sample with specified mean
      => 0.9386480625062965
-    > @sr.triangular(0, 2.5, 10) # Get triangular random sample with specified lower limit, mode, upper limit
+    > r.triangular(0, 2.5, 10)             # Get triangular random sample with specified lower limit, mode, upper limit
      => 3.1083306054169277
 ```
 
@@ -56,11 +56,12 @@ See [lib/simple-random.rb](lib/simple-random/simple_random.rb) for all available
 
 * Fork the project.
 * Make your feature addition or bug fix.
-* Add tests for it. This is important so I don't break it in a
-  future version unintentionally.
-* Commit, do not mess with rakefile, version, or history.
+* Add tests for it. This is important so I don't break it in a future version unintentionally.
+* Commit, but please do not mess with the gemspec, `Rakefile`, `VERSION`, `LICENSE`, or `.travis.yml`.
   (if you want to have your own version, that is fine but bump version in a commit by itself I can ignore when I pull)
-* Send me a pull request. Bonus points for topic branches.
+* Send me a pull request.
+
+HT: This list was modified from the default instructions that ship with [jeweler](https://github.com/technicalpickles/jeweler).
 
 ## Copyright
 
@@ -68,6 +69,14 @@ Distributed under the Code Project Open License, which is similar to MIT or BSD.
 
 ## History
 
+### 1.0.4 - 2026-02-18
+* Remove vulnerable development dependency declarations for `bundler` and `rdoc`
+* Use `https://rubygems.org` as the gem source in development
+
+### 1.0.3 - 2015-11-25
+* Attempt to reduce code complexity and improve readability
+* Change error handling somewhat to throw specific errors and improve messages
+
 ### 1.0.2 - 2015-11-24
 * Merge pull request from [cunchem](https://github.com/cunchem) to fix Laplace method
 

data/VERSION

@@ -1 +1 @@
-1.0.2
+1.0.4

data/lib/simple-random/multi_threaded_simple_random.rb

@@ -5,8 +5,7 @@ class MultiThreadedSimpleRandom < SimpleRandom
     @instances = nil
 
     def instance
-
-      unless @instances
+      unless instance_variable_defined?('@instances') && @instances
         extend MonitorMixin
 
         self.synchronize do
@@ -27,4 +26,4 @@ class MultiThreadedSimpleRandom < SimpleRandom
   end
 
   private_class_method :new
-end
+end

data/lib/simple-random/simple_random.rb

@@ -1,58 +1,67 @@
 class SimpleRandom
-  def initialize
-    @m_w = 521288629
-    @m_z = 362436069
-  end
+  class InvalidSeedArgument < StandardError; end
 
-  def set_seed(*args)
-    if args.size > 1
-      @m_w = args.first.to_i if args.first.to_i != 0
-      @m_z = args.last.to_i if args.last.to_i != 0
-    elsif args.first.is_a?(Numeric)
-      @m_z = args.first.to_i if args.first.to_i != 0
-    elsif args.first.is_a?(Time)
-      x = (args.first.to_f * 1000000).to_i
-      @m_w = x >> 16
-      @m_z = x % 4294967296 # 2 ** 32
+  I_32_BIT        = 4294967296
+  F_32_BIT        = 4294967296.0
+  DEFAULT_SEEDS   = [521288629, 362436069]
+
+  def initialize(*args)
+    if args.empty?
+      set_seed(*DEFAULT_SEEDS)
     else
-      x = (Time.now.to_f * 1000000).to_i
-      @m_w = x >> 16
-      @m_z = x % 4294967296 # 2 ** 32
+      set_seed(*args)
     end
+  end
+
+  def set_seed(*args)
+    validate_seeds!(*args)
 
-    @m_w %= 4294967296
-    @m_z %= 4294967296
+    @m_w, @m_z = determine_seeds(*args)
+
+    ensure_32bit_seeds!
+  end
+
+  def seeds=(value)
+    set_seed(*[value].flatten.compact)
+  end
+
+  def seeds
+    [@m_w, @m_z]
   end
 
   # Produce a uniform random sample from the open interval (lower, upper).
-  # The method will not return either end point.
   def uniform(lower = 0, upper = 1)
-    raise 'Invalid range' if upper <= lower
-    ((get_unsigned_int + 1) * (upper - lower) / 4294967296.0) + lower
+    fail ArgumentError, 'Upper bound must be greater than lower bound.' unless lower < upper
+
+    ((get_unsigned_int + 1) * (upper - lower) / F_32_BIT) + lower
   end
 
   # Sample normal distribution with given mean and standard deviation
   def normal(mean = 0.0, standard_deviation = 1.0)
-    raise 'Invalid standard deviation' if standard_deviation <= 0
+    fail ArgumentError, 'Standard deviation must be strictly positive' unless standard_deviation > 0
+
     mean + standard_deviation * ((-2.0 * Math.log(uniform)) ** 0.5) * Math.sin(2.0 * Math::PI * uniform)
   end
 
   # Get exponential random sample with specified mean
   def exponential(mean = 1)
-    raise 'Mean must be positive' if mean <= 0
+    fail ArgumentError, "Mean must be strictly positive" unless mean > 0
+
     -1.0 * mean * Math.log(uniform)
   end
 
   # Get triangular random sample with specified lower limit, mode, upper limit
   def triangular(lower, mode, upper)
-    raise 'Upper limit must be larger than lower limit' if upper < lower
-    raise 'Mode must lie between the upper and lower limits' if (mode < lower || mode > upper)
-    f_c = (mode - lower) / (upper - lower)
-    uniform_rand_num = uniform
-    if uniform_rand_num < f_c
-      lower + Math.sqrt(uniform_rand_num * (upper - lower) * (mode - lower))
+    fail ArgumentError, 'Upper bound must be greater than lower bound.' unless lower < upper
+    fail ArgumentError, 'Mode must lie between the upper and lower limits' if mode > upper || mode < lower
+
+    r = (upper - lower).to_f
+    u = uniform
+
+    if u < ((mode - lower) / r)
+      lower + Math.sqrt(u * r * (mode - lower))
     else
-      upper - Math.sqrt((1 - uniform_rand_num) * (upper - lower) * (upper - mode))
+      upper - Math.sqrt((1.0 - u) * r * (upper - mode))
     end
   end
 
@@ -60,28 +69,26 @@ class SimpleRandom
   # by George Marsaglia and Wai Wan Tsang.  ACM Transactions on Mathematical Software
   # Vol 26, No 3, September 2000, pages 363-372.
   def gamma(shape, scale)
-    if shape >= 1.0
-      d = shape - 1.0 / 3.0
-      c = 1 / ((9 * d) ** 0.5)
-      while true
-        v = 0.0
-        while v <= 0.0
-          x = normal
-          v = 1.0 + c * x
-        end
-        v = v ** 3
+    fail ArgumentError, 'Shape must be strictly positive' unless shape > 0
+    return scale * gamma(shape + 1.0, 1.0) * uniform ** -shape if shape < 1
+
+    d = shape - 1 / 3.0
+    c = (9 * d) ** -0.5
+
+    begin
+      z = normal
+
+      condition1 = z > (-1.0 / c)
+      condition2 = false
+
+      if condition1
         u = uniform
-        if u < (1.0 - 0.0331 * (x ** 4)) || Math.log(u) < (0.5 * (x ** 2) + d * (1.0 - v + Math.log(v)))
-          return scale * d * v
-        end
+        v = (1 + c * z) ** 3
+        condition2 = Math.log(u) < (0.5 * (z ** 2) + d * (1.0 - v + Math.log(v)))
       end
-    elsif shape <= 0.0
-      raise 'Shape must be positive'
-    else
-      g = gamma(shape + 1.0, 1.0)
-      w = uniform
-      return scale * g * (w ** (1.0 / shape))
-    end
+    end while !condition2
+
+    scale * d * v
   end
 
   def chi_square(degrees_of_freedom)
@@ -93,34 +100,36 @@ class SimpleRandom
   end
 
   def beta(a, b)
-    raise "Alpha and beta parameters must be positive. Received a = #{a} and b = #{b}." unless a > 0 && b > 0
+    fail ArgumentError, "Parameters must be strictly positive" unless a > 0 && b > 0
     u = gamma(a, 1)
     v = gamma(b, 1)
     u / (u + v)
   end
 
   def weibull(shape, scale)
-    raise 'Shape and scale must be positive' if shape <= 0.0 || scale <= 0.0
+    fail ArgumentError, 'Shape and scale must be positive' unless shape > 0 && scale > 0
 
     scale * ((-Math.log(uniform)) ** (1.0 / shape))
   end
 
   def cauchy(median, scale)
-    raise 'Scale must be positive' if scale <= 0
+    fail ArgumentError, 'Scale must be positive' unless scale > 0
 
     median + scale * Math.tan(Math::PI * (uniform - 0.5))
   end
 
   def student_t(degrees_of_freedom)
-    raise 'Degrees of freedom must be positive' if degrees_of_freedom <= 0
+    fail ArgumentError, 'Degrees of freedom must be strictly positive' unless degrees_of_freedom > 0
 
     normal / ((chi_square(degrees_of_freedom) / degrees_of_freedom) ** 0.5)
   end
 
   def laplace(mean, scale)
-    u_1 = uniform
+    u_1 = uniform(-0.5, 0.5)
     u_2 = uniform
-    mean + ((u_1 < 0.5 ? 1 : -1) * scale * Math.log( 1-2*(u_2-0.5).abs))
+
+    sign = u_1 / u_1.abs
+    mean + sign * scale * Math.log(1 - u_2)
   end
 
   def log_normal(mu, sigma)
@@ -144,65 +153,96 @@ class SimpleRandom
     ((@m_z << 16) + (@m_w & 65535)) % 4294967296
   end
 
-  def gamma_function(x)
-    g = [
-      1.0,
-      0.5772156649015329,
-      -0.6558780715202538,
-      -0.420026350340952e-1,
-      0.1665386113822915,
-      -0.421977345555443e-1,
-      -0.9621971527877e-2,
-      0.7218943246663e-2,
-      -0.11651675918591e-2,
-      -0.2152416741149e-3,
-      0.1280502823882e-3,
-      -0.201348547807e-4,
-      -0.12504934821e-5,
-      0.1133027232e-5,
-      -0.2056338417e-6,
-      0.6116095e-8,
-      0.50020075e-8,
-      -0.11812746e-8,
-      0.1043427e-9,
-      0.77823e-11,
-      -0.36968e-11,
-      0.51e-12,
-      -0.206e-13,
-      -0.54e-14,
-      0.14e-14
-    ]
-
-    r = 1.0
+  def validate_seeds!(*args)
+    return true if args.compact.empty?
+
+    unless args[0].to_f > 0
+      fail InvalidSeedArgument, 'Seeds must be strictly positive'
+    end
+
+    unless args[1].nil? || args[1].to_f > 0
+      fail InvalidSeedArgument, 'Seeds must be strictly positive'
+    end
 
+    true
+  end
+
+  def generate_temporal_seed(timestamp = Time.now)
+    x = (timestamp.to_f * 1000000).to_i
+
+    [x >> 16,  x % 4294967296]
+  end
+
+  def gamma_function(x)
     return 1e308 if x > 171.0
-    if x.is_a?(Fixnum) || x == x.to_i
-      if x > 0
-        ga = (2...x).inject(1.0) { |prod, i| prod * i }
-      else
-        1e308
-      end
+
+    if x.to_f == x.to_i
+      return unless x > 0
+      return 1 if x.to_i == 1
+
+      (1...x).inject(&:*)
     else
-      if x.abs > 1.0
-        r = (1..(x.abs.to_i)).inject(1.0) { |prod, i| prod * (x.abs - i) }
-        z = x.abs - x.abs.to_i
+      z = if x.abs > 1.0
+        x.abs - x.abs.to_i
       else
-        z = x
+        x
       end
 
-      gr = g[24]
-      23.downto(0).each do |i|
-        gr = gr * z + g[i]
+      gr = GAMMA_VALUES.inject(GAMMA_NAUGHT) do |sum, g|
+        sum * z + g
       end
-      ga = 1.0 / (gr * z)
-      if x.abs > 1
-        ga *= r
-        if x < 0
-          ga = -Math::PI / (x * ga * Math.sin(Math::PI * x))
-        end
+
+      r = if x.abs > 1
+        (1..(x.abs.to_i)).inject(1.0) { |prod, i| prod * (x.abs - i) }
+      else
+        1.0
       end
-    end
 
-    ga
+      if x < 0 && x.abs > 1
+        -Math::PI * gr * z / (x * r * Math.sin(Math::PI * x))
+      else
+        r / (gr * z)
+      end
+    end
   end
+
+  def ensure_32bit_seeds!
+    @m_w = @m_w.to_i % I_32_BIT
+    @m_z = @m_z.to_i % I_32_BIT
+  end
+
+  def determine_seeds(*args)
+    return generate_temporal_seed(args.first || Time.now) if args.empty? || args.first.respond_to?(:iso8601)
+    return [DEFAULT_SEEDS.first, args.first] if args.size < 2
+    args[0..1]
+  end
+
+  GAMMA_NAUGHT = 0.14e-14
+
+  GAMMA_VALUES = [
+    -5.4e-15,
+    -2.06e-14,
+    5.1e-13,
+    -3.6968e-12,
+    7.7823e-12,
+    1.043427e-10,
+    -1.1812746e-09,
+    5.0020075e-09,
+    6.116095e-09,
+    -2.056338417e-07,
+    1.133027232e-06,
+    -1.2504934821e-06,
+    -2.01348547807e-05,
+    0.0001280502823882,
+    -0.0002152416741149,
+    -0.0011651675918591,
+    0.007218943246663,
+    -0.009621971527877,
+    -0.0421977345555443,
+    0.1665386113822915,
+    -0.0420026350340952,
+    -0.6558780715202538,
+    0.5772156649015329,
+    1.0
+  ]
 end

data/lib/simple-random.rb

@@ -1,2 +1,2 @@
-require 'simple-random/simple_random'
-require 'simple-random/multi_threaded_simple_random'
+require File.join(File.dirname(__FILE__), 'simple-random', 'simple_random')
+require File.join(File.dirname(__FILE__), 'simple-random', 'multi_threaded_simple_random')

data/simple-random.gemspec

@@ -6,12 +6,12 @@
 
 Gem::Specification.new do |s|
   s.name = "simple-random"
-  s.version = "1.0.2"
+  s.version = "1.0.4"
 
   s.required_rubygems_version = Gem::Requirement.new(">= 0") if s.respond_to? :required_rubygems_version=
   s.require_paths = ["lib"]
   s.authors = ["John D. Cook", "Jason Adams"]
-  s.date = "2015-11-24"
+  s.date = "2026-02-18"
   s.description = "Simple Random Number Generator including Beta, Cauchy, Chi square, Exponential, Gamma, Inverse Gamma, Laplace (double exponential), Normal, Student t, Uniform, and Weibull.  Ported from John D. Cook's C# Code."
   s.email = "jasonmadams@gmail.com"
   s.extra_rdoc_files = [
@@ -22,7 +22,6 @@ Gem::Specification.new do |s|
     ".document",
     ".travis.yml",
     "Gemfile",
-    "Gemfile.lock",
     "LICENSE",
     "README.md",
     "Rakefile",
@@ -45,23 +44,17 @@ Gem::Specification.new do |s|
     if Gem::Version.new(Gem::VERSION) >= Gem::Version.new('1.2.0') then
       s.add_development_dependency(%q<minitest>, [">= 0"])
       s.add_development_dependency(%q<shoulda>, [">= 0"])
-      s.add_development_dependency(%q<rdoc>, ["~> 3.12"])
-      s.add_development_dependency(%q<bundler>, ["~> 1.0"])
       s.add_development_dependency(%q<jeweler>, ["~> 2.0.1"])
       s.add_development_dependency(%q<simplecov>, [">= 0"])
     else
       s.add_dependency(%q<minitest>, [">= 0"])
       s.add_dependency(%q<shoulda>, [">= 0"])
-      s.add_dependency(%q<rdoc>, ["~> 3.12"])
-      s.add_dependency(%q<bundler>, ["~> 1.0"])
       s.add_dependency(%q<jeweler>, ["~> 2.0.1"])
       s.add_dependency(%q<simplecov>, [">= 0"])
     end
   else
     s.add_dependency(%q<minitest>, [">= 0"])
     s.add_dependency(%q<shoulda>, [">= 0"])
-    s.add_dependency(%q<rdoc>, ["~> 3.12"])
-    s.add_dependency(%q<bundler>, ["~> 1.0"])
     s.add_dependency(%q<jeweler>, ["~> 2.0.1"])
     s.add_dependency(%q<simplecov>, [">= 0"])
   end

data/test/helper.rb

@@ -1,4 +1,5 @@
 require 'simplecov'
+require 'awesome_print'
 
 module SimpleCov::Configuration
   def clean_filters
@@ -33,4 +34,4 @@ require 'simple-random'
 class MiniTest::Test
 end
 
-MiniTest.autorun
+MiniTest.autorun

data/test/test_simple_random.rb

@@ -12,7 +12,7 @@ end
 class Array
   def mean
     if size > 0
-      inject(0.0) { |sum, i| sum + i } / size.to_f
+      inject(&:+) / size.to_f
     else
       0.0
     end
@@ -28,102 +28,72 @@ class Array
   end
 end
 
-class TestSimpleRandom < MiniTest::Test
-  context "A simple random number generator" do
-    setup do
-      @r = SimpleRandom.new
-    end
+def Time.now
+  new(2015, 11, 26, 12, 1, 15, '-05:00')
+end
 
-    should "generate random numbers from a uniform distribution in the interval (0, 1)" do
-      SAMPLE_SIZE.times do
-        u = @r.uniform
-        assert u < 1
-        assert u > 0
+class TestSimpleRandom < MiniTest::Test
+  context "Setting the seeds for a simple random number generator" do
+    context "on initialization" do
+      should "assign default seeds when none are specified" do
+        r = SimpleRandom.new
+        assert r.seeds == SimpleRandom::DEFAULT_SEEDS
       end
-    end
 
-    should "generate uniformly random numbers with mean approximately 0.5" do
-      numbers = generate_numbers(@r, :uniform)
-      epsilon = (0.5 - numbers.mean).abs
-
-      assert epsilon < MAXIMUM_EPSILON
-    end
-
-    should "generate random numbers from a normal distribution with mean approximately 0" do
-      numbers = generate_numbers(@r, :normal)
-      epsilon = (0.0 - numbers.mean).abs
-
-      assert epsilon < MAXIMUM_EPSILON
-    end
-
-    should "generate random numbers from a normal distribution with sample standard deviation approximately 1" do
-      numbers = generate_numbers(@r, :normal)
-      epsilon = (1.0 - numbers.standard_deviation).abs
-
-      assert epsilon < MAXIMUM_EPSILON
-    end
-
-    should "generate random numbers from an exponential distribution with mean approximately 1" do
-      numbers = generate_numbers(@r, :exponential)
-      epsilon = (1.0 - numbers.mean).abs
-
-      assert epsilon < MAXIMUM_EPSILON
-    end
+      should "assign the seeds specified in the initializer" do
+        r = SimpleRandom.new(1, 2)
+        assert r.seeds == [1, 2]
+      end
 
-    should "generate random numbers from triangular(0, 1, 1) in the range [0, 1]" do
-      SAMPLE_SIZE.times do
-        t = @r.triangular(0.0, 1.0, 1.0)
-        assert t <= 1.0
-        assert t >= 0.0
+      should "reject negative seed values" do
+        assert_raises SimpleRandom::InvalidSeedArgument do
+          SimpleRandom.new(-1, 3)
+        end
       end
     end
 
-    should "generate random numbers from triangular(0, 1, 1) with mean approximately 0.66" do
-      a = 0.0
-      c = 1.0
-      b = 1.0
-      numbers = generate_numbers(@r, :triangular, a, c, b)
-      mean = (a + b + c) / 3
-      epsilon = (mean - numbers.mean).abs
-
-      assert epsilon < MAXIMUM_EPSILON
-    end
+    context 'after initialization' do
+      setup do
+        @r = SimpleRandom.new
+      end
 
-    should "generate random numbers from triangular(0, 1, 1) with standard deviation approximately 0.23" do
-      a = 0.0
-      c = 1.0
-      b = 1.0
-      numbers = generate_numbers(@r, :triangular, a, c, b)
-      std_dev = Math.sqrt((a**2 + b**2 + c**2 - a*b - a*c - b*c) / 18)
-      epsilon = (std_dev - numbers.standard_deviation).abs
+      should "accept a single value and leave the first seed the same" do
+        @r.seeds = 1
+        assert @r.seeds == [521288629, 1]
+      end
 
-      assert epsilon < MAXIMUM_EPSILON
-    end
+      should "update the seeds when given an array of values" do
+        @r.seeds = [1, 2]
+        assert @r.seeds == [1, 2]
+      end
 
-    should "generate a random number sampled from a gamma distribution" do
-      assert @r.gamma(5, 2.3)
-    end
+      should "accept a timestamp instead of an numeric value" do
+        @r.seeds = Time.parse('2015-01-01T00:00:00.000-0500')
+        assert @r.seeds == [193992865, 413250560]
+      end
 
-    should "generate a random number sampled from an inverse gamma distribution" do
-      assert @r.inverse_gamma(5, 2.3)
+      should "use the current timestamp when nothing is specified" do
+        @r.seeds = nil
+        assert @r.seeds == [628393424, 2245012672]
+      end
     end
 
-    should "generate a random number sampled from a beta distribution" do
-      assert @r.beta(5, 2.3)
-    end
+    should "provide different results with different integer seeds" do
+      r1 = SimpleRandom.new
+      r1.set_seed(2)
+      r2 = SimpleRandom.new
+      r2.set_seed(1234512343214134)
 
-    should "generate a random number sampled from a chi-square distribution" do
-      assert @r.chi_square(10)
-    end
+      r1_randoms = 100.times.map { r1.uniform(0, 10).floor }
+      r2_randoms = 100.times.map { r2.uniform(0, 10).floor }
 
-    should "generate a random number using weibull" do
-      assert @r.weibull(5, 2.3)
+      assert r1_randoms != r2_randoms
     end
   end
 
-  context "A multi-threaded simple random number generator" do
+  context "A simple random number generator" do
     setup do
-      @r = MultiThreadedSimpleRandom.instance
+      @r = SimpleRandom.new
     end
 
     should "generate random numbers from a uniform distribution in the interval (0, 1)" do
@@ -159,7 +129,7 @@ class TestSimpleRandom < MiniTest::Test
       numbers = generate_numbers(@r, :exponential)
       epsilon = (1.0 - numbers.mean).abs
 
-      assert epsilon < MAXIMUM_EPSILON * 2
+      assert epsilon < MAXIMUM_EPSILON
     end
 
     should "generate random numbers from triangular(0, 1, 1) in the range [0, 1]" do
@@ -183,30 +153,37 @@ class TestSimpleRandom < MiniTest::Test
 
     should "generate random numbers from triangular(0, 1, 1) with standard deviation approximately 0.23" do
       a = 0.0
-      c = 1.0
       b = 1.0
-      numbers = generate_numbers(@r, :triangular, a, c, b)
-      std_dev = Math.sqrt((a**2 + b**2 + c**2 - a*b - a*c - b*c) / 18)
+      c = 1.0
+      numbers = generate_numbers(@r, :triangular, a, b, c)
+      std_dev = Math.sqrt((a ** 2 + b ** 2 + c ** 2 - a * b - a * c - b * c) / 18)
       epsilon = (std_dev - numbers.standard_deviation).abs
 
       assert epsilon < MAXIMUM_EPSILON
     end
 
-    should "generate random numbers from laplace(0, 1) with mean approximately 0" do
-      mean = 0.0
-      scale = 0.1
-      numbers = generate_numbers(@r, :laplace, mean, scale)
+    should "generate random numbers from triangular(0, 0.5, 1) with mean approximately 0.5" do
+      a = 0.0
+      b = 0.5
+      c = 1.0
+
+      numbers = generate_numbers(@r, :triangular, a, b, c)
+      mean = (a + b + c) / 3
       epsilon = (mean - numbers.mean).abs
 
       assert epsilon < MAXIMUM_EPSILON
     end
-    
+
     should "generate a random number sampled from a gamma distribution" do
       assert @r.gamma(5, 2.3)
+      assert @r.gamma(5.3, 2.7)
+      assert @r.gamma(2.3, 2)
     end
 
     should "generate a random number sampled from an inverse gamma distribution" do
       assert @r.inverse_gamma(5, 2.3)
+      assert @r.inverse_gamma(5.7, 2.8)
+      assert @r.inverse_gamma(3.2, 2)
     end
 
     should "generate a random number sampled from a beta distribution" do
@@ -221,6 +198,25 @@ class TestSimpleRandom < MiniTest::Test
       assert @r.weibull(5, 2.3)
     end
 
+    should "generate random number from a dirichlet distribution" do
+      assert @r.dirichlet(5.3, 2.7)
+    end
+
+    should "generate random numbers from laplace(0, 1) with mean approximately 0" do
+      mean = 0.0
+      scale = 0.1
+      numbers = generate_numbers(@r, :laplace, mean, scale)
+      epsilon = (mean - numbers.mean).abs
+
+      assert epsilon < MAXIMUM_EPSILON
+    end
+  end
+
+  context "A multi-threaded simple random number generator" do
+    setup do
+      @r = MultiThreadedSimpleRandom.instance
+    end
+
     should "work independently in every thread" do
       sample_count = 10
       thread_count = 10
@@ -241,17 +237,5 @@ class TestSimpleRandom < MiniTest::Test
       assert samples.size == thread_count
       assert samples.uniq.size == 1
     end
-
-    should "provide different results with different integer seeds" do
-      r1 = SimpleRandom.new
-      r1.set_seed(2)
-      r2 = SimpleRandom.new
-      r2.set_seed(1234512343214134)
-
-      r1_randoms = 100.times.map { r1.uniform(0, 10).floor }
-      r2_randoms = 100.times.map { r2.uniform(0, 10).floor }
-
-      assert r1_randoms != r2_randoms
-    end
   end
 end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: simple-random
 version: !ruby/object:Gem::Version
-  version: 1.0.2
+  version: 1.0.4
 platform: ruby
 authors:
 - John D. Cook
@@ -9,7 +9,7 @@ authors:
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2015-11-24 00:00:00.000000000 Z
+date: 2026-02-18 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: minitest
@@ -39,34 +39,6 @@ dependencies:
     - - ">="
       - !ruby/object:Gem::Version
         version: '0'
-- !ruby/object:Gem::Dependency
-  name: rdoc
-  requirement: !ruby/object:Gem::Requirement
-    requirements:
-    - - "~>"
-      - !ruby/object:Gem::Version
-        version: '3.12'
-  type: :development
-  prerelease: false
-  version_requirements: !ruby/object:Gem::Requirement
-    requirements:
-    - - "~>"
-      - !ruby/object:Gem::Version
-        version: '3.12'
-- !ruby/object:Gem::Dependency
-  name: bundler
-  requirement: !ruby/object:Gem::Requirement
-    requirements:
-    - - "~>"
-      - !ruby/object:Gem::Version
-        version: '1.0'
-  type: :development
-  prerelease: false
-  version_requirements: !ruby/object:Gem::Requirement
-    requirements:
-    - - "~>"
-      - !ruby/object:Gem::Version
-        version: '1.0'
 - !ruby/object:Gem::Dependency
   name: jeweler
   requirement: !ruby/object:Gem::Requirement
@@ -108,7 +80,6 @@ files:
 - ".document"
 - ".travis.yml"
 - Gemfile
-- Gemfile.lock
 - LICENSE
 - README.md
 - Rakefile
@@ -138,8 +109,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubyforge_project: 
-rubygems_version: 2.4.5.1
+rubygems_version: 3.0.3.1
 signing_key: 
 specification_version: 4
 summary: Simple Random Number Generator

data/Gemfile.lock

@@ -1,81 +0,0 @@
-GEM
-  remote: http://rubygems.org/
-  specs:
-    activesupport (4.1.4)
-      i18n (~> 0.6, >= 0.6.9)
-      json (~> 1.7, >= 1.7.7)
-      minitest (~> 5.1)
-      thread_safe (~> 0.1)
-      tzinfo (~> 1.1)
-    addressable (2.3.6)
-    builder (3.2.2)
-    descendants_tracker (0.0.4)
-      thread_safe (~> 0.3, >= 0.3.1)
-    docile (1.1.5)
-    faraday (0.9.0)
-      multipart-post (>= 1.2, < 3)
-    git (1.2.7)
-    github_api (0.11.3)
-      addressable (~> 2.3)
-      descendants_tracker (~> 0.0.1)
-      faraday (~> 0.8, < 0.10)
-      hashie (>= 1.2)
-      multi_json (>= 1.7.5, < 2.0)
-      nokogiri (~> 1.6.0)
-      oauth2
-    hashie (3.1.0)
-    highline (1.6.21)
-    i18n (0.6.9)
-    jeweler (2.0.1)
-      builder
-      bundler (>= 1.0)
-      git (>= 1.2.5)
-      github_api
-      highline (>= 1.6.15)
-      nokogiri (>= 1.5.10)
-      rake
-      rdoc
-    json (1.8.1)
-    jwt (1.0.0)
-    mini_portile (0.6.0)
-    minitest (5.4.0)
-    multi_json (1.10.1)
-    multi_xml (0.5.5)
-    multipart-post (2.0.0)
-    nokogiri (1.6.2.1)
-      mini_portile (= 0.6.0)
-    oauth2 (0.9.4)
-      faraday (>= 0.8, < 0.10)
-      jwt (~> 1.0)
-      multi_json (~> 1.3)
-      multi_xml (~> 0.5)
-      rack (~> 1.2)
-    rack (1.5.2)
-    rake (10.3.2)
-    rdoc (3.12.2)
-      json (~> 1.4)
-    shoulda (3.5.0)
-      shoulda-context (~> 1.0, >= 1.0.1)
-      shoulda-matchers (>= 1.4.1, < 3.0)
-    shoulda-context (1.2.1)
-    shoulda-matchers (2.6.1)
-      activesupport (>= 3.0.0)
-    simplecov (0.8.2)
-      docile (~> 1.1.0)
-      multi_json
-      simplecov-html (~> 0.8.0)
-    simplecov-html (0.8.0)
-    thread_safe (0.3.4)
-    tzinfo (1.2.1)
-      thread_safe (~> 0.1)
-
-PLATFORMS
-  ruby
-
-DEPENDENCIES
-  bundler (~> 1.0)
-  jeweler (~> 2.0.1)
-  minitest
-  rdoc (~> 3.12)
-  shoulda
-  simplecov