[ruby-dev:4788] Re: compare to NaN (Re: Re: call for finite/isnan testers)

ごとけんです

In message "[ruby-dev:4780] Re: compare to NaN (Re:  Re: call for finite/isnan testers)"
    on 99/02/06, Yukihiro Matsumoto <[e-mail]> writes:

>まつもと ゆきひろです
>
>In message "[ruby-dev:4760] compare to NaN (Re:  Re: call for finite/isnan testers)"
>    on 99/02/06, EGUCHI Osamu <[e-mail]> writes:
>
>|えぐち＠エスアンドイー です。
>
>| * Float::INFINITY の追加
>| * Float::NaN の追加
>| * NaN との大小判別例外 NaNCompareError クラスの追加
>| * FLoat#infinite? の追加(Inf:正数, -Inf:負数、otherwise false)
>| * FLoat#nan? の追加(また x == Float::NaN は常に偽)
>| * FLoat#finte? の追加(Inf, -Inf or NaN:false, otherwise true)
>| * FreeBSD 環境での浮動小数点関連のトラップの抑止(main.c)
>
>面白いですね．採用に関して前向きなのですが，気になる点がいく
>つかあります．
>
>仕様について
>
>  -INFINITYを示す定数は不要か(-Float::INFINITYで十分か)

あった方がよいように思いますが名前が難しいですね。
NEG_INFINITY かな。

>  infinite?, finite?の両方が必要か

IEEE 754 的には、これらは単に逆というわけではありません。

		1.0 +-Inf NaN
   finite?       T    F    F
   infinite?     F    T    F
   nan?          F    F    T

確かにどれか2つがあれば他の1つが作れますが。
def infinite; !x.finite? && !x.nan?; end; みたく。

>  Float::infinite?は真の時に正数，負数を返す仕様は妥当か
>
>  このFloat::INFINITYやFloat::NaNのシステムはごとけんさんの
>  Float抽象クラス化と衝突しないか

確かに実装上は double だから Float なんですけど、
実際は Numeric::INFINITY というのがあって、それを
比較のときは適当なのとくらべる方がすっきりするような
気もします。Complex だと実部か虚部かどっちかが
無限大だったら無限大というように、無限大自体は
数値ではないけど普通数値と言うと無限大を想定するのは
わりと自然です。

またこの場合、Integer::INFINITY, Float::INFINITY,
SingleFloat::INFINITY, DoubleFloat::INFINITY 
というのがあった方がよいかも知れません。NaN も同様です。

もう一つの考え方は、NaN, +-Inf それぞれ一つで、
一括して Numeric の定数だとする考え方です。

それから、いずれにしても定数に関して演算の規則を
規定しなければなりませんね。 IEEE 754 のアヤシイ訳を
読んで検討してみます。

>実装について
>
>  main.cでFreeBSD対応しているが，Init_Numeric()あたりでやる
>  べきでは？ いや，そうするとRubyをlinkしたアプリケーション
>  すべてに強制することになるか…．

ぼくもやるなら Init_Numeric() の方がよいと思います。

>  isinf(),isnan()の移植性をどうするか．ごとけんさんの調査待
>  ちか．

isnan() 相当はどこにでもあって、finite() もしくは isinf() 
ならどこにでもある模様です。 
BeOS, NeXT, それから 64 bit の Alpha とかはよく分かりませんが。

>  1.0/0.0でInfを0.0/0.0でNaNを作っているが移植性は？

木村さんの(ぢゃなくてもいいけど VC の)結果を待ってます。

>  static double infinity, nan, zero = 0.0;のstaticはいらない
>  んじゃない？

ですよね。

>  flo_infinite_p()で+1, -1を返してるが，これはINT2FIX()でく
>  くるべきだと思う．

# なるほどお。Float の抽象化難航中 (汗;

>数値に強い方の登場をお待ちしています．

>数値に強いってのと数学に強いってのはまた違うんだよな．私はどっ
>ちにも弱いけど．

僕だってどっちも強いと言うほどは知らないんだけど、
まつもとさんが凄いとこは、知らなくてもいろいろ実装
できちゃうとこですね。

-- gotoken