結果の合成(sumtally)

  • モンテカルロで別々に実行した計算結果は加算できる.

設定

  • 実行時のファイルはそのまま流用すればよい.

  • 他の条件はすべて同じにする必要あり.

  • [parameters] セクションにて、icntrl = 13 を設定

  • sumtallyしたいtally設定ファイルをコピー (sumtally_phits.inpなどへ)

  • sumtallyしたいtallyセクション中に以下を追記.

sumtally start
isumtally = 2
nfile = 2
    sum/file_1.out 1.0
    sum/file_2.out 1.0
sfile = sum/file_s.out
sumfactor = 1.0
sumtally end
  • 通常通り、phitsを実行

実行例

  • 同じファイルを2個 1:1で読み込ませたとき.

    • 平均値に変化はない

    • 相対誤差が 0.707倍 (1/root(2) ) 倍になる.

    • 統計数が2倍 ⇒ 不確かさが減少.