導入: なぜ想定形状配列が必要なのか？
数値計算において、サブルーチンに配列を渡す際、配列サイズを毎回引数として指定するのは非常に面倒でミスも起きやすいものです。特に、呼び出し元で確保した配列のサイズをサブルーチン側で正確に把握できないと、メモリ外アクセスなどの重大なバグにつながります。この課題を解決するのが「想定形状配列（Assumed-shape Array）」です。これを使うことで、配列の次元情報を安全に引き継ぎ、直感的なコードを書けるようになります。
基礎知識: 想定形状配列とは？
想定形状配列とは、サブルーチンの引数宣言において、サイズを数値で指定する代わりに「:（コロン）」を使って次元数だけを定義する手法です。Fortranの古い規格（F77）では、配列サイズを明示的に渡す必要がありましたが、現代のFortranでは、配列自体が自身のサイズ情報（メタデータ）を保持しています。想定形状配列は、その情報を呼び出し先で自動的に受け取る仕組みです。これにより、サブルーチン内でも配列の形状を意識せず、安全にスライシングや全体演算が可能になります。
実装/解決策: 基本的な使い方
実装は非常にシンプルです。サブルーチンの引数リストで、配列の次元分だけ「:」を記述するだけです。このとき、呼び出し元のプログラム側でインターフェースが明確である必要があります（モジュール内にサブルーチンを記述するのが最も簡単で推奨される方法です）。
サンプルプログラム
以下のコードをコピーして、コンパイル・実行してみてください。配列のサイズを気にせず、行列の和を計算する例です。
program main
implicit none
real, dimension(2, 3) :: data_a, data_b, result
! 配列に値を代入
data_a = 1.0
data_b = 2.0
! サブルーチンを呼び出す（サイズを渡す必要はありません）
call add_arrays(data_a, data_b, result)
print , “計算結果の要素(1,1): “, result(1,1)
contains
! 想定形状配列を使用したサブルーチン
subroutine add_arrays(a, b, res)
real, intent(in) :: a(:,:), b(:,:) ! 形状を自動で引き継ぐ
real, intent(out) :: res(:,:) ! 結果用配列
! 配列全体演算が可能
res = a + b
! スライシング（一部だけ取り出す）も安全に行える
! 例: res(1, 🙂 = a(1, 🙂 + b(1, 🙂
end subroutine add_arrays
end program main
応用・注意点: 現場で役立つアドバイス
1. インターフェースの重要性: 想定形状配列を利用するには、呼び出し元がそのサブルーチンの引数情報を知っている必要があります。必ず `module` 内にサブルーチンを記述するか、`interface` ブロックを使用してください。これを怠ると、コンパイル時にエラーになったり、予期せぬ挙動を引き起こしたりします。
2. 形状の確認: サブルーチン内で配列のサイズを知りたい場合は、`size(arr, 1)` や `ubound(arr, 1)` といった組み込み関数を使用してください。これにより、引数で渡された配列のサイズを動的に取得できます。
3. パフォーマンス: 現代のコンパイラは非常に最適化されており、想定形状配列を使用したからといって計算が遅くなることはほとんどありません。むしろ、サイズチェックが自動化されるため、バグを防ぐメリットの方が圧倒的に大きいです。
まずは既存のプログラムの引数を「:」に書き換えて、コードがどれだけスッキリするか試してみてください。
コメント