結果の合成(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倍 ⇒ 不確かさが減少.