ruby-dev

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

>>> In message [ruby-dev:5842] Re: Precision in Mathematica (Re: Bigfloat)
    On Mon, 1 Mar 1999 19:34:40 +0900, gotoken@math.sci.hokudai.ac.jp (GOTO Kentaro) said:

gotoken> ごとけんです
gotoken> 
gotoken> In message "[ruby-dev:5839] Re: Precision in Mathematica (Re: Bigfloat)"
gotoken>     on 99/03/01, EGUCHI Osamu <eguchi@shizuokanet.ne.jp> writes:
gotoken> 
...
gotoken> >keiju> >1. についてですが、これは『関数は引数の精度と同等の精度を
gotoken> >keiju> >得られる事が期待できるアルゴリズムを使用せよ』
gotoken> >keiju> >と読み取って良いのでしょうか？（物によってはメチャ厳しい ^^;;）
gotoken> >keiju> 
gotoken> >keiju> そうじゃなくて, 大変なんであきらめるってことです(^^;;;
gotoken> >
gotoken> >永久に返って来ないより、なんぼかましですね ^^;;
gotoken> >
gotoken> >既に Bignum ** Bignum のように、時間に対する問題に関しては
gotoken> >なにがなんでも値を求めるのでなく、実行時間配慮する方向で
gotoken> >態度が決まってると考えています。
gotoken> >＃どのへんで、あきらめるかは問題ですが、
gotoken> 
gotoken> いや，Bignum は時間がかかっても計算させたほうが
gotoken> いいんではないかと思ってます．実装の仕方によるけど
gotoken> そんなに時間かからないんじゃないかなぁ．

Bignum ** Bignum の一番小さな組合せは、

  (2 ** 30) ** (2 ** 30)

だと思いますが、これでも時間の問題と記憶容量の問題で
実際に評価すると 'failed to allocate memory' でしょう。
、、、例が悪かったなあ ^^;;;

gotoken> 四則演算は最大誤差を別個に評価して，仮の結果と最大誤差を
gotoken> 足して変化しない桁を調べればいいんではないかと．
gotoken> 関数と言っても，四則演算の級数か，ニュートン法でであらわ
gotoken> すんだから，引数の精度のうちで最小のものを使って必要な
gotoken> 項や収束条件は求まるから問題ないような．
gotoken> あ，ニュートン法は危険かも．

これも時間が問題？
正面突破がニュートン法より早い関数もあるだろし、、
＃libm のソース読も ^^;; <- きっとビット操作の塊だ！

gotoken> >keiju> >案外 3. が上の解釈の 1. そのものだったりして、、、
gotoken> >keiju> 
gotoken> >keiju> うーん. どうなんでしょう? >> ごとけんさん?
gotoken> 
gotoken> ぼくは，1. をいまいち理解してなかったけど，
gotoken> とりあえずは Math モジュールの関数は prec(Float) を呼ぶと
gotoken> 言うことでいいんではないでしょうか??

その精度で満足できない人は、自分で書きましょう ^^)l
でもって、順次公開しましょう。 (^^)

gotoken> >keiju> expの精度がどのように変化するのか興味ありますね.
gotoken> >
gotoken> >exp がそうだかわかりませんが、機械の浮動小数点で同じ関数を
gotoken> >実行して、第一近似を得て適当な精度になるまでトライする
gotoken> >という手法が使えそうです。
gotoken> 
gotoken> exp って，指数関数のことですか??

そぅです(ょね ^^;)

	えぐち

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@...>