導入
数値計算の現場において、計算結果のサイズが実行時まで分からない、あるいは条件によって変動するケースは非常に多いです。Fortranの「割付配列（Allocatable Array）」は、こうした動的なメモリ管理に欠かせない機能です。本記事では、サブルーチン内でメモリを確保し、その実体を呼び出し元へ安全かつ効率的に引き継ぐ「引き継ぎ割付け」のテクニックを解説します。この手法を習得することで、データ生成ファクトリのような柔軟なモジュール設計が可能になります。
基礎知識
Fortranにおける割付配列とは、プログラム実行時にヒープメモリを動的に確保・解放できる配列のことです。通常、`allocatable`属性を付与して宣言し、`allocate`文でメモリを確保します。
重要な点は、サブルーチンに渡す際の`intent`属性です。`intent(out)`を指定すると、サブルーチンに入る直前に配列のメモリが解放され、サブルーチン内で確保したメモリがそのまま呼び出し元へ引き継がれます。これにより、呼び出し元でサイズを事前に知らなくても、計算結果を受け取ることが可能になります。
実装/解決策
メモリの引き継ぎを実現する鍵は、サブルーチンの引数に`allocatable`属性を付与することです。これにより、単なる配列のポインタ渡しではなく、メモリ管理の所有権そのものをサブルーチンに受け渡しできます。
1. 引数に`intent(out)`を指定し、`allocatable`属性を付与します。
2. サブルーチン内で`allocate`文を実行します。
3. 呼び出し元では、配列を宣言するだけでメモリ確保は不要です。
サンプルプログラム
以下のコードは、指定されたサイズで配列を生成し、呼び出し元へデータを引き継ぐ実用的な例です。
! 割付配列の引き継ぎサンプル
program main
implicit none
! 配列を宣言するが、ここではメモリ確保は不要
real, allocatable :: data(:)
integer :: n
n = 10
! サブルーチンでメモリを確保し、結果を受け取る
call generate_array(data, n)
print , "生成された配列のサイズ:", size(data)
print , "先頭の値:", data(1)
! 使用後は明示的に解放することが推奨される
deallocate(data)
end program main
subroutine generate_array(a, n)
! intent(out)により、呼び出し元のメモリを一旦クリアし、
! この中で確保したメモリを呼び出し元に引き継ぐ
real, allocatable, intent(out) :: a(:)
integer, intent(in) :: n
! ここでメモリを確保。この実体が呼び出し元に返る
allocate(a(n))
a = 1.0 ! 初期値を代入
end subroutine generate_array
応用・注意点
現場での開発において注意すべき点がいくつかあります。
まず、メモリリークの防止です。`intent(out)`を使用した場合、サブルーチンに渡す前に配列が既に確保されていたとしても、古いメモリは自動的に解放されます。しかし、プログラム終了時には呼び出し元で必ず`deallocate`を実行する習慣をつけましょう。
また、意図しない解放を避けるため、既存のデータに追記したい場合は`intent(inout)`を使用してください。`intent(out)`は「古いデータを完全に破棄して新しく作り直す」場合に適しています。
最後に、大規模な計算では、頻繁な`allocate`/`deallocate`は断片化の原因となるため、再利用可能なバッファを用意するなどの工夫を組み合わせると、より堅牢な数値計算環境を構築できます。
コメント