導入
大規模な数値計算において、配列の受け渡しは避けて通れない処理です。しかし、Fortranなどの言語で配列を関数やサブルーチンに渡す際、「配列全体を渡す」場合と「部分配列（スライス）を渡す」場合では、コンパイラの振る舞いが大きく異なります。この違いを理解していないと、意図しないメモリコピーが発生し、計算速度の低下やメモリ不足を招く原因となります。本記事では、高性能な数値計算を実現するための配列の受け渡し戦略を解説します。
基礎知識
配列を引数として渡す際、コンパイラは「記述子（Descriptor）」を介してメモリ上のデータにアクセスします。
配列全体（A）を渡す場合、元の配列のメモリレイアウトを直接参照する記述子が生成されます。
一方、部分配列（A(1:N)）を渡す場合、コンパイラは「指定範囲のみを指す新しい記述子」を生成します。もし、この部分配列がメモリ上で連続していない（Non-contiguous）と判断されると、コンパイラは一時的なコピーを作成して連続領域に詰め直してからサブルーチンへ渡すというコストの高い処理を実行します。これが「スライス生成のパフォーマンス影響」の正体です。
実装/解決策
一時的なメモリコピーを回避するためには、サブルーチン側の引数定義で「メモリの連続性」をコンパイラに明示することが極めて重要です。具体的には、`contiguous`属性を付与することで、コンパイラに対して「この引数はメモリ上で連続したデータとして受け取る」という最適化のヒントを与えます。また、`intent`属性を併用することで、メモリの保護と最適化の精度を両立させます。
サンプルプログラム
以下は、連続性を保証した状態で配列を受け取るためのテンプレートコードです。
module array_util
implicit none
contains
! intent(inout)でデータの修正を許可し、contiguousでメモリ連続性を保証
subroutine process_data(arr)
! メモリ上で連続していることをコンパイラに伝える
real(8), intent(inout), contiguous :: arr(:)
integer :: i
! ループのベクトル化が促進されやすくなる
do i = 1, size(arr)
arr(i) = arr(i) 2.0d0
end do
end subroutine process_data
end module array_util
program main
use array_util
real(8), allocatable :: A(:)
allocate(A(100))
A = 1.0d0
! ケース1: 配列全体を渡す（最適）
call process_data(A)
! ケース2: 部分配列を渡す（contiguousが宣言されていればコピーを回避できる可能性がある）
call process_data(A(1:50))
end program main
応用・注意点
現場での開発において注意すべき点は、「データの断片化」です。例えば、多次元配列の列方向のスライスなどを渡す場合、物理的にメモリが連続していないケースが多々あります。このような状況で`contiguous`を強制すると、実行時エラーになる可能性があるため注意が必要です。
また、古いコンパイラや一部の最適化レベルでは、`contiguous`を指定していても内部でコピーが発生することがあります。重要な関数や計算のボトルネックとなる箇所では、コンパイラの最適化レポート（Intel Fortranであれば `-qopt-report` など）を確認し、「Temporary array created」という警告が出ていないかを必ずチェックしてください。メモリレイアウトの設計を最初から「列優先（Fortranの場合）」に最適化しておくことが、最も確実なパフォーマンス向上策となります。
コメント