From 6f3aeed89a76a6a128c5e5d8fd207ea00c2b0863 Mon Sep 17 00:00:00 2001 From: liangjg Date: Wed, 23 Jan 2019 16:53:31 -0500 Subject: [PATCH] fix test error --- docs/source/pythonapi/deplete.rst | 2 ++ openmc/deplete/integrator/leqi.py | 5 ++--- 2 files changed, 4 insertions(+), 3 deletions(-) diff --git a/docs/source/pythonapi/deplete.rst b/docs/source/pythonapi/deplete.rst index e38c00e96..7b24326d3 100644 --- a/docs/source/pythonapi/deplete.rst +++ b/docs/source/pythonapi/deplete.rst @@ -16,6 +16,8 @@ algorithms for depletion calculations. integrator.predictor integrator.cecm + integrator.celi + integrator.leqi integrator.cf4 integrator.epc_rk4 integrator.si_celi diff --git a/openmc/deplete/integrator/leqi.py b/openmc/deplete/integrator/leqi.py index 3c6d26367..04415dc21 100644 --- a/openmc/deplete/integrator/leqi.py +++ b/openmc/deplete/integrator/leqi.py @@ -117,12 +117,12 @@ def leqi(operator, timesteps, power=None, power_density=None, print_out=True): if i == 0: # Save results for the last step dt_l = dt - vec, t, op_res_last = celi_inner(operator, vec, p, i, i_res, + x_new, t, op_res_last = celi_inner(operator, vec, p, i, i_res, t, dt, print_out) continue # Perform remaining LE/QI - x = [copy.deepcopy(vec)] + x = [copy.deepcopy(x_new)] op_results = [operator(x[0], p)] inputs = list(zip(op_res_last.rates, op_results[0].rates, @@ -145,7 +145,6 @@ def leqi(operator, timesteps, power=None, power_density=None, print_out=True): Results.save(operator, x, op_results, [t, t+dt], p, i_res+i) # update results - x = [x_new] op_res_last = copy.deepcopy(op_results[0]) t += dt dt_l = dt