【Fortran学習｜豆知識】Fortranで堅牢なAPIを作る！PRESENT関数による引数省略のスマートな制御術

1. 導入：なぜ「引数省略」の判定が重要なのか

数値計算の現場で、似ているけれど少しだけ挙動が異なるサブルーチンをいくつも作っていませんか？例えば「境界条件ありの計算」と「なしの計算」を別々のルーチンにすると、コードの重複が増え、メンテナンスが困難になります。Fortranの PRESENT 関数を使うと、一つのサブルーチン内で引数の有無を判定し、必要な処理だけを実行する「スマートで堅牢なAPI」を設計できます。これにより、不要な計算を避けつつ、NULL参照のような予期せぬエラーを未然に防ぐことが可能になります。

2. 基礎知識：optional属性とは

Fortranにおいて、サブルーチンの引数に optional 属性を付与すると、呼び出し側はその引数を省略できるようになります。しかし、省略された引数をそのままサブルーチン内で使用しようとすると、プログラムはクラッシュします。ここで登場するのが PRESENT(引数名) です。これは、該当する引数が実際に呼び出し元から渡されているかどうかを論理型（.true. または .false.）で返す組み込み関数です。

3. 実装・解決策：条件分岐による柔軟な制御

実装の基本は、サブルーチンの冒頭で if (present(引数)) を使用し、引数が渡された場合のみその変数にアクセスする構造にすることです。これにより、引数が渡されなかった場合には、その処理自体をスキップしたり、あらかじめ決めておいたデフォルト値を適用したりといった柔軟な制御が可能になります。

4. サンプルプログラム

以下のコードは、重力計算の有無をオプションで切り替える例です。コピー＆ペーストして動作を確認してみてください。

subroutine compute_physics(velocity, gravity) implicit none ! velocityは必須引数、gravityは省略可能な引数 real, intent(inout) :: velocity real, intent(in), optional :: gravity real :: g_const


    ! 1. gravityが渡されているか判定

    if (present(gravity)) then

        g_const = gravity

        print , "重力あり: 加速度を設定します"

    else

        g_const = 0.0

        print , "重力なし: 重力の影響を無視します"

    end if
    ! 2. 判定結果に基づいた計算

    velocity = velocity + g_const

end subroutine
! --- 呼び出し側の例 ---

program test_optional

    real :: v = 10.0
    ! ケースA: gravityを渡す

    call compute_physics(v, 9.8)

! ケースB: gravityを省略する call compute_physics(v) end program

5. 応用・注意点：現場で役立つアドバイス

PRESENT を使う際、特に注意すべきは「省略された引数を、さらに別のサブルーチンへそのまま渡す」ケースです。引数が存在しない状態で別のサブルーチンに渡すとエラーになるため、必ず中継先のサブルーチンでも optional 属性が定義されていることを確認してください。また、引数にデフォルト値を持たせたい場合は、上記の例のように if (present()) ブロック内で変数を初期化する手法が最も安全でバグを誘発しにくい設計です。慣れてくると、一つのサブルーチンで多様な物理設定を切り替えられるようになり、コードの可読性が劇的に向上します。