Skip to content

add continue feature for depletion #3272

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
paulromano merged 17 commits into from
Mar 6, 2025
Merged

add continue feature for depletion #3272

Changes from 1 commit
Commits
Show all changes
17 commits
Select commit Hold shift + click to select a range
1730c21
add continue feature for depletion and test new code
lewisgross1296 Jan 22, 2025
5ed5a59
added flags to SIIntegrator, convereted np.equal to np.array_equal, u…
lewisgross1296 Jan 28, 2025
48f8478
implement get_source_rates
lewisgross1296 Jan 28, 2025
136289c
split failure cases into two tests and added match param
lewisgross1296 Jan 28, 2025
46ab8b4
source_rate does not have index (scalar variable)
lewisgross1296 Jan 28, 2025
4fdf333
fix off by 1
lewisgross1296 Feb 19, 2025
3f96f1a
split checks, need to compute dts, need timestep_units in get_times call
lewisgross1296 Feb 19, 2025
f3358e4
make sure that continue runs with no previous results do not break
lewisgross1296 Feb 24, 2025
09df7c1
fix indent issue, add some description
lewisgross1296 Mar 4, 2025
020aff2
addressing some of Patrick's feedback
lewisgross1296 Mar 4, 2025
00e575d
remove comment in favor of docstring
lewisgross1296 Mar 4, 2025
220f1a9
imrove docstring
lewisgross1296 Mar 4, 2025
1915e6d
move completed_source_rates declaration after success of timestep check
Mar 5, 2025
384e99c
Style changes to test file
pshriwise Mar 5, 2025
471e5b6
Revert change to pugixml submodule
paulromano Mar 5, 2025
b27c9df
Small stylistic fixes
paulromano Mar 5, 2025
e2d0ad1
Simplify tests for continue run
paulromano Mar 5, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
added flags to SIIntegrator, convereted np.equal to np.array_equal, u… 
…pdated doc strings
  • Loading branch information
@lewisgross1296
lewisgross1296 committed Mar 4, 2025
commit 5ed5a59107349fefa2038ac555ceb2502a43e369
28 changes: 26 additions & 2 deletions openmc/deplete/abc.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -572,6 +572,18 @@ class Integrator(ABC):
:attr:`solver`.

.. versionadded:: 0.12
continue_timesteps : bool, optional
Whether or not to treat the current solve as a continuation of a
previous simulation. Defaults to `False`. If `True`, the timesteps

@pshriwise pshriwise Mar 4, 2025

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Let's explain the intended behavior for a False value as well here.

provided to the `Integrator` must match exactly those that exist
in the `prev_results` passed to the `Opereator`. The `power`,
Comment thread
lewisgross1296 marked this conversation as resolved.
Resolved
`power_density`, or `source_rates` must match as well. It
is the user's responsibility to make sure that the continue
solve uses the same method of specifying `power`, `power_density`,
or `source_rates`.
Comment thread
lewisgross1296 marked this conversation as resolved.
Resolved

.. versionadded:: 0.15.1

Attributes
----------
operator : openmc.deplete.abc.TransportOperator
Expand Down Expand Up @@ -933,6 +945,17 @@ class SIIntegrator(Integrator):
:attr:`solver`.

.. versionadded:: 0.12
continue_timesteps : bool, optional
Whether or not to treat the current solve as a continuation of a
previous simulation. Defaults to `False`. If `True`, the timesteps
provided to the `Integrator` must match exactly those that exist
in the `prev_results` passed to the `Opereator`. The `power`,
`power_density`, or `source_rates` must match as well. It
is the user's responsibility to make sure that the continue
solve uses the same method of specifying `power`, `power_density`,
or `source_rates`.

.. versionadded:: 0.15.1

Attributes
----------
Expand Down Expand Up @@ -974,13 +997,14 @@ def __init__(
source_rates: Optional[Sequence[float]] = None,
timestep_units: str = 's',
n_steps: int = 10,
solver: str = "cram48"
solver: str = "cram48",
continue_timesteps: bool = False,
):
check_type("n_steps", n_steps, Integral)
check_greater_than("n_steps", n_steps, 0)
super().__init__(
operator, timesteps, power, power_density, source_rates,
timestep_units=timestep_units, solver=solver)
timestep_units=timestep_units, solver=solver, continue_timesteps=continue_timesteps)
self.n_steps = n_steps

def _get_bos_data_from_operator(self, step_index, step_power, n_bos):
Expand Down