rubylabs 0.5.4

data/LICENSE

@@ -0,0 +1,20 @@
+Copyright (c) 2009 John S. Conery
+
+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.rdoc

@@ -0,0 +1,9 @@
+= rubylabs
+
+The modules in this Gem have methods and data for lab projects in the textbook
+<em>The Science of Computing: A Problem Solving Approach</em>.  There is one module
+for each chapter, plus a set of common methods shared by more than one module.
+
+== Copyright
+
+Copyright (c) 2009 John S. Conery. See LICENSE for details.

data/Rakefile

@@ -0,0 +1,63 @@
+require 'rubygems'
+require 'rake'
+
+begin
+  require 'jeweler'
+  Jeweler::Tasks.new do |gem|
+    # gem is a Gem::Specification... see http://www.rubygems.org/read/chapter/20 for additional settings
+    gem.name = "rubylabs"
+    gem.summary = %Q{Software and data for lab projects in the Science of Computing text.}
+    gem.description = %Q{A set of modules for interactive experiments in an introductory computer science class.}
+    gem.email = "conery@cs.uoregon.edu"
+    gem.homepage = "http://github.com/conery/rubylabs"
+    gem.authors = ["conery"]
+    gem.files = FileList['[A-Z]*', 'lib/**/*.rb', 'data/*', 'bin/**/*.rb']
+    gem.test_files = FileList['test/**/*.rb']
+    gem.rubyforge_project = "rubylabs"
+  end
+
+  Jeweler::RubyforgeTasks.new  
+  Jeweler::GemcutterTasks.new
+
+rescue LoadError
+  puts "Jeweler (or a dependency) not available. Install it with: sudo gem install jeweler"
+end
+
+require 'rake/testtask'
+Rake::TestTask.new(:test) do |test|
+  test.libs << 'lib' << 'test'
+  test.pattern = 'test/**/*_test.rb'
+  test.verbose = true
+end
+
+begin
+  require 'rcov/rcovtask'
+  Rcov::RcovTask.new do |test|
+    test.libs << 'test'
+    test.pattern = 'test/**/*_test.rb'
+    # test.verbose = true
+  end
+rescue LoadError
+  task :rcov do
+    abort "RCov is not available. In order to run rcov, you must: sudo gem install spicycode-rcov"
+  end
+end
+
+
+task :default => :test
+
+require 'rake/rdoctask'
+Rake::RDocTask.new do |rdoc|
+  if File.exist?('VERSION.yml')
+    config = YAML.load(File.read('VERSION.yml'))
+    version = "#{config[:major]}.#{config[:minor]}.#{config[:patch]}"
+  else
+    version = ""
+  end
+
+  rdoc.rdoc_dir = 'rdoc'
+  rdoc.title = "rubylabs #{version}"
+  rdoc.rdoc_files.include('README*')
+  rdoc.rdoc_files.include('lib/**/*.rb')
+end
+

data/VERSION

@@ -0,0 +1 @@
+0.5.4

data/bin/bb.rb

@@ -0,0 +1,12 @@
+#! /usr/bin/ruby -rubygems
+
+# From:  http://www.rubyagent.com/2007/03/starting-multi-agent-system-environment
+
+require 'drb/drb'
+require 'rinda/tuplespace'
+
+DRb.start_service("druby://localhost:53783", Rinda::TupleSpace.new)
+
+puts 'The blackboard is running...'
+
+DRb.thread.join

data/bin/statistics2-0.53/ext/extconf.rb

@@ -0,0 +1,11 @@
+require "mkmf"
+create_makefile("statistics2")
+open("Makefile", "a") do |f|
+  f << "test:\n\truby -I./lib:. t.rb\n"
+  f << "\n"
+  f << "uninstall:\n"
+  f << "\t@echo rm $(RUBYLIBDIR)/statistics2.rb\n"
+  f << "\t@$(RM) $(RUBYLIBDIR)/statistics2.rb\n"
+  f << "\t@echo rm $(RUBYARCHDIR)/statistics2.{o,so}\n"
+  f << "\t@$(RM) $(RUBYARCHDIR)/statistics2.{o,so}\n"
+end

data/bin/statistics2-0.53/ext/show.rb

@@ -0,0 +1,11 @@
+#!/usr/bin/env ruby
+require "statistics2.so"
+if $0 == __FILE__
+  if ARGV.empty?
+    puts "Example:"
+    puts "        #$0 normaldist 0.01"
+    puts "        #$0 pf_x 2 3 0.01"
+    exit
+  end
+  p Statistics2.send(ARGV[0], *ARGV[1..-1].map{|x| eval(x)})
+end

data/bin/statistics2-0.53/ext/t.rb

@@ -0,0 +1,46 @@
+require 'test/unit'
+$test = true
+
+$LOAD_PATH.unshift ".."
+require 'sample-tbl'
+
+eval(File.read('../statistics2.rb').sub(/\bStatistics2\b/, 'Statistics20'))
+$mod0 = Statistics20
+
+require 'statistics2.so'
+$mod = Statistics2
+
+class T_Statistics2 < Test::Unit::TestCase
+
+  def test_normal
+    norm_tbl(0.0, 3.1) do |x|
+      a, b = $mod.normal___x(x), $mod0.normal___x(x)
+      assert_in_delta a, b, 0.000001
+    end
+  end
+
+  def test_chi
+    chi2_tbl() do |n, x|
+      a, b = $mod.pchi2_x(n, x), $mod0.pchi2_x(n, x)
+      assert_in_delta a/b, 1.0, 0.01
+    end
+  end
+
+  def test_t
+    t_tbl() do |n, x|
+      a, b = $mod.ptx__x(n, x), $mod.ptx__x(n, x)
+      assert_in_delta a, b, 0.001
+    end
+  end
+
+  def test_f
+    [0.01, 0.025, 0.05].each do |opt|
+      f_tbl(opt) do |n1, n2, y|
+        a, b = $mod.pf_x(n1, n2, y), $mod0.pf_x(n1, n2, y)
+        assert_in_delta a/b, 1.0, 0.01
+      end
+    end
+  end
+end
+
+require "t-inv.rb"

data/bin/statistics2-0.53/mklist.rb

@@ -0,0 +1,26 @@
+#!/usr/bin/env ruby
+cmd = "ruby sample-tbl.rb"
+param_tbl = [["norm", "tbl-nomral.tbl"],
+           ["chi2", "tbl-chi2.tbl"],
+           ["t", "tbl-t.tbl"],
+           ["f 0.05", "tbl-F50.tbl"],
+           ["f 0.025", "tbl-F25.tbl"],
+           ["f 0.01", "tbl-F10.tbl"]]
+
+unless dir = ARGV.shift
+  puts "Usage: $0 DIR"
+  puts "-- This script creates some tables of distributions:"
+  param_tbl.each do |mod, tbl|
+    puts "    " + mod + " => " + tbl
+  end
+  exit
+end
+
+unless File.directory? dir
+  system("mkdir #{dir}")
+end
+
+param_tbl.each do |pa, tbl|
+  file = File.join(dir, tbl)
+  system("#{cmd} #{pa} > #{file}")
+end

data/bin/statistics2-0.53/sample-tbl.rb

@@ -0,0 +1,129 @@
+#!/usr/bin/env ruby
+
+if $test
+  def nop(*a); ""; end
+  alias _sprintf nop
+  alias _printf nop
+  alias _puts nop
+else
+  alias _sprintf sprintf
+  alias _printf printf
+  alias _puts puts
+end
+
+def norm_line(x, n)
+  s = sprintf("%1.1f|", x)
+  (0.00).step(0.09, 0.01) do |y|
+    s << _sprintf(" %1.#{n}f", yield(x + y))
+  end
+  s
+end
+
+def norm_tbl(s, e, ln = nil, tn = nil, &b)
+  n = 4
+  unless ln
+    _printf("    " + (" "*(n-1) + "%1.2f")*10 + "\n", * (0..9).map {|x| x*0.01})
+  end
+  i = 0
+  (s.to_f).step(e.to_f, 0.1) do |x|
+    next if ln && ln.to_i != i+1
+    _puts norm_line(x, n, &b)
+    i += 1
+  end
+end
+
+def chi2_tbl(ln = nil, tn = nil)
+  pers = [0.995, 0.99, 0.975, 0.95, 0.05, 0.025, 0.01, 0.005]
+  arbi = (1..30).to_a + [40, 60, 80, 100]
+  form = "  %7.5f"
+  unless ln
+    _printf("     "); pers.each do |a|; _printf(form, a); end;  _puts
+  end
+  arbi.each_with_index do |n, i|
+    next if ln && ln.to_i != i+1
+    _printf("%4d|", n) unless tn
+    pers.each_with_index do |a, j|
+      next if tn && tn.to_i != j+1
+      form = case n
+             when 1; a >= 0.95 ? " %.4e" : " %6.3f"
+             when 2..5; a >= 0.95 ? "  %7.5f" : "  %7.3f"
+             when 6..24; "  %7.3f"
+             else 26..100; "  %7.2f"
+             end
+      _printf(form, yield(n, a))
+    end
+    _puts
+  end
+end
+
+def t_tbl(ln = nil, tn = nil)
+  pers = [0.5, 0.4, 0.3, 0.2, 0.1, 0.05, 0.02, 0.01, 0.001]
+  arbi = (1..30).to_a + [40, 60, 120]#, 9999]
+  form = " %7.3f"
+  unless ln
+    _printf("     "); pers.each do |a|; _printf(form, a); end;  _puts
+  end
+  arbi.each_with_index do |n, i|
+    next if ln && ln.to_i != i+1
+    _printf("%4d|", n) unless tn
+    pers.each_with_index do |a, j|
+      next if tn && tn.to_i != j+1
+      _printf(form, yield(n, a))
+    end
+    _puts
+  end
+end
+
+def f_tbl(a, k1 = nil, k2 = nil)
+  arbi1 = (1..10).to_a + [12, 15, 20, 24, 30, 40, 60, 120]#, 9999]
+  arbi2 = (1..30).to_a + [40, 60, 120]#, 9999]
+  unless k1
+    _printf("     "); arbi1.each do |n1|; _printf(" %4d", n1); end;  _puts
+  end
+  arbi2.each do |n2|
+    form = n2 == 1 ? " %4d" : n2 == 2 ? " %4.1f" : " %4.2f"
+    next if k2 && k2.to_i != n2
+    _printf("%4d|", n2) unless k2
+    arbi1.each do |n1|
+      next if k1 && k1.to_i != n1
+      _printf(form, yield(n1, n2, a.to_f))
+    end
+    _puts
+  end
+end
+
+def show_tbl(mod, dist, *opts)
+  case dist
+  when nil, "norm"
+    norm_tbl(0.0, 3.1, *opts) do |x|
+      mod.normal___x(x)
+    end
+  when "chi", "chi2"
+    chi2_tbl(*opts) do |n, x|
+      mod.pchi2_x(n, x)
+    end
+  when "t"
+    t_tbl(*opts) do |n, x|
+      mod.ptx__x(n, x)
+    end
+  when "f"
+    f_tbl(*opts) do |n1, n2, a|
+      mod.pf_x(n1, n2, a)
+    end
+  end
+end
+
+
+if $0 == __FILE__
+  if ARGV.empty?
+    puts "-- This script makes tables of distributions"
+    puts "Example:"
+    puts "      #$0 norm"
+    puts "      #$0 chi2"
+    puts "      #$0 t"
+    puts "      #$0 f 0.01"
+    exit
+  end
+  require "statistics2"
+  show_tbl(Statistics2, *ARGV)
+end