symbo 0.1.0pre

data/lib/symbo/sqrt.rb

@@ -0,0 +1,16 @@
+# frozen_string_literal: true
+
+require 'symbo/fraction'
+require 'symbo/power'
+
+module Symbo
+  class Sqrt
+    def self.[](x) # rubocop:disable Naming/MethodParameterName
+      Symbo::Power[x, Symbo::Fraction[1, 2]]
+    end
+  end
+
+  def √(x) # rubocop:disable Naming/MethodName, Naming/BinaryOperatorParameterName, Naming/MethodParameterName
+    Sqrt[x]
+  end
+end

data/lib/symbo/sum.rb

@@ -0,0 +1,211 @@
+# frozen_string_literal: true
+
+require 'symbo/expression'
+require 'symbo/factorial'
+require 'symbo/function'
+require 'symbo/mergeable'
+require 'symbo/product'
+
+module Symbo
+  # シンボリックな和
+  class Sum < Expression
+    include Mergeable
+
+    using Symbo
+
+    # :section: Power Transformation Methods
+
+    # べき乗の低
+    #
+    #   (:x + :y).base # => :x + :y
+    def base
+      dup
+    end
+
+    # べき指数
+    #
+    #   (:x + :y).exponent # => 1
+    def exponent
+      1
+    end
+
+    # :section: Basic Distributive Transformation Methods
+
+    # 同類項の項部分
+    #
+    #   (:x + :y).term # => Product[:x + :y]
+    def term
+      Product[self]
+    end
+
+    # 同類項の定数部分
+    #
+    #   (:x + :y).const # => 1
+    def const
+      1
+    end
+
+    # :section: Order Relation Methods
+
+    # 交換法則によるオペランド並べ替えに使う順序関係
+    #
+    # - 相手が和の場合
+    # 最右のオペランドから順に compare していき、異なるものがあればそれで順序を決定する。
+    # どちらかのオペランドがなくなれば、短いほうが左側。
+    #
+    #   (:a + :b).compare(:a + :c) # => true
+    #   Sum[:a, :c, :d].compare(Sum[:b, :c, :d]) # => true
+    #   (:c + :d).compare(Sum[:b, :c, :d]) # => true
+    #
+    # - 階乗、関数、シンボルの場合
+    # 相手を単項の和にして比較
+    #
+    #   (1 + :x).compare(:y) # => true
+    #
+    # - それ以外の場合
+    # 次のルールで比較
+    #
+    #   !other.compare(self)
+    #
+    # rubocop:disable Metrics/CyclomaticComplexity
+    def compare(other)
+      case other
+      when Sum
+        return @operands.last.compare(other.operands.last) if @operands.last != other.operands.last
+
+        m = length
+        n = other.length
+        if [m, n].min >= 2
+          1.upto([m, n].min) do |j|
+            return operand(m - j).compare(other.operand(n - j)) if operand(m - j) != other.operand(n - j)
+          end
+        end
+        m < n
+      when Factorial, Function, Symbol
+        compare Sum[other]
+      else
+        !other.compare(self)
+      end
+    end
+    # rubocop:enable Metrics/CyclomaticComplexity
+
+    # :section:
+
+    # rubocop:disable Metrics/PerceivedComplexity
+    # rubocop:disable Metrics/CyclomaticComplexity
+    def evaluate
+      v = @operands[0].evaluate
+      w = @operands[1].evaluate
+
+      if (v.integer? || v.complex?) && (w.integer? || w.complex?)
+        v.plus w
+      elsif v.constant? && w.constant?
+        Fraction[Sum[Product[v.numerator, w.denominator].evaluate, Product[w.numerator, v.denominator].evaluate].evaluate,
+                 Product[v.denominator, w.denominator].evaluate].evaluate
+      else
+        Sum[v.simplify, w.simplify]
+      end
+    end
+    # rubocop:enable Metrics/PerceivedComplexity
+    # rubocop:enable Metrics/CyclomaticComplexity
+
+    def simplify_rational_number
+      self
+    end
+
+    def to_s
+      elements = @operands.map do |each|
+        case each
+        when Sum
+          if each.length > 1
+            [' + ', "(#{each})"]
+          else
+            [' + ', each.to_s]
+          end
+        when Product
+          if each.length > 1 && each.operand(0) == -1
+            [' - ', Product[*each.operands[1..-1]].to_s]
+          else
+            [' + ', each.to_s]
+          end
+        else
+          [' + ', each.to_s]
+        end
+      end
+      elements.flatten[1..-1].join
+    end
+
+    protected
+
+    def _simplify
+      return UNDEFINED if @operands.include?(UNDEFINED)
+      return operand(0) if length == 1
+
+      v = simplify_rec(@operands)
+      if v.size == 1
+        v[0]
+      elsif v.size > 1
+        Sum[*v]
+      else
+        0
+      end
+    end
+
+    private
+
+    # rubocop:disable Metrics/PerceivedComplexity
+    # rubocop:disable Metrics/CyclomaticComplexity
+    def simplify_rec(l) # rubocop:disable Naming/MethodParameterName
+      case l
+      in Constant, Constant
+        p = Sum[*l].simplify_rne
+        p.zero? ? [] : [p]
+      in u1, u2 if u1.zero? && !u2.sum?
+        [u2]
+      in u1, u2 if u2.zero? && !u1.sum?
+        [u1]
+      in u1, u2 if l.none?(:sum?) && u1.term == u2.term
+        s = Sum[u1.const, u2.const].simplify
+        p = Product[u1.term, s].simplify
+        p.zero? ? [] : [p]
+      in u1, u2 if l.none?(&:sum?) && u2.compare(u1)
+        [u2, u1]
+      in _, _ if l.none?(&:sum?)
+        l
+      in Sum => u1, Sum => u2
+        merge u1.operands, u2.operands
+      in Sum => u1, u2
+        merge u1.operands, [u2]
+      in u1, Sum => u2
+        merge [u1], u2.operands
+      in Sum => u1, *rest
+        merge u1.operands, simplify_rec(rest)
+      in u1, *rest
+        merge [u1], simplify_rec(rest)
+      end
+    end
+    # rubocop:enable Metrics/PerceivedComplexity
+    # rubocop:enable Metrics/CyclomaticComplexity
+
+    # rubocop:disable Metrics/PerceivedComplexity
+    def simplify_rne_rec
+      if length == 1
+        v = operand(0).simplify_rne_rec
+        if v.undefined?
+          UNDEFINED
+        else
+          v
+        end
+      elsif length == 2
+        v = operand(0).simplify_rne_rec
+        w = operand(1).simplify_rne_rec
+        if v.undefined? || w.undefined?
+          UNDEFINED
+        else
+          Sum[v, w].evaluate
+        end
+      end
+    end
+    # rubocop:enable Metrics/PerceivedComplexity
+  end
+end

data/lib/symbo/symbol.rb

@@ -0,0 +1,102 @@
+# frozen_string_literal: true
+
+require 'symbo/algebraic_operators'
+require 'symbo/expression_type'
+require 'symbo/relational_operators'
+
+module Symbo
+  refine Symbol do
+    include AlgebraicOperators
+    include ExpressionType
+    include RelationalOperators
+
+    # :section: Simplification Methods
+
+    def simplify
+      self
+    end
+
+    # :section: Power Transformation Methods
+
+    # べき乗の低
+    #
+    #   :x.base # => :x
+    def base
+      self
+    end
+
+    # べき指数
+    #
+    #   :x.exponent # => 1
+    def exponent
+      1
+    end
+
+    # :section: Basic Distributive Transformation Methods
+
+    # 同類項の項部分
+    #
+    #   :x.term # => Product[:x]
+    def term
+      Product[self]
+    end
+
+    # 同類項の定数部分
+    #
+    #   :x.const # => 1
+    def const
+      1
+    end
+
+    # :section: Order Relation Methods
+
+    # 交換法則によるオペランド並べ替えに使う順序関係
+    #
+    # - 相手がシンボルの場合
+    # 辞書順で順序を決定
+    #
+    #   :a.compare(:b) # => true
+    #   :A.compare(:a) # => true
+    #   :v1.compare(:v2) # => true
+    #   :x1.compare(:xa) # => true
+    #
+    # - それ以外の場合
+    #   :x.compare(:x**2) # => true
+    #   :x.compare(Function(:x, :t)) # => true
+    #   :x.compare(Function(:y, :t)) # => true
+    def compare(other)
+      if other.is_a?(Symbol)
+        if self == :π
+          -1
+        else
+          self < other
+        end
+      else
+        !other.compare(self)
+      end
+    end
+
+    # :section: Relational Operator Methods
+
+    def ==(other)
+      other.is_a?(Symbol) && to_s == other.to_s
+    end
+
+    # :category: Simplification Methods
+
+    def simplify_rational_number
+      raise unless self == UNDEFINED
+
+      UNDEFINED
+    end
+  end
+end
+
+# Matrix などの中で使われる Symbol#+ などをハイジャック
+class Symbol
+  include Symbo::ExpressionType
+
+  def +(other)
+    Symbo::Sum[self, other]
+  end
+end

data/lib/symbo/tensor_product.rb

@@ -0,0 +1,15 @@
+# frozen_string_literal: true
+
+module Symbo
+  class TensorProduct
+    def self.[](matrix0, matrix1)
+      col_ms = (0...matrix0.column_size).map do |col|
+        ms = (0...matrix0.row_size).map do |row|
+          matrix0[row, col] * matrix1
+        end
+        Matrix.vstack(*ms)
+      end
+      Matrix.hstack(*col_ms)
+    end
+  end
+end

data/lib/symbo/trigonometric_function.rb

@@ -0,0 +1,26 @@
+# frozen_string_literal: true
+
+require 'symbo/function'
+
+module Symbo
+  class TrigonometricFunction < Function
+    using Symbo
+
+    def x
+      @operands[1]
+    end
+
+    def simplify
+      self.class.new(x.simplify)._simplify
+    end
+
+    private
+
+    # rubocop:disable Metrics/AbcSize
+    def kn_pi?
+      x.product? && x.length == 2 && x.operand(0).constant? && x.operand(1) == PI &&
+        [1, 2, 3, 4, 6].include?(x.operand(0).denominator) && x.operand(0).numerator.integer?
+    end
+    # rubocop:enable Metrics/AbcSize
+  end
+end

data/lib/symbo/version.rb

@@ -0,0 +1,5 @@
+# frozen_string_literal: true
+
+module Symbo
+  VERSION = '0.1.0pre'
+end

data/symbo.gemspec

@@ -0,0 +1,33 @@
+# frozen_string_literal: true
+
+lib = File.expand_path('lib', __dir__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+
+require 'symbo/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = 'symbo'
+  spec.version       = Symbo::VERSION
+  spec.authors       = ['Yasuhito Takamiya']
+  spec.email         = ['yasuhito@gmail.com']
+
+  spec.summary       = 'A simple CAS (computer algebra system) in Ruby'
+  spec.homepage      = 'https://github.com/yasuhito/symbo'
+  spec.license       = 'MIT'
+
+  spec.metadata['homepage_uri'] = spec.homepage
+  spec.metadata['source_code_uri'] = 'https://github.com/yasuhito/symbo'
+  spec.metadata['changelog_uri'] = 'https://github.com/yasuhito/symbo/blob/develop/CHANGELOG.md'
+
+  # Specify which files should be added to the gem when it is released.
+  # The `git ls-files -z` loads the files in the RubyGem that have been added into git.
+  spec.files = Dir.chdir(File.expand_path(__dir__)) do
+    `git ls-files -z`.split("\x0").reject { |f| f.match(%r{^(test|spec|features)/}) }
+  end
+  spec.bindir        = 'exe'
+  spec.executables   = spec.files.grep(%r{^exe/}) { |f| File.basename(f) }
+  spec.require_paths = ['lib']
+
+  spec.add_development_dependency 'bundler'
+  spec.add_development_dependency 'rake'
+end

metadata

@@ -0,0 +1,116 @@
+--- !ruby/object:Gem::Specification
+name: symbo
+version: !ruby/object:Gem::Version
+  version: 0.1.0pre
+platform: ruby
+authors:
+- Yasuhito Takamiya
+autorequire: 
+bindir: exe
+cert_chain: []
+date: 2020-03-03 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  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'
+- !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: 
+email:
+- yasuhito@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- ".github/workflows/ruby.yml"
+- ".gitignore"
+- ".rubocop.yml"
+- ".ruby-version"
+- CHANGELOG.md
+- Gemfile
+- Gemfile.lock
+- LICENSE
+- README.md
+- Rakefile
+- bin/console
+- bin/setup
+- lib/symbo.rb
+- lib/symbo/algebraic_operators.rb
+- lib/symbo/complex.rb
+- lib/symbo/constant.rb
+- lib/symbo/cos.rb
+- lib/symbo/diff.rb
+- lib/symbo/e.rb
+- lib/symbo/expression.rb
+- lib/symbo/expression_type.rb
+- lib/symbo/factorial.rb
+- lib/symbo/float.rb
+- lib/symbo/fraction.rb
+- lib/symbo/function.rb
+- lib/symbo/integer.rb
+- lib/symbo/matrix.rb
+- lib/symbo/mergeable.rb
+- lib/symbo/pi.rb
+- lib/symbo/power.rb
+- lib/symbo/product.rb
+- lib/symbo/qubit.rb
+- lib/symbo/quot.rb
+- lib/symbo/relational_operators.rb
+- lib/symbo/sin.rb
+- lib/symbo/sqrt.rb
+- lib/symbo/sum.rb
+- lib/symbo/symbol.rb
+- lib/symbo/tensor_product.rb
+- lib/symbo/trigonometric_function.rb
+- lib/symbo/version.rb
+- symbo.gemspec
+homepage: https://github.com/yasuhito/symbo
+licenses:
+- MIT
+metadata:
+  homepage_uri: https://github.com/yasuhito/symbo
+  source_code_uri: https://github.com/yasuhito/symbo
+  changelog_uri: https://github.com/yasuhito/symbo/blob/develop/CHANGELOG.md
+post_install_message: 
+rdoc_options: []
+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: 1.3.1
+requirements: []
+rubygems_version: 3.1.2
+signing_key: 
+specification_version: 4
+summary: A simple CAS (computer algebra system) in Ruby
+test_files: []