堀之内です。

> > (2)はやばいと思います。NArray 等の拡張ライブラリーを使う場合や、
> > Fortran や C での計算と比較する場合と、丸め誤差が一致しなくなり
> > ます。カオスなどに代表される、初期誤差が増大するケースでは浮動
> > 小数点方式がかわると結果が大きく変わります。そういう問題なんだ
> > からしょうがないという考えかたもありえますが、数値計算自体の検
> > 証可能性が確保できなくなります。
> 
> Fortran や C では、コード上での演算順序と
> コンパイル結果の演算順序が同じになる保証がないですよね。
> この問題に対しては、最適化を極力抑制 (-O0 とか) して
> 演算の順序が保存されることを期待して乗り切るんでしょうか。

Fortran や C も上げたのは勇み足でした。仰るとおり、
そこを気にする場合はオプション等で指定することになると
思います。

＃ 私は実際の研究ではそういう点でかなりアバウトな計算しか
   してなかったりします...。(気象畑です。計算量多すぎ。)



--
ML: yarv-dev quickml.atdot.net
Info: http://www.atdot.net/~ko1/quickml