ruby-dev

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

>>> In message [ruby-dev:6044] Re: Bigfloat [compare to NaN] 
    On Sat, 6 Mar 1999 18:43:15 +0900, keiju@rational.com (石塚圭樹 ) said:

keiju> けいじゅ＠日本ラショナルソフトウェアです.
keiju> 
keiju> In [ruby-dev :5966 ] the message: "[ruby-dev:5966] Re: Bigfloat
keiju> [compare to NaN] ", on Mar/04 14:05(JST) EGUCHI Osamu writes:
keiju> 
keiju> >えぐち＠エスアンドイー です。
keiju> 
keiju> >keiju> だんだんメイルが溜ってきた(^^;;;
keiju> >
keiju> >フロー制御の閾値の選択が、むつかしいですね ^^;;;
keiju> 
keiju> ですね. えぐちさんのメイルだいぶ捨ててしまいました(__;;;

そういうことわたしも（よく）あります、^^;;
＃しろやぎさんたらよまずにたべたさっきのメールのごようじなーに

keiju> >Bigfloat でも加算などで、結果の有効数字を指定できると、
keiju> >計算の中間精度を妥当なところで、見切ってしまえるので、
keiju> >速度面で（特定の場面ではメモリ使用量でも）有効だと思います。
keiju> 
keiju> いちおう, coerceで精度の低い方に合わせるのでもおなじ効果がありますが...

その時だけってケースもあるのでは？？？
あ、coerce を特異メソッドにすればいいのですね ^^;;;

keiju> >ただ、演算子の引数を増やす術がないので、（ですよね）
keiju> >
keiju> >	Bigfloat#add_with_prec(right, prec = @prec)
keiju> >
keiju> >のようなメソッドを、用意する必要があるのが難点です。
keiju> 
keiju> いや定義できるんですが.

え゛

keiju> class Foo
keiju>   def +(other, prec)
keiju>     ...
keiju>   end
keiju> end
keiju> 
keiju> は定義できるんですが, 呼び出しがsendしかできないという何点があります
keiju> (^^;;;
keiju> 
keiju>   foo.send :+, other, prec

をぉ、これは何かの時に使えそうだ。

keiju> >	low_prec(x, y)
keiju> >
keiju> >の様な引数で精度の小さい方の精度に、高い方の精度を
keiju> >合わせる関数を用意するのも手ですが、使いにくそ ^^;;;;
keiju> 
keiju> それならcoerceで自動的にやってくれた方が...

そうなんだけど、組み込みのクラスが coerce を使ってない
場合もあるので、 Fixnum#coerce を再定義しても

	1+1.1 => 2

には出来なかったりするのでした。
＃この例では、無限精度の Fixnum に合わせていて、逆ですが、、
(Fixnum#+ を再定義すればいいのだけれど、、、)

fix_plus(x, y) の

    if (TYPE(y) == T_FLOAT) {
        return rb_float_new((double)FIX2LONG(x) + RFLOAT(y)->value);
    }

の様な、「言語組み込みの数値の型変換」はなくすべきかも知れません。
＃効率の問題があるのですが ^^)l

keiju> >keiju> mpって四則演算だけですよね? それだったらあまり採用してもね...
keiju> >
keiju> >そうですね、四則の他は 2 ^ n, abs, sqrt や比較程度なので、
keiju> >スタートラインとしては、かなり厳しいかも知れない ^^;;;
keiju> 
keiju> ですね.

関数ライブラリの多くは、
  テーブルルックアップ法で第一近似を得てから
  ニュートン法などで必要な精度を得る
戦略で実装されていますが、可変精度でこれが使えるのか、
考えています。
sqrt に対する ** に様な、十分精度の出そうな質の良い逆関数が
あれば問題ないのですが、sin() などでも、目的の精度が出たかの判断に
難ありですね。
(atan と三平方で判定すると、それ自体の精度に疑問が出そう、、)

  def cos(t) sqrt(1 - sin(t) ** 2) end

でバランスさせる等の工夫も必要そうだし、、
＃ Float を越える精度で sin() を得たい人がいるかは疑問 ^^)l

	えぐち

Thread

Prev Next

In This Thread

Prev Next

[#5809] mod_ruby-0.0.7 — Shugo Maeda <shugo@...>

[#5868] Re: Precision in Mathematica (Re: Bigfloat) — "豊福" <toyofuku@...>

[#5895] [REQ] 中置の def — EGUCHI Osamu <eguchi@...>

[#5906] Re: [REQ] 中置の def — Inaba Hiroto <inaba@...>

[#5931] Re: [REQ] 中置の def — Inaba Hiroto <inaba@...>

[#5934] Re: Subseq of Array — Inaba Hiroto <inaba@...>

[#5945] \xff doesn't work fine? — Yasuhiro Fukuma <yasuf@...>

[#5948] Re: Precision in Mathematica (Re: Bigfloat) — "豊福" <toyofuku@...>

[#5955] Re: a genericity — toyofuku@...

[#5965] Re: Precision in Mathematica (Re: Bigfloat) — keiju@... (石塚圭樹 )

[#5970] [REQ] variable of regex — keiju@... (Keiju ISHITSUKA)

[#5978] [BUG] Rational#initialize w/ !Integer === den — EGUCHI Osamu <eguchi@...>

[#5982] Re: [REQ] 中置の def — "Takehiko Ueki" <takehiko@...>

[#5991] Re: Subseq of Array — Inaba Hiroto <inaba@...>

[#5995] Re: [REQ] 中置の def — "Takehiko Ueki" <takehiko@...>

[#6009] Re: [REQ] 中置の def — "Takehiko Ueki" <takehiko@...>

[#6016] [Req] Thread#name — Shugo Maeda <shugo@...>

[#6034] Re: [REQ] 中置の def — "Takehiko Ueki" <takehiko@...>

[#6039] Forward: [ruby-talk:00243] win32 lib — matz@... (Yukihiro Matsumoto)

[#6059] [PATCH] str_index — IKARASHI Akira <ikarashi@...>

[#6063] [PATCH] io_write — IKARASHI Akira <ikarashi@...>

[#6064] ruby-1.3.1-990225 — IKARASHI Akira <ikarashi@...>

[#6089] Thread.current.join — Shugo Maeda <shugo@...>

[#6090] Thread#kill — Shugo Maeda <shugo@...>

[#6093] Re: sieve.rb — Inaba Hiroto <inaba@...>

[#6100] [PATCH] my_setenv() for cygwin — WATANABE Hirofumi <watanabe@...>

[#6120] Re: [REQ] 中置の def — "Takehiko Ueki" <takehiko@...>

[#6121] Re: [REQ] 中置の def — "Takehiko Ueki" <takehiko@...>

[#6122] Re: [REQ] 中置の def — "Takehiko Ueki" <takehiko@...>

[#6141] Retry of mod_ruby — moochan <moochan@...>

[#6148] String#index — IKARASHI Akira <ikarashi@...>

[#6175] unpack — IKARASHI Akira <ikarashi@...>

[#6176] Re: [ruby-list:12717] String#unpack — matz@... (Yukihiro Matsumoto)

[#6181] Making of Interval class. — moochan <moochan@...>

[#6206] [bug?] ftplib.rb — ttate@...

[#6213] Ruby 1.3.1-990311 — Yukihiro Matsumoto <matz@...>

[#6220] Re: [ruby-list:12721] __FILE__,__LINE__ への代入 (Re: [REQ] rlex.rb) — keiju@... (石塚圭樹 )

[#6229] exception re-raise in Thread#join — Shugo Maeda <shugo@...>

[#6237] [REQ] p -10 == p(-10) — Koji Arai <JCA02266@...>

[#6242] "a\n".gsub(/$/, '-') — Wakou Aoyama <wakou@...>

[#6257] 米国からコンタクトレンズの個人輸入のご案内 !!! — DM <dmservice@...>

[#6259] str += — WATANABE Tetsuya <tetsu@...>

[#6276] Re: [ruby-list:12784] String#index, String#unpack — IKARASHI Akira <ikarashi@...>

[#6278] ruby-mode.el — IKARASHI Akira <ikarashi@...>

[#6281] Re: "a\n".gsub(/$/, '-') — Inaba Hiroto <inaba@...>

[#6282] Re: "a\n".gsub(/$/, '-') — Inaba Hiroto <inaba@...>

[#6288] [REQ] interh.h — Shugo Maeda <shugo@...>

[#6313] [PATCH] nkf.c — WATANABE Hirofumi <eban@...>

[#6317] Re: Bigfloat [compare to NaN] — "豊福" <toyofuku@...>

[#6319] Re: Bigfloat [compare to NaN] — "豊福" <toyofuku@...>

[#6323] [REQ] mkmf.rb — Yasuhiro Fukuma <yasuf@...>

[#6328] Re: Ruby 1.3.1-990315 — WATANABE Hirofumi <watanabe@...>

[#6330] sprintf — Tadayoshi Funaba <tadf@...>

[#6332] Re: Ruby 1.3.1-990315 — ttate@...

[#6334] Re: Ruby 1.3.1-990315 — Koji Arai <JCA02266@...>

[#6335] Re: Ruby 1.3.1-990315 — Koji Arai <JCA02266@...>

[#6336] security level — ttate@...

[#6337] Re: Ruby 1.3.1-990315 — WATANABE Tetsuya <tetsu@...>

[#6351] patch to ruby-1.3.1-990315, suppress warnings — Ryo HAYASAKA <hayasaka@...>

[#6356] [PATCH] ext/extmk.rb — WATANABE Hirofumi <watanabe@...>

[#6357] Re: Ruby 1.3.1-990315 — Yoshida Masato <yoshidam@...>

[#6358] Re: Bigfloat [compare to NaN] — toyofuku@...

[#6359] Re: array[last..first] (Re: Subseq of Array)" — toyofuku@...

[#6391] value of NODE_DEFN — Shugo Maeda <shugo@...>

[#6402] class_eval{} SEGVs in 1.3.1-990315 — Inaba Hiroto <inaba@...>

[#6403] [BUG?] bignum "&" and "or" — IKARASHI Akira <ikarashi@...>

[#6404] [BUG] select — Wakou Aoyama <wakou@...>

[#6408] Re: class_eval{} SEGVs in 1.3.1-990315 — Inaba Hiroto <inaba@...>

[#6411] [PATCH] END, at_exit (Re: profiler) — WATANABE Hirofumi <watanabe@...>

[#6415] Kconv.kconv("", Kconv::EUC, Kconv::JIS) — Wakou Aoyama <wakou@...>

[#6417] nonzero? should return nil — Inaba Hiroto <inaba@...>

[#6425] [BUG] unsaved context — ttate@...

[#6426] wrong line number — Koji Arai <JCA02266@...>

[#6427] invalid argument check in specify_eval() — Koji Arai <JCA02266@...>

[#6434] Time#+ — Koji Arai <JCA02266@...>

[#6446] sub("string", "string") — Wakou Aoyama <wakou@...>

[#6449] Re: Ruby 1.3.1-990315 — KIMURA Koichi <kimura@...>

[#6450] Re: [REQ] Object#+@ as identity (Re: [REQ] p -10 == p(-10)) — Inaba Hiroto <inaba@...>

[#6451] AC_C_CONST (Re: Ruby 1.3.1-990315) — EGUCHI Osamu <eguchi@...>

[#6467] Re: regexp — Inaba Hiroto <inaba@...>

[#6220] Re: [ruby-list:12721] FILE,LINE への代入 (Re: [REQ] rlex.rb) — keiju@... (石塚圭樹 )

[#6521] how to control hash in C — 助田　雅紀 <masaki.suketa@...>

[#6522] RE: how to control hash in C — 助田　雅紀 <masaki.suketa@...>

[#6526] Re: how to control hash in C — 助田　雅紀 <masaki.suketa@...>

[#6535] Re: how to control hash in C — 助田　雅紀 <masaki.suketa@...>

[#6536] Re: how to control hash in C — 助田　雅紀 <masaki.suketa@...>