- 1 名前：デフォルトの名無しさん [2009/07/20(月) 12:36:39 ]
- C++用 無限精度実数ライブラリRealLibについて語るスレです。
本家
The RealLib Project
ttp://www.brics.dk/~barnie/RealLib/
[XR] RealLib ソースコード勉強会 (1:終)
ttp://d.hatena.ne.jp/hzkr/20080915
[XR] RealLib ではじめる誤差ゼロ実数計算
ttp://d.hatena.ne.jp/hzkr/20080909
- 32 名前：デフォルトの名無しさん mailto:sage [2009/07/24(金) 20:26:15 ]
- お前の類似性判断基準はライブラリかそうでないかしかないのかw
- 33 名前：デフォルトの名無しさん mailto:sage [2009/07/24(金) 20:31:08 ]
- >>32
別にそれだけとは言わねぇけどさ。
- 34 名前：デフォルトの名無しさん mailto:sage [2009/07/24(金) 23:49:10 ]
- けどなんなんだよ
- 35 名前：デフォルトの名無しさん mailto:sage [2009/07/25(土) 00:11:15 ]
- まあそう煽るな。
ライブラリでないものはライブラリでないと見抜ける人でないと（C++を使うのは）難しい
ってことを伝えたかっただけだ。
- 36 名前：デフォルトの名無しさん mailto:sage [2009/07/25(土) 00:14:57 ]
- シンプレックス法とかの吐き出し法アルゴは
誤差に弱いし、超越関数使ってないし、効果すぐ出そうだよね。
ちょっと試してみようかな。
、、と思ったけど、入力パラメータが実数な時点で
既に誤差ってるんだった。
FPUとSIMD-FPの代わりに、こいつをHWインプリメントとか
出来んもんかね。
- 37 名前：デフォルトの名無しさん mailto:sage [2009/07/25(土) 04:37:56 ]
- Real(0.1) → 誤差あり
Real("0.1") → 誤差なし
か。なるほど。
- 38 名前：デフォルトの名無しさん mailto:sage [2009/07/25(土) 09:52:09 ]
- 「アルゴ」と略すと頭悪そうに見える件
- 39 名前：デフォルトの名無しさん mailto:sage [2009/07/25(土) 10:09:15 ]
- >>37
Real(0.1) → 誤差あり
Real("0.1") → 誤差なし
Real(1.5) → 誤差なし
Real("1.5") → 誤差なし
ふーん。なるほど。
- 40 名前：デフォルトの名無しさん [2009/07/25(土) 10:17:40 ]
- rとr'が完全に等しい時、
比較演算子==で比較すると無限ループするという仕様は
ちょっと悲しい。
妥協する方法はないのだろうか？
- 41 名前：デフォルトの名無しさん mailto:sage [2009/07/25(土) 10:20:55 ]
- アルゴ君ネタはもういいよ
- 42 名前：デフォルトの名無しさん mailto:sage [2009/07/25(土) 10:21:44 ]
- >>40
解析的に等しいかどうかのチェックはやっぱ難しいのかね
- 43 名前：デフォルトの名無しさん mailto:sage [2009/07/25(土) 11:32:38 ]
- しかし時間は有限なので
結局処理に充てられる時間の都合上有限精度
無限精度であることは唯の一度も無い！
- 44 名前：デフォルトの名無しさん mailto:sage [2009/07/25(土) 11:36:43 ]
- >>19 expression templateみたいなものか。
- 45 名前：デフォルトの名無しさん [2009/07/25(土) 12:05:15 ]
- >>43
そんな誰もが解り切っていることを言わなくても。
そもそも>>6によれば実装上の限界が定まっているらしいじゃないか。
- 46 名前：デフォルトの名無しさん mailto:sage [2009/07/25(土) 12:22:30 ]
- 1/2^2147483648 まで表現できれば十分・・・とは言えない場合ってどんなんだろう
- 47 名前：デフォルトの名無しさん [2009/07/25(土) 12:50:14 ]
- sqrt(2)とかを求めたらメモリーを使い切ってcoredumpってこと？
- 48 名前：デフォルトの名無しさん mailto:sage [2009/07/25(土) 12:52:33 ]
- sqrt(2) を sqrt(2) のまま扱う
有効桁数を指定して表示させることもできる
という Mathematica と同じようなことをしてると思われ
- 49 名前：デフォルトの名無しさん [2009/07/25(土) 12:58:01 ]
- じゃあ対数を多用する計算でも誤差伝搬を無くせるの？
- 50 名前：デフォルトの名無しさん mailto:sage [2009/07/25(土) 13:18:34 ]
- 複素数は無理なのか
- 51 名前：デフォルトの名無しさん mailto:sage [2009/07/25(土) 13:19:27 ]
- >>47
お前はこのスレで何を読んできたんだ？ｗ
- 52 名前：デフォルトの名無しさん mailto:sage [2009/07/25(土) 13:22:28 ]
- しかし sqrt(2) == sqrt(2) は無限ループ！
- 53 名前：デフォルトの名無しさん mailto:sage [2009/07/25(土) 15:00:40 ]
- 数桁比較して等しければ
解析的に等しいかどうかチェックした方がいいんだろうけど、
Mathematica でもこのあたり完璧にできるの？
等価な式を引き算して simplify したら 0 になるとか
- 54 名前：デフォルトの名無しさん [2009/07/26(日) 00:26:40 ]
- Mathematicaと区別つかないやつ大杉
- 55 名前：デフォルトの名無しさん mailto:sage [2009/07/26(日) 00:40:15 ]
- じゃあどこが違うんだよ
ライブラリかどうかとか言ったら嗤うぞw
- 56 名前：デフォルトの名無しさん mailto:sage [2009/07/26(日) 10:04:35 ]
- >>55
1. ライブラリかどうかも「当然」違うけど
2. 文字処理なども違うだろ？
3. 制作者も違えば
4. ライセンスも違う。
5. 普及率も違う。
むしろどこが同じなの？
- 57 名前：デフォルトの名無しさん mailto:sage [2009/07/26(日) 10:17:28 ]
- つまらん
- 58 名前：デフォルトの名無しさん mailto:sage [2009/07/26(日) 10:55:35 ]
- >>57
良いから答えろw
- 59 名前：デフォルトの名無しさん mailto:sage [2009/07/27(月) 17:35:54 ]
- 「アキレスと亀」をデモしてくれるんなら、使ってみる。
- 60 名前：デフォルトの名無しさん mailto:sage [2009/07/27(月) 19:19:29 ]
- アキレスと亀って、要はフラクタルでしょ
マンデルブローの拡大デモにアキレスと亀のスプライトを入れれば完成
- 61 名前：デフォルトの名無しさん mailto:sage [2009/07/27(月) 19:32:07 ]
- パラドックスライブラリでも作るか
- 62 名前：デフォルトの名無しさん [2009/07/27(月) 21:48:30 ]
- 無限精度なんてあるんだ。
面白いな。
いつ使うのか知らないけど。
- 63 名前：デフォルトの名無しさん mailto:sage [2009/07/27(月) 22:32:50 ]
- もちろんベンチマークする時に決まってるじゃん
- 64 名前：デフォルトの名無しさん [2009/07/28(火) 07:06:35 ]
- アキレスの亀の本質は
負ける前に時間を停めればその時刻までは負けてない
- 65 名前：デフォルトの名無しさん [2009/07/28(火) 10:35:06 ]
- そんなくだらない本質を語らなくても。
RealLibについて語るしかなかろう。
- 66 名前：デフォルトの名無しさん mailto:sage [2009/07/28(火) 10:54:25 ]
- おれの用途にはこれくらいの精度で十二分だなあ。
速度はどの程度でるんだろうか。
ちなみに「無限」とか言わずに、精度10^x倍ライブラリ、とした方が
関心持つやつが増えると思うよ＞開発者
- 67 名前：デフォルトの名無しさん mailto:sage [2009/07/28(火) 11:24:20 ]
- .dk で開発されてるライブラリについてここでアドバイスしても
届くわけねーだろ
- 68 名前：デフォルトの名無しさん mailto:sage [2009/07/28(火) 11:29:29 ]
- .dkって
デンマーク？
- 69 名前：デフォルトの名無しさん mailto:sage [2009/07/28(火) 11:46:26 ]
- ドナルド・クヌース
- 70 名前：デフォルトの名無しさん mailto:sage [2009/07/29(水) 13:30:58 ]
- 諸君、話題は無いのかね？
- 71 名前：デフォルトの名無しさん [2009/07/31(金) 18:56:38 ]
- 無いっすねぇ(>_<)
- 72 名前：デフォルトの名無しさん [2009/08/02(日) 23:54:59 ]
- 誰かタメになる話題をplz
- 73 名前：名無しさん [2009/08/04(火) 04:45:52 ]
- 「精度保証付き数値計算」でググれ
- 74 名前：デフォルトの名無しさん mailto:sage [2009/08/09(日) 11:08:28 ]
- ググってみた。
面白そうだが、
精度保証付き数値計算 C++ライブラリ
でググってもライブラリが見つからん。
- 75 名前：デフォルトの名無しさん [2009/08/09(日) 12:11:36 ]
- boostの有理数ライブラリと比べてどうかな？
あっちは除算によるまるめ誤差は出ないんだっけか？
- 76 名前：74 mailto:sage [2009/08/09(日) 12:15:56 ]
- >>75
有理数ライブラリは有理数演算用で、
小数計算には役立たない。
例えばπとかには完全に無力。
- 77 名前：デフォルトの名無しさん [2009/08/11(火) 16:00:41 ]
- 保守
- 78 名前：名無しさん [2009/08/12(水) 18:49:46 ]
- >>74
特定のライブラリではない汎用の方法がある。
条件は：
（１）IEEE-754をサポートしている
（２）対象とする計算が線形である
アセンブラかC/C++のインラインアセンブリで、丸め方法を捜査する。
（答え） ＝ ｛（全部切り上げで計算した結果）＋（全部切り下げで計算した結果）｝ / ２
（精度） ＝ ABS｛（全部切り上げで計算した結果）ー（全部切り下げで計算した結果）｝
- 79 名前：デフォルトの名無しさん mailto:sage [2009/08/12(水) 19:02:39 ]
- >>78
ありがとう。
でも線形であるって相当厳しい要求じゃないか？
sinどころか累乗が出てきただけで死ぬっしょ？
- 80 名前：デフォルトの名無しさん mailto:sage [2009/08/12(水) 19:18:28 ]
- 適当に線形化したら…
やっぱだめか
- 81 名前：名無しさん [2009/08/12(水) 21:29:30 ]
- CORDICとかはあかんね
でも、累乗は何とかなるかも
Matlabは使ったことないが、Octave風に書くと
X = [1, 2, 3, ..., N]'
X2 = X .* X;
X3 = X2 .* X;
（ダッシュ’は転置。 .* 内席ではなく要素同士の掛け算）
とすると
X3 == [1^3, 2^3, 3^3, 4^3, 5^3, 6^3, ... , N^2]'
- 82 名前：デフォルトの名無しさん [2009/08/15(土) 13:16:19 ]
- あげあげ
- 83 名前：デフォルトの名無しさん [2009/08/21(金) 00:31:20 ]
- 保守
- 84 名前：名無しさん＠そうだ選挙に行こう [2009/08/30(日) 07:50:30 ]
- 保守
|
 |
