ruby-dev

ごとけんです

In message "[ruby-dev:5839] Re: Precision in Mathematica (Re: Bigfloat)"
    on 99/03/01, EGUCHI Osamu <eguchi@shizuokanet.ne.jp> writes:

>keiju> >keiju> 1. 四則演算は厳密に精度を計算し, 関数系は引数の精度までを引き継ぐ
>keiju> >keiju> 2. 今迄通り精度は大きい方に合わせていく(float系の場合)
>keiju> >keiju> 3. Mathematicaに頑張って合わせる
>keiju> >keiju> 
>keiju> >keiju> ぐらいしか思い浮かばないですが...
>keiju> >
>keiju> >1. についてですが、これは『関数は引数の精度と同等の精度を
>keiju> >得られる事が期待できるアルゴリズムを使用せよ』
>keiju> >と読み取って良いのでしょうか？（物によってはメチャ厳しい ^^;;）
>keiju> 
>keiju> そうじゃなくて, 大変なんであきらめるってことです(^^;;;
>
>永久に返って来ないより、なんぼかましですね ^^;;
>
>既に Bignum ** Bignum のように、時間に対する問題に関しては
>なにがなんでも値を求めるのでなく、実行時間配慮する方向で
>態度が決まってると考えています。
>＃どのへんで、あきらめるかは問題ですが、

いや，Bignum は時間がかかっても計算させたほうが
いいんではないかと思ってます．実装の仕方によるけど
そんなに時間かからないんじゃないかなぁ．
四則演算は最大誤差を別個に評価して，仮の結果と最大誤差を
足して変化しない桁を調べればいいんではないかと．
関数と言っても，四則演算の級数か，ニュートン法でであらわ
すんだから，引数の精度のうちで最小のものを使って必要な
項や収束条件は求まるから問題ないような．
あ，ニュートン法は危険かも．

>keiju> >案外 3. が上の解釈の 1. そのものだったりして、、、
>keiju> 
>keiju> うーん. どうなんでしょう? >> ごとけんさん?

ぼくは，1. をいまいち理解してなかったけど，
とりあえずは Math モジュールの関数は prec(Float) を呼ぶと
言うことでいいんではないでしょうか??

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

exp って，指数関数のことですか??

-- gotoken

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