1. 導入：なぜ配列の「線形化」が必要なのか
数値計算の現場では、巨大な行列を一括で初期化したり、ファイルへ高速に書き出したりする場面が多々あります。Fortranのレガシーな機能である「EQUIVALENCE」は、メモリ上の同じ領域を複数の変数名で参照させるための仕組みです。かつては、2次元配列を1次元配列として扱うことで、ループ回数を減らしたり、メモリ操作を高速化したりするために重宝されていました。本記事では、この仕組みを理解し、現代のプログラミングにおいてどのように安全に書き換えるべきかを解説します。
2. 基礎知識：メモリの「列優先」と線形化の仕組み
Fortranの配列は、メモリ上で「列優先（Column-Major）」というルールで配置されています。例えば 2×2 の配列 A(2,2) がある場合、メモリ上では A(1,1), A(2,1), A(1,2), A(2,2) の順に並びます。
EQUIVALENCE文を使うと、このメモリ領域の先頭を「1次元配列 B(4)」と重ね合わせることができます。これにより、行列のすべての要素を、インデックス計算を介さずに一気に操作することが可能になります。
3. 実装：EQUIVALENCEの具体的な使い方
EQUIVALENCE (A(1,1), B(1)) と記述することで、2次元配列 A と 1次元配列 B が同じメモリ領域を共有します。これにより、Aの成分を操作すると、自動的にBの成分も書き換わるという性質を利用します。
4. サンプルプログラム：EQUIVALENCEによる線形化
以下は、2次元配列を1次元配列として一括初期化するコード例です。
program equivalence_example
implicit none
! 2次元配列と、それを線形化するための1次元配列を定義
integer, parameter :: NX = 2, NY = 3
integer :: A(NX, NY)
integer :: B(NX NY)
! EQUIVALENCEによりメモリ領域を共有させる
equivalence (A(1, 1), B(1))
! Bを操作することで、間接的にAを初期化する
! 全要素を0にセット（高速な一括処理の例）
B = 0
! 確認：Aの値が表示される
print , "A(1,1) =", A(1,1)
print , "A(2,3) =", A(2,3)
end program equivalence_example
5. 応用・注意点：現代的なアプローチへの移行
EQUIVALENCEは非常に強力ですが、メモリレイアウトの前提が崩れると意図しないバグ（セグメンテーション違反など）を引き起こすリスクがあります。
現代的な解決策：
最新のFortran（Fortran 90以降）では、EQUIVALENCEに頼る必要はありません。
1. RESHAPE関数： 配列の形状を再定義して別の配列へコピーする際に便利です。
2. ポインタ（POINTER）： ターゲットとなる配列を指し示すことで、メモリを共有したまま安全に形状を変換できます。
3. 型付きメモリコピー： C言語とのインターフェース（ISO_C_BINDING）を利用して、より厳密なメモリ操作を行うのが標準的です。
古いコードを保守する際はEQUIVALENCEの仕組みを理解しておくことは必須ですが、新規にコードを書く際は、コンパイラが最適化を行いやすい「モダンな書き方」を選択するようにしましょう。
コメント