1. 導入：なぜこの問題が重要なのか
数値計算の現場で古くから使われているFORTRANコードを保守していると、稀に「なぜかサブルーチンの計算結果が呼び出しごとに変わってしまう」「並列化すると値が壊れる」といった不可解なバグに遭遇することがあります。その原因の多くは、メモリの番地を共有するEQUIVALENCE文と、変数の寿命を固定するSAVE属性の予期せぬ干渉にあります。本稿では、この「隠れた波及効果」がなぜ危険なのか、そのメカニズムと対策を解説します。
2. 基礎知識：EQUIVALENCEとSAVE属性の役割
EQUIVALENCE文は、異なる名前の変数を同一のメモリ領域に割り当てる機能です。メモリを節約するために多用されていた過去の遺産ですが、現代のプログラミングではメモリ上の意図しない書き換えを招く原因となります。
一方、SAVE属性は、サブルーチンや関数が終了しても、その中の変数の値を保持し続けるための属性です。通常、ローカル変数は関数終了時に破棄されますが、SAVEを付けることで「静的変数」として扱われます。問題は、この2つを組み合わせたときに発生します。
3. 実装と論理的な解説
EQUIVALENCE文で接続された変数は、物理的に同じメモリアドレスを指しています。そのため、片方にSAVE属性を付与すると、たとえもう片方にSAVEを付けていなくても、そのメモリ領域全体がSAVE属性を持つことになります。
これにより、本来は「呼び出しのたびに初期化されるはずの変数（B）」が、意図せず「前回の値を保持する変数（A）」として振る舞うようになります。これがサブルーチンの「再入可能性（リエントラント性）」を破壊し、マルチスレッド環境での競合や、反復計算での不具合を誘発するのです。
4. サンプルプログラム
以下のコードは、意図せずSAVE属性が波及してしまう例です。
SUBROUTINE CALCULATE_VAL()
IMPLICIT NONE
REAL :: A, B
! AにSAVE属性を付与
SAVE A
! AとBを同じメモリ番地に割り当てる
EQUIVALENCE (A, B)
! 本来、Bはローカル変数として初期化されてほしいが、
! AがSAVEされているため、Bの値も保持されてしまう
PRINT , “現在のBの値:”, B
B = B + 1.0
PRINT , “更新後のBの値:”, B
END SUBROUTINE
このプログラムを複数回呼び出すと、Bの値は毎回0から始まるのではなく、前回の計算結果に1を加算し続ける挙動を示します。
5. 応用・注意点：現場で陥りやすいバグの回避策
現場でのトラブルを避けるためには、以下の原則を守ることが重要です。
・EQUIVALENCEの全廃：現代のFortran（Fortran 90以降）では、EQUIVALENCEを使う必要性はほぼありません。構造体（TYPE）やポインタ、あるいはALLOCATABLE配列を利用することで、メモリ管理を安全に行えます。
・SAVEの明示：SAVE属性を使う場合は、影響範囲が明確になるよう、特定の変数だけでなくサブルーチン全体にSAVEを付けるか、あるいはモジュール変数としてスコープを明確にする設計に移行してください。
・再入可能性の意識：並列計算（OpenMP等）を行う場合、SAVE属性を持つ変数は共有メモリとして扱われ、データ競合の直接の原因となります。可能な限りローカル変数を使用し、値を保持する必要がある場合は引数として受け渡す設計を推奨します。
レガシーコードの修正を行う際は、まずはEQUIVALENCE文を検索し、その変数がSAVE属性と絡んでいないかを確認することから始めてみてください。
コメント