Fortranの「PURE」こそ、並列計算の守護神である理由
こんにちは。宇宙航空の現場で何十年も数値計算の荒波を乗り越えてきた、Fortranエンジニアです。
C言語やPythonから来た皆さんが、Fortranを触って最初に驚くのは「なぜこんなに制約が多いのか？」という点かもしれません。特に、今回解説する`PURE`というキーワード。これは単なる規約ではなく、「コンパイラに極限まで最適化を許し、かつ並列実行時の悪夢（デッドロックやデータ競合）を未然に防ぐための契約書」なんです。
今日は、なぜ`PURE`が重要なのか、そしてどう使いこなすべきかを、現場の泥臭い経験を交えてお話しします。
—
1. なぜ「副作用」が計算科学の敵なのか
Pythonでグローバル変数にフラグを立てて、複数のスレッドから適当に値を書き換える……そんなコードを書いて、結果が実行のたびに微妙にズレる「再現不可能なバグ」に苦しんだことはありませんか？
Fortranの`PURE`は、そんな惨劇を言語仕様レベルで封印します。`PURE`を宣言した手続き（関数やサブルーチン）の中では、以下の行為が物理的に不可能になります。
- グローバル変数（`COMMON`ブロックや`MODULE`内の変数）の変更
- I/O処理（`print`や`write`など、外部への出力）
- ポインタを介した意図しないメモリ領域への干渉
つまり、「入力に対して出力が確定する」という数学的な関数そのものになるわけです。
—
2. 実際に書いてみよう：PUREの作法
まずは、最も基本的な`PURE`関数の書き方です。ここでは「ある粒子にかかる力を計算する」ような関数を想定してみましょう。
module physics_engine
implicit none
! PUREを付与することで、この手続きが「クリーン」であることを保証します
! 変数の状態を壊さないため、コンパイラは安心して並列化の判断を下せます
pure function calculate_force(mass, acceleration) result(force)
real(8), intent(in) :: mass, acceleration
real(8) :: force
! サイドエフェクト（外部変数の書き換えなど）を一切行わない
force = mass acceleration
end function calculate_force
end module physics_engine
もしここで「ついでにログを出力しよう」と思って`print , force`と書いた瞬間に、コンパイラは「お前、契約違反だぞ！」とコンパイルエラーを吐きます。最初は窮屈に感じるかもしれませんが、これこそが「並列化の安全性」を担保する最強の防御壁なのです。
—
3. DO CONCURRENTとPUREの蜜月関係
モダンFortranの醍醐味である`DO CONCURRENT`は、ループの順序を問いません。つまり、ループの全イテレーションを同時に（並列に）走らせることができます。
ここで、呼び出す手続きが`PURE`であれば、コンパイラは「このループの各処理は互いに干渉しない」と確信できるため、CPUのベクトル演算ユニットをフル活用した爆速バイナリを生成してくれます。
! 高速化の鉄則：PURE手続きをDO CONCURRENTで呼び出す
do concurrent (i = 1:1000000)
forces(i) = calculate_force(masses(i), accels(i))
end do
もし`calculate_force`が`PURE`でなかったら、コンパイラは「もしかしたらループの途中でグローバル変数を書き換えて、次のループに影響を与えるかもしれない」と疑心暗鬼になり、安全のために最適化を諦めてしまいます。`PURE`の有無だけで、実行速度が数倍変わることはザラにあるのです。
—
4. 現場の知見：コンパイラを味方にする設定
最後に、皆さんが数値計算の現場で戦うための、コンパイラフラグのヒントを共有します。Intel Fortran (ifort/ifx) や gfortran を使う際は、以下のオプションを意識してください。
- `-qopenmp` (Intel) / `-fopenmp` (gfortran): 並列化を有効にする。
- `-O3 -march=native`: PURE手続きが揃っていれば、これだけでSIMD（ベクトル化）が劇的に効きます。
- `-check all` (開発中のみ推奨): `PURE`の規約を破っていないか、実行時にも厳密にチェックしてくれます。
—
最後に：なぜ「泥臭い」必要があるのか
皆さんが書くコードが数千、数万のコアで動くスーパーコンピュータ上の計算なら、わずかな「副作用」が数日間の計算を無に帰すバグの温床になります。
`PURE`をつけることは、単なるコーディング規約ではなく、「この計算は数学的に正しい」ということをコンパイラと共有するコミュニケーションです。
最初は「`print`もできないなんて！」と憤慨するかもしれませんが、一度この作法に慣れると、デバッグの時間が劇的に減ることに気づくはずです。さあ、まずは小さな関数から`pure`を付けて、その堅牢な世界を体験してみてください！
何か分からないことがあれば、いつでも聞いてくださいね。現場の壁を乗り越えるための知恵は、まだまだたくさんありますから。
コメント