導入
数値計算の現場において、レガシーなFortranコードをメンテナンスしていると頻繁に遭遇するのが、引数なしの「SAVE文」です。これは手続き内のローカル変数の値を、呼び出し終了後も保持させる強力な機能ですが、現代の並列計算環境（OpenMPなど）では、意図しないデータ競合を引き起こす重大なリスク源となります。本記事では、SAVE文の仕組みと、現代的な開発において安全なコードに書き換えるための考え方を解説します。
基礎知識
Fortranにおいて、手続き内のローカル変数は原則としてその手続きの実行終了とともに解放されます（スタック領域への配置）。しかし、SAVE文を指定すると、その変数は「静的領域（Static Memory）」に割り当てられます。
引数なしのSAVE文を記述すると、そのスコープ内のすべてのローカル変数が静的領域に配置されます。これはF77時代にはメモリ節約や計算速度向上のために多用されましたが、現代ではマルチスレッド実行や再帰呼び出しが前提となっており、複数のスレッドから同一の静的変数を書き換えることで、計算結果が壊れるという深刻な課題を抱えています。
実装/解決策
レガシーコードをモダナイズする際は、SAVE文の「一括指定」をやめ、必要な変数のみを個別に指定する形にリファクタリングすることを推奨します。また、OpenMP環境では、該当する変数をスレッドごとに独立させる「THREADPRIVATE」属性や、そもそも静的変数に頼らない設計（引数として変数を渡す手法）への変更が求められます。
サンプルプログラム
以下のコードは、SAVE文が引き起こすスレッド不安全の典型例と、その改善の方向性を示したものです。
SUBROUTINE CALCULATE_TOTAL(INPUT_VAL)
! 従来のSAVE文：全変数が静的になり、並列実行時に値が競合する
! SAVE
INTEGER, INTENT(IN) :: INPUT_VAL
! 改善策：必要な変数のみに限定し、可能であれば引数で管理する
INTEGER, SAVE :: ACCUMULATOR = 0
! スレッド並列時にここが競合して計算が破綻するリスクがある
ACCUMULATOR = ACCUMULATOR + INPUT_VAL
PRINT , “現在の合計値:”, ACCUMULATOR
END SUBROUTINE
! 現代的な推奨実装：状態を保持せず引数で完結させる
SUBROUTINE CALCULATE_TOTAL_SAFE(INPUT_VAL, CURRENT_SUM)
INTEGER, INTENT(IN) :: INPUT_VAL
INTEGER, INTENT(INOUT) :: CURRENT_SUM
CURRENT_SUM = CURRENT_SUM + INPUT_VAL
END SUBROUTINE
応用・注意点
現場でSAVE文を削除または修正する際は、以下の点に注意してください。
1. 初期化のタイミング:
SAVE文を外すと、変数は呼び出しごとに初期化されます。以前のコードが「前回の計算結果を保持すること」を前提としたアルゴリズムの場合、プログラムの論理構造自体を修正する必要があります。
2. モジュール変数との混同:
モジュール内に定義された変数は、デフォルトでSAVE属性を持ちます。これらも同様に並列計算時のデータ競合の温床となります。可能であれば、計算に必要な状態は「Derived Type（構造体）」にまとめ、サブルーチン引数として渡す設計に切り替えるのが、最も堅牢でテストしやすいコードとなります。
3. コンパイラオプションの活用:
古いFortranコードを移植する際は、コンパイラの警告オプション（例：gfortranの -Wsurprising や -fcheck=all など）を有効にし、予期せぬ変数の静的化が行われていないか定期的にチェックすることをお勧めします。
コメント