template <int dim, int spacedim>
ParticleHandler<dim, spacedim>::ParticleHandler(
const Triangulation<dim, spacedim> &triangulation,
- const Mapping<dim, spacedim> &mapping,
+ const Mapping<dim, spacedim> & mapping,
const unsigned int n_properties)
: triangulation(&triangulation, typeid(*this).name())
, mapping(&mapping, typeid(*this).name())
void
ParticleHandler<dim, spacedim>::initialize(
const Triangulation<dim, spacedim> &new_triangulation,
- const Mapping<dim, spacedim> &new_mapping,
+ const Mapping<dim, spacedim> & new_mapping,
const unsigned int n_properties)
{
clear();
template <int dim, int spacedim>
typename ParticleHandler<dim, spacedim>::particle_iterator
ParticleHandler<dim, spacedim>::insert_particle(
- const Particle<dim, spacedim> &particle,
+ const Particle<dim, spacedim> & particle,
const typename Triangulation<dim, spacedim>::active_cell_iterator &cell)
{
return insert_particle(particle.get_location(),
template <int dim, int spacedim>
typename ParticleHandler<dim, spacedim>::particle_iterator
ParticleHandler<dim, spacedim>::insert_particle(
- const void *&data,
+ const void *& data,
const typename Triangulation<dim, spacedim>::active_cell_iterator &cell)
{
Assert(triangulation != nullptr, ExcInternalError());
template <int dim, int spacedim>
typename ParticleHandler<dim, spacedim>::particle_iterator
ParticleHandler<dim, spacedim>::insert_particle(
- const Point<spacedim> &position,
- const Point<dim> &reference_position,
+ const Point<spacedim> & position,
+ const Point<dim> & reference_position,
const types::particle_index particle_index,
const typename Triangulation<dim, spacedim>::active_cell_iterator &cell,
const ArrayView<const double> &properties)
ParticleHandler<dim, spacedim>::insert_global_particles(
const std::vector<Point<spacedim>> &positions,
const std::vector<std::vector<BoundingBox<spacedim>>>
- &global_bounding_boxes,
- const std::vector<std::vector<double>> &properties,
+ & global_bounding_boxes,
+ const std::vector<std::vector<double>> & properties,
const std::vector<types::particle_index> &ids)
{
if (!properties.empty())
compare_particle_association(
const unsigned int a,
const unsigned int b,
- const Tensor<1, dim> &particle_direction,
+ const Tensor<1, dim> & particle_direction,
const std::vector<Tensor<1, dim>> ¢er_directions)
{
const double scalar_product_a = center_directions[a] * particle_direction;
std::set<unsigned int> cell_to_neighbor_subdomain;
for (const unsigned int v : cell->vertex_indices())
{
- const auto vertex_ghost_neighbors = vertices_with_ghost_neighbors.find(cell->vertex_index(v));
+ const auto vertex_ghost_neighbors =
+ vertices_with_ghost_neighbors.find(cell->vertex_index(v));
if (vertex_ghost_neighbors !=
vertices_with_ghost_neighbors.end())
{
const std::map<
types::subdomain_id,
std::vector<typename Triangulation<dim, spacedim>::active_cell_iterator>>
- &send_cells,
+ & send_cells,
const bool build_cache)
{
Assert(triangulation != nullptr, ExcInternalError());
template <int dim, int spacedim>
void
ParticleHandler<dim, spacedim>::register_additional_store_load_functions(
- const std::function<std::size_t()> &size_callb,
+ const std::function<std::size_t()> & size_callb,
const std::function<void *(const particle_iterator &, void *)> &store_callb,
const std::function<const void *(const particle_iterator &, const void *)>
&load_callb)
const auto callback_function =
[this](
const typename Triangulation<dim, spacedim>::cell_iterator
- &cell_iterator,
+ & cell_iterator,
const typename Triangulation<dim, spacedim>::CellStatus cell_status) {
return this->pack_callback(cell_iterator, cell_status);
};
template <int dim, int spacedim>
void
ParticleHandler<dim, spacedim>::unpack_callback(
- const typename Triangulation<dim, spacedim>::cell_iterator &cell,
+ const typename Triangulation<dim, spacedim>::cell_iterator & cell,
const typename Triangulation<dim, spacedim>::CellStatus status,
const boost::iterator_range<std::vector<char>::const_iterator> &data_range)
{
// now update particle storage location and properties if necessary
switch (status)
{
- case parallel::TriangulationBase<dim, spacedim>::CELL_PERSIST: {
+ case parallel::TriangulationBase<dim, spacedim>::CELL_PERSIST:
+ {
// all particles are correctly inserted
}
break;
- case parallel::TriangulationBase<dim, spacedim>::CELL_COARSEN: {
+ case parallel::TriangulationBase<dim, spacedim>::CELL_COARSEN:
+ {
// all particles are in correct cell, but their reference location
// has changed
for (auto &particle : loaded_particles_on_cell)
}
break;
- case parallel::TriangulationBase<dim, spacedim>::CELL_REFINE: {
+ case parallel::TriangulationBase<dim, spacedim>::CELL_REFINE:
+ {
// we need to find the correct child to store the particles and
// their reference location has changed
typename particle_container::iterator &cache =