crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-27 05:35:49 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

Add ThermalScattering::has_nuclide method

Browse source
This commit is contained in:
Paul Romano 2018-08-15 10:54:02 -05:00
parent aab7b7ca41
commit 7aa6db0c7e
2 changed files with 15 additions and 10 deletions
Download patch file Download diff file Expand all files Collapse all files

22
src/thermal.cpp
View file

@ -1,6 +1,6 @@
#include "thermal.h"
#include <algorithm> // for sort, move, min, max
#include <algorithm> // for sort, move, min, max, find
#include <cmath> // for round, sqrt, fabs
#include <sstream> // for stringstream
@ -148,7 +148,7 @@ ThermalScattering::ThermalScattering(hid_t group, const std::vector<double>& tem
void
ThermalScattering::calculate_xs(double E, double sqrtkT, int* i_temp,
double* elastic, double* inelastic)
double* elastic, double* inelastic) const
{
// Determine temperature for S(a,b) table
double kT = sqrtkT*sqrtkT;
@ -180,12 +180,9 @@ ThermalScattering::calculate_xs(double E, double sqrtkT, int* i_temp,
auto& sab = data_[i];
// Get index and interpolation factor for inelastic grid
int i_grid;
double f;
if (E < sab.inelastic_e_in_.front()) {
i_grid = 0;
f = 0.0;
} else {
int i_grid = 0;
double f = 0.0;
if (E >= sab.inelastic_e_in_.front()) {
auto& E_in = sab.inelastic_e_in_;
i_grid = lower_bound_index(E_in.begin(), E_in.end(), E);
f = (E - E_in[i_grid]) / (E_in[i_grid+1] - E_in[i_grid]);
@ -221,7 +218,7 @@ ThermalScattering::calculate_xs(double E, double sqrtkT, int* i_temp,
}
// Get interpolation factor for elastic grid
f = (E - E_in[i_grid])/(E_in[i_grid+1] - E_in[i_grid]);
f = (E - E_in[i_grid]) / (E_in[i_grid+1] - E_in[i_grid]);
// Calculate S(a,b) elastic scattering cross section
auto& xs = sab.elastic_P_;
@ -233,6 +230,13 @@ ThermalScattering::calculate_xs(double E, double sqrtkT, int* i_temp,
}
}
bool
ThermalScattering::has_nuclide(const char* name) const
{
std::string nuc {name};
return std::find(nuclides_.begin(), nuclides_.end(), nuc) != nuclides_.end();
}
//==============================================================================
// ThermalData implementation
//==============================================================================

3
src/thermal.h
View file

@ -78,7 +78,8 @@ public:
double tolerance, const double* minmax);
void calculate_xs(double E, double sqrtkT, int* i_temp, double* elastic,
double* inelastic);
double* inelastic) const ;
bool has_nuclide(const char* name) const;
std::string name_; // name of table, e.g. "c_H_in_H2O"
double awr_; // weight of nucleus in neutron masses