Gecko [ www.corporacionhdm.com.pe ]

Name	Size	Permission	Date
bigdecimal	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
cgi	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
date	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
digest	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
dl	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
drb	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
io	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
irb	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
net	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
openssl	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
optparse	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
racc	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
rdoc	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
rexml	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
rinda	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
rss	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
runit	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
shell	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
soap	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
test	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
uri	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
webrick	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
wsdl	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
x86_64-linux	[ DIR ]	drwxr-xr-x	2024-03-03 22:50
xmlrpc	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
xsd	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
yaml	[ DIR ]	drwxr-xr-x	2024-03-03 22:48
English.rb	5.6 KB	-rw-r--r--	2007-02-12 23:01
Env.rb	274 B	-rw-r--r--	2007-02-12 23:01
abbrev.rb	2.5 KB	-rw-r--r--	2007-02-12 23:01
base64.rb	3.37 KB	-rw-r--r--	2007-07-16 15:47
benchmark.rb	17.73 KB	-rw-r--r--	2008-02-10 15:24
cgi-lib.rb	6.89 KB	-rw-r--r--	2007-02-12 23:01
cgi.rb	73.74 KB	-rw-r--r--	2009-12-14 02:40
complex.rb	12.84 KB	-rw-r--r--	2009-08-03 05:59
csv.rb	24.46 KB	-rw-r--r--	2007-02-12 23:01
date.rb	53.02 KB	-rw-r--r--	2010-06-08 04:45
date2.rb	128 B	-rw-r--r--	2007-02-12 23:01
debug.rb	20.61 KB	-rw-r--r--	2007-02-12 23:01
delegate.rb	8.81 KB	-rw-r--r--	2009-07-16 00:35
digest.rb	1.12 KB	-rw-r--r--	2007-07-28 00:40
drb.rb	19 B	-rw-r--r--	2007-02-12 23:01
e2mmap.rb	4.04 KB	-rw-r--r--	2007-02-12 23:01
erb.rb	21.38 KB	-rw-r--r--	2009-02-23 17:44
eregex.rb	487 B	-rw-r--r--	2007-02-12 23:01
expect.rb	633 B	-rw-r--r--	2008-02-18 01:17
fileutils.rb	42.23 KB	-rw-r--r--	2011-05-20 22:29
finalize.rb	5.38 KB	-rw-r--r--	2007-02-12 23:01
find.rb	1.84 KB	-rw-r--r--	2007-02-12 23:01
forwardable.rb	6.16 KB	-rw-r--r--	2008-06-06 08:05
ftools.rb	6.17 KB	-rw-r--r--	2007-02-12 23:01
generator.rb	8.1 KB	-rw-r--r--	2008-04-10 10:52
getoptlong.rb	14.88 KB	-rw-r--r--	2007-02-12 23:01
getopts.rb	2.25 KB	-rw-r--r--	2007-02-12 23:01
gserver.rb	6.43 KB	-rw-r--r--	2007-02-12 23:01
importenv.rb	590 B	-rw-r--r--	2007-02-12 23:01
ipaddr.rb	21.96 KB	-rw-r--r--	2008-07-12 15:08
irb.rb	7.43 KB	-rw-r--r--	2009-08-09 08:44
jcode.rb	4.3 KB	-rw-r--r--	2007-02-12 23:01
kconv.rb	8.12 KB	-rw-r--r--	2007-02-12 23:01
logger.rb	17.59 KB	-rw-r--r--	2011-05-30 02:08
mailread.rb	1.28 KB	-rw-r--r--	2007-02-12 23:01
mathn.rb	5.42 KB	-rw-r--r--	2007-02-12 23:01
matrix.rb	27.21 KB	-rw-r--r--	2009-01-24 15:02
md5.rb	411 B	-rw-r--r--	2007-03-06 10:09
mkmf.rb	50.65 KB	-rw-r--r--	2010-12-04 06:34
monitor.rb	7.93 KB	-rw-r--r--	2009-11-25 07:45
mutex_m.rb	2.07 KB	-rw-r--r--	2007-02-12 23:01
observer.rb	5.15 KB	-rw-r--r--	2007-02-12 23:01
open-uri.rb	20.49 KB	-rw-r--r--	2007-10-08 11:16
open3.rb	2.1 KB	-rw-r--r--	2009-12-14 04:28
openssl.rb	575 B	-rw-r--r--	2010-11-22 07:21
optparse.rb	47.12 KB	-rw-r--r--	2009-02-20 11:43
ostruct.rb	3.35 KB	-rw-r--r--	2009-05-26 12:06
parsearg.rb	1.55 KB	-rw-r--r--	2007-02-12 23:01
parsedate.rb	1.33 KB	-rw-r--r--	2007-02-12 23:01
pathname.rb	29.39 KB	-rw-r--r--	2010-11-23 08:21
ping.rb	1.48 KB	-rw-r--r--	2007-02-12 23:01
pp.rb	15.97 KB	-rw-r--r--	2007-06-07 10:06
prettyprint.rb	18.33 KB	-rw-r--r--	2007-02-16 19:53
profile.rb	90 B	-rw-r--r--	2007-02-12 23:01
profiler.rb	1.59 KB	-rw-r--r--	2007-02-12 23:01
pstore.rb	11.15 KB	-rw-r--r--	2008-06-06 08:05
rational.rb	12.05 KB	-rw-r--r--	2010-06-08 05:02
readbytes.rb	835 B	-rw-r--r--	2007-02-12 23:01
resolv-replace.rb	1.55 KB	-rw-r--r--	2008-08-08 01:58
resolv.rb	56.83 KB	-rw-r--r--	2010-12-23 03:22
rss.rb	504 B	-rw-r--r--	2007-10-21 12:19
rubyunit.rb	180 B	-rw-r--r--	2007-02-12 23:01
scanf.rb	20.63 KB	-rw-r--r--	2009-01-20 03:23
securerandom.rb	4.27 KB	-rw-r--r--	2011-12-10 12:17
set.rb	27.08 KB	-rw-r--r--	2008-06-09 09:20
sha1.rb	418 B	-rw-r--r--	2007-03-06 10:09
shell.rb	4.66 KB	-rw-r--r--	2007-02-12 23:01
shellwords.rb	3.99 KB	-rw-r--r--	2008-03-21 12:16
singleton.rb	8.08 KB	-rw-r--r--	2008-07-03 11:14
sync.rb	6.09 KB	-rw-r--r--	2009-02-19 16:41
tempfile.rb	4.86 KB	-rw-r--r--	2008-04-21 09:43
thread.rb	104 B	-rw-r--r--	2007-02-13 19:39
thwait.rb	4.32 KB	-rw-r--r--	2010-06-08 07:08
time.rb	31.58 KB	-rw-r--r--	2008-06-06 08:05
timeout.rb	3 KB	-rw-r--r--	2010-06-08 06:24
tmpdir.rb	3.69 KB	-rw-r--r--	2009-01-26 02:12
tracer.rb	2.73 KB	-rw-r--r--	2007-02-12 23:01
tsort.rb	7.99 KB	-rw-r--r--	2007-02-12 23:01
un.rb	4.54 KB	-rw-r--r--	2007-02-12 23:01
uri.rb	710 B	-rw-r--r--	2008-04-15 09:41
weakref.rb	2.68 KB	-rw-r--r--	2007-02-12 23:01
webrick.rb	811 B	-rw-r--r--	2007-02-12 23:01
yaml.rb	12.36 KB	-rw-r--r--	2008-04-19 11:45

Rename

#
#   mathn.rb - 
#   	$Release Version: 0.5 $
#   	$Revision: 1.1.1.1.4.1 $
#   	$Date: 1998/01/16 12:36:05 $
#   	by Keiju ISHITSUKA(SHL Japan Inc.)
#
# --
#
#   
#

require "complex.rb"
require "rational.rb"
require "matrix.rb"

class Integer

def gcd2(int)
    a = self.abs
    b = int.abs
    a, b = b, a if a < b
    
    pd_a = a.prime_division
    pd_b = b.prime_division
    
    gcd = 1
    for pair in pd_a
      as = pd_b.assoc(pair[0])
      if as
	gcd *= as[0] ** [as[1], pair[1]].min
      end
    end
    return gcd
  end
  
  def Integer.from_prime_division(pd)
    value = 1
    for prime, index in pd
      value *= prime**index
    end
    value
  end
  
  def prime_division
    raise ZeroDivisionError if self == 0
    ps = Prime.new
    value = self
    pv = []
    for prime in ps
      count = 0
      while (value1, mod = value.divmod(prime)
	     mod) == 0
	value = value1
	count += 1
      end
      if count != 0
	pv.push [prime, count]
      end
      break if prime * prime  >= value
    end
    if value > 1
      pv.push [value, 1]
    end
    return pv
  end
end
  
class Prime
  include Enumerable

def initialize
    @seed = 1
    @primes = []
    @counts = []
  end
  
  def succ
    i = -1
    size = @primes.size
    while i < size
      if i == -1
	@seed += 1
	i += 1
      else
	while @seed > @counts[i]
	  @counts[i] += @primes[i]
	end
	if @seed != @counts[i]
	  i += 1
	else
	  i = -1
	end
      end
    end
    @primes.push @seed
    @counts.push @seed + @seed
    return @seed
  end
  alias next succ

def each
    loop do
      yield succ
    end
  end
end

class Fixnum
  alias / quo
end

class Bignum
  alias / quo
end

class Rational
  Unify = true

def inspect
    format "%s/%s", numerator.inspect, denominator.inspect
  end

alias power! **

def ** (other)
    if other.kind_of?(Rational)
      other2 = other
      if self < 0
	return Complex.new!(self, 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.new!(num, den)
    elsif other.kind_of?(Float)
      Float(self) ** other
    else
      x , y = other.coerce(self)
      x ** y
    end
  end

def power2(other)
    if other.kind_of?(Rational)
      if self < 0
	return Complex(self, 0) ** other
      elsif other == 0
	return Rational(1,1)
      elsif self == 0
	return Rational(0,1)
      elsif self == 1
	return Rational(1,1)
      end
      
      dem = nil
      x = self.denominator.to_f.to_i
      neard = self.denominator.to_f ** (1.0/other.denominator.to_f)
      loop do
	if (neard**other.denominator == self.denominator)
	  dem = neaed
	  break
	end
      end
      nearn = self.numerator.to_f ** (1.0/other.denominator.to_f)
      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.new!(num, den)
    elsif other.kind_of?(Float)
      Float(self) ** other
    else
      x , y = other.coerce(self)
      x ** y
    end
  end
end

module Math
  def sqrt(a)
    if a.kind_of?(Complex)
      abs = sqrt(a.real*a.real + a.image*a.image)
#      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.image >= 0 
	Complex(x, y)
      else
	Complex(x, -y)
      end
    elsif a >= 0
      rsqrt(a)
    else
      Complex(0,rsqrt(-a))
    end
  end
  
  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

module_function :sqrt
  module_function :rsqrt
end

class Complex
  Unify = true
end