- 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
- 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 ]
- 保守
|
 |
