浮動小数点例外（NaN）の検出と捕捉

Fortranのレガシーコードや数値計算プログラムにおいて、NaN（非数）の出力は計算の異常を示唆する重要な手がかりです。NaNが発生した時点でプログラムを停止し、コアダンプを生成することで原因箇所を特定できます。Intelコンパイラ（ifort）とGNUコンパイラ（gfortran）の両方において、IEEE算術標準モジュールを使用した検出手法と、適切なコンパイルフラグの設定が重要です。

program ieee_except_demo
  use, intrinsic :: ieee_arithmetic
  use, intrinsic :: iso_fortran_env
  implicit none

  real(real64) :: variable_a, variable_b, result_val

  ! 変数をQuiet NaNで初期化
  variable_a = ieee_value(variable_a, ieee_quiet_nan)
  variable_b = 42.0_real64

  ! NaNを含む演算の実行
  result_val = variable_a + variable_b
  print *, "Addition: ", result_val

  result_val = variable_b / variable_a
  print *, "Division: ", result_val

  result_val = variable_a * 0.0_real64
  print *, "Multiplication: ", result_val

end program ieee_except_demo

以下に、浮動小数点例外を捕捉するための主要なコンパイルオプションを示します。gfortranでは-ffpe-trapを、ifortでは-fpe0を使用することで、無効演算やゼロ除算時に即座にトラップを発生させることが可能です。

• gfortran: -ffpe-trap=invalid,zero,overflow
• ifort: -fpe0 (全ての浮動小数点例外をトラップ)

デバッグ支援のためのコンパイラオプション比較

バグの早期発見には、コンパイラの診断機能を最大限に活用することが不可欠です。以下の表は、配列境界のチェックや未初期化変数の検出など、主要なデバッグオプションの比較です。

メモリリークの診断と監視

長時間実行されるFortranアプリケーションにおいて、メモリ使用量が単調増加し続ける場合、メモリリークの可能性が高いです。Fortranではallocatable属性を持つ配列はスコープを抜けると自動的に解放されますが、pointer属性を使用している場合は手動でdeallocateを行う必要があります。特に派生型（derived type）内にポインタメンバが存在する場合は、ファイナライザ（finalプロシージャ）を適切に定義し、解放漏れを防ぐ必要があります。

メモリの推移を監視するには、/procファイルシステムを利用したシェルスクリプトが有効です。

#!/bin/bash
# メモリ監視スクリプト (使用例: ./monitor_mem.sh <PID>)
TARGET_PID=$1
LOG_FILE="memory_usage.log"

while true; do
    if ps -p $TARGET_PID > /dev/null; then
        # VmSize (仮想メモリ) と VmRSS (物理メモリ) を取得
        VmSize=$(grep VmSize /proc/$TARGET_PID/status | awk '{print $2}')
        VmRSS=$(grep VmRSS /proc/$TARGET_PID/status | awk '{print $2}')
        echo "$(date +%s) $VmSize $VmRSS" >> $LOG_FILE
        sleep 1
    else
        break
    fi
done

アサーションと事前・事後条件のチェック

複雑な数値計算ルーチンでは、プログラマの意図しない状態に入ることを防ぐためにアサーション（表明）が有効です。計算の各フェーズにおいて、事前条件（入力パラメータの妥当性）と事後条件（計算結果の整合性）を検証するコードを挿入することで、エラーの発生箇所を迅速に特定できます。

module assert_utils
  implicit none
contains
  subroutine check_condition(is_valid, error_msg)
    logical, intent(in) :: is_valid
    character(len=*), intent(in) :: error_msg
    if (.not. is_valid) then
      print *, "Runtime Assertion Failed: ", trim(error_msg)
      error stop 1
    end if
  end subroutine check_condition
end module assert_utils

subroutine compute_divide(a, b, result)
  use assert_utils
  real, intent(in) :: a, b
  real, intent(out) :: result

  ! 事前条件チェック：ゼロ除算の防止
  call check_condition(abs(b) > 1.0e-12, "Divisor is too close to zero")

  result = a / b

  ! 事後条件チェック：結果がNaNでないこと
  call check_condition(result == result, "Result is NaN")
end subroutine compute_divide

整数オーバーフローとスタックオーバーフロー対策

大規模な解析（有限要素法など）において、整数オーバーフローは深刻な問題を引き起こします。特に配列のインデックス計算において、デフォルトの32ビット整数（integer(4)）では約20億（2GB）の要素しか扱えません。これを超える配列サイズを扱う場合は、integer(int64)やinteger(8)を明示的に使用する必要があります。インデックス計算の結果が負の値になる現象は、典型的なオーバーフローの症状です。

また、スタックオーバーフローは、関数内で大きなサイズの自動配列（静的配列）を宣言した際に発生します。この問題を回避するには、スタックサイズの制限（ulimit -s）を増やすか、大規模な配列をヒープ領域に確保するためにallocatable属性を使用する必要があります。

subroutine process_large_data(n_elements)
  use, intrinsic :: iso_fortran_env, only: int64
  integer, intent(in) :: n_elements
  ! 大規模配列を扱うためには64ビット整数を使用
  integer(int64) :: i, idx
  real, allocatable, dimension(:) :: large_array

  ! ヒープ領域に動的確保することでスタックオーバーフローを防止
  allocate(large_array(n_elements))

  do i = 1, int(n_elements, kind=i)
    idx = i * 2_int64
    ! 計算処理...
  end do

  deallocate(large_array)
end subroutine process_large_data