【入門編】手続きポインタ(PROCEDURE POINTER)と抽象インターフェース – モダンFortran言語仕様と実践実践マスター

こんにちは。長年、数値計算の最前線でスパコンのノードを焼き切る寸前まで追い込み、Fortranの冷徹なまでの最適化性能と付き合ってきた者です。

CやPythonから来た皆さんがFortranに触れるとき、「なぜこんなに古い言語で、関数ポインタのような抽象的な概念を使うのか？」と疑問に思うかもしれません。答えはシンプルです。「究極のパフォーマンスと、エンジニアリングの柔軟性を両立させるため」です。

今日は、Fortranの「手続きポインタ（Procedure Pointer）」と「抽象インターフェース（Abstract Interface）」を使って、コンパイル時の型安全性を維持しながら、実行時にアルゴリズムを華麗に切り替える技を伝授します。

—

C言語の関数ポインタは強力ですが、型チェックがザルになりがちで、セグメンテーションフォールトの温床になります。一方、Fortranの手続きポインタは、「インターフェース」という契約書を事前に交わすことで、型安全性を担保します。

まずは「抽象インターフェース」を定義しましょう。これは「こういう引数と戻り値を持つ関数であれば、中身が何であっても受け入れるよ」という約束事です。

! 抽象インターフェースの定義：関数の設計図です
abstract interface
function solver_func(x) result(y)
real, intent(in) :: x
real :: y
end function solver_func
end interface

これで `solver_func` という型が完成しました。これに適合する関数なら、後からいくらでもすり替えることが可能です。

—

次に、このインターフェースを利用するポインタ変数を宣言し、実行時に切り替えるコードを書いてみましょう。

module solver_module
implicit none

! 抽象インターフェース（再掲）
abstract interface
function op_type(x) result(y)
real, intent(in) :: x
real :: y
end function op_type
end interface

contains
! 具体的なアルゴリズムA
function linear_op(x) result(y)
real, intent(in) :: x
y = 2.0 x
end function linear_op

! 具体的なアルゴリズムB
function quadratic_op(x) result(y)
real, intent(in) :: x
y = x2
end function quadratic_op
end module solver_module

program main
use solver_module
implicit none

! 手続きポインタの宣言
procedure(op_type), pointer :: my_solver => null()
real :: input = 5.0

! ここで動的に切り替える！
my_solver => linear_op
print , “Linear: “, my_solver(input)

my_solver => quadratic_op
print , “Quadratic: “, my_solver(input)
end program main

—

なぜわざわざ `abstract interface` を挟むのでしょうか？

C++の仮想関数テーブル（vtable）のような動的ディスパッチは、キャッシュミスを誘発しやすく、超大規模計算では微小なオーバーヘッドが積もり積もって数時間のロスになります。

Fortranの手続きポインタは、コンパイラにとって「このポインタが指す先は、このインターフェースを守っている」と保証されているため、インライン展開などの最適化が効きやすい構造になっています。

プロシージャの属性: 渡す関数には `pure` や `elemental` を付与することを検討してください。これにより、コンパイラは「この関数は副作用がない」と判断し、ベクトル演算（SIMD）の最適化を積極的に適用してくれます。
列優先順位の意識: 関数内で配列を扱う場合、Fortranは列優先（Column-major）です。ポインタ経由で関数を呼び出す際も、渡す配列のメモリレイアウトがアクセスパターンに合致しているか、常に意識してください。

—

まずは、皆さんの手元のコードで、if文の分岐が多発している箇所を探してみてください。

! こういう書き方から…
if (mode == 1) then
y = calc_a(x)
else
y = calc_b(x)
end if

これを手続きポインタを使って、「初期化時にポインタをセットし、ループ内ではポインタを叩くだけ」の設計に書き換えてみてください。コードの可読性が上がるだけでなく、条件分岐の予測失敗（Branch Misprediction）が減り、驚くほど計算速度が改善することがあります。

Fortranは古い言語ですが、その設計思想は極めてモダンです。手続きポインタを使いこなせれば、大規模な数値シミュレーションの基盤を、柔軟かつ堅牢に構築できるはずです。

さあ、皆さんもコンパイラを回して、この「Fortran流の柔軟性」をその目で確かめてみてください。応援しています！