メモリ断片化との決別：ALLOCATABLEによる動的ヒープ制御の「極致」
多くの研究者がF77時代の名残である「巨大な静的COMMONブロック」や、スタック領域を食いつぶす「自動配列」の亡霊に憑りつかれている。数万コアを動員するHPC環境において、スタックサイズ制限（`ulimit -s`）を気にしながらコードを書くなど、現代の数値計算アーキテクトとしてはナンセンス極まりない。
今回は、Fortran 2003以降の洗練された`ALLOCATABLE`管理と、それが現代のメモリハイアラキーに与える深い影響について、現場の知見を叩き込む。
—
1. 「暗黙の再割当て」という甘美な罠
Fortran 2003から導入された「自動再割当て（Automatic Reallocation）」は、左辺の配列サイズが右辺と異なるとき、コンパイラが自動で`DEALLOCATE`と`ALLOCATE`を代行してくれる強力な機能だ。しかし、HPCの現場において、ループ内部でこれを不用意に叩くことは、ヒープ領域の断片化（Fragmentation）を招き、最悪の場合、メモリリークではなく「メモリの断片化によるアロケーション失敗」を引き起こす。
数千億要素のグリッドを扱う際、メモリを細かく確保・解放する挙動は、OSのカーネルレベルでのロック競合を引き起こす。
極意：
ループ内での自動再割当ては厳禁だ。コンパイルオプションで `-fno-realloc-lhs` を検討するほどシビアな設計をせよ。配列の形状が確定しているなら、最初に必要十分な量を確保し、それを再利用（Dirty reuse）するのが鉄則である。
—
2. キャッシュラインを意識したメモリレイアウト
Fortranは伝統的に「列優先（Column-major）」である。この仕様を理解せず、`ALLOCATABLE`で確保したメモリにアクセスするのは、CPUキャッシュを捨てる行為に等しい。
以下の例を見てほしい。
! 非推奨：キャッシュミスを誘発するループ順序
do j = 1, ny
do i = 1, nx
A(i, j) = B(i, j) C(i, j)
end do
end do
これに対し、`A(i, j)`を`i`方向（列方向）に走査する際は、キャッシュライン（通常64バイト）を最大限活用できるよう、メモリ上の連続性を保証しなければならない。
高度な最適化手法：
大規模行列演算を行う際、`ALLOCATE`する配列の境界をページサイズ（4KB）やキャッシュラインサイズでアライメントする工夫が必要になる場合がある。Intelコンパイラなら `!DIR$ ATTRIBUTES ALIGN:64 :: A` といった指示文を併用し、ベクトル化（AVX-512等）の効率を最大化せよ。
—
3. オブジェクトのライフサイクルと「メモリの墓場」
`ALLOCATABLE`な変数を`TYPE`の中に埋め込み、それを`MPI_COMM_WORLD`の数千プロセスで保持する場合、デストラクターの呼び出し（`DEALLOCATE`）を怠ると、クラスター全体が数時間でメモリ不足に陥る。
特に、Fortran 2008以降の `MOVE_ALLOC` を活用せよ。これはメモリのコピーを発生させずにポインタの所有権を移動させる極めて強力な組み込みサブルーチンだ。
! メモリコピーなしで配列の所有権を移譲する例
subroutine resize_buffer(old_arr, new_size)
real(8), allocatable, intent(inout) :: old_arr(:)
real(8), allocatable :: temp_arr(:)
integer, intent(in) :: new_size
allocate(temp_arr(new_size))
! 必要であればデータをコピー
temp_arr(1:size(old_arr)) = old_arr
! 所有権を移動し、古い領域を自動的にDEALLOCATE
call move_alloc(from=temp_arr, to=old_arr)
end subroutine
—
4. VTune/Scalascaを用いたボトルネックの可視化
「感覚」でメモリ管理を行ってはならない。`Intel VTune Profiler`を用いれば、メモリ帯域幅の飽和や、NUMAノードを跨いだアクセス（Remote Access）によるレイテンシの増大をミリ秒単位で特定できる。
- NUMA最適化: `ALLOCATE`したスレッドと、実際にそのメモリを計算で使うスレッドが同じノードに配置されているか？ `numactl –interleave=all` 等で誤魔化すのではなく、`OpenMP`の`proc_bind`と`places`を駆使し、メモリ局所性（Memory Locality）を物理的に固定せよ。
まとめ：アーキテクトとしての矜持
モダンFortranは、もはや「古い言語」ではない。`ALLOCATABLE`や`MOVE_ALLOC`、そして`ISO_C_BINDING`を駆使すれば、C++に劣らないメモリ制御能力を、Fortran特有の最適化のしやすさと共に享受できる。
コードを書くとき、常に頭の中に「メモリコントローラからCPUキャッシュまでのデータ経路」を描け。コンパイラを信じるな、コンパイラが吐き出したバイナリをプロファイラで糾弾せよ。それが、スパコンの計算資源を極限まで引き出す、我々エンジニアの責務である。
次回の講義では、MPIの非同期通信（Non-blocking communication）とキャッシュ汚染の関係について、より深いレイヤへ踏み込むとしよう。現場からは以上だ。
コメント