ruby-dev

原です。

Tadayoshi Funaba さんは書きました:
>> $ ruby -I../nurat-0.0.4/lib:../nurat-0.0.3/ext bench.rb nurat-0.0.4 >>
>> benchmark.log
> 
> これだと 0.0.3 をつかってしまいそうですが。ま、あまり変りないかもしれ
> ませんが。

あ、失敗。

> 原さんのはかけ算が圧倒的に速いですね。ただ、加減算などになるとあまり違
> いがなくなってきて、特に数が大きくなると、nurat のほうが速くなったりす
> るようです。いいところをうまく取りこめるといいですが。

> --- hara-bench.rb	Tue Feb 19 06:49:34 2008
> +++ hara-bench-n.rb	Tue Feb 19 06:59:31 2008

> --- hara-bench.rb	Tue Feb 19 06:49:34 2008
> +++ hara-bench-n2.rb	Tue Feb 19 07:00:08 2008

確かに加算の rational-1.19 はいまいちだなあ。

> --- hara-bench.rb	Tue Feb 19 06:49:34 2008
> +++ hara-bench-n3.rb	Tue Feb 19 07:01:04 2008
> @@ -22,13 +22,13 @@
>  end
>  
>  def prod1
> -  n = 10**4
> +  n = 10**5
>    a = sa('a', 1, 100, n).zip(sa('b', 1, 100, n)).map{|x, y| Rational(x, y)}
>  
>    bm("product-small", $libname) do
> -    s = Rational(1)
> +    s = Rational(2)
>      a.each do |x|
> -      s *= x
> +      s **= x
>      end
>    end
>  end

これは、掛け算ではなく **(aRational) の問題みたい。nurat では、

  N = Float::MAX.to_i * 2
  p Rational(N + 2, N + 1).to_f

とすると NaN ですが、rational-1.19 では複雑なことをしていて、1.0 にな
ります。しかし、工夫の余地があるということだなあ。

-----^ bm.rb
#!/usr/bin/ruby
require "rational.rb"

$libname = ARGV.shift || "(unknown)"

def sa(name, from, to, size)
  fname = sprintf("INT%s_%d_%d_%d", name, from, to, size)
  unless File.exist?(fname)
    a = (1..size).map{Rational(rand(to - from + 1) + from)}
    open(fname, "w") {|f| f << a.join("\n")}
  end
  File.read(fname).split.map{|x| x.to_i}
end

def bm(name, routine = nil, n = 1)
  start = Time.now
  n.times do
    yield
  end
  stop = Time.now
  printf("= %-20s %-20s %f sec.\n",  name, "(#{$libname})", stop - start)#-
end

def rats(a, b, n)
  sa(*(a+[n])).zip(sa(*(b+[n]))).map{|x, y| Rational(x, y)}
end

def plus0(title, x, y, n)
  a = rats(['a', *x], ['b', *y], n)
  bm(title, $libname) do
    s = Rational(0)
    a.each do |x|
      s += x
    end
  end
end

def mul0(title, x, y, n)
  a = rats(['a', *x], ['b', *y], n)
  bm(title, $libname) do
    s = Rational(1)
    a.each do |x|
      s *= x
    end
  end
end

def rat_pow0(title, x, y, n)
  a = rats(['a', *x], ['b', *y], n)
  bm(title, $libname) do
    s = Rational(2)
    a.each do |x|
      s **= x
    end
  end
end


if $0 == __FILE__
  plus0("addtion-small",  [-100, 100], [1, 100], 10**3)
  mul0("product-small", [1, 100], [1, 100], 10**3)
  plus0("addition-big", [10**30, 10**40], [10**30, 10**40], 50)
  mul0("product-big", [10**30, 10**40], [10**30, 10**40], 50)
  rat_pow0("exponential-rat-big", [1, 100], [1, 100], 10**5)
end
-----$ bm.rb