fortran66のブログ

fortran について書きます。

Kahan和

Intel Fortran の intrinsic function の SUM は、大きな配列の和をとる時、配列の並びで大きく値が変わる。(浮動小数点では交換則が成り立っていないが、その誤差が大きく見えてくる。)これは、倍精度変数をアキュームレータにするか、いわゆる Kahan 和をとることで解消される。