64bit OSで動かすと32bitの時と処理結果が異なるプログラムで悩んだメモです。

調べていくとfloatを使った浮動小数点演算で誤差が出ることがわかりました。分からなかったのが、IEEE754でfloatの演算精度は32bitと規定されているのになぜ同じソースから作ったプログラムで結果が違うのかということです。

x86では浮動小数点演算にx87を使うのですが、この演算精度は32bitや64bitではなく、80bitで行われるそうです。この辺の話は 浮動小数点演算ではまった話 – bkブログとかBK通信に詳しく載ってます。要約すると32bitの精度で誤差が生じる計算が含まれていても”たまたま”うまく計算できることがある ということです。

たまたま精度よく計算できるプログラムをx86_64に持っていってコンパイルしようとします。このとき、x86_64では必ずSSEが使えるため、コンパイラはデフォルトでSSEを使うようにコンパイルします。SSEで浮動小数点演算を行うと演算精度がx87の80bitよりも悪くなるため、本来生じるはずだった演算誤差が表面化することになります。

高林さんがBK通信で書いているMac OSで演算結果が違うというのも、Mac OSが64bitOSであることが関係している気がします。

バイナリーハックス買おうかな。

Binary Hacks ―ハッカー秘伝のテクニック100選

posted with amazlet at 08.12.09

高林 哲 鵜飼 文敏 佐藤 祐介 浜地 慎一郎 首藤 一幸
オライリー・ジャパン
売り上げランキング: 23306

おすすめ度の平均:

ハードコア？なソフトウエア
大工さんにおける電動工具の紹介本
当然教科書ではない。でも、とても参考になります。
バイナリアンの基本

Amazon.co.jp で詳細を見る

オレンジニュースでlog5j – Google Codeが紹介されていたので、log4XXって各言語にあるのかなと気になったので調べてみました。

• log4c: Log4c : Logging for C Library(C)
• Log for C++(C++)
• Log4cplus Project Page(C++)
• log4cxx – About(C++)
• Apache log4j 1.2 – log4j 1.2(Java)
• log5j – Google Code(Java)
• Apache log4net – Apache log4net: Home(.NET)
• log4php – Apache log4php(PHP)
• Log4r – Manual(Ruby)
• SourceForge.net: Log4Py(Python)
• perl.com: Retire your debugger, log smartly with Log::Log4perl!(Perl)
• Log4js – The Logging Framework for JavaScript(JavaScript)
• Log4E – Home(Eiffel)

思いつく言語で見つかったのだとこんな感じかな? 他にもご存知の方がいましたら教えて下さい。