Constructs a grid with the given minimum grid coordinates and grid length.
If you need periodic boundaries you have to use this constructor to set the correct length to use for wrapping particles around the borders.
122 :
length_(min_and_length.second) {
123 const auto min_position = min_and_length.first;
124 const SizeType particle_count = particles.size();
132 cells_.emplace_back(particles.copy_to_vector());
143 const int max_cells = numeric_cast<int>(
145 ? std::floor(std::cbrt(particle_count))
146 : std::max(2., std::floor(std::cbrt(particle_count))));
151 std::array<double, 3> index_factor = {1. / max_interaction_length,
152 1. / max_interaction_length,
153 1. / max_interaction_length};
158 std::numeric_limits<int>::max() / index_factor[i])) {
159 throw std::overflow_error(
160 "An integer overflow would occur constructing the system grid.\n"
161 "Impossible to (further) simulate the provided system using "
163 "Refer to the user guide for more information (see the Modi "
167 static_cast<int>(std::floor(
length_[i] * index_factor[i]));
179 std::string error_box_too_small =
180 "Input error: With the chosen time step (Delta_Time), your box is\n"
181 "too small for the grid. Using the provided time step, the minimal\n"
182 "length of the box should be " +
183 std::to_string(2 * max_interaction_length) +
184 "fm. Using a smaller time step\n"
185 "will reduce the minimal needed box size. The use of test particles\n"
186 "also helps reducing the minimum needed size. Have a look to the\n"
187 "user guide (e.g. box modus page) for further information.\n"
188 "Please, adjust your config file and run SMASH again.";
189 throw std::runtime_error(error_box_too_small);
197 if (
length_[i] >= max_interaction_length) {
202 index_factor[i] = std::nextafter(index_factor[i], 0.);
217 " cells. Therefore the Grid falls back to a single cell / "
221 if (include_unformed_particles) {
224 cells_.emplace_back(particles.copy_to_vector());
228 cells_.front().reserve(particles.size());
229 std::copy_if(particles.begin(), particles.end(),
230 std::back_inserter(
cells_.front()),
231 [&](
const ParticleData &
p) {
232 return p.xsec_scaling_factor(timestep_duration) > 0.0;
245 "\nindex_factor: ", index_factor);
256 auto &&cell_index_for = [&](
const ParticleData &
p) {
258 numeric_cast<SizeType>(std::floor(
259 (
p.position()[1] - min_position[0]) * index_factor[0])),
260 numeric_cast<SizeType>(std::floor(
261 (
p.position()[2] - min_position[1]) * index_factor[1])),
262 numeric_cast<SizeType>(std::floor(
263 (
p.position()[3] - min_position[2]) * index_factor[2])));
265 for (
const auto &
p : particles) {
266 if (!include_unformed_particles &&
267 (
p.xsec_scaling_factor(timestep_duration) <= 0.0)) {
270 const auto idx = cell_index_for(
p);
275 "\nan out-of-bounds access would be necessary for the "
278 "\nfor a grid with the following parameters:\nmin: ", min_position,
280 "\nindex_factor: ", index_factor,
"\ncells_.size: ",
cells_.size(),
281 "\nrequested index: ", idx);
282 throw std::runtime_error(
"out-of-bounds grid access on construction");
int SizeType
A type to store the sizes.
const std::array< double, 3 > length_
The 3 lengths of the complete grid. Used for periodic boundary wrapping.
std::array< int, 3 > number_of_cells_
The number of cells in x, y, and z direction.
double cell_volume_
The volume of a single cell.
std::vector< ParticleList > cells_
The cell storage.
SizeType make_index(SizeType x, SizeType y, SizeType z) const
#define SMASH_SOURCE_LOCATION
Hackery that is required to output the location in the source code where the log statement occurs.
std::array< einhard::Logger<>, std::tuple_size< LogArea::AreaTuple >::value > logg
An array that stores all pre-configured Logger objects.
#define unlikely(x)
Tell the branch predictor that this expression is likely false.
@ ParticleNumber
Limit the number of cells to the number of particles.
bool all_of(Container &&c, UnaryPredicate &&p)
Convenience wrapper for std::all_of that operates on a complete container.
@ Normal
Without ghost cells.
@ PeriodicBoundaries
With ghost cells for periodic boundaries.
static constexpr int LGrid
@ Largest
Make cells as large as possible.
Public Types inherited from smash::GridBase