導入：なぜ初期化の書き方が重要なのか？
数値計算プログラムにおいて、動的メモリ確保（ALLOCATE）は頻繁に行う処理の一つです。通常、メモリを確保した後に「DOループ」などで値を代入していく書き方が一般的ですが、これには「コードが長くなる」「処理ステップ数が増える」という欠点があります。今回紹介する「SOURCEキーワード」を使った初期化手法を用いることで、メモリ確保と値の代入を1行で完了させ、プログラムの可読性と実行速度を同時に向上させることができます。
基礎知識：動的メモリ確保と初期化の仕組み
FORTRANなどの言語における動的メモリ確保とは、プログラムの実行中に必要な分だけメモリ領域を確保する仕組みです。通常、ALLOCATE文で領域を確保した直後の配列は「不定値」となっており、そのまま計算に使うとバグの原因になります。そのため、従来は確保後に値を書き込む必要がありました。しかし、モダンなFORTRAN（Fortran 2003以降）では、確保と同時に値を流し込む機能が標準化されています。これが今回解説するSOURCEオプションです。
実装：SOURCEを用いた一括初期化
実装は非常にシンプルです。ALLOCATE文の括弧内に「SOURCE = [値のリスト]」と記述するだけです。この記述により、コンパイラは「この配列は確保直後にこの値で埋まる」と判断し、内部的に最適化を施します。具体的には、メモリ確保と同時に値をコピーするルーチンが呼び出されるため、手動でループを回すよりも効率的なメモリ操作が行われます。
サンプルプログラム：そのまま動かせるコード例
以下のコードをコピーして、コンパイラで実行してみてください。
program init_example
implicit none
! 動的配列の宣言
real, allocatable :: vec(:)
integer :: i
! allocate文で領域確保と同時に初期値を代入
! SOURCEを指定することで、ループを書かずに初期化が完了します
allocate(vec(3), source=[1.0, 2.0, 3.0])
! 結果の確認
print , “初期化された配列の内容:”
do i = 1, size(vec)
print , “vec(“, i, “) =”, vec(i)
end do
! 使い終わったら忘れずに解放
deallocate(vec)
end program init_example
応用・注意点：現場で陥りやすい罠
この手法を用いる際、いくつか注意すべきポイントがあります。
1. データ型の不一致に注意
SOURCEに指定する値の型（整数や実数）と、宣言した配列の型は一致させてください。型が異なると、予期せぬ型変換が発生して性能が低下したり、値が正しく代入されなかったりすることがあります。
2. 配列サイズの整合性
確保するサイズとSOURCEの要素数は必ず一致させてください。例えば `allocate(vec(5), source=[1, 2, 3])` のように書くとエラーや予期せぬ動作の原因となります。
3. 大規模データへの適用
定数配列を直接書く手法は、要素数が少ない場合には非常に有効ですが、数万個の要素を初期化する場合は、ソースコードが肥大化してしまいます。その場合は、別のファイルから読み込むか、別の初期化関数を併用するなど、適材適所での利用を心がけましょう。
この一括初期化をマスターすることで、コードが一段と引き締まり、意図が明確なプログラムを作成できるようになります。ぜひ日々の開発に取り入れてみてください。
コメント