mathpack 0.4.9 → 0.5.0

checksums.yaml

@@ -1,7 +1,15 @@
 ---
-SHA1:
-  metadata.gz: 3d83a43c2f840fff63926eec5d9548e0f5051fcd
-  data.tar.gz: 949637138199c144e04ce0ed085d8b249d40a4dc
+!binary "U0hBMQ==":
+  metadata.gz: !binary |-
+    ZDY4OTU2ZDNlZGQwMTAyM2U2M2I3NjA4NzU1MDQ5ZGM4OGY5OGRmOA==
+  data.tar.gz: !binary |-
+    ZjRiMjFkYTdmMTc5YjJlZTRhMTk5MzUyZWJkMWJkYmY4OWU0MmVjMg==
 SHA512:
-  metadata.gz: 257fe6c5dac7b567ead47d8546aebd2f945a6d32a4cd87f7fb4e14ad4420eadf026a2139dc8770bf8caee294cb946beb250bb80136ce8b9947954ce8f6f2c5d9
-  data.tar.gz: a18a58ed4da327ffdd05f92d88c5f11940346273a16cc2b3d5b33a53f5b9c973ff9d57eee1cb0fc466ce08f46aa610fc7186a901cca741d92476981201ca01a7
+  metadata.gz: !binary |-
+    OGJhMzY1N2ViM2ZmZWM3MWM3ZDU2YmVlMDllN2M5NGJiMGQ4YmM4NTdjOTUx
+    MDFjZjc0NWYyNDM4ZGU4ZjA4ZWQ5OTM1NGM1Zjk5ZGNiNzdiZWFjYjFiOGEx
+    OWQ2ZmMzZjNlMWE4ZmI3YzAzOWUzZmI5ZDBmODMxYTM4YjdkNjY=
+  data.tar.gz: !binary |-
+    MjczOGM0NTJjZTEzMmIzYTUyZjg5OWMxY2IzOTM3ODE4ZjQ2N2MyNTc1MWUx
+    NDQ1OGFmM2M1NTBhNGRkOWU5OTZlYTMzYzEwYjQ1YjVmZGRiNjJmMzg0ODgx
+    MzI0ODMxM2I2ZWZjZGE0ZjY3M2NmZGNmZmJiOWI4ZTM5MWQ3ZWU=

data/.travis.yml

@@ -1,10 +1,20 @@
-language: "ruby"
+language: ruby
 rvm:
-  - 2.0.0
-  - 2.1.2
-  - 2.1.3
-  - 2.2.0
-  - 2.3.0
-script: "bundle exec rspec"
+- 2.0.0
+- 2.1.2
+- 2.1.3
+- 2.2.0
+- 2.3.0
+script: bundle exec rspec
+before_install:
+- gem update bundler
+deploy:
+  provider: rubygems
+  api_key:
+    secure: O7YI3MNFgWH9xrjpTaYu6zvfhgpgQrMQullQebWR9B7sHXUw69XRW4mSy6PdsAHhFddymmfIWBW1rqS5HhiO7I50r/o/S90Oqy2vftQ3Fm8EXhKywc60r6o60EUOiqw/QObnqT/yspX0PpJ4uGlkDpc6XoDY/01QvhvoK13TnS4=
+  gem: mathpack
+  on:
+    branch: master
+    repo: maxmilan/mathpack
 notifications:
   email: false

data/README.md

@@ -20,16 +20,16 @@ Or install it yourself as:
 
 ## Information
 `Mathpack` includes following modules:
-- **SLE**. Solves system of linear equations
-- **Statistics**. Provides methods to analyze data samples
-- **Functions**. Collects mathematical functions
-- **Approximation**. Allows to approximate table and analytical functions by polynom
-- **NonlinearEquations**. Solves unlinear mathematical equations
-- **IntegralEquations**. Solves integral second order Fredholm and Volter equations
-- **DifferentialEquations**. Solves system of differential equations with left border initial conditions
-- **Integration**. Integrates functions
-- **IO**. Prints data
-- **Functional**. Lambda functions
+* [SLE](#sle). Solves system of linear equations
+* [Statistics](#statistics). Provides methods to analyze data samples
+* [Functions](#functions). Collects mathematical functions
+* [Approximation](#approximation). Allows to approximate table and analytical functions by polynom
+* [NonlinearEquations](#nonlinearequations). Solves nonlinear mathematical equations
+* [IntegralEquations](#integralequations). Solves integral second order Fredholm and Volter equations
+* [DifferentialEquations](#differentialequations). Solves system of differential equations with left border initial conditions
+* [Integration](#integration). Integrates functions
+* [IO](#io). Prints data
+* [Functional](#functional). Lambda functions
 
 ## Statistics
 `Statistics` class have following methods
@@ -98,14 +98,14 @@ stat.trend(polynom_power: 1) #=> 1.7999999999999996*x - 0.9999999999999987
 #### solve(params = {})
 returns solution of nonlinear equation. 
 ##### Parameters
-- ***start*** - point to start iteration process
-- ***eps*** - calculations accuraccy
+* ***start*** - point to start iteration process
+* ***eps*** - calculations accuraccy
 
 ### Usage
 Now you have no problems solving nonlinear equations. If you want, for example, to solve equation ![equation](http://latex.codecogs.com/gif.latex?x%5E%7B2%7D%20%3D%20%5Csin%28%7Bx+1%7D%29)
 
 You need to complete the following steps:
-- Equation should look like ![equation](http://latex.codecogs.com/gif.latex?%5Ctiny%20f%28x%29%20%3D%200)
+- Equation should look like ![equation](http://latex.codecogs.com/gif.latex?f%28x%29%3D0)
 - For our equation ![equation](http://latex.codecogs.com/gif.latex?%5Csmall%20f%28x%29%20%3D%20x%5E%7B2%7D%20-%20%5Csin%28x+1%29)
 - Choose the calculations accurracy. For example ![equation](http://latex.codecogs.com/gif.latex?%5Csmall%200.00001)
 - Choose some point near the expected root of equation. For example ![equation](http://latex.codecogs.com/gif.latex?%5Csmall%200)
@@ -124,10 +124,10 @@ Mathpack::NonlinearEquations.solve(start: 0.01, eps: 0.00001){|x| Math.exp(x-2)
 #### solve_system(params = {})
 returns solution of system of nonlinear equations by *Newton method*
 ##### Parameters
-- ***start*** - vector to start iteration process
-- ***eps*** - calculations accuraccy
-- ***f*** - vector of right part lambdas
-- ***w_matrix*** - matrix *W* in Newton method
+* ***start*** - vector to start iteration process
+* ***eps*** - calculations accuraccy
+* ***f*** - vector of right part lambdas
+* ***w_matrix*** - matrix *W* in Newton method
 
 ### Usage
 If you have system of equations ![equation](http://latex.codecogs.com/gif.latex?f_%7Bi%7D%28x%29%20%3D%200%2C%20i%20%3D%200%2C%201%2C%20...%2C%20N-1)
@@ -155,12 +155,12 @@ Mathpack::NonlinearEquations.solve_system(start: [1, 5], eps: 1e-4, f: f, w_matr
 #### solve_fredholm_2(params={})
 returns solution of integral equation as a hash of nodes and values arrays.
 ##### Parameters
-- ***from*** - left border
-- ***to*** - right border
-- ***lambda*** - *lambda* parameter
-- ***k*** - kernel function (2 arguements)
-- ***f*** - inhomogeneity function (1 arguement)
-- ***eps*** - accuracy
+* ***from*** - left border
+* ***to*** - right border
+* ***lambda*** - *lambda* parameter
+* ***k*** - kernel function (2 arguements)
+* ***f*** - inhomogeneity function (1 arguement)
+* ***eps*** - accuracy
 
 ### Usage
 Let we have the following integral equation
@@ -178,12 +178,12 @@ Mathpack::IntegralEquations.solve_fredholm_2(from: 1.0, to: 2.0, lambda: 0.5, ep
 #### solve_volter_2(params{})
 returns solution of integral equation as a hash of nodes and values arrays.
 ##### Parameters
-- ***from*** - left border
-- ***to*** - right border
-- ***lambda*** - *lambda* parameter
-- ***k*** - kernel function (2 arguements)
-- ***f*** - inhomogeneity function (1 arguement)
-- ***eps*** - accuracy
+* ***from*** - left border
+* ***to*** - right border
+* ***lambda*** - *lambda* parameter
+* ***k*** - kernel function (2 arguements)
+* ***f*** - inhomogeneity function (1 arguement)
+* ***eps*** - accuracy
 
 ### Usage
 Let we have the following integral equation
@@ -202,11 +202,11 @@ Mathpack::IntegralEquations.solve_volter_2(from: 0.0, to: 1.0, lambda: 0.5, eps:
 #### solve_cauchie_system(params={})
 returns solution of differential equations system as a hash of nodes and values arrays for each function. For example, { x: [], u1: [], u2: [], ...  }
 ##### Parameters
-- ***from*** - left border
-- ***to*** - right border
-- ***system*** - array of lambdas representing each row of system
-- ***y0*** - array of values of derivatives on left border. Starts with first derivative
-- ***eps*** - accuracy
+* ***from*** - left border
+* ***to*** - right border
+* ***system*** - array of lambdas representing each row of system
+* ***y0*** - array of values of derivatives on left border. Starts with first derivative
+* ***eps*** - accuracy
 
 ### Usage
 Let we have following system of differential equations
@@ -228,8 +228,8 @@ Mathpack::DifferentialEquations.solve_cauchie_system(from: 1.0, to: 2.0, eps: 1e
 #### solve(params = {})
 returns solution of system of linear equations.
 ##### Parameters
-- ***matrix*** - system matrix
-- ***f*** - right part vector
+* ***matrix*** - system matrix
+* ***f*** - right part vector
 
 ### Usage
 Let's solve some system of linear equations. It can be written as
@@ -256,16 +256,16 @@ Mathpack::SLE.solve(matrix: a, f: b) #=> Matrix[[-1.0, 2.0, 4.0]]
 #### approximate_by_polynom(params = {})
 returns array of coefficients of polynom, approximating given function on [from, to] segment.
 ##### Parameters
-- ***x*** - array of approximation nodes
-- ***polynom_power*** - power of approximation polynom
-- ***f*** - functions values in *x* if you have table function
+* ***x*** - array of approximation nodes
+* ***polynom_power*** - power of approximation polynom
+* ***f*** - functions values in *x* if you have table function
 
 #### generate_nodes(params = {}) 
 returns nodes for approximation with some step.
 ##### Parameters
-- ***from*** - first node
-- ***to*** - last node
-- ***step*** - step between nodes
+* ***from*** - first node
+* ***to*** - last node
+* ***step*** - step between nodes
 
 #### print_polynom(coefficients)
 returns a string representing polynom with given coefficients.
@@ -297,8 +297,8 @@ Mathpack::Approximation.print_polynom(result) #=> x^2
 #### integrate(params = {})
 returns integral value.
 ##### Parameters
-- ***from*** - start of integration
-- ***to*** - end of integration
+* ***from*** - start of integration
+* ***to*** - end of integration
 
 ### Usage
 Let you have the following integral:
@@ -321,10 +321,10 @@ Mathpack::Integration.integrate(from: -Float::INFINITY, to: Float::INFINITY){ |x
 #### print_table_function(params = {})
 writes table function values to file
 ##### Parameters
-- ***filename*** - name of output file
-- ***x*** - arguements array
-- ***y*** - function values array
-- ***labels*** - hash of labels names for *x* and *y* column
+* ***filename*** - name of output file
+* ***x*** - arguements array
+* ***y*** - function values array
+* ***labels*** - hash of labels names for *x* and *y* column
 
 #### read_table_function(filename)
 returns table function values hash, written to **filename**

data/lib/mathpack/functions.rb

@@ -21,7 +21,7 @@ module Mathpack
     end
 
     def self.heaviside(x)
-      x <= 0 ? 0 : 1
+      x <= 0.0 ? 0 : 1
     end
   end
 end

data/lib/mathpack/statistics.rb

@@ -10,19 +10,19 @@ module Mathpack
     end
 
     def mean
-      raw_moment(1)
+      @mean ||= raw_moment(1)
     end
 
     def variance
-      central_moment(2)
+      @variance ||= central_moment(2)
     end
 
     def skewness
-      central_moment(3) / variance**1.5
+      @skewness ||= central_moment(3) / variance**1.5
     end
 
     def kurtosis
-      central_moment(4) / variance**2 - 3.0
+      @kurtosis ||= central_moment(4) / variance**2 - 3.0
     end
 
     def max
@@ -35,25 +35,22 @@ module Mathpack
 
     def raw_moment(power)
       raw_moment = 0.0
-      @data_set.each_index do |i|
-        raw_moment += @frequency[i] * @data_set[i]**power
+      @data_set.zip(@frequency).each do |value, frequency|
+        raw_moment += frequency * value**power
       end
       raw_moment / number
     end
 
     def central_moment(power)
       central_moment = 0.0
-      m = mean
-      @data_set.each_index do |i|
-        central_moment += @frequency[i] * (@data_set[i] - m)**power
+      @data_set.zip(@frequency).each do |value, frequency|
+        central_moment += frequency * (value - mean)**power
       end
       central_moment / number
     end
 
     def empirical_cdf(x)
-      result = 0.0
-      @series.each { |val| result += Mathpack::Functions.heaviside(x - val) }
-      result / number
+      1.0 / number * @series.count { |value| Mathpack::Functions.heaviside(x - value) > 0 }
     end
 
     def print_empirical_cdf(filename)
@@ -65,9 +62,7 @@ module Mathpack
 
     def empirical_pdf(x)
       h = variance**0.5 * number**(-1.0 / 6)
-      result = 0.0
-      @series.each { |val| result += (Mathpack::Functions.heaviside(x - val + h) - Mathpack::Functions.heaviside(x - val - h)) / (2 * h) }
-      result / number
+      1.0 / number * @series.inject(0) { |sum, value| sum + (Mathpack::Functions.heaviside(x - value + h) - Mathpack::Functions.heaviside(x - value - h)) / (2 * h) }
     end
 
     def print_empirical_pdf(filename)

data/lib/mathpack/version.rb

@@ -1,3 +1,3 @@
 module Mathpack
-  VERSION = '0.4.9'
+  VERSION = '0.5.0'
 end

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 'bundler', '~> 1.11.2'
   spec.add_development_dependency 'rake', '~>0'
 end

data/spec/sle_spec.rb

@@ -13,15 +13,15 @@ describe 'SLE' do
     end
 
     it 'raises error when matrix is singular' do
-      expect { Mathpack::SLE.solve(matrix: unsolved, f: b) }.to raise_error
+      expect { Mathpack::SLE.solve(matrix: unsolved, f: b) }.to raise_error(RuntimeError)
     end
 
     it 'raises error when matrix has incorrect size' do
       b << 1
-      expect { Mathpack::SLE.solve(matrix: a, f: b) }.to raise_error
+      expect { Mathpack::SLE.solve(matrix: a, f: b) }.to raise_error(RuntimeError)
       b.pop
       a << [1, 2, 3, 4]
-      expect { Mathpack::SLE.solve(matrix: a, f: b) }.to raise_error
+      expect { Mathpack::SLE.solve(matrix: a, f: b) }.to raise_error(RuntimeError)
     end
 
     it 'returns vector if matrix class is given' do

metadata

@@ -1,41 +1,41 @@
 --- !ruby/object:Gem::Specification
 name: mathpack
 version: !ruby/object:Gem::Version
-  version: 0.4.9
+  version: 0.5.0
 platform: ruby
 authors:
 - maxmilan
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2016-01-21 00:00:00.000000000 Z
+date: 2016-03-11 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - "~>"
+    - - ~>
       - !ruby/object:Gem::Version
-        version: '1.6'
+        version: 1.11.2
   type: :development
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - "~>"
+    - - ~>
       - !ruby/object:Gem::Version
-        version: '1.6'
+        version: 1.11.2
 - !ruby/object:Gem::Dependency
   name: rake
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - "~>"
+    - - ~>
       - !ruby/object:Gem::Version
         version: '0'
   type: :development
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - "~>"
+    - - ~>
       - !ruby/object:Gem::Version
         version: '0'
 description: Includes collection of mathematical methods
@@ -45,8 +45,8 @@ executables: []
 extensions: []
 extra_rdoc_files: []
 files:
-- ".gitignore"
-- ".travis.yml"
+- .gitignore
+- .travis.yml
 - Gemfile
 - LICENSE.txt
 - README.md
@@ -84,17 +84,17 @@ require_paths:
 - lib
 required_ruby_version: !ruby/object:Gem::Requirement
   requirements:
-  - - ">="
+  - - ! '>='
     - !ruby/object:Gem::Version
       version: '0'
 required_rubygems_version: !ruby/object:Gem::Requirement
   requirements:
-  - - ">="
+  - - ! '>='
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
 rubyforge_project: 
-rubygems_version: 2.4.8
+rubygems_version: 2.5.1
 signing_key: 
 specification_version: 4
 summary: Summary