rsl 0.0.1

data/CHANGELOG

@@ -0,0 +1,6 @@
+CHANGELOG
+
+== 0.0.1
+
+* Initial release.
+* Simple least squares fit

data/README

@@ -0,0 +1,11 @@
+== Ruby Scientific Library
+ 
+Native Ruby scientific/numeric programs, available for use in any application.
+
+For more information see http://www.sbrew.com/rsl
+
+RSL is distributed under the Ruby License.
+
+== WARNING
+
+This is alpha-quality software.  It is still very young.

data/Rakefile

@@ -0,0 +1,198 @@
+require 'rubygems'
+require 'rake'
+require 'rake/testtask'
+require 'rake/rdoctask'
+require 'rake/packagetask'
+require 'rake/gempackagetask'
+require 'rake/contrib/rubyforgepublisher'
+
+$:.unshift(File.dirname(__FILE__) + "/lib")
+require 'rsl'
+
+PKG_NAME      = 'rsl'
+PKG_VERSION   = Rsl::VERSION
+PKG_FILE_NAME = "#{PKG_NAME}-#{PKG_VERSION}"
+
+RELEASE_NAME  = "REL #{PKG_VERSION}"
+
+RUBY_FORGE_PROJECT = "rsl"
+RUBY_FORGE_USER    = "timbreitkreutz"
+
+desc "Default Task"
+task :default => [ :clean, :test ]
+
+desc "Clean images generated by tests"
+task :clean do
+  rm FileList['test/output/*.png']
+  rm_rf 'pkg'
+  rm_rf 'doc'
+end
+
+# Run the unit tests
+Rake::TestTask.new { |t|
+  t.libs << "test"
+  t.pattern = 'test/*_test.rb'
+  t.verbose = true
+}
+
+
+# Genereate the RDoc documentation
+Rake::RDocTask.new { |rdoc|
+  rdoc.rdoc_dir = 'doc'
+  rdoc.title    = "RSL -- Ruby Scientific Library"
+#  rdoc.options << '--line-numbers --inline-source --main README --accessor adv_attr_accessor=M'
+  rdoc.template = "#{ENV['template']}.rb" if ENV['template']
+  rdoc.rdoc_files.include('README', 'CHANGELOG')
+  rdoc.rdoc_files.include('lib/rsl.rb')
+  rdoc.rdoc_files.include('lib/rsl/*.rb')
+}
+
+
+# Create compressed packages
+spec = Gem::Specification.new do |s|
+  s.platform = Gem::Platform::RUBY
+  s.name = PKG_NAME
+  s.summary = "A native Ruby Scientific/Numerical library."
+  s.description = %q{Find and contribute numeric/scientific programs in native Ruby.}
+  s.version = PKG_VERSION
+
+  s.author = "Tim Breitkreutz"
+  s.email = "tim@sbrew.com"
+  s.rubyforge_project = RUBY_FORGE_PROJECT
+  s.homepage = "http://www.sbrew.com/rsl"
+
+  s.has_rdoc = true
+  s.requirements << 'none'
+  s.require_path = 'lib'
+  s.autorequire = 'rsl'
+
+  s.files = [ "Rakefile", "README", "CHANGELOG" ]
+  s.files = s.files + Dir.glob( "lib/**/*" ).delete_if { |item| item.include?( "\.svn" ) }
+  s.files = s.files + Dir.glob( "test/**/*" ).delete_if { |item| item.include?( "\.svn" ) }
+end
+  
+Rake::GemPackageTask.new(spec) do |p|
+  p.gem_spec = spec
+  p.need_tar = true
+  p.need_zip = true
+end
+
+desc "Publish the API documentation"
+task :pgem => [:package] do 
+  Rake::SshFilePublisher.new("tim@sbrew.com", "public_html/gems/gems", "pkg", "#{PKG_FILE_NAME}.gem").upload
+end
+
+desc "Publish the release files to RubyForge."
+task :release => [:package] do
+  files = ["gem", "tgz", "zip"].map { |ext| "pkg/#{PKG_FILE_NAME}.#{ext}" }
+
+  if RUBY_FORGE_PROJECT then
+    require 'net/http'
+    require 'open-uri'
+
+    project_uri = "http://rubyforge.org/projects/#{RUBY_FORGE_PROJECT}/"
+    project_data = open(project_uri) { |data| data.read }
+    group_id = project_data[/[?&]group_id=(\d+)/, 1]
+    raise "Couldn't get group id" unless group_id
+
+    # This echos password to shell which is a bit sucky
+    if ENV["RUBY_FORGE_PASSWORD"]
+      password = ENV["RUBY_FORGE_PASSWORD"]
+    else
+      print "#{RUBY_FORGE_USER}@rubyforge.org's password: "
+      password = STDIN.gets.chomp
+    end
+
+    login_response = Net::HTTP.start("rubyforge.org", 80) do |http|
+      data = [
+        "login=1",
+        "form_loginname=#{RUBY_FORGE_USER}",
+        "form_pw=#{password}"
+      ].join("&")
+      http.post("/account/login.php", data)
+    end
+
+    cookie = login_response["set-cookie"]
+    raise "Login failed" unless cookie
+    headers = { "Cookie" => cookie }
+
+    release_uri = "http://rubyforge.org/frs/admin/?group_id=#{group_id}"
+    release_data = open(release_uri, headers) { |data| data.read }
+    package_id = release_data[/[?&]package_id=(\d+)/, 1]
+    raise "Couldn't get package id" unless package_id
+
+    first_file = true
+    release_id = ""
+
+    files.each do |filename|
+      basename  = File.basename(filename)
+      file_ext  = File.extname(filename)
+      file_data = File.open(filename, "rb") { |file| file.read }
+
+      puts "Releasing #{basename}..."
+
+      release_response = Net::HTTP.start("rubyforge.org", 80) do |http|
+        release_date = Time.now.strftime("%Y-%m-%d %H:%M")
+        type_map = {
+          ".zip"    => "3000",
+          ".tgz"    => "3110",
+          ".gz"     => "3110",
+          ".gem"    => "1400"
+        }; type_map.default = "9999"
+        type = type_map[file_ext]
+        boundary = "rubyqMY6QN9bp6e4kS21H4y0zxcvoor"
+
+        query_hash = if first_file then
+          {
+            "group_id" => group_id,
+            "package_id" => package_id,
+            "release_name" => RELEASE_NAME,
+            "release_date" => release_date,
+            "type_id" => type,
+            "processor_id" => "8000", # Any
+            "release_notes" => "",
+            "release_changes" => "",
+            "preformatted" => "1",
+            "submit" => "1"
+          }
+        else
+          {
+            "group_id" => group_id,
+            "release_id" => release_id,
+            "package_id" => package_id,
+            "step2" => "1",
+            "type_id" => type,
+            "processor_id" => "8000", # Any
+            "submit" => "Add This File"
+          }
+        end
+
+        query = "?" + query_hash.map do |(name, value)|
+          [name, URI.encode(value)].join("=")
+        end.join("&")
+
+        data = [
+          "--" + boundary,
+          "Content-Disposition: form-data; name=\"userfile\"; filename=\"#{basename}\"",
+          "Content-Type: application/octet-stream",
+          "Content-Transfer-Encoding: binary",
+          "", file_data, ""
+          ].join("\x0D\x0A")
+
+        release_headers = headers.merge(
+          "Content-Type" => "multipart/form-data; boundary=#{boundary}"
+        )
+
+        target = first_file ? "/frs/admin/qrs.php" : "/frs/admin/editrelease.php"
+        http.post(target + query, data, release_headers)
+      end
+
+      if first_file then
+        release_id = release_response.body[/release_id=(\d+)/, 1]
+        raise("Couldn't get release id") unless release_id
+      end
+
+      first_file = false
+    end
+  end
+end

data/lib/rsl/base.rb

@@ -0,0 +1,3 @@
+module Rsl
+  VERSION = '0.0.1'
+end

data/lib/rsl/least_squares_fit.rb

@@ -0,0 +1,43 @@
+
+
+module Rsl
+  # Linear fit--least squares regression
+  # Fits the data in xes and ys to Y = b + a X
+  # returns (a, b)
+  def Rsl.least_squares_fit( xes, ys )
+    raise "input vectors must not be nil" if xes == nil || ys == nil
+    raise "input vectors must be larger" if xes.size < 3 || ys.size < 3
+    raise "input vectors must be of equal size" if xes.size != ys.size
+
+    run_x = 0.0
+    run_y = 0.0
+
+    count = 1.0
+
+    xes.to_a.each { |x| 
+      raise "X vector contains a nil" if x == nil
+      run_x += ( x - run_x ) / count 
+      count += 1.0 
+    }
+    count = 1.0
+    ys.to_a.each { |y| 
+      raise "Y vector contains a nil" if y == nil
+      run_y += ( ( y ) - run_y ) / count
+      count += 1.0 
+    }
+
+    dx2, dxdy = 0.0, 0.0
+
+    (1..xes.size).each do |count|
+      dx = xes[count-1] - run_x
+      dy = ( ys[count-1] ) - run_y
+      dx2 += ( dx * dx - dx2 ) / count.to_f
+      dxdy += ( dx * dy - dxdy ) / count.to_f
+    end
+  
+    b = dxdy / dx2 if dx2 != 0.0
+    a = run_y - run_x * b
+   
+    return a, b
+  end
+end

data/lib/rsl.rb

@@ -0,0 +1,2 @@
+require File.dirname(__FILE__) + '/rsl/base'
+require File.dirname(__FILE__) + '/rsl/least_squares_fit'

data/test/least_squares_test.rb

@@ -0,0 +1,63 @@
+#!/usr/bin/env ruby
+
+require 'test/unit'
+require 'rsl'
+require 'matrix'
+require 'mathn'
+
+class LinearFitTest < Test::Unit::TestCase
+
+  def test_fit
+    a, b = Rsl.least_squares_fit([1, 2, 3], [10, 11, 12 ])
+    assert_equal(sprintf("%.2f", a), "9.00")
+    assert_equal(sprintf("%.2f", b), "1.00")
+
+    a, b = Rsl.least_squares_fit([1, 2, 3], [10, 10.5, 12 ])
+    assert_equal(sprintf("%.2f", a), "8.83")
+    assert_equal(sprintf("%.2f", b), "1.00")
+
+    a, b = Rsl.least_squares_fit([1, 2, 3, 8, 11, 22], [10, 10.5, 12, 12.55, 13.55, 15.01 ])
+    assert_equal(sprintf("%.2f", a), "10.53")
+    assert_equal(sprintf("%.2f", b), "0.22")
+
+    # Some simple boundaries, zeros and negative inputs 
+
+    a, b = Rsl.least_squares_fit([1, 0, 3, 8, 11, 22], [10, 10.5, 12, 12.55, 13.55, 15.01 ])
+    assert_equal(sprintf("%.2f", a), "10.65")
+    assert_equal(sprintf("%.2f", b), "0.22")
+
+    a, b = Rsl.least_squares_fit([1, 8, 3, 8, 11, 22], [10, 0, 12, 12.55, 13.55, 15.01 ])
+    assert_equal(sprintf("%.2f", a), "8.33")
+    assert_equal(sprintf("%.2f", b), "0.25")
+
+    a, b = Rsl.least_squares_fit([-100, 8, 3, 8, 11, 22], [10, 0, -12, 12.55, 13.55, 15.01 ])
+    assert_equal(sprintf("%.2f", a), "6.41")
+    assert_equal(sprintf("%.2f", b), "-0.01")
+  end
+
+  def test_nil_vector
+    assert_raise(RuntimeError) { a, b = Rsl.least_squares_fit(nil, nil) }
+  end
+
+  def test_nil_member
+    assert_raise(RuntimeError) { a, b = Rsl.least_squares_fit([nil, 1, 2, 3], [2, 3, 4, 5]) }
+  end
+
+  def test_wrong_length
+    assert_raise(RuntimeError) { a, b = Rsl.least_squares_fit([9, 1, 2, 3], [2, 3, 4, 5, 5]) }
+  end
+
+  def test_too_short
+    assert_raise(RuntimeError) { a, b = Rsl.least_squares_fit([9], [2]) }
+    assert_raise(RuntimeError) { a, b = Rsl.least_squares_fit([nil], [nil]) }
+    assert_raise(RuntimeError) { a, b = Rsl.least_squares_fit([2,2], [3,3]) }
+  end
+
+  def test_vectors
+    x = Vector[1, 2, 3]
+    y = Vector[10, 11, 12]
+    a, b = Rsl.least_squares_fit(x, y)
+    assert_equal(sprintf("%.2f", a), "9.00")
+    assert_equal(sprintf("%.2f", b), "1.00")
+  end
+end

metadata

@@ -0,0 +1,52 @@
+--- !ruby/object:Gem::Specification 
+rubygems_version: 0.8.11
+specification_version: 1
+name: rsl
+version: !ruby/object:Gem::Version 
+  version: 0.0.1
+date: 2006-03-30 00:00:00 -07:00
+summary: A native Ruby Scientific/Numerical library.
+require_paths: 
+- lib
+email: tim@sbrew.com
+homepage: http://www.sbrew.com/rsl
+rubyforge_project: rsl
+description: Find and contribute numeric/scientific programs in native Ruby.
+autorequire: rsl
+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: 
+authors: 
+- Tim Breitkreutz
+files: 
+- Rakefile
+- README
+- CHANGELOG
+- lib/rsl
+- lib/rsl.rb
+- lib/rsl/base.rb
+- lib/rsl/least_squares_fit.rb
+- test/least_squares_test.rb
+test_files: []
+
+rdoc_options: []
+
+extra_rdoc_files: []
+
+executables: []
+
+extensions: []
+
+requirements: 
+- none
+dependencies: []
+