導入: なぜインデックスの再定義が重要なのか
数値計算の現場では、物理モデルに基づいた配列のインデックスを 0 から開始したい場合と、計算ライブラリや既存コードの仕様で 1 から開始したい場合が混在します。通常、データ構造が異なると「新しい配列を作成して値をコピーする」という手間が発生し、メモリ消費量と計算速度の両面で損失が生じます。Fortranの「想定形状配列（Assumed-shape array）」の引数宣言を活用すれば、メモリコピーを一切行わずに、呼び出し側の配列を「別のインデックス体系」として解釈させることが可能です。
基礎知識: 想定形状配列とは
Fortranにおける想定形状配列とは、サブルーチンの引数宣言において配列のサイズを明示せず、コロン（:）を用いる手法です。通常は `real :: a(:)` のように宣言しますが、ここに `0:` や `100:` といった下限値を指定することで、そのサブルーチン内でのみ「配列の開始番号」を強制的に再定義できます。これは物理シミュレーションにおける境界条件の管理や、数学的な添字とプログラム上の添字を一致させたい際に非常に強力な機能です。
実装/解決策: コピーなしの座標変換
この手法の核は、呼び出し側のメモリ配置をそのまま利用し、ポインタの解釈だけを書き換える点にあります。例えば、1から10まで格納された配列 `A(1:10)` を `sub(A)` と渡す際、引数側で `v(0:)` と受け取れば、サブルーチン内での `v(0)` は呼び出し側の `A(1)` を指すことになります。これにより、物理モデル上の「0点」をコード上で直感的に扱えるようになります。
サンプルプログラム
以下のコードは、1から始まる配列を、0から始まるインデックス体系として受け取り、計算を行う例です。
program main
implicit none
real :: data(1:5) = [1.0, 2.0, 3.0, 4.0, 5.0]
! 通常の呼び出しだが、サブルーチン側でインデックスが再定義される
call process_array(data)
end program main
subroutine process_array(v)
! v(0:) と指定することで、サブルーチン内では0番目から開始される
real, intent(inout) :: v(0:)
integer :: i
! v(0)は呼び出し側のdata(1)を指す
! v(4)は呼び出し側のdata(5)を指す
do i = 0, 4
print , "インデックス", i, "の値:", v(i)
end do
end subroutine process_array
応用・注意点: 現場で役立つ回避策
この手法を利用する際に注意すべき点は、配列の境界外参照（Out of Bounds）です。例えば、呼び出し側が `A(1:10)` で、サブルーチン側が `v(0:)` と定義した場合、`v(10)` にアクセスすると、実際のメモリ領域を超えてアクセスしてしまい、セグメンテーションフォールトを引き起こします。
また、コンパイラの設定によっては、想定形状配列の受け渡し時に一時的なメモリ領域（テンポラリ）が生成される場合があります。これを避けるためには、必ずインターフェースブロックを明示するか、モジュール内にサブルーチンを記述して、明示的なインターフェースを用意してください。これにより、コンパイラは「コピーなしの参照渡し」を正しく最適化できるようになります。大規模な行列計算や流体シミュレーションにおいて、このメモリ節約術は非常に有効ですので、ぜひ活用してください。
コメント