dydx 0.1.31 → 0.1.314

data/lib/dydx/helper.rb

@@ -141,5 +141,10 @@ module Dydx
     def ==(x)
       to_s == x.to_s
     end
+
+    # Refactor
+    def **(x)
+      self ^ (x)
+    end
   end
 end

data/lib/dydx/integrand.rb

@@ -7,13 +7,14 @@ module Dydx
     end
 
     def [](a, b, n = 1000)
-      raise ArgumentError, 'b should be greater than a' if a > b
-      f = function
       # HOT FIX: should implement Infinity class
       a = - 1000 if a == - Float::INFINITY
       b = 1000 if b == Float::INFINITY
 
       a, b = [a, b].map(&:to_f)
+      raise ArgumentError, 'b should be greater than a' if a > b
+      $int_f = function
+
       n = [n, (b - a) * 2].max
       n += 1 if n.to_i.odd?
       h = (b - a) / n
@@ -24,6 +25,10 @@ module Dydx
       round_8( (h / 3) * (f(a) + f(b) + 2 * even_sum + 4 * odd_sum) )
     end
 
+    def f(vars)
+      int_f(vars)
+    end
+
     def round_8(num)
       return num if num.abs == Float::INFINITY
       (num * 10 ** 8).round * 10.0 ** (-8)

data/lib/dydx/version.rb

@@ -1,3 +1,3 @@
 module Dydx
-  VERSION = "0.1.31"
+  VERSION = "0.1.314"
 end

data/spec/dydx_spec.rb

@@ -5,35 +5,16 @@ describe Dydx do
     expect(Dydx::VERSION).not_to be nil
   end
 
-  context 'ex1' do
-    $a = (:x ^ :n)
-    let(:d1){ da/dx }
-    let(:d2){ d/dx($a) }
-    it{ expect(d1.to_s).to eq('( n * ( x ^ ( n - 1 ) ) )') }
-    it{ expect(d2.to_s).to eq('( n * ( x ^ ( n - 1 ) ) )') }
-  end
-
-  context 'ex2' do
-    $b = (:x ^ (:x * 2))
-    let(:d1){ db/dx }
-    let(:d2){ d/dx($b) }
-    it{ expect(d1.to_s).to eq('( ( 2 * x ) * ( x ^ ( ( 2 * x ) - 1 ) ) )') }
-    it{ expect(d2.to_s).to eq('( ( 2 * x ) * ( x ^ ( ( 2 * x ) - 1 ) ) )') }
-  end
-
-  context 'ex3' do
-    $c = (:t ^ 2) / 2
-    let(:d1){ dc/dt }
-    let(:d2){ d/dt($c) }
-    it{ expect(d1.to_s).to eq('t') }
-    it{ expect(d2.to_s).to eq('t') }
-  end
+  it 'demo' do
+    f(x, y) <= x + x*y + y
+    expect(f(x, y)).to eq(x * (1 + y) + y)
+    expect(f(a, 2)).to eq(3*a + 2)
+    expect(f(1, a + b)).to eq(1 + 2 * ( a + b ))
+    expect(f(1, 2)).to eq(5)
+    expect(d/dx(f(x, y))).to eq(1 + y)
 
-  context 'ex4' do
-    $i = 2 * (e ^ (2 * :z))
-    let(:d1){ di/dz }
-    let(:d2){ d/dz($i) }
-    it{ expect(d1.to_s).to eq('( 4 * ( e ^ ( 2 * z ) ) )') }
-    it{ expect(d2.to_s).to eq('( 4 * ( e ^ ( 2 * z ) ) )') }
+    g(x) <= sin(x)
+    expect(d/dx(g(x))).to eq(cos(x))
+    expect(S(g(x), dx)[0, pi/2]).to eq(1.0)
   end
 end

data/spec/lib/algebra/formula_spec.rb

@@ -28,6 +28,12 @@ describe Dydx::Algebra::Formula do
     it{ expect( (addition * multiplication).to_s ).to eq('( ( x + y ) * ( x * y ) )') }
   end
 
+  describe '#subst' do
+    it{ expect((x + y).subst(x: 3, y: 3)).to eq(6) }
+    it{ expect((x + y).subst(x: 3)).to eq(3 + y) }
+    it{ expect((x + y + pi).subst(x: 3, y: 3).to_f).to eq(Math::PI + 6) }
+  end
+
   describe '#differentiate' do
     it{ expect(addition.d(:x)).to eq(1) }
     it{ expect(addition.d(:y)).to eq(1) }

data/spec/lib/algebra/set_spec.rb

@@ -0,0 +1,227 @@
+require 'spec_helper'
+# HOTFIX: Refactor by using context or describe
+describe Dydx::Algebra::Set do
+  it{ expect(_(1)).to eq(_(1)) }
+  it{ expect(_(-1)).to eq(_(-1)) }
+  it{ expect(e0).to eq(e0) }
+  it{ expect(e1).to eq(e1) }
+
+  it{ expect(e).to eq(e) }
+  it{ expect(pi).to eq(pi) }
+
+  it { expect(log10(1)).to eq(e0) }
+  it { expect(log10(10)).to eq(e1) }
+  it { expect(log10(10 ^ n)).to eq(n) }
+  it { expect(log10(3 ^ n)).to eq(n * log10( 3 )) }
+
+  it { expect(log2(1)).to eq(e0) }
+  it { expect(log2(2)).to eq(e1) }
+  it { expect(log2(2 ^ n)).to eq(n) }
+  it { expect(log2(3 ^ n)).to eq(n * log2( 3 )) }
+
+  it{ expect(log(1)).to eq(0) }
+  it{ expect(log(e)).to eq(1) }
+  it{ expect(log(e ^ n)).to eq(n) }
+
+  it{ expect(sin(pi)).to eq(0) }
+
+  it{ expect(cos(0)).to eq(1) }
+  it{ expect(cos(pi)).to eq(-1) }
+  it{ expect(cos(2 * pi)).to eq(1) }
+
+  describe '#to_s' do
+    it{ expect(e.to_s).to eq('e') }
+
+    it{ expect(1.to_s).to eq('1') }
+
+    it{ expect(1.0.to_s).to eq('1.0') }
+
+    it{ expect(e1.to_s).to eq('1') }
+
+    it{ expect(pi.to_s).to eq('pi') }
+
+    it{ expect(cos(x).to_s).to eq('cos( x )') }
+    it{ expect(sin(x).to_s).to eq('sin( x )') }
+    it{ expect(tan(x).to_s).to eq('tan( x )') }
+
+    it{ expect(x.to_s).to eq('x') }
+  end
+
+  describe '#==' do
+    it{ expect(_(1) == _(1)).to be_true }
+  end
+
+  describe '#subst' do
+    it{ expect(1.subst).to eq(1) }
+    it{ expect(1.0.subst).to eq(1.0) }
+    it{ expect(e.subst).to eq(e) }
+    it{ expect(pi.subst).to eq(pi) }
+    it{ expect(sin(x).subst(x: 3)).to eq(sin(3)) }
+    it{ expect(cos(x).subst(x: pi)).to eq(-1) }
+    it{ expect(tan(0).subst(x: pi)).to eq(0) }
+    it{ expect(log(x).subst(x: e)).to eq(1) }
+    it{ expect(log10(x).subst(x: 7)).to eq(log10(7)) }
+    it{ expect(log2(2).subst(x: 2)).to eq(1) }
+    it{ expect(x.subst(x: 2)).to eq(2) }
+    it{ expect(x.subst(y: 2)).to eq(x) }
+  end
+
+  describe '#differentiate' do
+    it { expect(e.d).to eq(0) }
+    it { expect((e ^ x).d).to eq(e ^ x) }
+    it { expect((e ^ (x + y)).d).to eq(e ^ ( x + y )) }
+
+    it { expect(pi.d(x)).to eq(0) }
+
+    it { expect(1.d).to eq(0) }
+    it { expect(3.d).to eq(0) }
+    it { expect(3.0.d).to eq(0) }
+
+    it { expect(sin(x).d).to eq(cos(x)) }
+    it { expect(cos(x).d).to eq(- sin(x)) }
+    it { expect(tan(x).d).to eq(1 / cos(x) ^ 2) }
+
+    it { expect(log(x).d).to eq(1 / x) }
+    it { expect(log10(x).d).to eq(1 / ( x * log( 10 ) )) }
+    it { expect(log2(x).d).to eq(1 / ( x * log( 2 ) )) }
+
+    it{ expect(x.d(x)).to eq(1) }
+  end
+
+  describe 'Calculate' do
+    context 'E with Fixnum' do
+      it{ expect(e + 0).to eq(e) }
+      it{ expect(e - 0).to eq(e) }
+      it{ expect((e * 0).to_s).to eq('0') }
+      it{ expect(e * 1).to eq(e) }
+      it{ expect{(e / 0).to_s}.to raise_error(ZeroDivisionError) }
+      it{ expect(e / 1).to eq(e) }
+      it{ expect((e ^ 0).to_s).to eq('1') }
+    end
+
+    context 'Fixnum with Formula' do
+      let(:formula) { (:x + :y) }
+      it{ expect((0 + formula).to_s).to eq(formula.to_s) }
+      it{ expect((0 - formula).to_s).to eq('( - ( x + y ) )') }
+      it{ expect((0 * formula).to_s).to eq('0') }
+      it{ expect((1 * formula).to_s).to eq(formula.to_s) }
+      it{ expect((0 / formula).to_s).to eq('0') }
+      it{ expect((1 / formula).to_s).to eq('( 1 / ( x + y ) )') }
+      it{ expect((0 ^ formula).to_s).to eq('0') }
+      it{ expect((1 ^ formula).to_s).to eq('1') }
+    end
+
+    context 'Fixnum with Symbol' do
+      it{ expect(0 + :x).to eq(:x) }
+      it{ expect((0 - :x).to_s).to eq('( - x )') }
+      it{ expect((0 * :x).to_s).to eq('0') }
+      it{ expect(1 * :x).to eq(:x) }
+      it{ expect((0 / :x).to_s).to eq('0') }
+      it{ expect((1 / :x).to_s).to eq('( 1 / x )') }
+      it{ expect((0 ^ :x).to_s).to eq('0') }
+      it{ expect((1 ^ :x).to_s).to eq('1') }
+    end
+
+    context 'Fixnum with Fixnum' do
+      it{ expect(0 + 3).to eq(3) }
+      it{ expect(3 + 0).to eq(3) }
+      it{ expect(2 + 3).to eq(5) }
+
+      it{ expect(0 - 3).to eq(-3) }
+      it{ expect(3 - 0).to eq(3) }
+      it{ expect(2 - 3).to eq(-1) }
+
+      it{ expect(0 * 3).to eq(0) }
+      it{ expect(3 * 0).to eq(0) }
+      it{ expect(1 * 3).to eq(3) }
+      it{ expect(3 * 1).to eq(3) }
+      it{ expect(3 * 2).to eq(6) }
+
+      it{ expect((0 / 3).to_s).to eq('0') }
+      it{ expect{(3 / 0).to_s}.to raise_error(ZeroDivisionError) }
+      it{ expect((3 / 1).to_s).to eq('3') }
+      # TODO:
+      it{ expect((2 / 3).to_s).to eq('0') }
+
+
+      it{ expect((0 ^ 3).to_s).to eq('0') }
+      it{ expect((3 ^ 0).to_s).to eq('1') }
+      it{ expect((1 ^ 3).to_s).to eq('1') }
+      it{ expect((3 ^ 1).to_s).to eq('3') }
+      it{ expect((3 ^ 2).to_s).to eq('9') }
+    end
+
+    context 'Float with Formula' do
+      let(:formula) { (:x + :y) }
+      it{ expect((0.0 + formula).to_s).to eq(formula.to_s) }
+      it{ expect((0.0 - formula).to_s).to eq('( - ( x + y ) )') }
+      it{ expect((0.0 * formula).to_s).to eq('0') }
+      it{ expect((1.0 * formula).to_s).to eq(formula.to_s) }
+      it{ expect((0.0 / formula).to_s).to eq('0') }
+      it{ expect((1.0 / formula).to_s).to eq('( 1 / ( x + y ) )') }
+      it{ expect((0.0 ^ formula).to_s).to eq('0') }
+      it{ expect((1.0 ^ formula).to_s).to eq('1') }
+    end
+
+    context 'Float with Symbol' do
+      it{ expect(0.0 + :x).to eq(:x) }
+      it{ expect((0.0 - :x).to_s).to eq('( - x )') }
+      it{ expect((0.0 * :x).to_s).to eq('0') }
+      it{ expect(1.0 * :x).to eq(:x) }
+      it{ expect((0.0 / :x).to_s).to eq('0') }
+      it{ expect((1.0 / :x).to_s).to eq('( 1 / x )') }
+      it{ expect((0.0 ^ :x).to_s).to eq('0') }
+      it{ expect((1.0 ^ :x).to_s).to eq('1') }
+    end
+
+    context 'Float with Float' do
+      it{ expect(0.0 + 3.0).to eq(3.0) }
+      it{ expect(3.0 + 0.0).to eq(3.0) }
+      it{ expect(2.0 + 3.0).to eq(5.0) }
+
+      it{ expect(0.0 - 3.0).to eq(-3.0) }
+      it{ expect(3.0 - 0.0).to eq(3.0) }
+      it{ expect(2.0 - 3.0).to eq(-1.0) }
+
+      it{ expect(0.0 * 3.0).to eq(0.0) }
+      it{ expect(3.0 * 0.0).to eq(0.0) }
+      it{ expect(1.0 * 3.0).to eq(3.0) }
+      it{ expect(3.0 * 1.0).to eq(3.0) }
+      it{ expect(3.0 * 2.0).to eq(6.0) }
+
+      it{ expect(0.0 / 3.0).to eq(0.0) }
+      it{ expect(3.0 / 0.0).to eq(oo) }
+      it{ expect(3.0 / 1.0).to eq(3.0) }
+      # TODO:
+      it{ expect(2.0 / 3.0).to eq(0.6666666666666666) }
+
+
+      it{ expect(0.0 ^ 3.0).to eq(0.0) }
+      it{ expect(3.0 ^ 0.0).to eq(1.0) }
+      it{ expect(1.0 ^ 3.0).to eq(1.0) }
+      it{ expect(3.0 ^ 1.0).to eq(3.0) }
+      it{ expect(3.0 ^ 2.0).to eq(9.0) }
+    end
+
+    context 'Pi with Fixnum' do
+      it{ expect(pi + 0).to eq(pi) }
+      it{ expect(pi - 0).to eq(pi) }
+      it{ expect((pi * 0).to_s).to eq('0') }
+      it{ expect(pi * 1).to eq(pi) }
+      it{ expect{(pi / 0).to_s}.to raise_error(ZeroDivisionError) }
+      it{ expect(pi / 1).to eq(pi) }
+      it{ expect((pi ^ 0).to_s).to eq('1') }
+    end
+
+    context 'Symbol with Fixnum' do
+      it{ expect(:x + 0).to eq(:x) }
+      it{ expect(:x - 0).to eq(:x) }
+      it{ expect((:x * 0).to_s).to eq('0') }
+      it{ expect(:x * 1).to eq(:x) }
+      it{ expect{(:x / 0).to_s}.to raise_error(ZeroDivisionError) }
+      it{ expect(:x / 1).to eq(:x) }
+      it{ expect((:x ^ 0).to_s).to eq('1') }
+      it{ expect(:x ^ 1).to eq(:x) }
+    end
+  end
+end

data/spec/lib/delta_spec.rb

@@ -8,4 +8,30 @@ describe Dydx:Delta do
   it{ expect(dx(y).var).to eq(:x) }
   it{ expect(dx(y).function).to eq(:y) }
   it{ expect{dxy}.to raise_error(NameError) }
+
+  before(:each) do
+    $y = example
+  end
+
+  context 'ex1' do
+    let(:example) { x ^ n }
+    it { expect(dy/dx).to eq( n * ( x ^ ( n - 1 ) ) ) }
+    it { expect(d/dx($y)).to eq( n * ( x ^ ( n - 1 ) ) ) }
+  end
+
+  context 'ex2' do
+    let(:example) { x ^ (x * 2) }
+    it{ expect(dy/dx).to eq(( 2 * x ) * ( x ^ ( ( 2 * x ) - 1 ) )) }
+    it{ expect(d/dx($y)).to eq(( 2 * x ) * ( x ^ ( ( 2 * x ) - 1 ) )) }
+  end
+
+  context 'ex3' do
+    let(:example) { (t ^ 2) / 2 }
+    it{ expect(dy/dt).to eq(t) }
+  end
+
+  context 'ex4' do
+    let(:example) { 2 * (e ^ (2 * z)) }
+    it{ expect(dy/dz).to eq(4 * ( e ^ ( 2 * z ) )) }
+  end
 end

data/spec/lib/function_spec.rb

@@ -1,39 +1,67 @@
 require 'spec_helper'
 
 describe Dydx:Function do
-  # TODO: refactor
-  it{ expect(f(x, y)).to eq(f(x, y)) }
-  it{ expect(f(x, y)).to eq($f) }
-  it{ expect(f(x, y).algebra).to be_nil }
-  it{ expect(f(x, y).vars).to eq([:x, :y]) }
-  it{ expect{f(x, y, z)}.to raise_error(ArgumentError) }
-  it{ expect{f(x)}.to raise_error(ArgumentError) }
-  it{ expect(f(x, y) <= x * y ).to eq(f(x, y)) }
-  it{ expect(f(x, y)).to eq(f(x, y)) }
-  it{ expect(f(x, y)).to eq($f) }
-  it{ expect(f(x, y).algebra).to eq(x * y) }
-  it{ expect(f(x, y)).to eq(x * y) }
-  it{ expect(eval(f(x, y).to_s)).to eq(f(x, y)) }
-  it{ expect(f(x, y)).to eq(eval(f(x, y).to_s)) }
-  it{ expect(f(a, b)).to eq(a * b) }
-  it{ expect(f(2, 3)).to eq(6) }
-  it{ expect(f(a + b, c)).to eq((a + b) * c) }
-  it{ expect(d/dx(f(x, y))).to eq(y) }
-  it{ expect(d/dy(f(x, y))).to eq(x) }
-  it{ expect(d/dz(f(x, y))).to eq(0) }
-
-  it{ expect(g(a, b) <= f(a + b, b)).to eq(g(a, b)) }
-  it{ expect(g(a, b)).to eq($g) }
-  it{ expect(g(2, 3)).to eq(15) }
-
-  it{ expect(h(a, b, c) <= d/db(g(a, b))).to eq(h(a, b, c)) }
-  it{ expect(h(a, b, c) <= d/db(g(a, b))).to eq($h) }
-  it{ expect(h(a, b, c)).to eq(( a + ( 2 * b ) )) }
-  it{ expect(h(a, b, c).algebra).to eq(( a + ( 2 * b ) )) }
-
-  it 'ex.4' do
-    $f = nil
+  before(:each) do
+    reset
+  end
+
+  it 'ex1' do
+    expect(f(x, y)).to eq(f(x, y))
+    expect(f(x, y)).to eq($f)
+    expect(f(x, y).algebra).to be_nil
+    expect(f(x, y).vars).to eq([:x, :y])
+    expect{f(x, y, z)}.to raise_error(ArgumentError)
+    expect{f(x)}.to raise_error(ArgumentError)
+  end
+
+  it 'ex2' do
+    expect(f(x, y) <= x * y ).to eq(f(x, y))
+    expect(f(x, y)).to eq(f(x, y))
+    expect(f(x, y)).to eq($f)
+    expect(f(x, y).algebra).to eq(x * y)
+    expect(f(x, y)).to eq(x * y)
+    expect(f(x, y)).to eq(eval(f(x, y).to_s))
+    expect(f(a, b)).to eq(a * b)
+    expect(f(2, x)).to eq(2 * x)
+    expect(f(2, 3)).to eq(6)
+    expect(f(a + b, c)).to eq((a + b) * c)
+    expect(d/dx(f(x, y))).to eq(y)
+    expect(d/dy(f(x, y))).to eq(x)
+    expect(d/dz(f(x, y))).to eq(0)
+  end
+
+  it 'ex3' do
+    f(x, y) <= x * y
+    g(a, b) <= f(a + b, b)
+    expect(g(a, b)).to eq((a + b) * b)
+    expect(g(a, 3)).to eq(9 + 3 * a)
+    expect(g(2, 3)).to eq(15)
+  end
+
+  it 'ex4' do
+    f(a, b) <= (a + b) * b
+    h(a, b, c) <= d/db(f(a, b))
+    expect(h(a, b, c)).to eq(( a + ( 2 * b ) ))
+    expect(h(a, b, c).algebra).to eq(( a + ( 2 * b ) ))
+  end
+
+  it 'ex5' do
     f(x) <= log(x)
-    expect(f(a)).to eq(log(a))
+    expect(f(e)).to eq(1)
+    expect(f(0)).to eq(-oo)
+    expect(f(y)).to eq(log(y))
+
+    g(x) <= d/dx(f(x))
+    expect(g(1)).to eq(1)
+  end
+
+  it 'ex6' do
+    f(x) <= sin(x)
+    expect(f(pi)).to eq(0)
+    expect(f(pi / 2)).to eq(1)
+    expect(f(y)).to eq(sin(y))
+
+    g(x) <= d/dx(f(x))
+    expect(g(pi)).to eq(-1)
   end
 end