mathpack 0.4.5 → 0.4.6

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: a8ce50434e6550fb905aa1cb215106fc208af310
-  data.tar.gz: 2bce06f9c4b43424bf26923d0ecf7c91789253d3
+  metadata.gz: 1e848a2ca014612e04f7d621ff4667bc4dad2603
+  data.tar.gz: 85a146ea4bea8cb8720f4a7542c4582849f42f11
 SHA512:
-  metadata.gz: 8568cd24de9e9404a584e7f2723c62f604cfd5a9e8235f393f9f6438770b865dd0e569d8ce57d736f567b76c94edc7c1182efd484c17cc178bd6fdf01f2fcdc0
-  data.tar.gz: 5eb2691333cc2af964d925a4f7c3c6c75598c08c6229b6d321a6011e9003e7ceb89e463c13d99e47a7af149f2e4fbb44cf10bcb83ec4c9fe3ad05fa06dee9eff
+  metadata.gz: e7f5d93c396d007c6b634e3b68874eef3e64b29e7d3e8397de39f0325d3f24b303391d4f8155e3afd106340d2fd527372c586696d2fc4077b4357dbf38861e5b
+  data.tar.gz: 1f55b2454e8271adaa34f955cf161cd9e52807ba3090d0c5f925c53568de2a2eb04b5db01c4b8318d70b939794c4273bde6b31e914b0e7c6a1541da5a1068371

data/lib/mathpack.rb

@@ -7,3 +7,4 @@ require 'mathpack/integration'
 require 'mathpack/io'
 require 'mathpack/functions'
 require 'mathpack/functional'
+require 'mathpack/integral_equations'

data/lib/mathpack/integral_equations.rb

@@ -0,0 +1,59 @@
+module Mathpack
+  module IntegralEquations
+    def self.solve_fredholm_2(params = {})
+      step = (params[:to] - params[:from]) / 2.0
+      y_next = solve_freedholm_with_step(step, params)
+      loop do
+        y_prev = y_next
+        step /= 2.0
+        y_next = solve_freedholm_with_step(step, params)
+        break if Mathpack::IO.count_diff(y_next, y_prev) / (2**4 - 1) <= params[:eps]
+      end
+      { x: Mathpack::Approximation.generate_nodes(from: params[:from], to: params[:to], step: step), f: y_next }
+    end
+
+    def self.solve_freedholm_with_step(step, params = {})
+      nodes = Mathpack::Approximation.generate_nodes(from: params[:from], to: params[:to], step: step)
+      a = Array.new(nodes.length) do |i|
+        i % 2 == 1 ? 4.0 * step / 3.0 : 2.0 * step / 3.0
+      end
+      a[0], a[-1] = 1.0 / 3.0 * step, 1.0 / 3.0 * step
+      matrix = Array.new(nodes.length) { Array.new nodes.length }
+      for i in 0...matrix.length do
+        for j in 0...matrix[i].length do
+          matrix[i][j] = i == j ? 1.0 - params[:lambda] * a[i] * params[:k].call(nodes[i], nodes[i]) : - params[:lambda] * a[j] * params[:k].call(nodes[i], nodes[j])
+        end
+      end
+      f = Array.new(nodes.length){ |i| params[:f].call(nodes[i]) }
+      Mathpack::SLE.solve(matrix: matrix, f: f)
+    end
+
+    def self.solve_volter_2(params = {})
+      step = (params[:to] - params[:from]) / 2.0
+      y_next = solve_volter_with_step(step, params)
+      loop do
+        y_prev = y_next
+        step /= 2.0
+        y_next = solve_volter_with_step(step, params)
+        break if Mathpack::IO.count_diff(y_next, y_prev) / (2**2 - 1) <= params[:eps]
+      end
+      { x: Mathpack::Approximation.generate_nodes(from: params[:from], to: params[:to], step: step), f: y_next }
+    end
+
+    def self.solve_volter_with_step(step, params = {})
+      y = []
+      nodes = Mathpack::Approximation.generate_nodes(from: params[:from], to: params[:to], step: step)
+      y[0] = params[:f].call(nodes[0])
+      y[1] = (params[:f].call(nodes[1]) + 0.5 * params[:lambda] * step * params[:k].call(nodes[1], nodes[0]) * y[0]) / (1.0 - 0.5 * params[:lambda] * step * params[:k].call(nodes[1], nodes[1]))
+      for i in 2...nodes.length do
+        sum = 0.0
+        for j in 1..(i-1) do
+          sum += params[:k].call(nodes[i], nodes[j]) * y[j]
+        end
+        sum *= step
+        y[i] = (0.5 * params[:lambda] * step * params[:k].call(nodes[i], nodes[0]) * y[0] + params[:lambda] * sum + params[:f].call(nodes[i])) / (1.0 - 0.5 * params[:lambda] * step * params[:k].call(nodes[i], nodes[i]))
+      end
+      y
+    end
+  end
+end

data/lib/mathpack/version.rb

@@ -1,3 +1,3 @@
 module Mathpack
-  VERSION = '0.4.5'
+  VERSION = '0.4.6'
 end

data/spec/integral_equations_spec.rb

@@ -0,0 +1,27 @@
+describe 'IntegralEquation' do
+  require 'mathpack/integral_equations'
+
+  describe '#solve_fredholm_2' do
+    let(:eps) { 1e-5 }
+
+    it 'solves fredholm 2 equations correctly' do
+      k = -> x, t { 1.0 / (x + t**2)**0.5 }
+      f = -> x { 2*x + (x + 1)**0.5 - (x + 4)**0.5 }
+      result = Mathpack::IntegralEquations.solve_fredholm_2(from: 1, to: 2, lambda: 0.5, eps: eps, k: k, f: f)[:f]
+      solution = Mathpack::Approximation.fill_f(Mathpack::Approximation.generate_nodes(from: 1, to: 2, step: 1.0 / (result.count - 1))) { |x| 2*x }
+      expect(Mathpack::IO.count_diff(result, solution) <= eps).to eq(true)
+    end
+  end
+
+  describe '#solve_volter_2' do
+    let(:eps) { 1e-5 }
+
+    it 'solves volter 2 equations correctly' do
+      k_volter = -> x, t { Math.cos(x - t) }
+      f_volter = -> x { 0.75 * Math.exp(x) + 0.25 * Math.cos(x) - 0.25 * Math.sin(x) }
+      result = Mathpack::IntegralEquations.solve_volter_2(from: 0, to: 1, lambda: 0.5, eps: eps, k: k_volter, f: f_volter)[:f]
+      solution = Mathpack::Approximation.fill_f(Mathpack::Approximation.generate_nodes(from: 0, to: 1, step: 1.0 / (result.count - 1))) { |x| Math.exp(x) }
+      expect(Mathpack::IO.count_diff(result, solution) <= eps).to eq(true)
+    end
+  end
+end

data/spec/nonlinear_equations_spec.rb

@@ -1,4 +1,4 @@
-describe 'Equation' do
+describe 'NonlinearEquation' do
   require 'mathpack/nonlinear_equations'
 
   describe '#solve' do

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: mathpack
 version: !ruby/object:Gem::Version
-  version: 0.4.5
+  version: 0.4.6
 platform: ruby
 authors:
 - maxmilan
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2015-04-05 00:00:00.000000000 Z
+date: 2015-05-08 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
@@ -55,6 +55,7 @@ files:
 - lib/mathpack/approximation.rb
 - lib/mathpack/functional.rb
 - lib/mathpack/functions.rb
+- lib/mathpack/integral_equations.rb
 - lib/mathpack/integration.rb
 - lib/mathpack/io.rb
 - lib/mathpack/nonlinear_equations.rb
@@ -65,6 +66,7 @@ files:
 - spec/approximation_spec.rb
 - spec/functional_spec.rb
 - spec/functions_spec.rb
+- spec/integral_equations_spec.rb
 - spec/integration_spec.rb
 - spec/io_spec.rb
 - spec/nonlinear_equations_spec.rb
@@ -90,7 +92,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
       version: '0'
 requirements: []
 rubyforge_project: 
-rubygems_version: 2.4.6
+rubygems_version: 2.4.3
 signing_key: 
 specification_version: 4
 summary: Summary