ふなばです。
rational/rationalize の話をしたので、原さんの rational 1.00 の
Rational(aFloat) と比べてみました。
原さんの Rational 1.00 ではこうなりますが、
Rational(0.3) => 2702159776422297/9007199254740992
いくつかの Common Lisp や Haskell 処理系などで、浮動小数点数の書式が同
じ場合、rational や toRational の結果はちゃんと一致しているようです。
というか、そうじゃないものはみつけられませんでした。
(rational 0.3f0) => 5033165/16777216
(rational 0.3d0) => 5404319552844595/18014398509481984
toRational (0.3::Float) => 5033165 % 16777216
toRational (0.3::Double) => 5404319552844595 % 18014398509481984
(inexact->exact 0.3) => 5404319552844595/18014398509481984
0.3 asFraction => (5404319552844595/18014398509481984)
など。
とりあえず、原さんの Rational 1.00 も rb_rat_flo_to_rat() をつぎのよう
にすると、同様になりそうです。
static VALUE
rb_rat_flo_to_rat(klass, x)
double x;
{
double f;
int e;
f = frexp(x, &e);
#if 1 /* dbl */
f = ldexp(f, DBL_MANT_DIG);
e -= DBL_MANT_DIG;
#else /* flt */
f = (float)ldexp(f, FLT_MANT_DIG);
e -= FLT_MANT_DIG;
#endif
return RT_S_NEW2(klass, RT_MUL(rb_dbl2big(f),
RT_POW(rb_int2big(2), INT2FIX(e))), INT2FIX(1));
}
# コード上の数字 2 は、FLT_RADIX にすべきかも。
CL の仕様書によれば「浮動小数点数は (数学的には) s*f*b^(e-p) の形式の
有理数である」らしいので。これに従えば、問題なく中身が取り出せているな
らば、結果も定まるように思います。
今のところ、どうすべきか、ということはないのですが、なぜ違うのかはちょっ
と疑問です。
rationalize のほうは、モノにもよりますが、rational の10倍前後遅いみた
いですね。なんか難しいことをしているのでしょうか。これを to_r にするの
は少し厳しいのかも。とりあえず、僕がみたかぎりでは、rational のほうが
普通っぽいようです。