1. 導入：なぜ「コピー」が大規模計算のボトルネックになるのか
数値計算の現場において、処理速度を低下させる隠れた要因の一つが、手続き（サブルーチンや関数）への引数受け渡し時に発生する「一時的なメモリコピー」です。特に、多次元配列の一部（配列セクション）を切り出して渡す際、コンパイラはメモリ上の整合性を保つために、一度連続したメモリ領域へ値をコピー（コピー・イン）し、計算終了後に元の場所へ書き戻す（コピー・アウト）という処理を自動で行うことがあります。
データサイズが数ギガバイト規模になると、このコピー処理だけで膨大な時間とメモリ帯域が浪費されます。本記事では、FortranのCONTIGUOUS属性を活用し、この不要なオーバーヘッドを排除する方法を解説します。
2. 基礎知識：なぜコピーが発生するのか
Fortranにおいて、配列セクション（例: A(1:100:2) といったストライド指定）は、メモリ上で飛び飛びの配置になることがあります。手続き側が「連続したメモリ空間」を期待している場合、コンパイラは引数の不整合を防ぐために、一時的な作業領域を作成して値を詰め替えます。
これを回避するためのキーワードがCONTIGUOUS属性です。これを宣言することで、コンパイラに対して「この引数は必ずメモリ上で連続している」という保証を与え、無駄なコピー処理の生成を抑止できます。
3. 実装・解決策
解決策は以下の2点です。
1. 手続き側の仮引数に CONTIGUOUS 属性を明示する。
2. 呼び出し側で、配列が連続していることを保証する（または連続した配列を渡す）。
手続き側でCONTIGUOUSを指定すると、コンパイラは「非連続なデータが渡された場合、コンパイルエラー（あるいは実行時エラー）にする」という制約を課すようになります。これにより、開発者は意図しないコピーが発生していないことを確信できるようになります。
4. サンプルプログラム
以下は、行列計算の一部を処理するサブルーチンの例です。
! 手続き側：CONTIGUOUS属性を付与することでコピーを抑制する
subroutine process_column(vec)
! 連続したメモリ配置を期待する宣言
real(8), contiguous, intent(inout) :: vec(:)
integer :: i
! 連続メモリとコンパイラが判断できるため、最適化が効きやすくなる
do i = 1, size(vec)
vec(i) = vec(i) 2.0d0
end do
end subroutine process_column
program main
real(8), target :: matrix(100, 100)
! 呼び出し側：連続している列を渡すことでコピーを回避
! 1列目はメモリ上で連続しているため、コピーは発生しない
call process_column(matrix(:, 1))
end program main
5. 応用・注意点
現場で活用する際の重要な注意点を挙げます。
1. 実行時エラーのリスク
CONTIGUOUSを指定した引数に対して、非連続な配列セクション（例: matrix(1:100, 1:100:2) など）を渡すと、プログラムが異常終了する可能性があります。呼び出し側の配列が必ず連続しているか、事前にis_contiguous()組み込み関数でチェックする設計にしてください。
2. コンパイラの最適化オプションとの併用
コピー・イン/アウトは、コンパイラの最適化レベル（-O3など）に依存します。CONTIGUOUSを指定しても、コードの構造次第ではSIMDベクトル化が阻害される場合があります。プロファイラを使用して、実際にコピーが発生していないか確認することが、大規模計算の現場では不可欠です。
3. 現代的なFortranの活用
もし新しいコードを書くのであれば、配列の受け渡しをポインタではなく、形状引き継ぎ配列（Assumed-shape array）とCONTIGUOUSを組み合わせるのが、現代的なFortranにおける最も高速かつ安全な手法です。
コメント