From 931bc405cf357f3e647528f0150c9658f683828f Mon Sep 17 00:00:00 2001 From: Jingang Liang Date: Sat, 24 Aug 2019 08:40:22 +0000 Subject: [PATCH] use finer energy piece decomposition for faster fitting --- openmc/data/multipole.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/openmc/data/multipole.py b/openmc/data/multipole.py index f7a09d0799..c5c35a32fc 100644 --- a/openmc/data/multipole.py +++ b/openmc/data/multipole.py @@ -469,7 +469,7 @@ def _vectfit_nuclide(endf_file, njoy_error=5e-4, vf_error=1e-3, vf_pieces=None, mts = [2, 27] if log: - print(" MTs: {})".format(mts)) + print(" MTs: {}".format(mts)) print(" Energy range: {:.3e} to {:.3e} eV ({} points)".format( E_min, E_max, n_points)) @@ -480,8 +480,8 @@ def _vectfit_nuclide(endf_file, njoy_error=5e-4, vf_error=1e-3, vf_pieces=None, # divide into pieces for complex nuclides peaks, _ = find_peaks(total_xs) n_peaks = peaks.size - if n_peaks > 300 or n_points > 50000 or n_peaks * n_points > 100*30000: - vf_pieces = max(5, n_peaks // 80, n_points // 5000) + if n_peaks > 200 or n_points > 30000 or n_peaks * n_points > 100*10000: + vf_pieces = max(5, n_peaks // 50, n_points // 2000) else: vf_pieces = 1 piece_width = (sqrt(E_max) - sqrt(E_min)) / vf_pieces