ベクトル化を殺す「見えない鎖」：Fortranにおけるエイリアシング問題の克服
大規模な流体解析や構造解析のコードをプロファイリングすると、驚くほど単純なループでベクトル化（SIMD）が阻害されているケースに出くわす。原因の多くは、コンパイラが「この配列とあの配列は、メモリ上で重なっているかもしれない（エイリアシング）」と疑念を抱き、安全側（保守的）なコード生成を選択してしまうことに起因する。
Fortranは本来、C言語のポインタとは異なり、配列のメモリアクセスに関して強力な最適化のヒントをコンパイラに与えられる言語だ。今回は、この「エイリアシングの壁」を突き破り、ハードウェアの理論性能を限界まで引き出すための設計指針を伝授する。
—
1. なぜ「重なり」が最適化を止めるのか
コンパイラは、ループ内の `A(i) = B(i) + C(i)` をベクトル化する際、もし `A` と `B` がメモリ上で重なっていたら、代入の順序によって計算結果が変わってしまう可能性を考慮しなければならない。
! コンパイラが「エイリアシング」を疑う典型的なケース
subroutine update_field(a, b, n)
integer, intent(in) :: n
real(8), intent(inout) :: a(n), b(n)
integer :: i
do i = 1, n
a(i) = a(i) + b(i) ! もし a(1) と b(2) が同じアドレスを指していたら？
end do
end subroutine
この疑心暗鬼を取り除くのが、モダンFortranの重要な役割だ。
—
2. `CONTIGUOUS` 属性：最適化への強力なパスポート
Fortran 2008から導入された `CONTIGUOUS` 属性は、極めて強力だ。これはコンパイラに対し、「この配列（あるいはポインタ）はメモリ上で連続しており、かつ他のメモリ領域と重複していない」という強力な保証を与える。
特に、サブルーチン引数として配列を受け取る際、スライス（`A(1:N:2)` のような飛び飛びのアクセス）ではなく、連続したメモリブロックであることを明示することで、コンパイラは迷いなくSIMD命令を注入する。
実装例：最適化を最大化する設計
subroutine compute_kernel(n, vec_a, vec_b)
integer, intent(in) :: n
! CONTIGUOUS属性により、配列が連続メモリにあるとコンパイラに確約する
! これにより、コンパイラはポインタの重複チェックを省略できる
real(8), intent(inout), contiguous :: vec_a(:)
real(8), intent(in), contiguous :: vec_b(:)
integer :: i
! $OMP SIMD (明示的なベクトル化指示を添えるとなお良い)
do i = 1, n
vec_a(i) = vec_a(i) vec_b(i)
end do
end subroutine
—
3. `RESTRICT` 相当の挙動をどう導き出すか
C言語の `restrict` キーワードに相当するものは、Fortranにおいては「配列が引数として渡される際、それらは重ならない（Alias-freeである）」という標準の仕様ルールがある。しかし、ポインタやターゲット属性を持つ変数を使用すると、この保証が崩れる。
避けるべきアンチパターンと修正案
! 【悪い例】ターゲット属性の多用はエイリアシングの温床
real(8), pointer, target :: p1(:), p2(:)
! 【良い設計】
! サブルーチンに渡す際は、可能な限り生の配列を渡し、
! 内部で仮引数にCONTIGUOUSを付与する。
! どうしてもポインタを使うなら、C_LOCやC_F_POINTERの適切な管理が不可欠。
—
4. 現場で効くコンパイラフラグの調整
どれだけコードを綺麗に書いても、コンパイラの設定が「安全第一」では宝の持ち腐れだ。以下は、Intel Fortran (ifx/ifort) や gfortran で必須となるオプションである。
- Intel Fortran (ifx/ifort):
`-O3 -xHost -qopt-report=5 -qopt-report-phase=vec`
- `-qopt-report` を必ず活用せよ。どのループがなぜベクトル化できなかったか、コンパイラが律儀に教えてくれる。
- gfortran:
`-O3 -march=native -ffast-math -fopt-info-vec-optimized`
- `-ffast-math` は浮動小数点の結合法則を無視することを許容するオプションだが、ベクトル化においては劇的な性能向上をもたらすことが多い。
—
まとめ：堅牢かつ高速なコードのために
1. 引数には必ず `CONTIGUOUS` を付与せよ：これによりコンパイラはエイリアシングの懸念を捨て、アグレッシブなループアンロールとSIMD命令を展開する。
2. `intent(in)` を活用せよ：コンパイラは `intent(in)` が指定された配列に対し、他の配列からの書き込みがないという追加の最適化ヒントを得る。
3. プロファイリングを怠るな：`opt-report` を読み解く習慣を持つ者が、結局のところ一番速いコードを書く。
数値計算における「10%の高速化」は、一週間分の計算時間が数時間短縮されることを意味する。エイリアシングという「見えない鎖」を断ち切り、演算器を休ませないコードを書くこと。それが我々エンジニアの矜持だ。
コメント