hpsqrt 1.0.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: 44cd55553d14006a807f7f752e7a9ff5d72f3059da918af90c3a4b78d00bcfcb
+  data.tar.gz: 024a8bc0f8b197334342f78b6f91d761d2eb29a247c135c2b7f1ff45a39ce4c6
+SHA512:
+  metadata.gz: dc828d520080593730770dee6ad5977209c6e37705551500942ccf7f5e09ef572e60f47d196c2d47fe37d50227f69b373a90a23e2ff58ea5c1b4fdb808856fa8
+  data.tar.gz: 55802d1412f4f89513578a0a8f5fc1da994873238f6d5efd6177d44dea8195f5cead352a62e5f882a874c9823a7f49a48d98c4a3bcccb021b89535e7b17f3f53

data/.gitignore

@@ -0,0 +1,8 @@
+/.bundle/
+/.yardoc
+/_yardoc/
+/coverage/
+/doc/
+/pkg/
+/spec/reports/
+/tmp/

data/.travis.yml

@@ -0,0 +1,6 @@
+---
+language: ruby
+cache: bundler
+rvm:
+  - 2.7.0
+before_install: gem install bundler -v 2.1.0.pre.1

data/Gemfile

@@ -0,0 +1,4 @@
+source "https://rubygems.org"
+
+# Specify your gem's dependencies in hpsqrt.gemspec
+gemspec

data/LICENSE.txt

@@ -0,0 +1,21 @@
+The MIT License (MIT)
+
+Copyright (c) 2019 yoshida
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in
+all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
+THE SOFTWARE.

data/README.md

@@ -0,0 +1,49 @@
+# HpSqrt
+
+High precision square root library for Ruby.
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'hpsqrt'
+```
+
+And then execute:
+
+    $ bundle install
+
+Or install it yourself as:
+
+    $ gem install hpsqrt
+
+## Usage
+
+    require 'hpsqrt/core_ext'
+    
+    p Sqrt(2)
+    # => #<HpSqrt:0x00003ff699526258 value=(1.4142135623730951+0.0i) expr=(√2)> 
+    
+    p Sqrt(2).to_f
+    #  => 1.4142135623730951
+    
+    p (Sqrt(2) ** 2).to_f
+    #  => 2.0 
+    
+    p (Sqrt(3) * Sqrt(5) * Sqrt(15)).to_f
+    #  => 15.0 
+    
+    p (Sqrt(1i) ** 4).to_c
+    #  => (-1.0+0.0i) 
+    
+    p ((Sqrt(7) + Sqrt(11)) * (Sqrt(7) - Sqrt(11))).to_f
+    #  => -4.0 
+
+## Contributing
+
+Bug reports and pull requests are welcome on GitHub at https://github.com/yoshida-eth0/ruby-sqrt.
+
+## License
+
+The gem is available as open source under the terms of the [MIT License](https://opensource.org/licenses/MIT).

data/Rakefile

@@ -0,0 +1,10 @@
+require "bundler/gem_tasks"
+require "rake/testtask"
+
+Rake::TestTask.new(:test) do |t|
+  t.libs << "test"
+  t.libs << "lib"
+  t.test_files = FileList["test/**/*_test.rb"]
+end
+
+task :default => :test

data/bin/console

@@ -0,0 +1,14 @@
+#!/usr/bin/env ruby
+
+require "bundler/setup"
+require "hpsqrt"
+
+# You can add fixtures and/or initialization code here to make experimenting
+# with your gem easier. You can also use a different console, if you like.
+
+# (If you use this, don't forget to add pry to your Gemfile!)
+# require "pry"
+# Pry.start
+
+require "irb"
+IRB.start(__FILE__)

data/bin/setup

@@ -0,0 +1,8 @@
+#!/usr/bin/env bash
+set -euo pipefail
+IFS=$'\n\t'
+set -vx
+
+bundle install
+
+# Do any other automated setup that you need to do here

data/hpsqrt.gemspec

@@ -0,0 +1,32 @@
+lib = File.expand_path("lib", __dir__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require "hpsqrt/version"
+
+Gem::Specification.new do |spec|
+  spec.name          = "hpsqrt"
+  spec.version       = HpSqrt::VERSION
+  spec.authors       = ["Yoshida Tetsuya"]
+  spec.email         = ["yoshida.eth0@gmail.com"]
+
+  spec.summary       = %q{High precision square root library}
+  spec.description   = %q{High precision square root library}
+  spec.homepage      = "https://github.com/yoshida-eth0/ruby-sqrt"
+  spec.license       = "MIT"
+
+  spec.metadata["homepage_uri"] = spec.homepage
+  spec.metadata["source_code_uri"] = spec.homepage
+  spec.metadata["changelog_uri"] = spec.homepage
+
+  # 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('..', __FILE__)) 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", "~> 2.1.a"
+  spec.add_development_dependency "rake", "~> 12.0"
+  spec.add_development_dependency "minitest", "~> 5.0"
+end

data/lib/hpsqrt.rb

@@ -0,0 +1,242 @@
+require 'hpsqrt/inspect_mode'
+require 'hpsqrt/value'
+require 'hpsqrt/version'
+
+
+class HpSqrt < Numeric
+
+  @@inspect_mode = INSPECT_MODE::DUMP
+
+
+  attr_reader :values
+
+  def initialize(values)
+    @values = values.freeze
+    freeze
+  end
+
+  def -@
+    values = self.values.map{|v,c| [v, -c]}.to_h
+    self.class.new(values)
+  end
+
+  def +(other)
+    other = self.class.create(other)
+
+    values = @values.merge(other.values) {|v, c1, c2|
+      c1 + c2
+    }
+    values.delete_if {|v,c| c==0}
+
+    self.class.new(values)
+  end
+
+  def -(other)
+    other = self.class.create(other)
+
+    other_values = other.values.map{|v,c| [v, -c]}.to_h
+    values = @values.merge(other_values) {|v, c1, c2|
+      c1 + c2
+    }
+    values.delete_if {|v,c| v==0}
+
+    self.class.new(values)
+  end
+
+  def *(other)
+    other = self.class.create(other)
+
+    values = {}
+    @values.each {|v1, c1|
+      other.values.each {|v2, c2|
+        v = v1 * v2
+        c = c1 * c2
+
+        values[v] ||= 0
+        values[v] += c
+      }
+    }
+    values.delete_if {|v,c| c==0}
+
+    self.class.new(values)
+  end
+
+  def /(other)
+    other = self.class.create(other)
+    other_inv = Value.new(number: Rational(1, other.value))
+
+    values = {}
+    @values.each {|v, c|
+      v *= other_inv
+
+      values[v] ||= 0
+      values[v] += c
+    }
+    values.delete_if {|v,c| c==0}
+
+    self.class.new(values)
+  end
+
+  def **(other)
+    other = self.class.create(other)
+
+    if other.int?
+      result = self.class.create(1)
+      other_i = other.real.to_i
+      other_i.abs.times {|i|
+        result *= self
+      }
+      if other_i<0
+        result = Rational(1, result)
+      end
+      result
+    else
+      self.class.number(self.value ** other.value)
+    end
+  end
+
+  def coerce(other)
+    [self.class.create(other), self]
+  end
+
+  def ==(other)
+    if self.class==other
+      self.value==other.value
+    elsif Numeric===other
+      self.value==other
+    else
+      super.==(other)
+    end
+  end
+
+  def value
+    @values.map {|v, c|
+      v.number * Math.sqrt(Complex(v.sqrt)) * c
+    }.sum
+  end
+
+  def real
+    value.real
+  end
+
+  def imag
+    Complex(value).imag
+  end
+
+  def to_i
+    c = to_c
+    if c.imag.zero?
+      c.real.to_i
+    else
+      raise RangeError, "can't convert %s into Integer" % c
+    end
+  end
+
+  def to_f
+    c = to_c
+    if c.imag.zero?
+      c.real.to_f
+    else
+      raise RangeError, "can't convert %s into Float" % c
+    end
+  end
+
+  def to_c
+    value.to_c
+  end
+
+  def to_r
+    value.to_r
+  end
+
+  def expr
+    value_to_s = -> (v) {
+      if Complex===v && v.imag.zero?
+        v = v.real
+      elsif Rational===v
+        v = v.to_s.sub(/\/1$/, "")
+      end
+      v = v.to_s
+      if v !~ /^[\d\.]+$/
+        v = "(%s)" % v
+      end
+      v
+    }
+
+    result = @values.map {|v, c|
+      n = v.number * c
+      s = v.sqrt
+
+      if s!=1
+        if n==1
+          "\u221A%s" % value_to_s[s]
+        elsif 0<n.real
+          "%s\u221A%s" % [value_to_s[n], value_to_s[s]]
+        elsif n==-1
+          "(-\u221A%s)" % value_to_s[s]
+        else
+          "(%s\u221A%s)" % [value_to_s[n], value_to_s[s]]
+        end
+      else
+        value_to_s[n]
+      end
+    }
+
+    if 0<result.length
+      result.join(" + ")
+    else
+      "0"
+    end
+  end
+
+  def inspect
+    to_s
+  end
+
+  def to_s
+    case @@inspect_mode
+    when INSPECT_MODE::VALUE
+      value.to_s
+    when INSPECT_MODE::EXPR
+      expr
+    when INSPECT_MODE::DUMP
+      "#<%s:0x%016x value=(%s) expr=(%s)>" % [self.class.name, self.object_id, value, expr]
+    end
+  end
+
+  def real?
+    imag.zero?
+  end
+
+  def imag?
+    !real?
+  end
+
+  def int?
+    v = value
+    is_imag = Complex===v && !v.imag.zero?
+    !is_imag && v.real==v.real.to_i
+  end
+
+  def float?
+    v = value
+    is_imag = Complex===v && !v.imag.zero?
+    !is_imag && Float===v.real && v.real!=v.real.to_i
+  end
+
+  def self.create(v)
+    if self===v
+      v
+    else
+      number(v)
+    end
+  end
+
+  def self.number(v)
+    new({Value.new(number: v) => 1})
+  end
+
+  def self.sqrt(v)
+    new({Value.new(sqrt: v) => 1})
+  end
+end

data/lib/hpsqrt/core_ext.rb

@@ -0,0 +1,8 @@
+require 'hpsqrt'
+
+
+module Kernel
+  def Sqrt(v)
+    HpSqrt.sqrt(v)
+  end
+end

data/lib/hpsqrt/inspect_mode.rb

@@ -0,0 +1,10 @@
+
+class HpSqrt < Numeric
+
+  module INSPECT_MODE
+    VALUE = :value
+    EXPR = :expr
+    DUMP = :dump
+  end
+
+end

data/lib/hpsqrt/value.rb

@@ -0,0 +1,36 @@
+class HpSqrt < Numeric
+
+  class Value
+    attr_reader :number
+    attr_reader :sqrt
+
+    def initialize(number: 1, sqrt: 1)
+      @number = number
+      @sqrt = sqrt
+      freeze
+    end
+
+    def *(other)
+      if self.class===other
+        n = @number * other.number
+        s = 1
+        if @sqrt==other.sqrt
+          n *= @sqrt
+        else
+          s = @sqrt * other.sqrt
+        end
+
+        self.class.new(number: n, sqrt: s)
+      end
+    end
+
+    def eql?(other)
+      self.class===other && @number==other.number && @sqrt==other.sqrt
+    end
+
+    def hash
+      [@number, @sqrt].hash
+    end
+  end
+
+end

data/lib/hpsqrt/version.rb

@@ -0,0 +1,3 @@
+class HpSqrt < Numeric
+  VERSION = "1.0.0"
+end

metadata

@@ -0,0 +1,102 @@
+--- !ruby/object:Gem::Specification
+name: hpsqrt
+version: !ruby/object:Gem::Version
+  version: 1.0.0
+platform: ruby
+authors:
+- Yoshida Tetsuya
+autorequire: 
+bindir: exe
+cert_chain: []
+date: 2019-06-11 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 2.1.a
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 2.1.a
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '12.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '12.0'
+- !ruby/object:Gem::Dependency
+  name: minitest
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '5.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '5.0'
+description: High precision square root library
+email:
+- yoshida.eth0@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- ".travis.yml"
+- Gemfile
+- LICENSE.txt
+- README.md
+- Rakefile
+- bin/console
+- bin/setup
+- hpsqrt.gemspec
+- lib/hpsqrt.rb
+- lib/hpsqrt/core_ext.rb
+- lib/hpsqrt/inspect_mode.rb
+- lib/hpsqrt/value.rb
+- lib/hpsqrt/version.rb
+homepage: https://github.com/yoshida-eth0/ruby-sqrt
+licenses:
+- MIT
+metadata:
+  homepage_uri: https://github.com/yoshida-eth0/ruby-sqrt
+  source_code_uri: https://github.com/yoshida-eth0/ruby-sqrt
+  changelog_uri: https://github.com/yoshida-eth0/ruby-sqrt
+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: '0'
+requirements: []
+rubygems_version: 3.0.3
+signing_key: 
+specification_version: 4
+summary: High precision square root library
+test_files: []