Fortranの「名前の再利用術」：MODULE PROCEDUREでコードを美しく磨き上げる
皆さん、こんにちは。宇宙航空の現場で数値計算の最適化を追い求めてきた身として、これからFortranの世界に足を踏み入れる皆さんに一つだけ伝えたいことがあります。
「Fortranは古臭い」なんていう先入観は、今日この瞬間に捨ててください。今のFortranは、大規模な並列計算を支える非常に強力で、かつ人間味のある言語です。今日は、皆さんがC言語やPythonで慣れ親しんだ「オーバーロード（多重定義）」を、Fortranでいかに美しく、かつ爆速に実装するか、その極意を伝授します。
—
なぜ「MODULE PROCEDURE」なのか？
C++なら`overload`、Pythonなら`functools.singledispatch`など、引数の型や形状に合わせて処理を切り替える仕組みがありますよね。Fortranでは`INTERFACE`ブロックと`MODULE PROCEDURE`を組み合わせることで、これと同じ……いや、Fortranのメモリ配置を理解した者にとっては、それ以上のパフォーマンスを発揮する構造が作れます。
例えば、「ベクトルの和」を取る関数を作るとき、`add_int`、`add_real`、`add_double`……と関数名を変えていたら、コードはすぐにスパゲッティ状態になります。これを`add`という一つの名前で呼び出せたら、どんなに幸せでしょうか。
—
ステップ1：インターフェースという「窓口」を作る
まずは、本体となる関数を`MODULE`の中に閉じ込めます。そして、外部から呼び出すための「窓口（INTERFACE）」を用意しましょう。これが、Fortranにおけるオーバーロードの基本形です。
module arithmetic_ops
implicit none
private
public :: add ! 外部から見えるのは’add’という名前だけにする
! 窓口となるインターフェース定義
interface add
module procedure add_real_scalar, add_real_array
end interface
contains
! 実装1：スカラ値の加算
function add_real_scalar(a, b) result(res)
real, intent(in) :: a, b
real :: res
res = a + b
end function
! 実装2：配列の加算（Fortranの真骨頂！）
function add_real_array(a, b) result(res)
real, intent(in) :: a(:), b(:)
real :: res(size(a))
! ここでループを書く必要はありません！
! Fortranの配列演算は、コンパイラがSIMD命令をフル活用して展開します
res = a + b
end function
end module
なぜこれが「現場の武器」になるのか
ここで鋭い皆さんは、「これ、内部で分岐させているだけじゃないの？」と思うかもしれません。しかし、Fortranのインターフェース解決はコンパイル時に行われます。
C言語の`void`を使った汎用関数や、Pythonの動的型付けと違い、Fortranはコンパイラが「どの関数を呼ぶか」を確定させてからバイナリを生成します。つまり、実行時のオーバーヘッドがゼロなのです。
さらに、Fortranは「列優先（Column-major）」というメモリ配置を持っています。配列を処理する際、メモリの連続性を意識して書けば、キャッシュミスが激減します。`MODULE PROCEDURE`で型ごとに最適化された処理を分けておけば、コンパイラは「この型ならこの最適化ができる」という判断を迷わず下せるわけです。
—
実践：ビルド時の最適化フラグを忘れずに
どれだけコードを綺麗に書いても、コンパイル設定で足を引っ張っては意味がありません。例えば、Intel Fortran (ifort/ifx) を使う場合、以下のようなフラグを意識してみてください。
現場でよく使う最適化の基本セット
-O3: 最強の最適化
-xHost: 今使っているCPUの命令セットをフル活用する
-ipo: モジュール間最適化（MODULEを跨いでインライン展開を狙う）
ifx -O3 -xHost -ipo -o simulation main.f90
`-ipo`（Interprocedural Optimization）を付与することで、`add`というインターフェースを介して呼ばれた先の実装関数が、呼び出し元のコードに完全にインライン展開されます。結果として、抽象化の代償である「関数呼び出しコスト」すら消滅するのです。
—
さあ、次はあなたの番です
最初は「`MODULE`の定義」や「`INTERFACE`の構文」に戸惑うかもしれません。しかし、一度この「型ごとに窓口を分ける」という書き方に慣れると、複雑な物理シミュレーションのコードが見違えるように整理されます。
まずは、皆さんが今書いているコードの関数を一つだけ、`MODULE`の中に移動させてみてください。そして、`INTERFACE`という名の「窓口」を作ってみる。それだけで、あなたの書くコードは「ただ動くだけのプログラム」から「プロの計算機科学者が書くソフトウェア」へと一歩進化します。
分からないことがあれば、いつでも聞いてください。Fortranという深淵で、皆さんと肩を並べてデバッグできる日を楽しみにしていますよ！
コメント