distribution 0.7.3 → 0.8.0

data/spec/uniform_spec.rb

@@ -0,0 +1,154 @@
+require File.expand_path(File.dirname(__FILE__) + "/spec_helper.rb")
+require 'distribution/uniform'
+
+describe Distribution::Uniform do
+
+  shared_examples_for "uniform engine" do
+
+    it ".rng should generate sequences with the right mean & variance" do
+      skip("This method is very innacurrate due to the low convergence rate")
+      # samples = 100_000
+      # sum = 0
+      # ss = 0
+      # lower = 0
+      # upper = 20
+
+      # # Expectations
+      # exp_mean = (upper + lower) / 2
+      # exp_variance = ((upper - lower) ** 2) / 12
+      # rng = @engine.rng(lower, upper)
+
+      # samples.times do
+      #   v = rng.call
+      #   sum += v
+      #   ss += (v - exp_mean) ** 2
+      # end
+
+      # mean = sum.to_f / samples
+      # variance = ss.to_f / samples
+      # mean.should be_within(1e-5).of(exp_mean)
+      # variance.should be_within(1e-5).of(exp_variance)
+    end
+
+    it ".rng with a specified seed should be reproducible" do
+      seed = Random.new_seed.modulo 100000007
+      gen_a = @engine.rng(0, 1, seed)
+      gen_b = @engine.rng(0, 1, seed)
+
+      expect((gen_a.call)).to eq(gen_b.call)
+    end
+
+    it ".pdf should return correct pdf for values within the defined range" do
+      if @engine.respond_to? :pdf
+        10.times do
+          low, width = rand, rand
+          x = low + rand * width
+          epdf = 1.0 / width
+          expect(@engine.pdf(x, low, low + width)).to be_within(1e-10).of(epdf)
+        end
+
+      else
+        pending("No #{@engine}.pdf")
+      end
+    end
+
+    it ".pdf should return 0 for values outside the defined range" do
+      if @engine.respond_to? :pdf
+        10.times do
+          low, width = rand, rand
+          # x lies just outside of  where the pdf exists as a non-zero value
+          # A small amount (1e-10) is removed from bad_x to ensure no overlap
+          x = low - (1 + rand) * width - 1e-10
+          expect(@engine.pdf(x, low, low + width)).to be_within(1e-10).of(0.0)
+        end
+
+      else
+        pending("No #{@engine}.pdf")
+      end
+    end
+
+
+    it ".cdf should return 0 for values smaller than the lower bound" do
+      if @engine.respond_to? :cdf
+        low, width = rand, rand
+        x = low - rand * width
+        expect(@engine.cdf(x, low, low + width)).to be_within(1e-10).of(0.0)
+      else
+        pending("No #{@engine}.cdf")
+      end
+    end
+
+    it ".cdf should return correct cdf for x within defined range" do
+      if @engine.respond_to? :cdf
+        low, width = rand, rand
+        x = low + rand * width
+        ecdf = (x - low) / width
+        expect(@engine.cdf(x, low, low + width)).to be_within(1e-10).of(ecdf)
+      else
+        pending("No #{@engine}.cdf")
+      end
+    end
+
+    it ".cdf should return 1 for values greater than the upper bound" do
+      if @engine.respond_to? :cdf
+        low, width = rand, rand
+        x = low + (1 + rand) * (width)
+        expect(@engine.cdf(x, low, low + width)).to be_within(1e-10).of(1.0)
+      else
+        pending("No #{@engine}.cdf")
+      end
+    end
+    
+    it ".quantile should return correct inverse cdf" do
+      if @engine.respond_to? :quantile
+        low, width = rand, rand
+        scale = rand
+        x = low + scale * width
+        qn = (x - low) / width
+        expect(@engine.quantile(qn, low, low + width)).to be_within(1e-10).of(x)
+      else
+        pending("No #{@engine}.quantile")
+      end
+    end
+    
+    it ".p_value should return same result as .quantile" do
+      if @engine.respond_to? :p_value and @engine.respond_to? :quantile
+        low, width = rand, rand
+        scale = rand
+        x = low + scale * width
+        qn = (x - low) / width
+        
+        expect(@engine.quantile(qn, low, low + width)).to eq(@engine.p_value(qn, low, low + width))
+      else
+        pending("No #{@engine}.p_value")
+      end
+    end
+  end
+
+
+  describe "singleton" do
+    before do
+      @engine = Distribution::Uniform
+    end
+    it_should_behave_like "uniform engine"
+  end
+
+  describe Distribution::Uniform::Ruby_ do
+    before do
+      @engine = Distribution::Uniform::Ruby_
+    end
+    it_should_behave_like "uniform engine"
+
+  end
+
+  if Distribution.has_gsl?
+    describe Distribution::Uniform::GSL_ do
+      before do
+        @engine = Distribution::Uniform::GSL_
+      end
+      it_should_behave_like "uniform engine"
+    end
+  end
+
+
+end

data/spec/weibull_spec.rb

@@ -1,18 +1,17 @@
-require File.expand_path(File.dirname(__FILE__)+"/spec_helper.rb")
+require File.expand_path(File.dirname(__FILE__) + '/spec_helper.rb')
 
 include ExampleWithGSL
 
 describe Distribution::Weibull do
-
-	 #shared_examples for "Weibull engine" do
-	        it "should return correct pdf" do
-		   Distribution::Weibull::pdf(2, 1, 1).should be_within(1e-3).of(0.13533)
-		end
-	 	it "should return correct cdf" do 
-		   Distribution::Weibull::cdf(2, 1, 1).should be_within(1e-3).of(0.86466)
-		end
-		it "should return correct p-value" do 
-		   Distribution::Weibull::p_value(0.86466, 1, 1).should be_within(1e-3).of(2.0)
-		end
-	 #end
-end
+  # shared_examples for "Weibull engine" do
+  it 'should return correct pdf' do
+    expect(Distribution::Weibull.pdf(2, 1, 1)).to be_within(1e-3).of(0.13533)
+  end
+  it 'should return correct cdf' do
+    expect(Distribution::Weibull.cdf(2, 1, 1)).to be_within(1e-3).of(0.86466)
+  end
+  it 'should return correct p-value' do
+    expect(Distribution::Weibull.p_value(0.86466, 1, 1)).to be_within(1e-3).of(2.0)
+  end
+  # end
+end

metadata

@@ -1,14 +1,15 @@
 --- !ruby/object:Gem::Specification
 name: distribution
 version: !ruby/object:Gem::Version
-  version: 0.7.3
+  version: 0.8.0
 platform: ruby
 authors:
 - Claudio Bustos
-autorequire: 
+- Carlos Agarie
+autorequire:
 bindir: bin
 cert_chain: []
-date: 2015-05-28 00:00:00.000000000 Z
+date: 2020-07-05 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
@@ -52,9 +53,11 @@ dependencies:
     - - "~>"
       - !ruby/object:Gem::Version
         version: '3.2'
-description: Distribution
+description: Distribution is a gem with several probabilistic distributions. Pure
+  Ruby is used by default, C (GSL) or Java extensions are used if available
 email:
 - clbustos@gmail.com
+- carlos.agarie@gmail.com
 executables:
 - distribution
 extensions: []
@@ -62,6 +65,7 @@ extra_rdoc_files: []
 files:
 - ".gitignore"
 - ".travis.yml"
+- ".yardopts"
 - Gemfile
 - History.txt
 - LICENCE.md
@@ -100,6 +104,7 @@ files:
 - lib/distribution/chisquare/java.rb
 - lib/distribution/chisquare/ruby.rb
 - lib/distribution/chisquare/statistics2.rb
+- lib/distribution/distributable.rb
 - lib/distribution/exponential.rb
 - lib/distribution/exponential/gsl.rb
 - lib/distribution/exponential/ruby.rb
@@ -140,11 +145,15 @@ files:
 - lib/distribution/poisson/gsl.rb
 - lib/distribution/poisson/java.rb
 - lib/distribution/poisson/ruby.rb
+- lib/distribution/shorthand.rb
 - lib/distribution/t.rb
 - lib/distribution/t/gsl.rb
 - lib/distribution/t/java.rb
 - lib/distribution/t/ruby.rb
 - lib/distribution/t/statistics2.rb
+- lib/distribution/uniform.rb
+- lib/distribution/uniform/gsl.rb
+- lib/distribution/uniform/ruby.rb
 - lib/distribution/version.rb
 - lib/distribution/weibull.rb
 - lib/distribution/weibull/gsl.rb
@@ -167,12 +176,13 @@ files:
 - spec/spec.opts
 - spec/spec_helper.rb
 - spec/t_spec.rb
+- spec/uniform_spec.rb
 - spec/weibull_spec.rb
 - vendor/java/commons-math-2.2.jar
 homepage: https://github.com/sciruby/distribution
 licenses: []
 metadata: {}
-post_install_message: 
+post_install_message:
 rdoc_options: []
 require_paths:
 - lib
@@ -187,9 +197,8 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubyforge_project: 
-rubygems_version: 2.4.5
-signing_key: 
+rubygems_version: 3.1.2
+signing_key:
 specification_version: 4
 summary: Distribution
 test_files: []