[ruby-dev:5200] Bigfloat [compare to NaN]

けいじゅ＠日本ラショナルソフトウェアです.

In [ruby-dev :5167 ] the message: "[ruby-dev:5167] Re: compare to NaN
(Re: Re: call for finite/isnan testers) ", on Feb/16 22:55(JST) EGUCHI
Osamu writes:

>えぐち＠エスアンドイー です。

>keiju> >仮数部固定長(１０進３２桁)で Bigfloat を実装しました。
>keiju> >Ruby のコードです。  添付します。
>keiju> おおっ!! すばらしい.
>
>９８％の完成度の物を、８０％まで完成度を落したようで
>心苦しいのですが ^^;;;

全然そんなことないですよ(^^;;;

>keiju> >BASE はともかく、PREC を可変出来るように設計したつもりです。
>keiju> これがないとね(^^;;;
>
>『BASE はともかく』が、後を引く事になってしまいました ^^;;;;

(^^;;

>keiju> >あと、無限／非数／正規化への配慮もしました。
>keiju> さすが(^^;;
>まだ、ゼロの符合に関する配慮が欠けています。

+1 0 0 と -1 0 0 ってことですかね?

>keiju> > def inf?(n)
>keiju> >     n == n * 1.1
>keiju> > end
>keiju> >が良い様です。(1.01 でも 1.000000000001 でもいいけど、、)
>keiju> うーん. 微妙な...
>
>これ嘘です。
>
>	def inf?(n)
>	   n == n * 1.1 + 1.0
>	end
>
>でないと 0.0 が Inf と判断されちゃいます ;p

それはそうですね(^^;;;

>keiju> 些細なことですが, 幾つか気がついたことを.
>添削ありがとうございます。:-)

いえ. そういうつもりでは(^^;;;

>keiju> >require "e2mmap"
>keiju> これ^無理して使って下さらなくても良いですよ(^^;;;
>keiju> 
>keiju> >class Bigfloat < Numeric
>keiju> >  extend Exception2MessageMapper
>keiju> >  def_e2message(TypeError, "wrong argument type %s (expected %s)")
>keiju> 
>keiju> このあたり^ですね.
>
>いやぁ、カッコ良くて気に入ってるんですけど :)

(^^;;;

>keiju> >  attr :sign # plus(and zero):+1, minus:-1, +Inf:+Inf, -Inf:-Inf, NaN:NaN
>keiju> >  attr :prov # PREC'digits prov.
>keiju> >  attr :exp  # base BASE of exp.
>keiju> 
>keiju> 仮数とか指数ってのはどっちかというと実装上の都合ですので丸見えにするこ
>keiju> ともないような.
>
>隠蔽しても良いのですが、 def + とかで、楽しちゃったもので、、

なるなる. protectedにするか,

@sign, @prov, @exp = num.instace_eval{[@sign, @prov, @exp]}

を用いるというてもありますね.

>keiju> @provは, 常にPRECだけの桁数を持つBignumなんですね? 確かにその方がよい
>keiju> に違いない(^^;;;
>
>有効数字だけ仮数を充填しておく事で正規化してます。

ですね. 気がつかなかったなあ...

>あれ、無限や非数を正規化してるじゃん！
>また、指数を増す方の正規化に、丸めのへ配慮が欠けてました。

>keiju> >      @exp = exp10 - FLOAT_PREC
>keiju> >      @prov = Integer(num / (BASE ** @exp))
>keiju> 
>keiju> あと, これだと丸め誤差が気にならないでもないんですが... 確かごとけんさ
>keiju> んが固定小数点数を作られているみたいなので, BigfloatはFloatとの相性を
>keiju> 考えてBASEを2などの数にした方がよいかも. そうすれば, Math.frexpを使っ
>keiju> て誤差なしで変換できますし.
>
>２進数は考えましたが、 to_s があまりにややこしかったので挫折しました。;_;

そうです? 変わらないと思うんですが...

>私的には、２進化はＣ言語版になってからのの課題です。

やはり, 期待して待っていよう(^^;;;

>keiju> ここどういうつもりですか? って, 私のコードがこうなっていますね(^^;;; 
>keiju> 昔はこうでなくては駄目だったんですが, 今はこっちの方が良いでしょう:
>keiju> 
>keiju>   @sign, @prov, @exp = num.instace_eval{[@sign, @prov, @exp]}
>
>え、こんなのあり？

今はありなんですね(^^;;; 

>keiju> ここもBASEが仮定されているような...
>keiju> 
>keiju>     n, d = @prov.divmod(BASE ** PREC)
>keiju> 
>keiju> ??
>
>ここは、１０進数文字列の 10 です。
>
>	3.14
>
>の 3 を取り出しています。
>＃他がてきと〜だと、狼少年になっちゃうなぁ ^^)l

良く読んでいませんでした. すいませーん(__;;;

>ごろごろ出て来ますね。 @_@

(^^;;;

>keiju> >  def <=>(other)
>keiju> >  def +(other)
>keiju> >  def -(other)
>keiju> >  def *(other)
>keiju> >  def coerce(other)
>keiju> 
>keiju> ここなんですが, PRECの違うBigfloatについてまだ考慮がなされていませんね.

>そうなんです。
>実は、prov の大きい方に合わせるというアイディアを持っていてのですが、
>実装には移ってないです。

provですか??

そういえば,

  def +(other)
    case other
    when Bigfloat
      #break
    when Integer, Rational, Float
      other = Bigfloat(other)
                            ^^^^

ってなっていますよね. 有効桁数の低いものから高いものにしても意味がない
ので, 逆の方が良いと思っているんですが?

Bigfloat間でも PRECの低い方にあわせた方が良いと思っているんですが.



>keiju> あ. あと, クラス定数を使うと良いといいましたが駄目でした(__;;; 
>
>これ、しばらく考えたのですが。だめでした。
>で、降参してこのまま ML に投げたです ^^;;

(^^;;;

>keiju> # でも, モジュールにするとこれでは駄目でもうちょっと工夫が必要ですが...
>
>モジュールにすると、どんな使い道が出るのか今一つ判らないのですが、

抽象クラスだから, モジュールでもいいかなと思ったぐらいです. 無理に拘ら
なくても良いと思います.

>自由な精度(と基数？)の浮動小数点クラスが簡単に作れると良いですね。

これは, どちらでも可能です.

>keiju> PRECの違ったBigfloatどうしの演算は, 上記のようにPRECの違うBigfloatをサ
>keiju> ブクラス化するか否やかにもかかってきますね.
>keiju> 
>keiju> あと, 表に出すかどうかも...
>keiju> 
>keiju> また, インスタンスがprecをもつけれども, 演算ではprecの違いは考えず変換
>keiju> はcoerceで対応するというのもあるかなと... この方法だとクラスは1つにな
>keiju> ります.
>
>これだと
>	1.0/3.0
>もんだいの決着が必要ですね。

??

でも, 最後の案は, Presicionと相性がすごく悪いのでやめましょう.
Presicionはクラスが1つの確定した精度を持たないと困るみたいです.

# ですよね? >> ごとけんさん

__
................................石塚 圭樹＠日本ラショナルソフトェア...
----------------------------------->> e-mail: [e-mail] <<---