############################################################## 計算の再開方法について ############################################################## * PHITSでは、計算途中に、処理を中断しても、途中から :blue:`計算を再開` できる. * 1000粒子を100回投げる場合に、50回投げた後で、中断しても、51回目を実行して今までためてきた統計を鑑みれば良い. ========================================================= 計算結果の記録 ========================================================= * 計算開始後、 :blue:`出力ファイルが生成・逐次記録` が始まる ( 各タリー出力 ). * 残りバッチ数、次回再開の乱数シードなどが、 :blue:`batch.out` に記録・更新されていく. ========================================================= 計算の中止 ========================================================= * 正しい中止: batch.outの :red:`残りバッチ数を "0" に書き換える` こと. + おそらく、Ctrl+C で途中中止しても問題はない. + 記載が中途半端になったりする可能性もあり得るので、上記正しい手順でやるべき. * batch.out, 各タリー出力に、中断しても、そこまでの出力ファイルが残る. ========================================================= 計算の再開(設定) ========================================================= * istdev の値を :red:`負(-1/-2)` に設定.( -1: バッチ分散モード、-2: 全ヒストリ分散モード ) - 通常は、 :blue:`istdev=-2` でよい. * 各 tally の :blue:`resfile` を指定する. - resfileのデフォルトはfileとマニュアルにあるが、設定しないと、「resfileがない」、とおこられ、istdev=+2にされる. * maxbchで、指定したバッチ回数、計算し、それまでにタリーして記載されていた粒子と加算する. - (途中やめ、途中復活モードではない. :blue:`「5/10 でやめたら、6/10からはじまる」ではない` ことに注意) * 統計誤差(r.err)を観察すると、たしかに、くり返すごとに減っていく.( :blue:`統計がたまっていく` ) ========================================================= 計算の再開(実行) ========================================================= * 通常通り、実行する.( $ go__phits.py inp/main_phits.inp etc. )