three141 1.0.0

data/History.txt

@@ -0,0 +1,3 @@
+== 1.0.0 / 2007-06-24
+
+* First Version

data/Manifest.txt

@@ -0,0 +1,7 @@
+History.txt
+Manifest.txt
+README.txt
+Rakefile
+bin/three141
+lib/three141.rb
+test/test_three141.rb

data/README.txt

@@ -0,0 +1,50 @@
+three141
+    by Sandro Paganotti
+    www.314159265358979323846264338327950288419716939937510.net
+
+== DESCRIPTION:
+  
+The BBP algorithm is very powerful because it let you calculate the nth pi digit without the need to store the previous n-1th digits. In this website is avaiable ruby-bbp, a ruby porting of the bbp algorithm.
+
+== FEATURES/PROBLEMS:
+  
+It calculates the tenth pi digits starting at a given position
+
+== SYNOPSIS:
+
+a = Three141.new
+a.digits(2)
+#     => ["3", "F", "6", "A", "8", "8", "8", "5", "A", "3", "0", "8", "E", "0", "0", "0"]
+
+== REQUIREMENTS:
+
+Ruby 1.8.5 or above
+
+== INSTALL:
+
+sudo gem install three141
+
+== LICENSE:
+
+(The MIT License)
+
+Copyright (c) 2007 FIX
+
+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/Rakefile

@@ -0,0 +1,18 @@
+# -*- ruby -*-
+
+require 'rubygems'
+require 'hoe'
+require './lib/three141.rb'
+
+Hoe.new('three141', Three141::VERSION) do |p|
+  p.rubyforge_name = 'three141'
+  p.author = 'Sandro Paganotti'
+  p.email = 'info@314159265358979323846264338327950288419716939937510.net'
+  p.summary = 'A BBP Porting in Ruby'
+  p.remote_rdoc_dir = 'three141' # Release to root
+  p.description = p.paragraphs_of('README.txt', 2..5).join("\n\n")
+  p.url = p.paragraphs_of('README.txt', 0).first.split(/\n/)[1..-1]
+  p.changes = p.paragraphs_of('History.txt', 0..1).join("\n\n")
+end
+
+# vim: syntax=Ruby

data/lib/three141.rb

@@ -0,0 +1,123 @@
+class Three141
+  
+  VERSION = '1.0.0'
+  
+  def initialize
+    @tp = Array.new
+    @tp1 = 0
+  end
+  
+  
+  # It returns the first 10 pi digits (hex) starting from position = id
+  # 
+  # eg: Bbp.new.digits(2)
+  #     => ["3", "F", "6", "A", "8", "8", "8", "5", "A", "3", "0", "8", "E", "0", "0", "0"]
+  
+  def digits(id)
+    s1 = series(1, id).to_f;
+    s2 = series(4, id).to_f;
+    s3 = series(5, id).to_f;
+    s4 = series(6, id).to_f;
+    
+    pid = (4.0 * s1 - 2.0 * s2 - s3 - s4).to_f;
+    pid = pid - pid.to_i;
+    if (pid < 0.0) 
+      pid = pid + 1.0;
+    end
+    chx = ihex(pid, 16);
+    return chx[0..10]
+  end
+  
+  private
+  
+  def ihex(x,nhx)
+    hx = %w(0 1 2 3 4 5 6 7 8 9 A B C D E F)
+    y= x.abs.to_f
+    chx = Array.new
+    (0...nhx).to_a.each do |i|
+      y = 16.0 * (y - y.floor);
+      chx[i] = hx[y.to_i];
+    end
+    return chx
+  end 
+  
+   
+  def series(m,id)
+    s =0.0
+    ak=0.0  
+    p=0.0 
+    s=0.0 
+    t=0.0
+    eps = (10**(-17)).to_f
+    
+    (0...id).to_a.each do |k|
+      ak = 8 * k + m;
+      p = id - k;
+      t = expm(p, ak);
+      s = s + t / ak;
+      s = s - s.to_i;
+    end
+    
+    (id...(id +100)).to_a.each do |k|
+      ak = 8 * k + m;
+      t = (16.0 ** (id - k).to_f) / ak;
+      if (t < eps) 
+        break;
+      end
+      s = s + t;
+      s = s - s.to_i;
+    end
+    
+    return s
+  end 
+  
+  def expm(p,ak)
+    ntp = 25
+    exi = 0
+    p1 = 0.0 
+    pt = 0.0 
+    r = 0.0
+
+      if (@tp1 == 0) 
+        @tp1 = 1
+        @tp[0] = 1.0
+
+        (1...ntp).to_a.each do |i|
+          @tp[i] = 2.0 * @tp[i-1];
+        end
+        exi = ntp
+      end
+
+      if (ak == 1.0) 
+        return 0.0
+      end
+
+      exi = ntp
+      (1...ntp).to_a.each do |i|
+        exi = i
+        if (@tp[i] > p) 
+           break
+        end
+      end
+        
+      pt = @tp[exi-1].to_f;
+      p1 = p;
+      r = 1.0;
+
+      (1..exi).to_a.each do |j|
+        if (p1 >= pt)
+          r = 16.0 * r;
+          r = r - ((r / ak).to_i * ak).to_i
+          p1 = p1 - pt
+        end
+        pt = 0.5 * pt
+        if (pt >= 1.0)
+          r = r * r
+          r = r - ((r / ak).to_i * ak).to_i
+        end
+      end
+
+      return r;
+  end
+  
+end

metadata

@@ -0,0 +1,63 @@
+--- !ruby/object:Gem::Specification 
+rubygems_version: 0.9.2
+specification_version: 1
+name: three141
+version: !ruby/object:Gem::Version 
+  version: 1.0.0
+date: 2007-06-24 00:00:00 +02:00
+summary: A BBP Porting in Ruby
+require_paths: 
+- lib
+email: info@314159265358979323846264338327950288419716939937510.net
+homepage: "    by Sandro Paganotti"
+rubyforge_project: three141
+description: "== FEATURES/PROBLEMS:  It calculates the tenth pi digits starting at a given position  == SYNOPSIS:  a = Three141.new a.digits(2) #     => [\"3\", \"F\", \"6\", \"A\", \"8\", \"8\", \"8\", \"5\", \"A\", \"3\", \"0\", \"8\", \"E\", \"0\", \"0\", \"0\"]  == REQUIREMENTS:"
+autorequire: 
+default_executable: 
+bindir: bin
+has_rdoc: true
+required_ruby_version: !ruby/object:Gem::Version::Requirement 
+  requirements: 
+  - - ">"
+    - !ruby/object:Gem::Version 
+      version: 0.0.0
+  version: 
+platform: ruby
+signing_key: 
+cert_chain: 
+post_install_message: 
+authors: 
+- Sandro Paganotti
+files: 
+- History.txt
+- Manifest.txt
+- README.txt
+- Rakefile
+- bin/three141
+- lib/three141.rb
+- test/test_three141.rb
+test_files: 
+- test/test_three141.rb
+rdoc_options: 
+- --main
+- README.txt
+extra_rdoc_files: 
+- History.txt
+- Manifest.txt
+- README.txt
+executables: 
+- three141
+extensions: []
+
+requirements: []
+
+dependencies: 
+- !ruby/object:Gem::Dependency 
+  name: hoe
+  version_requirement: 
+  version_requirements: !ruby/object:Gem::Version::Requirement 
+    requirements: 
+    - - ">="
+      - !ruby/object:Gem::Version 
+        version: 1.2.1
+    version: