1

-O0 -g -traceback -fpe3デバッグ モード (フラグをオン)で Intel コンパイラ バージョン 13.1.3.192 を使用して巨大な Fortran コードを実行しています。次の出力メッセージが表示されました。

... ...
forrtl: warning (402): fort: (1): In call to MPI_ALLGATHER, an array temporary was created for argument #1

forrtl: error (65): floating invalid
Image              PC                Routine            Line        Source
arts               00000000016521D9  pentadiagonal_ser         385  pentadiagonal.f90
arts               0000000001644862  pentadiagonal_             62  pentadiagonal.f90
arts               00000000004DF167  implicit_solve_x_        1201  implicit.f90
arts               0000000000538079  scalar_bquick_inv         383  scalar_bquick.f90
arts               00000000004EFEAC  scalar_step_              190  scalar.f90
arts               0000000000401744  simulation_run_           182  simulation.f90
arts               0000000000401271  MAIN__                     10  main.f90
arts               0000000000400FAC  Unknown               Unknown  Unknown
arts               000000000420E444  Unknown               Unknown  Unknown
arts               0000000000400E81  Unknown               Unknown  Unknown

エラーの原因は、五対角行列を解くためのサブルーチン pentadiagonal_serial にあります。

subroutine pentadiagonal_serial(A,B,C,D,E,R,n,lot)
  use precision
  implicit none

  integer, intent(in) :: n,lot
  real(WP), dimension(lot,n) :: A     ! LOWER-2
  real(WP), dimension(lot,n) :: B     ! LOWER-1
  real(WP), dimension(lot,n) :: C     ! DIAGONAL
  real(WP), dimension(lot,n) :: D     ! UPPER+1
  real(WP), dimension(lot,n) :: E     ! UPPER+2
  real(WP), dimension(lot,n) :: R     ! RHS - RESULT
  real(WP), dimension(lot) :: const
  integer :: i

  if (n .eq. 1) then
    ! Solve 1x1 system
    R(:,1) = R(:,1)/C(:,1)
    return
  else if (n .eq. 2) then
    ! Solve 2x2 system
    const(:) = B(:,2)/C(:,1)
    C(:,2) = C(:,2) - D(:,1)*const(:)
    R(:,2) = R(:,2) - R(:,1)*const(:)
    R(:,2) = R(:,2)/C(:,2)
    R(:,1) = (R(:,1) - D(:,1)*R(:,2))/C(:,1)
    return
 end if

 ! Forward elimination
 do i=1,n-2
     ! Eliminate A(2,i+1)
     const(:) = B(:,i+1)/(C(:,i)+tiny(1.0_WP))
     C(:,i+1) = C(:,i+1) - D(:,i)*const(:)
     D(:,i+1) = D(:,i+1) - E(:,i)*const(:)
     R(:,i+1) = R(:,i+1) - R(:,i)*const(:)

その行で

  const(:) = B(:,i+1)/(C(:,i)+tiny(1.0_WP))

エラーが発生します。の値を出力しようとしましたが、値const(:)が存在することがわかりましたInfinity。しかし、なぜそれが無限大を生成できるのか理解できません。私が見る限り、分母がゼロになるのを避けるために、tiny(1.0_WP)が追加されC(:,i)、分母がゼロになることはほとんど不可能になりました...このサブルーチンが呼び出されたときに、すべてが初期化されるか、宣言後に値が与えられることも確認しました。そのため、どこが間違っているのかわかりませんでした。

4

1 に答える 1

1

（コメントの回答。回答のない質問を参照してください。ただし、コメントで問題が解決されています（またはチャットで拡張されています）。コメントには多くのチャットがあり、実際の回答を抽出するのは困難ですが、OPは解決済みを示しています. )

@SethMMorton は次のように書いています。

の値のいずれかCが前後-tinyにある場合でも、分母がゼロになる可能性があります。

の値を出力しましたCか? 小さな値は の無限大と相関していconstますか?

@Alexander Vogt は次のように書いています。

さて、無限大に何かを追加すると、浮動小数点例外が発生します...そして、ifort が無限大による除算がゼロであることを検出できるかどうかはわかりませんが、それは別の浮動小数点例外であると思います! C 内で無限値が発生しないようにする必要があります。

@ステファンは書いた：

回避策: ?gbsvMKL/Lapack のルーチンはどうですか? 次に、最後のヒント: でコンパイルしてみてくださいgfortran。コンパイラが異なれば、コンパイル時/実行時の警告やエラーも異なりgfortran、問題を早期に検出できる可能性があります。

于 2015-01-25T17:27:00.363 に答える