21887062786

Committed 10 Feb 2026 11:55PM UTC coverage: 81.899% (-0.1%) from 82.009%

Build # 21887062786

Build Type

Pull #3493

github

Committed by

web-flow

Commit Message

Merge 320b68711 into 3f20a5e22

Pull Request Pull Request #3493: Implement vector fitting to replace external `vectfit` package

Run Details

17361 of 24292 branches covered (71.47%)

Branch coverage included in aggregate %.

185 of 218 new or added lines in 3 files covered. (84.86%)

2770 existing lines in 69 files now uncovered.

56369 of 65733 relevant lines covered (85.75%)

51144917.16 hits per line

Source File
Press 'n' to go to next uncovered line, 'b' for previous

98.79

/src/finalize.cpp

#include "openmc/finalize.h"

#include "openmc/bank.h"
#include "openmc/capi.h"
#include "openmc/cmfd_solver.h"
#include "openmc/collision_track.h"
#include "openmc/constants.h"
#include "openmc/cross_sections.h"
#include "openmc/dagmc.h"
#include "openmc/eigenvalue.h"
#include "openmc/event.h"
#include "openmc/geometry.h"
#include "openmc/geometry_aux.h"
#include "openmc/material.h"
#include "openmc/mesh.h"
#include "openmc/message_passing.h"
#include "openmc/mgxs_interface.h"
#include "openmc/nuclide.h"
#include "openmc/photon.h"
#include "openmc/plot.h"
#include "openmc/random_lcg.h"
#include "openmc/random_ray/random_ray_simulation.h"
#include "openmc/settings.h"
#include "openmc/simulation.h"
#include "openmc/source.h"
#include "openmc/surface.h"
#include "openmc/tallies/tally.h"
#include "openmc/thermal.h"
#include "openmc/timer.h"
#include "openmc/volume_calc.h"
#include "openmc/weight_windows.h"

#include "xtensor/xview.hpp"

namespace openmc {

void free_memory()
{
  free_memory_geometry();
  free_memory_surfaces();
  free_memory_material();
  free_memory_volume();
  free_memory_simulation();
  free_memory_photon();
  free_memory_settings();
  free_memory_thermal();
  library_clear();
  nuclides_clear();
  free_memory_source();
  free_memory_mesh();
  free_memory_tally();
  free_memory_bank();
  free_memory_plot();
  free_memory_weight_windows();
  if (mpi::master) {
    free_memory_cmfd();
  }
  if (settings::event_based) {
    free_event_queues();
  }
}

} // namespace openmc

using namespace openmc;

int openmc_finalize()
{
  if (simulation::initialized)
    openmc_simulation_finalize();

  // Clear results
  openmc_reset();

  // Reset timers
  reset_timers();

  // Reset global variables
  settings::assume_separate = false;
  settings::check_overlaps = false;
  settings::collision_track_config = CollisionTrackConfig {};
  settings::confidence_intervals = false;
  settings::create_fission_neutrons = true;
  settings::create_delayed_neutrons = true;
  settings::electron_treatment = ElectronTreatment::LED;
  settings::delayed_photon_scaling = true;
  settings::energy_cutoff = {0.0, 1000.0, 0.0, 0.0};
  settings::time_cutoff = {INFTY, INFTY, INFTY, INFTY};
  settings::entropy_on = false;
  settings::event_based = false;
  settings::free_gas_threshold = 400.0;
  settings::gen_per_batch = 1;
  settings::legendre_to_tabular = true;
  settings::legendre_to_tabular_points = -1;
  settings::material_cell_offsets = true;
  settings::max_lost_particles = 10;
  settings::max_order = 0;
  settings::max_particles_in_flight = 100000;
  settings::max_secondaries = 10000;
  settings::max_particle_events = 1'000'000;
  settings::max_history_splits = 10'000'000;
  settings::max_tracks = 1000;
  settings::max_write_lost_particles = -1;
  settings::n_log_bins = 8000;
  settings::n_inactive = 0;
  settings::n_particles = -1;
  settings::output_summary = true;
  settings::output_tallies = true;
  settings::particle_restart_run = false;
  settings::path_cross_sections.clear();
  settings::path_input.clear();
  settings::path_output.clear();
  settings::path_particle_restart.clear();
  settings::path_sourcepoint.clear();
  settings::path_statepoint.clear();
  settings::photon_transport = false;
  settings::reduce_tallies = true;
  settings::rel_max_lost_particles = 1.0e-6;
  settings::res_scat_on = false;
  settings::res_scat_method = ResScatMethod::rvs;
  settings::res_scat_energy_min = 0.01;
  settings::res_scat_energy_max = 1000.0;
  settings::restart_run = false;
  settings::run_CE = true;
  settings::run_mode = RunMode::UNSET;
  settings::solver_type = SolverType::MONTE_CARLO;
  settings::source_latest = false;
  settings::source_rejection_fraction = 0.05;
  settings::source_separate = false;
  settings::source_write = true;
  settings::ssw_cell_id = C_NONE;
  settings::ssw_cell_type = SSWCellType::None;
  settings::ssw_max_particles = 0;
  settings::ssw_max_files = 1;
  settings::survival_biasing = false;
  settings::temperature_default = 293.6;
  settings::temperature_method = TemperatureMethod::NEAREST;
  settings::temperature_multipole = false;
  settings::temperature_range = {0.0, 0.0};
  settings::temperature_tolerance = 10.0;
  settings::trigger_on = false;
  settings::trigger_predict = false;
  settings::trigger_batch_interval = 1;
  settings::uniform_source_sampling = false;
  settings::ufs_on = false;
  settings::urr_ptables_on = true;
  settings::verbosity = -1;
  settings::weight_cutoff = 0.25;
  settings::weight_survive = 1.0;
  settings::weight_windows_file.clear();
  settings::weight_windows_on = false;
  settings::write_all_tracks = false;
  settings::write_initial_source = false;

  simulation::keff = 1.0;
  simulation::need_depletion_rx = false;
  simulation::ssw_current_file = 1;
  simulation::total_gen = 0;

  simulation::entropy_mesh = nullptr;
  simulation::ufs_mesh = nullptr;

  data::energy_max = {INFTY, INFTY, INFTY, INFTY};
  data::energy_min = {0.0, 0.0, 0.0, 0.0};
  data::temperature_min = 0.0;
  data::temperature_max = INFTY;
  data::mg = {};
  model::root_universe = -1;
  model::plotter_seed = 1;
  openmc::openmc_set_seed(DEFAULT_SEED);
  openmc::openmc_set_stride(DEFAULT_STRIDE);

  // Deallocate arrays
  free_memory();

#ifdef OPENMC_LIBMESH_ENABLED
  settings::libmesh_init.reset();
#endif

  // Free all MPI types
#ifdef OPENMC_MPI
  if (mpi::source_site != MPI_DATATYPE_NULL) {
    MPI_Type_free(&mpi::source_site);
  }
  if (mpi::collision_track_site != MPI_DATATYPE_NULL) {
    MPI_Type_free(&mpi::collision_track_site);
  }
#endif

  openmc_finalize_random_ray();

  return 0;
}

int openmc_reset()
{

  model::universe_level_counts.clear();

  for (auto& t : model::tallies) {
    t->reset();
  }

  // Reset global tallies
  simulation::n_realizations = 0;
  xt::view(simulation::global_tallies, xt::all()) = 0.0;

  simulation::k_col_abs = 0.0;
  simulation::k_col_tra = 0.0;
  simulation::k_abs_tra = 0.0;
  simulation::k_sum = {0.0, 0.0};
  simulation::satisfy_triggers = false;

  settings::cmfd_run = false;

  simulation::n_lost_particles = 0;

  return 0;
}

int openmc_reset_timers()
{
  reset_timers();
  return 0;
}

int openmc_hard_reset()
{
  // Reset all tallies and timers
  openmc_reset();
  reset_timers();

  // Reset total generations and keff guess
  simulation::keff = 1.0;
  simulation::total_gen = 0;

  // Reset the random number generator state
  openmc::openmc_set_seed(DEFAULT_SEED);
  openmc::openmc_set_stride(DEFAULT_STRIDE);
  return 0;
}

1	#include "openmc/finalize.h"
2
3	#include "openmc/bank.h"
4	#include "openmc/capi.h"
5	#include "openmc/cmfd_solver.h"
6	#include "openmc/collision_track.h"
7	#include "openmc/constants.h"
8	#include "openmc/cross_sections.h"
9	#include "openmc/dagmc.h"
10	#include "openmc/eigenvalue.h"
11	#include "openmc/event.h"
12	#include "openmc/geometry.h"
13	#include "openmc/geometry_aux.h"
14	#include "openmc/material.h"
15	#include "openmc/mesh.h"
16	#include "openmc/message_passing.h"
17	#include "openmc/mgxs_interface.h"
18	#include "openmc/nuclide.h"
19	#include "openmc/photon.h"
20	#include "openmc/plot.h"
21	#include "openmc/random_lcg.h"
22	#include "openmc/random_ray/random_ray_simulation.h"
23	#include "openmc/settings.h"
24	#include "openmc/simulation.h"
25	#include "openmc/source.h"
26	#include "openmc/surface.h"
27	#include "openmc/tallies/tally.h"
28	#include "openmc/thermal.h"
29	#include "openmc/timer.h"
30	#include "openmc/volume_calc.h"
31	#include "openmc/weight_windows.h"
32
33	#include "xtensor/xview.hpp"
34
35	namespace openmc {
36
37	void free_memory()	7,505✔
38	{
39	free_memory_geometry();	7,505✔
40	free_memory_surfaces();	7,505✔
41	free_memory_material();	7,505✔
42	free_memory_volume();	7,505✔
43	free_memory_simulation();	7,505✔
44	free_memory_photon();	7,505✔
45	free_memory_settings();	7,505✔
46	free_memory_thermal();	7,505✔
47	library_clear();	7,505✔
48	nuclides_clear();	7,505✔
49	free_memory_source();	7,505✔
50	free_memory_mesh();	7,505✔
51	free_memory_tally();	7,505✔
52	free_memory_bank();	7,505✔
53	free_memory_plot();	7,505✔
54	free_memory_weight_windows();	7,505✔
55	if (mpi::master) {	7,505✔
56	free_memory_cmfd();	6,581✔
57	}
58	if (settings::event_based) {	7,505!
UNCOV 59	free_event_queues();	×
60	}
61	}	7,505✔
62
63	} // namespace openmc
64
65	using namespace openmc;
66
67	int openmc_finalize()	7,505✔
68	{
69	if (simulation::initialized)	7,505✔
70	openmc_simulation_finalize();	3✔
71
72	// Clear results
73	openmc_reset();	7,505✔
74
75	// Reset timers
76	reset_timers();	7,505✔
77
78	// Reset global variables
79	settings::assume_separate = false;	7,505✔
80	settings::check_overlaps = false;	7,505✔
81	settings::collision_track_config = CollisionTrackConfig {};	7,505✔
82	settings::confidence_intervals = false;	7,505✔
83	settings::create_fission_neutrons = true;	7,505✔
84	settings::create_delayed_neutrons = true;	7,505✔
85	settings::electron_treatment = ElectronTreatment::LED;	7,505✔
86	settings::delayed_photon_scaling = true;	7,505✔
87	settings::energy_cutoff = {0.0, 1000.0, 0.0, 0.0};	7,505✔
88	settings::time_cutoff = {INFTY, INFTY, INFTY, INFTY};	7,505✔
89	settings::entropy_on = false;	7,505✔
90	settings::event_based = false;	7,505✔
91	settings::free_gas_threshold = 400.0;	7,505✔
92	settings::gen_per_batch = 1;	7,505✔
93	settings::legendre_to_tabular = true;	7,505✔
94	settings::legendre_to_tabular_points = -1;	7,505✔
95	settings::material_cell_offsets = true;	7,505✔
96	settings::max_lost_particles = 10;	7,505✔
97	settings::max_order = 0;	7,505✔
98	settings::max_particles_in_flight = 100000;	7,505✔
99	settings::max_secondaries = 10000;	7,505✔
100	settings::max_particle_events = 1'000'000;	7,505✔
101	settings::max_history_splits = 10'000'000;	7,505✔
102	settings::max_tracks = 1000;	7,505✔
103	settings::max_write_lost_particles = -1;	7,505✔
104	settings::n_log_bins = 8000;	7,505✔
105	settings::n_inactive = 0;	7,505✔
106	settings::n_particles = -1;	7,505✔
107	settings::output_summary = true;	7,505✔
108	settings::output_tallies = true;	7,505✔
109	settings::particle_restart_run = false;	7,505✔
110	settings::path_cross_sections.clear();	7,505✔
111	settings::path_input.clear();	7,505✔
112	settings::path_output.clear();	7,505✔
113	settings::path_particle_restart.clear();	7,505✔
114	settings::path_sourcepoint.clear();	7,505✔
115	settings::path_statepoint.clear();	7,505✔
116	settings::photon_transport = false;	7,505✔
117	settings::reduce_tallies = true;	7,505✔
118	settings::rel_max_lost_particles = 1.0e-6;	7,505✔
119	settings::res_scat_on = false;	7,505✔
120	settings::res_scat_method = ResScatMethod::rvs;	7,505✔
121	settings::res_scat_energy_min = 0.01;	7,505✔
122	settings::res_scat_energy_max = 1000.0;	7,505✔
123	settings::restart_run = false;	7,505✔
124	settings::run_CE = true;	7,505✔
125	settings::run_mode = RunMode::UNSET;	7,505✔
126	settings::solver_type = SolverType::MONTE_CARLO;	7,505✔
127	settings::source_latest = false;	7,505✔
128	settings::source_rejection_fraction = 0.05;	7,505✔
129	settings::source_separate = false;	7,505✔
130	settings::source_write = true;	7,505✔
131	settings::ssw_cell_id = C_NONE;	7,505✔
132	settings::ssw_cell_type = SSWCellType::None;	7,505✔
133	settings::ssw_max_particles = 0;	7,505✔
134	settings::ssw_max_files = 1;	7,505✔
135	settings::survival_biasing = false;	7,505✔
136	settings::temperature_default = 293.6;	7,505✔
137	settings::temperature_method = TemperatureMethod::NEAREST;	7,505✔
138	settings::temperature_multipole = false;	7,505✔
139	settings::temperature_range = {0.0, 0.0};	7,505✔
140	settings::temperature_tolerance = 10.0;	7,505✔
141	settings::trigger_on = false;	7,505✔
142	settings::trigger_predict = false;	7,505✔
143	settings::trigger_batch_interval = 1;	7,505✔
144	settings::uniform_source_sampling = false;	7,505✔
145	settings::ufs_on = false;	7,505✔
146	settings::urr_ptables_on = true;	7,505✔
147	settings::verbosity = -1;	7,505✔
148	settings::weight_cutoff = 0.25;	7,505✔
149	settings::weight_survive = 1.0;	7,505✔
150	settings::weight_windows_file.clear();	7,505✔
151	settings::weight_windows_on = false;	7,505✔
152	settings::write_all_tracks = false;	7,505✔
153	settings::write_initial_source = false;	7,505✔
154
155	simulation::keff = 1.0;	7,505✔
156	simulation::need_depletion_rx = false;	7,505✔
157	simulation::ssw_current_file = 1;	7,505✔
158	simulation::total_gen = 0;	7,505✔
159
160	simulation::entropy_mesh = nullptr;	7,505✔
161	simulation::ufs_mesh = nullptr;	7,505✔
162
163	data::energy_max = {INFTY, INFTY, INFTY, INFTY};	7,505✔
164	data::energy_min = {0.0, 0.0, 0.0, 0.0};	7,505✔
165	data::temperature_min = 0.0;	7,505✔
166	data::temperature_max = INFTY;	7,505✔
167	data::mg = {};	7,505✔
168	model::root_universe = -1;	7,505✔
169	model::plotter_seed = 1;	7,505✔
170	openmc::openmc_set_seed(DEFAULT_SEED);	7,505✔
171	openmc::openmc_set_stride(DEFAULT_STRIDE);	7,505✔
172
173	// Deallocate arrays
174	free_memory();	7,505✔
175
176	#ifdef OPENMC_LIBMESH_ENABLED
177	settings::libmesh_init.reset();	1,553✔
178	#endif
179
180	// Free all MPI types
181	#ifdef OPENMC_MPI
182	if (mpi::source_site != MPI_DATATYPE_NULL) {	3,598✔
183	MPI_Type_free(&mpi::source_site);	3,520✔
184	}
185	if (mpi::collision_track_site != MPI_DATATYPE_NULL) {	3,598✔
186	MPI_Type_free(&mpi::collision_track_site);	3,520✔
187	}
188	#endif
189
190	openmc_finalize_random_ray();	7,505✔
191
192	return 0;	7,505✔
193	}
194
195	int openmc_reset()	14,741✔
196	{
197
198	model::universe_level_counts.clear();	14,741✔
199
200	for (auto& t : model::tallies) {	64,345✔
201	t->reset();	49,604✔
202	}
203
204	// Reset global tallies
205	simulation::n_realizations = 0;	14,741✔
206	xt::view(simulation::global_tallies, xt::all()) = 0.0;	14,741✔
207
208	simulation::k_col_abs = 0.0;	14,741✔
209	simulation::k_col_tra = 0.0;	14,741✔
210	simulation::k_abs_tra = 0.0;	14,741✔
211	simulation::k_sum = {0.0, 0.0};	14,741✔
212	simulation::satisfy_triggers = false;	14,741✔
213
214	settings::cmfd_run = false;	14,741✔
215
216	simulation::n_lost_particles = 0;	14,741✔
217
218	return 0;	14,741✔
219	}
220
221	int openmc_reset_timers()	188✔
222	{
223	reset_timers();	188✔
224	return 0;	188✔
225	}
226
227	int openmc_hard_reset()	170✔
228	{
229	// Reset all tallies and timers
230	openmc_reset();	170✔
231	reset_timers();	170✔
232
233	// Reset total generations and keff guess
234	simulation::keff = 1.0;	170✔
235	simulation::total_gen = 0;	170✔
236
237	// Reset the random number generator state
238	openmc::openmc_set_seed(DEFAULT_SEED);	170✔
239	openmc::openmc_set_stride(DEFAULT_STRIDE);	170✔
240	return 0;	170✔
241	}

openmc-dev / openmc / 21887062786

Source File Press 'n' to go to next uncovered line, 'b' for previous

Source File
Press 'n' to go to next uncovered line, 'b' for previous