20278515727

Committed 16 Dec 2025 06:26PM UTC coverage: 81.822% (-0.1%) from 81.92%

Build # 20278515727

Build Type

Pull #3493

github

Committed by

web-flow

Commit Message

Merge 7f36869a3 into bbfa18d72

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

Coverage Stats

17013 of 23575 branches covered (72.17%)

Branch coverage included in aggregate %.

188 of 207 new or added lines in 3 files covered. (90.82%)

3101 existing lines in 56 files now uncovered.

54973 of 64404 relevant lines covered (85.36%)

41385524.25 hits per line

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

98.77

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

openmc-dev / openmc / 20278515727

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