【入門編】SIMDベクトル化を阻害する配列エイリアシング問題 – モダンFortran言語仕様と実践実践マスター

Fortranの「罠」を突破せよ：SIMDを殺す「エイリアシング」と戦う技術

こんにちは。長年、スパコンの冷却ファンの音を子守唄にして、数値計算の泥沼を泳いできた元研究機関のエンジニアです。

これからFortranを触る皆さんは、きっと「なぜFortranは計算が速いのか？」という伝説に惹かれていることでしょう。しかし、書き方次第でその伝説はあっさり崩壊します。今日は、現代のコンパイラが「本気」を出してSIMD（単一命令複数データ）演算を叩き込むのを、我々人間が「うっかり」邪魔してしまうエイリアシング問題についてお話しします。

1. なぜ「重なり」が速度を殺すのか

C言語やPythonの経験がある方なら、「ポインタ」の扱いで一度は頭を抱えたことがあるはずです。Fortranでも、サブルーチンに配列を渡すとき、コンパイラは心の中でこう呟いています。

「おい、この渡された2つの配列、実はメモリ上の同じ場所を指していないか？もしそうなら、一方の値を書き換えた瞬間に、もう一方の値も変わってしまう。SIMDで並列計算なんてしたら結果が壊れるぞ。……念のため、SIMDは諦めて逐次処理にしておこう」

これがエイリアシング（別名参照）問題です。コンパイラは安全側に倒して最適化を放棄します。しかし、我々エンジニアは「この配列とあの配列が重なるはずがない」と知っている。この「確信」をコンパイラに伝える方法を知るだけで、コードの速度が数倍から数十倍に跳ね上がります。

2. 現場で使える「最強の呪文」：CONTIGUOUS

現代のFortran（Fortran 2008以降）には、この問題を解決するための非常に強力な属性が用意されています。それが `CONTIGUOUS` です。

まずは、よくある「やってはいけない」例を見てみましょう。

subroutine compute_update(a, b, n)
integer, intent(in) :: n
real(8), intent(inout) :: a(n)
real(8), intent(in) :: b(n)
integer :: i

! コンパイラは「aとbが重なっているかも？」と疑い、
! SIMD化（ベクトル化）を躊躇します。
do i = 1, n
a(i) = a(i) + b(i)
end do
end subroutine

これに対し、`CONTIGUOUS`属性を追加して「メモリ上で重なりはなく、連続している」ことを保証します。

subroutine compute_update_optimized(a, b, n)
integer, intent(in) :: n
! ここで「メモリ上で連続している」とコンパイラに宣言する
real(8), intent(inout), contiguous :: a(n)
real(8), intent(in), contiguous :: b(n)
integer :: i

! これにより、コンパイラは迷いなくSIMD命令を生成します！
do i = 1, n
a(i) = a(i) + b(i)
end do
end subroutine

`CONTIGUOUS`を付けることで、コンパイラは「この配列はメモリ上できっちり隙間なく並んでいるので、安心してベクトルレジスタにまとめてロードしてくれ」という強力なサインを受け取ることができます。

3. さらに一歩先へ：コンパイラへの「強制力」

もしあなたがIntelコンパイラ（ifort/ifx）やgfortranを使っているなら、コードの修正に加えて、コンパイラの最適化オプションを見直すことも忘れないでください。

Intelコンパイラの場合: `-qopt-report` を付けてコンパイルすると、どのループがベクトル化できて、どのループがエイリアシングによって「阻害」されたかがログに全て出力されます。これがデバッグの最強の武器です。
gfortranの場合: `-fno-alias` や `-fargument-noalias` というオプションを試してみてください。これは「関数の引数同士は決して重ならない」という前提をコンパイラに強制するフラグです。

4. 若手エンジニアへのアドバイス：まずは「意識」から

「なぜ自分のコードが遅いのか？」と思ったとき、まずはそのループがSIMD化されているかを確認してください。

1. 配列のアクセス順序は正しいか？（Fortranは列優先。`a(i, j)` なら `i` を内側のループにするのが鉄則です）
2. エイリアシングを疑え。 `CONTIGUOUS` は使ったか？
3. コンパイラの最適化レポートを読め。 自分の書いたコードが、機械にどう「翻訳」されたかを確認する習慣こそが、一流への近道です。

Fortranは古い言語と言われますが、メモリレイアウトをこれほどまでに制御できる言語は他にありません。皆さんが書くその計算コードが、スパコンの演算ユニットをフル稼働させる喜びを、ぜひ体験してください。

もし行き詰まったら、またここへ戻ってきてください。次は「配列のストライドとキャッシュヒット率」という、さらに深い沼でお待ちしていますよ！