導入
現代のプログラミングでは動的メモリ確保が一般的ですが、数値計算の現場、特にレガシーなFortranコードや高性能計算（HPC）ライブラリを扱う際、「整合配列（Adjustable Array）」の概念は避けて通れません。整合配列は、手続きの呼び出し時に配列のサイズが確定する仕組みです。なぜ今さら古い技術を学ぶ必要があるのか。それは、BLASやLAPACKといった標準的な数値計算ライブラリが、メモリの連続性を前提としたこの形式をベースに設計されているためです。本記事では、整合配列の仕組みと、現代的な開発現場で安全に扱うためのポイントを解説します。
基礎知識
整合配列とは、関数の引数として渡される配列のサイズを、別の引数（NやMなど）で指定する手法を指します。
例えば、`subroutine sub(a, n)` という手続きの中で `real :: a(n)` と宣言すると、呼び出し側で指定されたサイズに従って配列が「形状」付けられます。
重要なのは、この手法が「メモリの物理的なコピー」を発生させない点です。呼び出し元のメモリ領域を直接参照するため、計算コストが極めて低く、キャッシュ効率を最大化できます。ただし、コンパイラによる次元チェックが実行時に依存するため、境界外アクセスには注意が必要です。
実装/解決策
整合配列を扱う際の鉄則は、「実引数のサイズを確実に渡すこと」です。多次元配列の場合、特に列優先（Column-major）のメモリ配置を意識する必要があります。
以下に、行列演算の基本となる二次元配列を渡す際の実装例を示します。
サンプルプログラム
以下のコードは、行列の各要素に値を加算する簡単な例です。
subroutine process_matrix(a, n, m)
! 引数n, mは行列の行数と列数
integer, intent(in) :: n, m
! 整合配列の宣言：n行m列の行列として扱う
real, intent(inout) :: a(n, m)
integer :: i, j
! 行列要素へのアクセス
do j = 1, m
do i = 1, n
! 整合配列はメモリ上で連続しているため、
! このようにインデックスで直接操作が可能
a(i, j) = a(i, j) + 1.0
end do
end do
end subroutine
! 呼び出し側の例
! program main
! real :: mat(10, 20)
! call process_matrix(mat, 10, 20)
! end program
応用・注意点
整合配列を使用する現場で最も注意すべきは、「形状の不一致によるセグメンテーションフォールト」です。呼び出し側で確保したメモリサイズよりも、手続き側で指定したサイズ（n, m）が大きい場合、メモリの境界を越えて不正な領域を読み書きしてしまいます。
これを回避するための現代的なTipsとして、以下の2点を推奨します。
1. インターフェースブロックの活用
明示的なインターフェース（module内で手続きを定義する等）を使用することで、コンパイラに引数のチェックを行わせることが可能です。これにより、誤ったサイズの配列を渡した際にビルド時や実行時に警告を出しやすくなります。
2. 配列ポインタやAllocatable配列との併用
レガシーなライブラリに渡す際は整合配列を利用しつつ、自身のコード内では動的配列（allocatable）を使用する「ハイブリッドな設計」を心がけてください。ポインタを経由することで、メモリ管理の安全性と、BLASライブラリとの親和性を両立できます。
整合配列は古い手法に見えますが、メモリレイアウトを直接制御できるという点で、数値計算エンジニアにとっては依然として強力な武器となります。仕組みを正しく理解し、安全に活用していきましょう。
コメント