1. 導入：なぜこの制約を知っておく必要があるのか
数値計算の現場では、メモリを効率的に扱うために「割付配列（Allocatable Array）」を多用します。Fortranの強力な機能の一つに、代入時に配列サイズが自動で調整される「自動再割付け」がありますが、実はこれには大きな落とし穴があります。それは「型の異なるデータは代入できない」という制約です。この仕様を知らないと、プログラムが予期せぬエラーで停止したり、バグの原因になったりします。今回は、この動的メモリ管理の限界と、正しく実装するためのコツを解説します。
2. 基礎知識：割付配列と自動再割付けの仕組み
Fortranの割付配列は、プログラムの実行中に必要なサイズだけメモリを確保できる便利な機能です。
「自動再割付け」とは、例えば `A = B` と記述した際、`A` のサイズが `B` と異なっていても、自動的に `A` のサイズを `B` に合わせて拡張・縮小してくれる機能を指します。
しかし、この機能はあくまで「同じ型（例：どちらも単精度実数）」であることを前提としています。単精度（REAL）の配列に倍精度（DOUBLE PRECISION）の配列を代入しようとすると、コンパイラは「型が一致しない」としてエラーを返します。これはメモリ安全性を担保するための重要な設計です。
3. 実装と解決策：型変換が必要な場合の正しい作法
型が異なる配列間で値をコピーしたい場合は、自動再割付けに頼るのではなく、明示的な手順を踏む必要があります。
1. 目的の型で配列を再確保（allocate）する。
2. 型変換を伴う代入を行う。
という手順が基本です。
4. サンプルプログラム：型安全な動的メモリ伸長
以下は、単精度のデータを倍精度の配列へ安全に転送する例です。そのままコピーしてコンパイルと実行を試してみてください。
プログラム：
program dynamic_allocation
implicit none
! 単精度の配列を定義
real, allocatable :: single_arr(:)
! 倍精度の配列を定義
double precision, allocatable :: double_arr(:)
integer :: i
! 1. 単精度配列にデータを格納
allocate(single_arr(3))
single_arr = [1.1, 2.2, 3.3]
! 2. 倍精度配列へ変換してコピーする準備
! 自動再割付けはできないため、事前に形状を合わせる
if (allocated(double_arr)) deallocate(double_arr)
allocate(double_arr(size(single_arr)))
! 3. 型変換を伴う代入を実行
double_arr = dble(single_arr)
! 結果の表示
print , “変換後の倍精度配列:”
do i = 1, size(double_arr)
print , double_arr(i)
end do
! 後片付け
deallocate(single_arr, double_arr)
end program dynamic_allocation
5. 応用・注意点：現場で陥りやすい罠
現場で最も注意すべきは、「暗黙の型変換」を期待してコードを書かないことです。
特に、関数の引数として配列を渡す際、仮引数と実引数の型が異なると、コンパイラオプションによっては警告が出ないまま誤った値が渡されることがあります。
また、`allocate` を繰り返す際に `deallocate` を忘れるとメモリリークの原因となります。必ず `if (allocated(arr)) deallocate(arr)` というイディオムを用いて、既存のメモリを解放してから再確保する習慣をつけましょう。数値計算の安定性は、こうしたメモリ管理の厳密さから生まれます。
コメント