コンパイラのバージョンの違いが引き起こす計算結果の大きな相違

しばらくの間，Fortranコンパイラのバージョンの違いによって結果が大きく変わる（古いコンパイラではうまく計算できるが，新しい方では発散する）現象に悩まされました．単精度計算をしていたのですが，倍精度で計算する必要があったというのが原因でした．新しいコンパイラはより単精度に忠実になったせいなのでしょうか．倍精度計算をデフォルトにし，その後単精度計算を適用するのが安全です．