導入
数値計算の現場、特に有限差分法（FDM）などで物理シミュレーションを行う際、計算領域の外側に「ゴーストセル（仮想境界）」を設けることは一般的です。しかし、配列のインデックスを 1 から開始するように設計すると、境界条件の適用時に `i-1` や `i+1` といった複雑なオフセット計算が頻発し、バグの温床となります。本稿では、Fortranの `ALLOCATE` 文で配列の下限（Lower Bound）を直接指定することで、論理座標と物理インデックスを一致させ、コードの可読性と堅牢性を劇的に向上させる手法を解説します。
基礎知識
Fortranにおける `ALLOCATE` は、実行時に必要なメモリサイズを確保する命令です。通常、`allocate(a(n))` と書くとインデックスは `1` から `n` までとなります。しかし、Fortranの配列宣言は `(下限:上限)` の形式をサポートしており、これを利用することで `0` から開始したり、あるいは負のインデックスを許容したりすることが可能です。これにより、物理的な座標系とプログラム上の添字を直感的にマッピングできます。
実装/解決策
解決策は極めてシンプルです。`ALLOCATE` 文の中で、各次元の範囲を `(下限:上限)` の形式で記述するだけです。これにより、配列の添字がそのまま物理的な位置関係を反映するようになります。例えば、左側に1点、右側に2点のゴーストセルが必要な場合、インデックスを `-1` から開始するように確保すれば、境界条件の計算式にオフセットを混ぜる必要がなくなり、コードが非常にクリーンになります。
サンプルプログラム
以下のコードは、ゴーストセルを含む2次元領域を動的に確保し、インデックスを直接操作する例です。
program ghost_cell_example
implicit none
integer :: nx, ny
! ポインタとして配列を定義（またはallocatable属性を使用）
real(8), allocatable :: q(:, 🙂
nx = 100
ny = 100
! 下限を-1、上限をnx+2と設定することで、
! q(-1, -1) から q(nx+2, ny+2) までの領域を確保する
allocate(q(-1:nx+2, -1:ny+2))
! 初期化例：境界以外を0にする
q = 0.0d0
! インデックスをそのまま物理座標として扱えるため、
! 境界条件の記述が直感的になる
! 例: 左端のゴーストセルへのアクセス
q(-1, 🙂 = 1.0d0
print , “領域の確保完了: 下限(-1,-1) 上限(“, nx+2, “,”, ny+2, “)”
deallocate(q)
end program ghost_cell_example
応用・注意点
この手法を用いる際に注意すべき点がいくつかあります。
1. lbound / ubound 関数の活用
サブルーチンに配列を渡す際、下限が 1 でないことを前提にコードを書くと汎用性が下がります。配列の範囲を問い合わせる `lbound(array, dim)` および `ubound(array, dim)` 関数を必ず使用し、下限が変動しても動作するように実装してください。
2. メモリレイアウトの考慮
Fortranは列優先（Column-major）言語です。多次元配列の下限を変更してもメモリ上の連続性は保たれますが、計算速度を重視する場合は、内側のループ（通常は第一添字）を回す際のアクセスパターンがキャッシュ効率に影響を与えないか確認してください。
3. デバッグ時の注意
`bounds-check` オプションを有効にしてコンパイルすることで、動的に確保した範囲外へのアクセスを即座に検知できます。開発中は必ずこのオプションを有効にしてください。
この手法を導入することで、境界条件の複雑なオフセット計算から解放され、物理モデルそのものの実装に集中できる環境を構築できます。ぜひ次回のプロジェクトから取り入れてみてください。
コメント