Gecko [ beriskala.com ]

Name	Size	Permission	Date
bigdecimal	[ DIR ]	drwxr-xr-x	2026-05-01 04:24
cgi	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
date	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
digest	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
dl	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
drb	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
fiddle	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
io	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
irb	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
json	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
matrix	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
minitest	[ DIR ]	drwxr-xr-x	2026-05-01 04:24
net	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
openssl	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
optparse	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
psych	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
racc	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
rake	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
rbconfig	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
rdoc	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
rexml	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
rinda	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
ripper	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
rss	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
rubygems	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
shell	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
syslog	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
test	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
tk	[ DIR ]	drwxr-xr-x	2026-05-01 04:23
tkextlib	[ DIR ]	drwxr-xr-x	2026-05-01 04:23
uri	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
webrick	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
x86_64-linux	[ DIR ]	drwxr-xr-x	2026-05-01 04:24
xmlrpc	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
yaml	[ DIR ]	drwxr-xr-x	2026-04-30 23:44
English.rb	6.44 KB	-rw-r--r--	2013-02-04 06:20
abbrev.rb	3.31 KB	-rw-r--r--	2013-02-21 21:05
base64.rb	2.63 KB	-rw-r--r--	2013-05-19 07:40
benchmark.rb	17.95 KB	-rw-r--r--	2013-09-20 20:35
cgi.rb	9.84 KB	-rw-r--r--	2013-07-18 18:20
cmath.rb	8.93 KB	-rw-r--r--	2015-07-03 13:42
complex.rb	380 B	-rw-r--r--	2009-08-16 20:04
csv.rb	81.68 KB	-rw-r--r--	2014-09-14 19:55
date.rb	946 B	-rw-r--r--	2023-07-26 17:52
debug.rb	29.08 KB	-rw-r--r--	2013-12-14 18:18
delegate.rb	11.13 KB	-rw-r--r--	2014-02-02 17:16
digest.rb	2.34 KB	-rw-r--r--	2023-07-26 17:52
dl.rb	280 B	-rw-r--r--	2023-07-26 17:52
drb.rb	19 B	-rw-r--r--	2009-10-02 14:15
e2mmap.rb	3.89 KB	-rw-r--r--	2013-12-17 14:50
erb.rb	26.07 KB	-rw-r--r--	2014-08-11 19:25
expect.rb	2.14 KB	-rw-r--r--	2023-07-26 17:52
fiddle.rb	1.65 KB	-rw-r--r--	2023-07-26 17:52
fileutils.rb	47.17 KB	-rw-r--r--	2015-04-28 09:46
find.rb	2.36 KB	-rw-r--r--	2014-10-15 19:01
forwardable.rb	7.86 KB	-rw-r--r--	2013-05-24 02:16
getoptlong.rb	15.38 KB	-rw-r--r--	2013-05-19 19:20
gserver.rb	8.86 KB	-rw-r--r--	2014-07-06 18:27
ipaddr.rb	17.05 KB	-rw-r--r--	2015-10-29 18:06
irb.rb	20.03 KB	-rw-r--r--	2016-03-25 12:52
json.rb	1.74 KB	-rw-r--r--	2023-07-26 17:52
kconv.rb	5.74 KB	-rw-r--r--	2023-07-26 17:52
logger.rb	22.37 KB	-rw-r--r--	2016-02-25 14:34
mathn.rb	6.52 KB	-rw-r--r--	2011-08-27 02:52
matrix.rb	45.43 KB	-rw-r--r--	2014-07-11 18:38
mkmf.rb	80.5 KB	-rw-r--r--	2023-07-26 17:49
monitor.rb	6.93 KB	-rw-r--r--	2012-11-16 20:25
multi-tk.rb	98.4 KB	-rw-r--r--	2023-07-26 17:52
mutex_m.rb	2 KB	-rw-r--r--	2013-02-20 06:21
observer.rb	5.81 KB	-rw-r--r--	2013-11-21 08:14
open-uri.rb	24.27 KB	-rw-r--r--	2014-02-22 12:24
open3.rb	20.37 KB	-rw-r--r--	2013-12-01 06:43
openssl.rb	528 B	-rw-r--r--	2023-07-26 17:52
optparse.rb	50.8 KB	-rw-r--r--	2013-11-28 11:04
ostruct.rb	7.7 KB	-rw-r--r--	2013-10-23 18:44
pathname.rb	15.3 KB	-rw-r--r--	2023-07-26 17:52
pp.rb	14.1 KB	-rw-r--r--	2013-10-22 12:59
prettyprint.rb	16.33 KB	-rw-r--r--	2013-05-19 07:40
prime.rb	13.2 KB	-rw-r--r--	2013-07-15 08:51
profile.rb	205 B	-rw-r--r--	2009-10-02 14:15
profiler.rb	4.51 KB	-rw-r--r--	2013-05-20 03:34
pstore.rb	14.85 KB	-rw-r--r--	2012-11-11 07:53
psych.rb	14.88 KB	-rw-r--r--	2023-07-26 17:52
rake.rb	2.12 KB	-rw-r--r--	2013-10-12 01:05
rational.rb	308 B	-rw-r--r--	2009-09-24 04:12
rdoc.rb	4.92 KB	-rw-r--r--	2013-12-22 13:55
remote-tk.rb	14.11 KB	-rw-r--r--	2023-07-26 17:52
resolv-replace.rb	1.73 KB	-rw-r--r--	2013-03-11 17:17
resolv.rb	72.31 KB	-rw-r--r--	2015-05-13 10:03
ripper.rb	2.53 KB	-rw-r--r--	2023-07-26 17:52
rss.rb	2.84 KB	-rw-r--r--	2011-05-11 14:52
rubygems.rb	30.97 KB	-rw-r--r--	2023-07-26 17:49
scanf.rb	23.54 KB	-rw-r--r--	2013-12-14 06:25
securerandom.rb	8.49 KB	-rw-r--r--	2013-11-07 21:34
set.rb	18.7 KB	-rw-r--r--	2013-11-23 03:20
shell.rb	10.3 KB	-rw-r--r--	2013-05-19 07:40
shellwords.rb	5.94 KB	-rw-r--r--	2013-12-14 21:56
singleton.rb	4.02 KB	-rw-r--r--	2011-05-18 18:39
socket.rb	25.6 KB	-rw-r--r--	2023-07-26 17:52
sync.rb	7.25 KB	-rw-r--r--	2013-05-19 07:40
tcltk.rb	10.06 KB	-rw-r--r--	2023-07-26 17:52
tempfile.rb	11.4 KB	-rw-r--r--	2013-11-21 12:58
thwait.rb	3.38 KB	-rw-r--r--	2013-11-07 20:32
time.rb	21.32 KB	-rw-r--r--	2015-08-17 12:11
timeout.rb	3.7 KB	-rw-r--r--	2015-08-17 12:25
tk.rb	163.72 KB	-rw-r--r--	2023-07-26 17:52
tkafter.rb	57 B	-rw-r--r--	2023-07-26 17:52
tkbgerror.rb	63 B	-rw-r--r--	2023-07-26 17:52
tkcanvas.rb	60 B	-rw-r--r--	2023-07-26 17:52
tkclass.rb	950 B	-rw-r--r--	2023-07-26 17:52
tkconsole.rb	63 B	-rw-r--r--	2023-07-26 17:52
tkdialog.rb	60 B	-rw-r--r--	2023-07-26 17:52
tkentry.rb	57 B	-rw-r--r--	2023-07-26 17:52
tkfont.rb	54 B	-rw-r--r--	2023-07-26 17:52
tkmacpkg.rb	60 B	-rw-r--r--	2023-07-26 17:52
tkmenubar.rb	63 B	-rw-r--r--	2023-07-26 17:52
tkmngfocus.rb	66 B	-rw-r--r--	2023-07-26 17:52
tkpalette.rb	63 B	-rw-r--r--	2023-07-26 17:52
tkscrollbox.rb	69 B	-rw-r--r--	2023-07-26 17:52
tktext.rb	54 B	-rw-r--r--	2023-07-26 17:52
tkvirtevent.rb	69 B	-rw-r--r--	2023-07-26 17:52
tkwinpkg.rb	60 B	-rw-r--r--	2023-07-26 17:52
tmpdir.rb	4.15 KB	-rw-r--r--	2013-05-19 07:40
tracer.rb	6.4 KB	-rw-r--r--	2013-07-18 18:20
tsort.rb	14.14 KB	-rw-r--r--	2013-10-17 19:29
ubygems.rb	268 B	-rw-r--r--	2009-10-02 14:15
un.rb	8.87 KB	-rw-r--r--	2013-07-05 18:13
uri.rb	3.07 KB	-rw-r--r--	2011-05-14 00:33
weakref.rb	3.23 KB	-rw-r--r--	2013-11-10 02:12
webrick.rb	6.69 KB	-rw-r--r--	2013-10-06 03:09
xmlrpc.rb	8.49 KB	-rw-r--r--	2013-12-12 06:39
yaml.rb	2.31 KB	-rw-r--r--	2013-08-12 08:19

Rename

#--
# $Release Version: 0.5 $
# $Revision: 1.1.1.1.4.1 $

##
# = mathn
#
# mathn is a library for changing the way Ruby does math.  If you need
# more precise rounding with multiple division or exponentiation
# operations, then mathn is the right tool.
#
# Without mathn:
#
#   3 / 2 => 1 # Integer
#
# With mathn:
#
#   3 / 2 => 3/2 # Rational
#
# mathn features late rounding and lacks truncation of intermediate results:
#
# Without mathn:
#
#   20 / 9 * 3 * 14 / 7 * 3 / 2 # => 18
#
# With mathn:
#
#   20 / 9 * 3 * 14 / 7 * 3 / 2 # => 20
#
#
# When you require 'mathn', the libraries for Prime, CMath, Matrix and Vector
# are also loaded.
#
# == Copyright
#
# Author: Keiju ISHITSUKA (SHL Japan Inc.)
#--
# class Numeric follows to make this documentation findable in a reasonable
# location

class Numeric; end

require "cmath.rb"
require "matrix.rb"
require "prime.rb"

require "mathn/rational"
require "mathn/complex"

unless defined?(Math.exp!)
  Object.instance_eval{remove_const :Math}
  Math = CMath # :nodoc:
end

##
# When mathn is required, Fixnum's division and exponentiation are enhanced to
# return more precise values from mathematical expressions.
#
#   2/3*3  # => 0
#   require 'mathn'
#   2/3*3  # => 2

class Fixnum
  remove_method :/

##
  # +/+ defines the Rational division for Fixnum.
  #
  #   1/3  # => (1/3)

alias / quo

alias power! ** unless method_defined? :power!

##
  # Exponentiate by +other+

def ** (other)
    if self < 0 && other.round != other
      Complex(self, 0.0) ** other
    else
      power!(other)
    end
  end

end

##
# When mathn is required Bignum's division and exponentiation are enhanced to
# return more precise values from mathematical expressions.

class Bignum
  remove_method :/

##
  # +/+ defines the Rational division for Bignum.
  #
  #   (2**72) / ((2**70) * 3)  # => 4/3

alias / quo

alias power! ** unless method_defined? :power!

##
  # Exponentiate by +other+

def ** (other)
    if self < 0 && other.round != other
      Complex(self, 0.0) ** other
    else
      power!(other)
    end
  end

end

##
# When mathn is required Rational is changed to simplify the use of Rational
# operations.
#
# Normal behaviour:
#
#   Rational.new!(1,3) ** 2 # => Rational(1, 9)
#   (1 / 3) ** 2            # => 0
#
# require 'mathn' behaviour:
#
#   (1 / 3) ** 2            # => 1/9

class Rational
  remove_method :**

##
  # Exponentiate by +other+
  #
  #   (1/3) ** 2 # => 1/9

def ** (other)
    if other.kind_of?(Rational)
      other2 = other
      if self < 0
        return Complex(self, 0.0) ** other
      elsif other == 0
        return Rational(1,1)
      elsif self == 0
        return Rational(0,1)
      elsif self == 1
        return Rational(1,1)
      end

npd = numerator.prime_division
      dpd = denominator.prime_division
      if other < 0
        other = -other
        npd, dpd = dpd, npd
      end

for elm in npd
        elm[1] = elm[1] * other
        if !elm[1].kind_of?(Integer) and elm[1].denominator != 1
          return Float(self) ** other2
        end
        elm[1] = elm[1].to_i
      end

for elm in dpd
        elm[1] = elm[1] * other
        if !elm[1].kind_of?(Integer) and elm[1].denominator != 1
          return Float(self) ** other2
        end
        elm[1] = elm[1].to_i
      end

num = Integer.from_prime_division(npd)
      den = Integer.from_prime_division(dpd)

Rational(num,den)

elsif other.kind_of?(Integer)
      if other > 0
        num = numerator ** other
        den = denominator ** other
      elsif other < 0
        num = denominator ** -other
        den = numerator ** -other
      elsif other == 0
        num = 1
        den = 1
      end
      Rational(num, den)
    elsif other.kind_of?(Float)
      Float(self) ** other
    else
      x , y = other.coerce(self)
      x ** y
    end
  end
end

##
# When mathn is required, the Math module changes as follows:
#
# Standard Math module behaviour:
#   Math.sqrt(4/9)     # => 0.0
#   Math.sqrt(4.0/9.0) # => 0.666666666666667
#   Math.sqrt(- 4/9)   # => Errno::EDOM: Numerical argument out of domain - sqrt
#
# After require 'mathn', this is changed to:
#
#   require 'mathn'
#   Math.sqrt(4/9)      # => 2/3
#   Math.sqrt(4.0/9.0)  # => 0.666666666666667
#   Math.sqrt(- 4/9)    # => Complex(0, 2/3)

module Math
  remove_method(:sqrt)

##
  # Computes the square root of +a+.  It makes use of Complex and
  # Rational to have no rounding errors if possible.
  #
  #   Math.sqrt(4/9)      # => 2/3
  #   Math.sqrt(- 4/9)    # => Complex(0, 2/3)
  #   Math.sqrt(4.0/9.0)  # => 0.666666666666667

def sqrt(a)
    if a.kind_of?(Complex)
      abs = sqrt(a.real*a.real + a.imag*a.imag)
#      if not abs.kind_of?(Rational)
#        return a**Rational(1,2)
#      end
      x = sqrt((a.real + abs)/Rational(2))
      y = sqrt((-a.real + abs)/Rational(2))
#      if !(x.kind_of?(Rational) and y.kind_of?(Rational))
#        return a**Rational(1,2)
#      end
      if a.imag >= 0
        Complex(x, y)
      else
        Complex(x, -y)
      end
    elsif a.respond_to?(:nan?) and a.nan?
      a
    elsif a >= 0
      rsqrt(a)
    else
      Complex(0,rsqrt(-a))
    end
  end

##
  # Compute square root of a non negative number. This method is
  # internally used by +Math.sqrt+.

def rsqrt(a)
    if a.kind_of?(Float)
      sqrt!(a)
    elsif a.kind_of?(Rational)
      rsqrt(a.numerator)/rsqrt(a.denominator)
    else
      src = a
      max = 2 ** 32
      byte_a = [src & 0xffffffff]
      # ruby's bug
      while (src >= max) and (src >>= 32)
        byte_a.unshift src & 0xffffffff
      end

answer = 0
      main = 0
      side = 0
      for elm in byte_a
        main = (main << 32) + elm
        side <<= 16
        if answer != 0
          if main * 4  < side * side
            applo = main.div(side)
          else
            applo = ((sqrt!(side * side + 4 * main) - side)/2.0).to_i + 1
          end
        else
          applo = sqrt!(main).to_i + 1
        end

while (x = (side + applo) * applo) > main
          applo -= 1
        end
        main -= x
        answer = (answer << 16) + applo
        side += applo * 2
      end
      if main == 0
        answer
      else
        sqrt!(a)
      end
    end
  end

class << self
    remove_method(:sqrt)
  end
  module_function :sqrt
  module_function :rsqrt
end

##
# When mathn is required, Float is changed to handle Complex numbers.

class Float
  alias power! **

##
  # Exponentiate by +other+

def ** (other)
    if self < 0 && other.round != other
      Complex(self, 0.0) ** other
    else
      power!(other)
    end
  end

end