random_bell 0.1.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: a642dbc691a26cf7c63cbeec3df54322c221e2e5
+  data.tar.gz: 6afe3bc00805d3dca7bf613e44feac36123a851f
+SHA512:
+  metadata.gz: 64283a0a0fd759d0ceb355a6bc4f1ddba4307b872a543b18447f46de7037df12205663be849e2e9d466e4e15e56be192728c4b2c568331c95f365d973034125c
+  data.tar.gz: 331fc647a8a1e5f71fcb78914e611517338abe7402b9402e03ee02b6a5b19edcf438a5ca2570ade3b4df370a42c6f90eceac4851973842f0e6c07d94d9bca6a6

data/.gitignore

@@ -0,0 +1,22 @@
+*.gem
+*.rbc
+.bundle
+.config
+.yardoc
+Gemfile.lock
+InstalledFiles
+_yardoc
+coverage
+doc/
+lib/bundler/man
+pkg
+rdoc
+spec/reports
+test/tmp
+test/version_tmp
+tmp
+*.bundle
+*.so
+*.o
+*.a
+mkmf.log

data/.rspec

@@ -0,0 +1 @@
+--color

data/.travis.yml

@@ -0,0 +1,4 @@
+language: ruby
+rvm:
+- 2.0.0
+- 2.1.0

data/Gemfile

@@ -0,0 +1,4 @@
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in random_bell.gemspec
+gemspec

data/LICENSE

@@ -0,0 +1,21 @@
+The MIT License (MIT)
+
+Copyright (c) 2014 OSA Shunsuke
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

data/LICENSE.txt

@@ -0,0 +1,22 @@
+Copyright (c) 2014 OSA Shunsuke
+
+MIT License
+
+Permission is hereby granted, free of charge, to any person obtaining
+a copy of this software and associated documentation files (the
+"Software"), to deal in the Software without restriction, including
+without limitation the rights to use, copy, modify, merge, publish,
+distribute, sublicense, and/or sell copies of the Software, and to
+permit persons to whom the Software is furnished to do so, subject to
+the following conditions:
+
+The above copyright notice and this permission notice shall be
+included in all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
+EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
+MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
+NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
+LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
+WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

data/README.md

@@ -0,0 +1,35 @@
+# RandomBell
+
+RandomBell is nomal(Gaussian) random number generator.
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'random_bell'
+```
+
+And then execute:
+
+```
+$ bundle
+```
+
+Or install it yourself as:
+
+```
+$ gem install random_bell
+```
+
+## Usage
+
+TODO: Write usage instructions here
+
+## Contributing
+
+1. Fork it ( https://github.com/s-osa/random_bell/fork )
+2. Create your feature branch (`git checkout -b my-new-feature`)
+3. Commit your changes (`git commit -am 'Add some feature'`)
+4. Push to the branch (`git push origin my-new-feature`)
+5. Create a new Pull Request

data/Rakefile

@@ -0,0 +1,13 @@
+require "bundler/gem_tasks"
+require "rspec/core/rake_task"
+
+# Exclude histogram specs by default
+RSpec::Core::RakeTask.new("spec") do |task|
+  task.rspec_opts = "--tag ~histogram"
+end
+
+RSpec::Core::RakeTask.new("spec:histogram")
+
+task default: [:spec]
+
+

data/lib/random_bell/version.rb

@@ -0,0 +1,3 @@
+class RandomBell
+  VERSION = "0.1.0"
+end

data/lib/random_bell.rb

@@ -0,0 +1,59 @@
+require "random_bell/version"
+
+class RandomBell
+  def initialize(mu: 0.5, sigma: 0.2, range: 0.0..1.0, method: :box_muller)
+    @mu     = mu
+    @sigma  = sigma
+    @range  = range
+    @method = method
+  end
+
+  def standard
+    send(@method)
+  end
+
+  def normal
+    @mu + standard * @sigma
+  end
+
+  def rand
+    loop do
+      num = normal
+      return num if @range.include?(num)
+    end
+  end
+
+  def to_histogram(rows: 20, size: 10 ** 4)
+    histogram = ""
+
+    numbers = (1..size).to_a.map{ rand }
+
+    klasses, klass_width = [], (@range.max - @range.min).to_f / rows
+    (0..rows).each_cons(2) do |min, max|
+      klasses << Range.new(klass_width * min, klass_width * max)
+    end
+
+    divider = (klasses.map{|klass| numbers.select{|num| klass.include?(num) }.size }.max / 50.0).ceil
+    klasses.each do |klass|
+      n = numbers.select{|num| klass.include?(num) }.size / divider
+      histogram << sprintf("%+3.3f", klass.max) << ": " << "*" * n << "\n"
+    end
+
+    histogram
+  end
+
+  private
+
+  def box_muller
+    theta  = 2 * Math::PI * Random.rand
+    r      = Math.sqrt( -2 * Math.log(Random.rand) )
+    z1, z2 = r * Math.cos(theta), r * Math.sin(theta)
+    [z1, z2].sample
+  end
+
+  def central_limit
+    arr = []
+    12.times{ arr << Random.rand }
+    arr.reduce(0.0, &:+) - 6
+  end
+end

data/random_bell.gemspec

@@ -0,0 +1,24 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'random_bell/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = "random_bell"
+  spec.version       = RandomBell::VERSION
+  spec.authors       = ["OSA Shunsuke"]
+  spec.email         = ["hhelibebcnofnenamg@gmail.com"]
+  spec.summary       = %q{Normal(Gaussian) random number generator.}
+  spec.description   = %q{Normal(Gaussian) random number generator.}
+  spec.homepage      = ""
+  spec.license       = "MIT"
+
+  spec.files         = `git ls-files -z`.split("\x0")
+  spec.executables   = spec.files.grep(%r{^bin/}) { |f| File.basename(f) }
+  spec.test_files    = spec.files.grep(%r{^(test|spec|features)/})
+  spec.require_paths = ["lib"]
+
+  spec.add_development_dependency "bundler"
+  spec.add_development_dependency "rake"
+  spec.add_development_dependency "rspec"
+end

data/spec/random_bell_spec.rb

@@ -0,0 +1,129 @@
+require 'spec_helper'
+
+describe RandomBell do
+  REPEAT = 10 ** 5
+  WITHIN = 0.1
+
+  describe "#standard" do
+    before :all do
+      @bell = RandomBell.new
+    end
+
+    describe "results" do
+      before :all do
+        @results = []
+        REPEAT.times{ @results << @bell.standard }
+      end
+
+      it "should be 0 as average" do
+        expect(@results.average).to be_within(WITHIN).of(0)
+      end
+
+      it "should be 1 as variance" do
+        expect(@results.variance).to be_within(WITHIN).of(1)
+      end
+
+      it "should be 1 as standard diviation" do
+        expect(@results.standard_diviation).to be_within(WITHIN).of(1)
+      end
+
+      it "drow histgram", histogram: true do
+        puts @results.to_histogram
+      end
+    end
+  end
+
+  describe "#normal" do
+    before :all do
+      @mu, @sigma = 1, 0.5
+      @bell = RandomBell.new(mu: @mu, sigma: @sigma)
+    end
+
+    describe "results" do
+      before :all do
+        @results = []
+        REPEAT.times{ @results << @bell.normal }
+      end
+
+      it "should be mu as average" do
+        expect(@results.average).to be_within(WITHIN).of(@mu)
+      end
+
+      it "should be square sigma as variance" do
+        expect(@results.variance).to be_within(WITHIN).of(@sigma ** 2)
+      end
+
+      it "should be sigma as normal diviation" do
+        expect(@results.standard_diviation).to be_within(WITHIN).of(@sigma)
+      end
+
+      it "drow histgram", histogram: true do
+        puts @results.to_histogram
+      end
+    end
+  end
+
+  describe "#rand" do
+    context "range not given" do
+      before :all do
+        @bell = RandomBell.new(mu: 1, sigma: 0.5)
+      end
+
+      it "should return float" do
+        expect(@bell.rand).to be_an_instance_of Float
+      end
+
+      describe "results" do
+        before :all do
+          @results = []
+          REPEAT.times{ @results << @bell.rand }
+        end
+
+        it "should be 1 as max" do
+          expect(@results.max).to be <= 1
+        end
+
+        it "should be 1 as min" do
+          expect(@results.min).to be >= 0
+        end
+
+        it "drow histgram", histogram: true do
+          puts @results.to_histogram
+        end
+      end
+    end
+
+    context "range given" do
+      before :all do
+        @mu, @sigma, @range = 50, 20, 40..100
+        @bell = RandomBell.new(mu: @mu, sigma: @sigma, range: @range)
+      end
+
+      describe "results" do
+        before :all do
+          @results = []
+          REPEAT.times{ @results << @bell.rand }
+        end
+
+        it "should be 100 as max" do
+          expect(@results.max).to be <= 100
+        end
+
+        it "should be 40 as min" do
+          expect(@results.min).to be >= 40
+        end
+
+        it "drow histgram", histogram: true do
+          puts @results.to_histogram
+        end
+      end
+    end
+  end
+
+  describe "#to_histogram" do
+    it "should return histogram string" do
+      bell = RandomBell.new(mu: 0.7)
+      expect(bell.to_histogram).to match(/\*+/)
+    end
+  end
+end

data/spec/spec_helper.rb

@@ -0,0 +1,19 @@
+project_root = File.join(File.dirname(__FILE__), '..')
+$: << project_root
+
+require 'lib/random_bell'
+
+class Array
+  def average
+    reduce(0.0, &:+) / size
+  end
+
+  def variance
+    avg = average
+    reduce(0.0){|sum, f| sum += ( f - avg ) ** 2 } / size
+  end
+
+  def standard_diviation
+    Math.sqrt(variance)
+  end
+end

metadata

@@ -0,0 +1,101 @@
+--- !ruby/object:Gem::Specification
+name: random_bell
+version: !ruby/object:Gem::Version
+  version: 0.1.0
+platform: ruby
+authors:
+- OSA Shunsuke
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2014-06-05 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: rspec
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+description: Normal(Gaussian) random number generator.
+email:
+- hhelibebcnofnenamg@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- ".rspec"
+- ".travis.yml"
+- Gemfile
+- LICENSE
+- LICENSE.txt
+- README.md
+- Rakefile
+- lib/random_bell.rb
+- lib/random_bell/version.rb
+- random_bell.gemspec
+- spec/random_bell_spec.rb
+- spec/spec_helper.rb
+homepage: ''
+licenses:
+- MIT
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.2.2
+signing_key: 
+specification_version: 4
+summary: Normal(Gaussian) random number generator.
+test_files:
+- spec/random_bell_spec.rb
+- spec/spec_helper.rb