malge 0.0.3 → 0.0.4

data/CHANGES

@@ -1,6 +1,9 @@
 = vasputils changelog
 
-== Master (for 0.0.4)
+== Master (for 0.0.5)
+
+== Version 0.0.4
+* Add Malge::ErrorFittedFunction and inherited classes.
 
 == Version 0.0.3
 * Add Malge::LeastSquare::a_exp_bx

data/Gemfile

@@ -10,5 +10,5 @@ group :development do
   gem "bundler", "~> 1.2.2"
   gem "jeweler", "~> 1.8.3"
   gem "simplecov", ">= 0"
-  #gem "psych", ">= 0"
+  gem "psych", ">= 0"
 end

data/VERSION

@@ -1 +1 @@
-0.0.3
+0.0.4

data/lib/malge.rb

@@ -1,6 +1,11 @@
 module Malge; end
 
-require "malge/simultaneousequations.rb"
+require "malge/errorfittedfunction.rb"
+require "malge/errorfittedfunction/aexpbx.rb"
+require "malge/errorfittedfunction/aexpbx32.rb"
+require "malge/errorfittedfunction/axinv.rb"
+require "malge/errorfittedfunction/axinv32.rb"
 require "malge/leastsquare.rb"
 require "malge/matrix.rb"
+require "malge/simultaneousequations.rb"
 require "malge/vector.rb"

data/lib/malge/errorfittedfunction.rb

@@ -0,0 +1,87 @@
+#! /usr/bin/env ruby
+#coding: utf-8
+
+#Class for function to which x[i] and y[i] fit.
+#This class include mathematical dealing for errors.
+#Not simple fitting;
+#this class assume true value with the expected most precise value and
+#use absolute value of difference between a data point and the assumed true value.
+class Malge::ErrorFittedFunction
+  attr_reader :coefficients
+
+  class TypeError < Exception; end
+  class SizeMismatchError < Exception; end
+  class NotImplementedError < Exception; end
+
+  #Generate error fitted function from x[i] and y[i].
+  #Errors are generated y[i] and y_true assumed by the expected most precise y[i].
+  def initialize(data_pairs)
+
+    raise TypeError unless data_pairs.class == Array
+    data_pairs.each do |pair|
+      raise TypeError unless pair.class == Array
+      raise SizeMismatchError unless pair.size == 2
+    end
+
+    @raw_pairs = data_pairs
+    #pp finest_y
+    @diff_abs_pairs = data_pairs.map { |pair| [pair[0], (pair[1] - finest_y).abs] }
+    @diff_abs_pairs.delete_at(-1)
+    fit
+  end
+
+  ##Return expected error 
+  #def expected_errors(x)
+  #  raise NotImplementedError, "Define #{__method__}() in the inherited class."
+  #end
+
+  #Return expected error at x, condition variable, on the fitted function.
+  #Note that this method does not return the error between actual and true value.
+  def expected_error(x)
+    raise NotImplementedError, "Define #{__method__}() in the inherited class."
+  end
+
+  #Return variance of distribution between each pair
+  #of expected error and y[i], actual data point.
+  #Note that this method does not return variance of difference
+  #between each pair of TRUE value and y[i].
+  #Ignore last value of x[i] and y[i].
+  def variance
+    sum = 0.0
+    #pp @diff_abs_pairs
+    @diff_abs_pairs.each do |pair|
+      #pp pair
+      sum += (pair[1] - expected_error(pair[0]) )**2
+    end
+    sum
+  end
+
+  #Return value of 'x', assumed condition value, correspond to extpected value of 'y'.
+  def x(y)
+    raise NotImplementedError, "Define #{__method__}() in the inherited class."
+  end
+
+  #Return the value of y[i] at which the most precise data is expected to be obtained.
+  def finest_y
+    raise NotImplementedError, "Define #{__method__}() in the inherited class."
+
+    #In the most case, it would be sufficient to select belows.
+    @raw_pairs.max_by { |pair| pair[0] }[1]
+    @raw_pairs.min_by { |pair| pair[0] }[1]
+  end
+
+  private
+
+  #Fit the data pairs to a certain function.
+  #Call from initalize().
+  #You need to design to set @coefficients.
+  def fit
+    #raise NotImplementedError, "Define fit() in the inherited class."
+    raise NotImplementedError, "Define #{__method__}() in the inherited class."
+  end
+
+  #Return the most precise y[i].
+  #This depends on each data and function.
+  #E.g., does a higher or lower value lead precise?
+end
+

data/lib/malge/errorfittedfunction/aexpbx.rb

@@ -0,0 +1,39 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+#Assumed y = a[0] * exp(a[1] *x)
+#
+#NOTE: @coefficients[0] might become negative value.
+# Need discussion for dealing?
+class Malge::ErrorFittedFunction::AExpBX < Malge::ErrorFittedFunction
+
+  def fit
+    inv_pairs =  @diff_abs_pairs.map {|pair|
+      x = pair[0]
+      y = Math::log(pair[1])
+      [x,y]
+    }
+    @coefficients = Malge::LeastSquare.least_square_1st_degree(inv_pairs)
+    @coefficients[0] = Math::exp @coefficients[0]
+  end
+
+  def expected_error(x)
+    @coefficients[0] * Math::exp( @coefficients[1] * x)
+  end
+
+  #y = a[0] * exp(a[1] *x)
+  #a[0] * exp(a[1] *x) = y
+  #exp(a[1] *x) = y/a[0]
+  #a[1] *x = log( y/a[0])
+  #x = log( y/a[0])/a[1]
+  def x(y)
+    return Math::log( y / @coefficients[0])/@coefficients[1]
+  end
+
+  def finest_y
+    @raw_pairs.max_by { |pair| pair[0] }[1]
+  end
+
+end
+
+

data/lib/malge/errorfittedfunction/aexpbx32.rb

@@ -0,0 +1,40 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+#Assumed y = a[0] * exp(a[1] *x)
+#
+#NOTE: @coefficients[0] might become negative value.
+# Need discussion for dealing?
+class Malge::ErrorFittedFunction::AExpBX32 < Malge::ErrorFittedFunction
+
+  def fit
+    inv_pairs =  @diff_abs_pairs.map {|pair|
+      x = pair[0] ** (3.0/2.0)
+      y = Math::log(pair[1])
+      [x,y]
+    }
+    @coefficients = Malge::LeastSquare.least_square_1st_degree(inv_pairs)
+    @coefficients[0] = Math::exp @coefficients[0]
+  end
+
+  def expected_error(x)
+    @coefficients[0] * Math::exp( @coefficients[1] * x **(3.0/2.0))
+  end
+
+  #y = a[0] * exp(a[1] *x^(3/2))
+  #a[0] * exp(a[1] *x^(3/2)) = y
+  #exp(a[1] *x^(3/2)) = y/a[0]
+  #a[1] *x^(3/2) = log( y/a[0])
+  #x^(3/2) = log( y/a[0])/a[1]
+  #x = (log( y/a[0])/a[1])^(2/3)
+  def x(y)
+    return (Math::log( y / @coefficients[0])/@coefficients[1]) **(2.0/3.0)
+  end
+
+  def finest_y
+    @raw_pairs.max_by { |pair| pair[0] }[1]
+  end
+
+end
+
+

data/lib/malge/errorfittedfunction/axinv.rb

@@ -0,0 +1,32 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+#Assumed y = a[0] + a[1]/x
+#
+#NOTE: @coefficients[0] might become negative value.
+# Need discussion for dealing?
+class Malge::ErrorFittedFunction::AXInv < Malge::ErrorFittedFunction
+
+  def fit
+    inv_pairs =  @diff_abs_pairs.map {|pair| [1.0/pair[0], pair[1]]}
+    @coefficients = Malge::LeastSquare.least_square_1st_degree(inv_pairs)
+  end
+
+  def expected_error(x)
+    #@coefficients[0] + @coefficients[1] * x
+    @coefficients[0] + @coefficients[1] * (1.0/x)
+  end
+
+  #        y = a[0] + a[1]/x
+  # y - a[0] = a[1]/x
+  #        x = a[1] / (y - a[0])
+  def x(y)
+    return @coefficients[1] / (y - @coefficients[0])
+  end
+
+  def finest_y
+    @raw_pairs.max_by { |pair| pair[0] }[1]
+  end
+
+end
+

data/lib/malge/errorfittedfunction/axinv32.rb

@@ -0,0 +1,32 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+#Assumed y = a[0] + a[1]/x^{3/2}
+#
+#NOTE: @coefficients[0] might become negative value.
+# Need discussion for dealing?
+class Malge::ErrorFittedFunction::AXInv32 < Malge::ErrorFittedFunction
+
+  def fit
+    inv32_pairs =
+      @diff_abs_pairs.map {|pair| [1.0/(pair[0] ** (3.0/2.0)), pair[1]]}
+    @coefficients = Malge::LeastSquare.least_square_1st_degree(inv32_pairs)
+  end
+
+  def expected_error(x)
+    @coefficients[0] + @coefficients[1] * (1.0/(x** (3.0/2.0)))
+  end
+
+  #        y = a[0] + a[1]/x
+  # y - a[0] = a[1]/x
+  #        x = a[1] / (y - a[0])
+  def x(y)
+    return (@coefficients[1] / (y - @coefficients[0])) ** (2.0/3.0)
+  end
+
+  def finest_y
+    @raw_pairs.max_by { |pair| pair[0] }[1]
+  end
+
+end
+

data/malge.gemspec

@@ -5,11 +5,11 @@
 
 Gem::Specification.new do |s|
   s.name = "malge"
-  s.version = "0.0.3"
+  s.version = "0.0.4"
 
   s.required_rubygems_version = Gem::Requirement.new(">= 0") if s.respond_to? :required_rubygems_version=
   s.authors = ["ippei94da"]
-  s.date = "2013-02-08"
+  s.date = "2013-02-14"
   s.description = "Mathematical library to deal with basic problems in algebla.\n  "
   s.email = "ippei94da@gmail.com"
   s.extra_rdoc_files = [
@@ -25,19 +25,29 @@ Gem::Specification.new do |s|
     "Rakefile",
     "VERSION",
     "lib/malge.rb",
+    "lib/malge/errorfittedfunction.rb",
+    "lib/malge/errorfittedfunction/aexpbx.rb",
+    "lib/malge/errorfittedfunction/aexpbx32.rb",
+    "lib/malge/errorfittedfunction/axinv.rb",
+    "lib/malge/errorfittedfunction/axinv32.rb",
     "lib/malge/leastsquare.rb",
     "lib/malge/matrix.rb",
     "lib/malge/simultaneousequations.rb",
     "lib/malge/vector.rb",
     "malge.gemspec",
     "test/helper.rb",
+    "test/test_errorfittedfunction.rb",
+    "test/test_errorfittedfunction_aexpbx.rb",
+    "test/test_errorfittedfunction_aexpbx32.rb",
+    "test/test_errorfittedfunction_axinv.rb",
+    "test/test_errorfittedfunction_axinv32.rb",
     "test/test_leastsquare.rb",
     "test/test_simultaneousequations.rb"
   ]
   s.homepage = "http://github.com/ippei94da/malge"
   s.licenses = ["MIT"]
   s.require_paths = ["lib"]
-  s.rubygems_version = "1.8.11"
+  s.rubygems_version = "1.8.23"
   s.summary = "Math ALGEbra library."
 
   if s.respond_to? :specification_version then
@@ -48,17 +58,20 @@ Gem::Specification.new do |s|
       s.add_development_dependency(%q<bundler>, ["~> 1.2.2"])
       s.add_development_dependency(%q<jeweler>, ["~> 1.8.3"])
       s.add_development_dependency(%q<simplecov>, [">= 0"])
+      s.add_development_dependency(%q<psych>, [">= 0"])
     else
       s.add_dependency(%q<rdoc>, ["~> 3.12"])
       s.add_dependency(%q<bundler>, ["~> 1.2.2"])
       s.add_dependency(%q<jeweler>, ["~> 1.8.3"])
       s.add_dependency(%q<simplecov>, [">= 0"])
+      s.add_dependency(%q<psych>, [">= 0"])
     end
   else
     s.add_dependency(%q<rdoc>, ["~> 3.12"])
     s.add_dependency(%q<bundler>, ["~> 1.2.2"])
     s.add_dependency(%q<jeweler>, ["~> 1.8.3"])
     s.add_dependency(%q<simplecov>, [">= 0"])
+    s.add_dependency(%q<psych>, [">= 0"])
   end
 end
 

data/test/test_errorfittedfunction.rb

@@ -0,0 +1,87 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+require "test/unit"
+require "malge.rb"
+require "malge/errorfittedfunction.rb"
+
+class Malge::ErrorFittedFunction::Dummy00 < Malge::ErrorFittedFunction
+  def fit; end
+end
+
+class Malge::ErrorFittedFunction::Dummy01 < Malge::ErrorFittedFunction
+  def expected_error(x)
+    @coefficients[0] + @coefficients[1] * x
+  end
+
+  #Assume: y = 1.0 + 2.0 * x
+  #x = 1/2( y - 1.0)
+  def x(y)
+    return 1.0/2.0 * (y - 1.0)
+  end
+
+  #Assume: y = 1.0 + 2.0 * x
+  def fit
+    @coefficients = [1.0, 2.0]
+  end
+
+  def finest_y
+    #max
+    finest_pair = @raw_pairs.max_by { |pair| pair[0] }
+    finest_pair[1]
+
+    #@raw_pairs[x][1]
+  end
+
+end
+
+class TC_ErrorFittedFunction < Test::Unit::TestCase
+  def setup
+    #@eff00 = Malge::ErrorFittedFunction::Dummy00.new([], [])
+    @eff01 = Malge::ErrorFittedFunction::Dummy01.new(
+      #2x+1 = 1, -1 
+      #2x+1 = 3, +2 
+      #2x+1 = 5, -1 
+      #[0.0, 1.0, 2.0],
+      #[0.0, 5.0, 4.0]
+      [
+        [0.0, 0.0],
+        [1.0, 5.0],
+        [2.0, 4.0],
+      ]
+    )
+  end
+
+  def test_variance
+    assert_equal(13.0, @eff01.variance)
+    #diff_abs = [4,1]
+    #expected = [1,3]
+  end
+
+  def test_expected_error
+    assert_equal(1.0, @eff01.expected_error(0.0))
+    assert_equal(3.0, @eff01.expected_error(1.0))
+    assert_equal(5.0, @eff01.expected_error(2.0))
+  end
+
+  def test_initialize
+    assert_raise(Malge::ErrorFittedFunction::NotImplementedError){
+      Malge::ErrorFittedFunction.new([])
+    }
+    assert_raise(Malge::ErrorFittedFunction::SizeMismatchError){
+      Malge::ErrorFittedFunction.new([[1,2], [1,3, 3]])
+    }
+    assert_raise(Malge::ErrorFittedFunction::TypeError){
+      Malge::ErrorFittedFunction.new([2, [1,2,3]])
+    }
+    assert_raise(Malge::ErrorFittedFunction::TypeError){
+      Malge::ErrorFittedFunction.new([[1,2], 1])
+    }
+  end
+
+  def test_finest_y
+    assert_equal(4.0, @eff01.finest_y)
+  end
+
+end
+

data/test/test_errorfittedfunction_aexpbx.rb

@@ -0,0 +1,73 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+require "test/unit"
+require "malge.rb"
+require "malge/errorfittedfunction.rb"
+
+class Malge::ErrorFittedFunction::AExpBX
+  public :fit
+end
+
+class TC_ErrorFittedFunction_AExpBX < Test::Unit::TestCase
+  TOLERANCE = 1.0E-10
+  def setup
+    @aebx00 = Malge::ErrorFittedFunction::AExpBX.new(
+      [
+        [ 0.0, 104.0],
+        [ 1.0, 102.0],
+        [ 2.0, 101.0],
+        [ 3.0, 100.0],
+      ]
+    )
+
+    @aebx01 = Malge::ErrorFittedFunction::AExpBX.new(
+      [
+        [ 0.0, 104.0],
+        [ 1.0,  98.0],
+        [ 2.0, 101.0],
+        [ 3.0, 100.0],
+      ]
+    )
+  end
+
+  #def test_initialize
+  #end
+
+  def test_fit
+    corrects = [4.0, - Math::log(2.0)]
+    results = @aebx00.coefficients
+    assert_equal(2, results.size)
+    assert_in_delta(corrects[0] , results[0], TOLERANCE)
+    assert_in_delta(corrects[1] , results[1], TOLERANCE)
+
+    results = @aebx01.coefficients
+    corrects = [4.0, - Math::log(2.0)]
+    assert_equal(2, results.size)
+    assert_in_delta(corrects[0] , results[0], TOLERANCE)
+    assert_in_delta(corrects[1] , results[1], TOLERANCE)
+  end
+
+  def test_expected_error
+    assert_in_delta(4.0, @aebx00.expected_error(0.0), TOLERANCE)
+    assert_in_delta(2.0, @aebx00.expected_error(1.0), TOLERANCE)
+    assert_in_delta(1.0, @aebx00.expected_error(2.0), TOLERANCE)
+    assert_in_delta(0.5, @aebx00.expected_error(3.0), TOLERANCE)
+  end
+
+  def test_finest_y
+    assert_equal(100.0, @aebx00.finest_y)
+  end
+
+  def test_variance
+    assert_in_delta( 0.0, @aebx00.variance, TOLERANCE)
+    #diff_abs = [4,1]
+    #expected = [1,3]
+  end
+
+  def test_x
+    assert_in_delta(1.0, @aebx00.x(2.0), TOLERANCE)
+  end
+
+end
+

data/test/test_errorfittedfunction_aexpbx32.rb

@@ -0,0 +1,92 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+require "test/unit"
+require "malge.rb"
+#require "malge/errorfittedfunction.rb"
+
+class Malge::ErrorFittedFunction::AExpBX32
+  public :fit
+end
+
+class TC_ErrorFittedFunction_AExpBX32 < Test::Unit::TestCase
+  TOLERANCE = 1.0E-09 #very large error in calculation at x()
+  def setup
+    @aebx00 = Malge::ErrorFittedFunction::AExpBX32.new(
+      [
+        [ 0.0, 100.0 + 3.0* 2.0**(-2.0 *  0    )],
+        [ 1.0, 100.0 + 3.0* 2.0**(-2.0 * 1.0  ) ],
+        [ 2.0, 100.0 + 3.0* 2.0**(-2.0 * (2.0**(3.0/2.0)))],
+        [ 3.0, 100.0           ],
+      ]
+    )
+
+    @aebx01 = Malge::ErrorFittedFunction::AExpBX32.new(
+      [
+        [ 0.0, 100.0 + 3.0* 2.0**(-2.0 *  0    )],
+        [ 1.0, 100.0 - 3.0* 2.0**(-2.0 * 1.0  ) ],
+        [ 2.0, 100.0 + 3.0* 2.0**(-2.0 * (2.0**(3.0/2.0)))],
+        [ 3.0, 100.0           ],
+      ]
+    )
+  end
+
+  #def test_initialize
+  #end
+
+  def test_fit
+    corrects = [3.0, - 2.0*Math::log(2.0)]
+    results = @aebx00.coefficients
+    assert_equal(2, results.size)
+    assert_in_delta(corrects[0] , results[0], TOLERANCE)
+    assert_in_delta(corrects[1] , results[1], TOLERANCE)
+
+    corrects = [3.0, - 2.0*Math::log(2.0)]
+    results = @aebx01.coefficients
+    assert_equal(2, results.size)
+    assert_in_delta(corrects[0] , results[0], TOLERANCE)
+    assert_in_delta(corrects[1] , results[1], TOLERANCE)
+  end
+
+  def test_expected_error
+    corrects =
+      [
+        3.0* 2.0**(-2.0 * 0.0)    ,
+        3.0* 2.0**(-2.0 * 1.0),
+        3.0* 2.0**(-2.0 * (2.0**(3.0/2.0)))
+      ]
+    assert_in_delta(corrects[0], @aebx00.expected_error(0.0), TOLERANCE)
+    assert_in_delta(corrects[1], @aebx00.expected_error(1.0), TOLERANCE)
+    assert_in_delta(corrects[2], @aebx00.expected_error(2.0), TOLERANCE)
+
+    corrects =
+      [
+        3.0* 2.0**(-2.0 * 0.0)    ,
+        3.0* 2.0**(-2.0 * 1.0),
+        3.0* 2.0**(-2.0 * (2.0**(3.0/2.0)))
+      ]
+    assert_in_delta(corrects[0], @aebx01.expected_error(0.0), TOLERANCE)
+    assert_in_delta(corrects[1], @aebx01.expected_error(1.0), TOLERANCE)
+    assert_in_delta(corrects[2], @aebx01.expected_error(2.0), TOLERANCE)
+
+  end
+
+  def test_finest_y
+    assert_equal(100.0, @aebx00.finest_y)
+    assert_equal(100.0, @aebx01.finest_y)
+  end
+
+  def test_variance
+    assert_in_delta( 0.0, @aebx00.variance, TOLERANCE)
+    assert_in_delta( 0.0, @aebx01.variance, TOLERANCE)
+  end
+
+  def test_x
+    assert_in_delta(0.0, @aebx00.x(3.00), TOLERANCE)
+    assert_in_delta(1.0, @aebx00.x(0.75), TOLERANCE)
+    assert_in_delta(0.0, @aebx01.x(3.00), TOLERANCE)
+    assert_in_delta(1.0, @aebx01.x(0.75), TOLERANCE)
+  end
+
+end
+

data/test/test_errorfittedfunction_axinv.rb

@@ -0,0 +1,57 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+require "test/unit"
+require "malge.rb"
+require "malge/errorfittedfunction.rb"
+
+#class Malge::ErrorFittedFunction::AXInv
+#  public :fit
+#end
+
+class TC_ErrorFittedFunction_AXInv < Test::Unit::TestCase
+  def setup
+    @axi01 = Malge::ErrorFittedFunction::AXInv.new(
+      #2x+1 = 1, -1 
+      #2x+1 = 3, +2 
+      #2x+1 = 5, -1 
+      #[0.0, 1.0, 2.0],
+      #[0.0, 5.0, 4.0]
+      [
+        [1.0, 8.0],
+        [2.0, 6.0],
+        [4.0, 5.0],
+        [8.0, 4.0],
+      ]
+    )
+  end
+
+  #def test_initialize
+  #end
+
+  #def test_fit
+  #end
+
+  def test_expected_error
+    assert_equal(4.0, @axi01.expected_error(1.0))
+    assert_equal(2.0, @axi01.expected_error(2.0))
+    assert_equal(1.0, @axi01.expected_error(4.0))
+    assert_equal(0.5, @axi01.expected_error(8.0))
+  end
+
+  def test_finest_y
+    assert_equal(4.0, @axi01.finest_y)
+  end
+
+  def test_variance
+    assert_equal( 0.0, @axi01.variance)
+    #diff_abs = [4,1]
+    #expected = [1,3]
+  end
+
+  def test_x
+    assert_equal(2.0, @axi01.x(2.0))
+  end
+
+end
+

data/test/test_errorfittedfunction_axinv32.rb

@@ -0,0 +1,71 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+require "test/unit"
+require "malge.rb"
+#require "malge/errorfittedfunction.rb"
+
+#class Malge::ErrorFittedFunction::AXInv32
+#  public :fit
+#end
+
+class TC_ErrorFittedFunction_AXInv32 < Test::Unit::TestCase
+  TOLERANCE = 1.0E-10
+
+  def setup
+    @axi00 = Malge::ErrorFittedFunction::AXInv32.new(
+      [
+        [ 1, 164.0],
+        [ 4, 108.0],
+        [16, 101.0],
+        [64, 100.0],
+      ]
+    )
+
+    @axi01 = Malge::ErrorFittedFunction::AXInv32.new(
+      [
+        [ 1, 164.0],
+        [ 4,  92.0],
+        [16, 101.0],
+        [64, 100.0],
+      ]
+    )
+
+  end
+
+  #def test_initialize
+  #end
+
+  def test_fit
+    corrects = [0.0, 64.0]
+    results = @axi00.coefficients
+    assert_in_delta(corrects[0], results[0], TOLERANCE)
+    assert_in_delta(corrects[1], results[1], TOLERANCE)
+  end
+
+  def test_expected_error
+    assert_equal(64.0, @axi00.expected_error(1.0))
+    assert_equal( 8.0, @axi00.expected_error(4.0))
+    assert_equal( 1.0, @axi00.expected_error(16.0))
+
+    assert_equal(64.0, @axi01.expected_error(1.0))
+    assert_equal( 8.0, @axi01.expected_error(4.0))
+    assert_equal( 1.0, @axi01.expected_error(16.0))
+  end
+
+  def test_finest_y
+    assert_equal(100.0, @axi00.finest_y)
+    assert_equal(100.0, @axi01.finest_y)
+  end
+
+  def test_variance
+    assert_equal( 0.0, @axi00.variance)
+    assert_equal( 0.0, @axi01.variance)
+  end
+
+  def test_x
+    assert_in_delta(16.0, @axi00.x(1.0), TOLERANCE)
+  end
+
+end
+

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: malge
 version: !ruby/object:Gem::Version
-  version: 0.0.3
+  version: 0.0.4
   prerelease: 
 platform: ruby
 authors:
@@ -9,11 +9,11 @@ authors:
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2013-02-08 00:00:00.000000000 Z
+date: 2013-02-14 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rdoc
-  requirement: &86230380 !ruby/object:Gem::Requirement
+  requirement: !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - ~>
@@ -21,10 +21,15 @@ dependencies:
         version: '3.12'
   type: :development
   prerelease: false
-  version_requirements: *86230380
+  version_requirements: !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ~>
+      - !ruby/object:Gem::Version
+        version: '3.12'
 - !ruby/object:Gem::Dependency
   name: bundler
-  requirement: &86229650 !ruby/object:Gem::Requirement
+  requirement: !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - ~>
@@ -32,10 +37,15 @@ dependencies:
         version: 1.2.2
   type: :development
   prerelease: false
-  version_requirements: *86229650
+  version_requirements: !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ~>
+      - !ruby/object:Gem::Version
+        version: 1.2.2
 - !ruby/object:Gem::Dependency
   name: jeweler
-  requirement: &86229110 !ruby/object:Gem::Requirement
+  requirement: !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - ~>
@@ -43,10 +53,15 @@ dependencies:
         version: 1.8.3
   type: :development
   prerelease: false
-  version_requirements: *86229110
+  version_requirements: !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ~>
+      - !ruby/object:Gem::Version
+        version: 1.8.3
 - !ruby/object:Gem::Dependency
   name: simplecov
-  requirement: &86228170 !ruby/object:Gem::Requirement
+  requirement: !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - ! '>='
@@ -54,7 +69,28 @@ dependencies:
         version: '0'
   type: :development
   prerelease: false
-  version_requirements: *86228170
+  version_requirements: !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: psych
+  requirement: !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '0'
 description: ! "Mathematical library to deal with basic problems in algebla.\n  "
 email: ippei94da@gmail.com
 executables: []
@@ -71,12 +107,22 @@ files:
 - Rakefile
 - VERSION
 - lib/malge.rb
+- lib/malge/errorfittedfunction.rb
+- lib/malge/errorfittedfunction/aexpbx.rb
+- lib/malge/errorfittedfunction/aexpbx32.rb
+- lib/malge/errorfittedfunction/axinv.rb
+- lib/malge/errorfittedfunction/axinv32.rb
 - lib/malge/leastsquare.rb
 - lib/malge/matrix.rb
 - lib/malge/simultaneousequations.rb
 - lib/malge/vector.rb
 - malge.gemspec
 - test/helper.rb
+- test/test_errorfittedfunction.rb
+- test/test_errorfittedfunction_aexpbx.rb
+- test/test_errorfittedfunction_aexpbx32.rb
+- test/test_errorfittedfunction_axinv.rb
+- test/test_errorfittedfunction_axinv32.rb
 - test/test_leastsquare.rb
 - test/test_simultaneousequations.rb
 homepage: http://github.com/ippei94da/malge
@@ -94,7 +140,7 @@ required_ruby_version: !ruby/object:Gem::Requirement
       version: '0'
       segments:
       - 0
-      hash: -142397791
+      hash: 215000701
 required_rubygems_version: !ruby/object:Gem::Requirement
   none: false
   requirements:
@@ -103,7 +149,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
       version: '0'
 requirements: []
 rubyforge_project: 
-rubygems_version: 1.8.11
+rubygems_version: 1.8.23
 signing_key: 
 specification_version: 3
 summary: Math ALGEbra library.