mathpack 0.2.0 → 0.3.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 8c14dffa0f11862781796118a019139c018cb2eb
-  data.tar.gz: b2be5e7d8b8bfdac4d40003a72109b6b0b9ea366
+  metadata.gz: e0f0d47957524843b49cef7c52ac64230e6f7d68
+  data.tar.gz: cb662a79d6977fb3f2d891f2b2eb2881b2c00cf9
 SHA512:
-  metadata.gz: a4dcfdfcd380f2f3e23587210cfb594e533c0cc341b6a5862ec4e7891b8d72d72c8e6f09f7c7156f03bf74d87c204ee54f4baa16f76895ac55e3b62a5e964d01
-  data.tar.gz: edd6fce95f19953b01d6fc92a89b645fa52a855579c4136306e8b2e7ba33480a622d6b07abf3bb92f6772c5287f9e8d86a30ee6e0aebd8d1363f4fb4611ca7fe
+  metadata.gz: 1d134a1000e9bf0068a36923484f373cf546492e8a6795539ee8616f06d09ec6e31886e096929f6c463a23a8d48dbcbf973b035309d3ccac92dceec6973f88e0
+  data.tar.gz: 0a8fc14ec5126df7b9afe59ca2141b52539cea8a710875c6ca6bb10e5aae7fc4238cca77a0f6b359be6bc6d0b215f85d220439f0f5b99475592616d8ed264e30

data/lib/mathpack/approximation.rb

@@ -0,0 +1,73 @@
+require 'mathpack'
+
+module Mathpack
+  module Approximation
+    def self.basic_fi(x, power)
+      x**power
+    end
+
+    def self.approximate_by_polynom(params = {}, &function)
+      f = block_given? ? fill_f(params[:x], &function) : params[:f]
+      mnk_matrix, mnk_vector = count_mnk_values(params[:x], f, params[:polynom_power])
+      Mathpack::SLE.solve(matrix: mnk_matrix, f: mnk_vector).reverse
+    end
+
+    def self.scalar_function_composition(first_vect, second_vect, first_power, second_power)
+      result = 0.0
+      first_vect.each_index do |i|
+        result += basic_fi(first_vect[i], first_power) * basic_fi(second_vect[i], second_power)
+      end
+      result
+    end
+
+    def self.count_mnk_values(x, f, polynom_power)
+      mnk_matrix = Array.new(polynom_power + 1) { Array.new(polynom_power + 1) }
+      mnk_vector = Array.new(polynom_power + 1)
+      0.upto(polynom_power) do |i|
+        0.upto(polynom_power) do |j|
+          mnk_matrix[i][j] = scalar_function_composition(x, x, i, j)
+        end
+        mnk_vector[i] = scalar_function_composition(x, f, i, 1)
+      end
+      [mnk_matrix, mnk_vector]
+    end
+
+    def self.print_polynom(coefficients)
+      polynom = ''
+      coefficients.each_index do |i|
+        next if coefficients[i] == 0
+        if i == 0
+          polynom += '-' if coefficients[i] < 0.0
+        else
+          polynom += coefficients[i] < 0.0 ? ' - ' : ' + '
+        end
+        polynom += coefficients[i].abs.to_s if coefficients[i].abs != 1 || i == coefficients.length - 1
+        polynom += '*' if i != coefficients.length - 1 && coefficients[i].abs != 1
+        case i
+          when coefficients.length - 1
+
+          when coefficients.length - 2
+            polynom += 'x'
+          else
+            polynom += coefficients.length - i >= 11 ? "x^(#{coefficients.length - 1 - i})" : "x^#{coefficients.length - 1 - i}"
+        end
+      end
+      polynom
+    end
+
+    def self.generate_nodes(params = {})
+      nodes = []
+      x = params[:start]
+      loop do
+        nodes << x
+        x += params[:step]
+        break if x > params[:end]
+      end
+      nodes
+    end
+
+    def self.fill_f(x, &function)
+      x.map { |val| function.call(val) }
+    end
+  end
+end

data/lib/mathpack/version.rb

@@ -1,3 +1,3 @@
 module Mathpack
-  VERSION = '0.2.0'
+  VERSION = '0.3.0'
 end

data/lib/mathpack.rb

@@ -2,3 +2,4 @@ require 'mathpack/version'
 require 'mathpack/statistics'
 require 'mathpack/equation'
 require 'mathpack/sle'
+require 'mathpack/approximation'

data/mathpack.gemspec

@@ -18,6 +18,6 @@ Gem::Specification.new do |spec|
   spec.test_files    = spec.files.grep(%r{^(test|spec|features)/$})
   spec.require_paths = ["lib"]
 
-  spec.add_development_dependency "bundler", "~> 1.6"
-  spec.add_development_dependency "rake"
+  spec.add_development_dependency 'bundler', '~> 1.6'
+  spec.add_development_dependency 'rake'
 end

data/spec/approximation_spec.rb

@@ -0,0 +1,38 @@
+describe 'Approximation' do
+  require 'mathpack/approximation'
+
+  context '#mnk' do
+    let(:x_start) { 0.0 }
+    let(:x_end) { 10.0 }
+    let(:step) { 0.25 }
+    let(:nodes) { Mathpack::Approximation.generate_nodes(start: x_start, end: x_end, step: step) }
+    let(:values) { Mathpack::Approximation.fill_f(nodes){ |val| val**2 } }
+
+    it 'should generate nodes' do
+      expect(nodes.length).to be_eql(41)
+      expect(nodes.first).to be_eql(x_start)
+      expect(nodes.last).to be_eql(x_end)
+    end
+
+    it 'should fill function table' do
+      expect(values.length).to be_eql(nodes.length)
+      expect(values.last).to be_eql(nodes.last**2)
+    end
+
+    it 'should print polynoms' do
+      expect(Mathpack::Approximation.print_polynom([5])).to be_eql('5')
+      expect(Mathpack::Approximation.print_polynom([2, 3])).to be_eql('2*x + 3')
+      expect(Mathpack::Approximation.print_polynom([1, -1, 1])).to be_eql('x^2 - x + 1')
+      expect(Mathpack::Approximation.print_polynom([1, 3, 5])).to be_eql('x^2 + 3*x + 5')
+      expect(Mathpack::Approximation.print_polynom([1, -2, 3, -4])).to be_eql('x^3 - 2*x^2 + 3*x - 4')
+      expect(Mathpack::Approximation.print_polynom([1, 0, 1])).to be_eql('x^2 + 1')
+    end
+
+    it 'should approximate by polynom' do
+      coefficients = Mathpack::Approximation.approximate_by_polynom(x: nodes, f: values, polynom_power: 2)
+      expect(coefficients).to be_eql([1.0, 0.0, 0.0])
+      coefficients = Mathpack::Approximation.approximate_by_polynom(x: nodes, f: values, polynom_power: 0)
+      expect(coefficients).to be_eql([33.75])
+    end
+  end
+end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: mathpack
 version: !ruby/object:Gem::Version
-  version: 0.2.0
+  version: 0.3.0
 platform: ruby
 authors:
 - maxmilan
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2014-12-31 00:00:00.000000000 Z
+date: 2015-01-25 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
@@ -52,11 +52,13 @@ files:
 - README.md
 - Rakefile
 - lib/mathpack.rb
+- lib/mathpack/approximation.rb
 - lib/mathpack/equation.rb
 - lib/mathpack/sle.rb
 - lib/mathpack/statistics.rb
 - lib/mathpack/version.rb
 - mathpack.gemspec
+- spec/approximation_spec.rb
 - spec/equation_spec.rb
 - spec/sle_spec.rb
 - spec/statistics_spec.rb