15#ifndef dealii_particles_particle_handler_h
16#define dealii_particles_particle_handler_h
37#include <boost/range/iterator_range.hpp>
38#include <boost/serialization/map.hpp>
39#include <boost/signals2.hpp>
63 template <
int dim,
int spacedim = dim>
98 const unsigned int n_properties = 0);
113 const unsigned int n_properties = 0);
167 reserve(
const std::size_t n_particles);
395 std::map<unsigned int, IndexSet>
399 &global_bounding_boxes,
400 const std::vector<std::vector<double>> &properties = {},
401 const std::vector<types::particle_index> &ids = {});
433 std::map<unsigned int, IndexSet>
435 const std::vector<Particle<dim, spacedim>> &
particles,
437 &global_bounding_boxes);
482 template <
typename VectorType>
484 std::is_convertible_v<VectorType *, Function<spacedim> *> ==
false>
486 const bool displace_particles =
true);
510 const bool displace_particles =
true);
533 const bool displace_particles =
true);
562 template <
typename VectorType>
565 const bool add_to_output_vector =
false);
583 const bool add_to_output_vector =
false);
699 PropertyPool<dim, spacedim> &
801 template <
class Archive>
803 serialize(Archive &ar,
const unsigned int version);
832 boost::signals2::signal<void(
1012#ifdef DEAL_II_WITH_MPI
1044 &new_cells_for_particles = std::map<
1048 const bool enable_cache =
false);
1064 &particles_to_send);
1136 const boost::iterator_range<std::vector<char>::const_iterator>
1143 typename particle_container::iterator
1150 typename particle_container::iterator
1157 typename particle_container::iterator
1164 typename particle_container::iterator
1173 template <
int dim,
int spacedim>
1182 template <
int dim,
int spacedim>
1193 template <
int dim,
int spacedim>
1202 template <
int dim,
int spacedim>
1211 template <
int dim,
int spacedim>
1220 template <
int dim,
int spacedim>
1231 template <
int dim,
int spacedim>
1240 template <
int dim,
int spacedim>
1249 template <
int dim,
int spacedim>
1256 typename particle_container::iterator begin =
1263 template <
int dim,
int spacedim>
1268 return owned_particles_end;
1273 template <
int dim,
int spacedim>
1278 typename particle_container::iterator begin = owned_particles_end;
1284 template <
int dim,
int spacedim>
1289 typename particle_container::iterator end =
1296 template <
int dim,
int spacedim>
1297 template <
class Archive>
1308 ar &global_number_of_particles &global_max_particles_per_cell
1309 &next_free_particle_index;
1314 template <
int dim,
int spacedim>
1315 template <
typename VectorType>
1316 inline std::enable_if_t<
1317 std::is_convertible_v<VectorType *, Function<spacedim> *> ==
false>
1319 const VectorType &input_vector,
1320 const bool displace_particles)
1323 get_next_free_particle_index() * spacedim);
1324 for (
auto &p : *
this)
1327 const auto id = p.get_id();
1329 if (displace_particles)
1330 for (
unsigned int i = 0; i < spacedim; ++i)
1331 position[i] += input_vector[
id * spacedim + i];
1333 for (
unsigned int i = 0; i < spacedim; ++i)
1334 position[i] = input_vector[
id * spacedim + i];
1336 sort_particles_into_subdomains_and_cells();
1341 template <
int dim,
int spacedim>
1342 template <
typename VectorType>
1345 VectorType &output_vector,
1346 const bool add_to_output_vector)
1349 get_next_free_particle_index() * spacedim);
1350 for (
const auto &p : *
this)
1352 auto point = p.get_location();
1353 const auto id = p.get_id();
1354 if (add_to_output_vector)
1355 for (
unsigned int i = 0; i < spacedim; ++i)
1356 output_vector[
id * spacedim + i] += point[i];
1358 for (
unsigned int i = 0; i < spacedim; ++i)
1359 output_vector[
id * spacedim + i] = point[i];
1361 if (add_to_output_vector)
Abstract base class for mapping classes.
std::list< ParticlesInCell > particle_container
std::function< void *(const particle_iterator &, void *)> store_callback
void register_additional_store_load_functions(const std::function< std::size_t()> &size_callback, const std::function< void *(const particle_iterator &, void *)> &store_callback, const std::function< const void *(const particle_iterator &, const void *)> &load_callback)
void exchange_ghost_particles(const bool enable_ghost_cache=false)
unsigned int tria_attached_data_index
types::particle_index n_global_particles() const
particle_iterator end_ghost() const
unsigned int global_max_particles_per_cell
void unpack_after_coarsening_and_refinement()
void register_data_attach()
internal::GhostParticlePartitioner< dim, spacedim > ghost_particles_cache
particle_container::iterator particle_container_ghost_begin() const
void prepare_for_serialization()
unsigned int n_properties_per_particle() const
types::particle_index global_number_of_particles
void get_particle_positions(VectorType &output_vector, const bool add_to_output_vector=false)
particle_container::iterator particle_container_owned_end() const
const double tolerance_inside_cell
void update_cached_numbers()
void prepare_for_coarsening_and_refinement()
particle_container::iterator particle_container_ghost_end() const
boost::iterator_range< particle_iterator > particle_iterator_range
void post_mesh_change_action()
void send_recv_particles(const std::map< types::subdomain_id, std::vector< particle_iterator > > &particles_to_send, const std::map< types::subdomain_id, std::vector< typename Triangulation< dim, spacedim >::active_cell_iterator > > &new_cells_for_particles=std::map< types::subdomain_id, std::vector< typename Triangulation< dim, spacedim >::active_cell_iterator > >(), const bool enable_cache=false)
ParticleIterator< dim, spacedim > particle_iterator
ObserverPointer< const Mapping< dim, spacedim >, ParticleHandler< dim, spacedim > > mapping
void send_recv_particles_properties_and_location(const std::map< types::subdomain_id, std::vector< particle_iterator > > &particles_to_send)
types::particle_index number_of_locally_owned_particles
particle_iterator begin_ghost() const
PropertyPool< dim, spacedim > & get_property_pool() const
void initialize(const Triangulation< dim, spacedim > &tria, const Mapping< dim, spacedim > &mapping, const unsigned int n_properties=0)
void serialize(Archive &ar, const unsigned int version)
std::unique_ptr< PropertyPool< dim, spacedim > > property_pool
types::particle_index get_max_local_particle_index() const
void reserve(const std::size_t n_particles)
particle_iterator begin() const
std::map< unsigned int, IndexSet > 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={}, const std::vector< types::particle_index > &ids={})
void notify_ready_to_unpack(const bool serialization)
void connect_to_triangulation_signals()
std::enable_if_t< std::is_convertible_v< VectorType *, Function< spacedim > * >==false > set_particle_positions(const VectorType &input_vector, const bool displace_particles=true)
std::function< std::size_t()> size_callback
std::vector< char > pack_callback(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const CellStatus status) const
void insert_particles(const std::multimap< typename Triangulation< dim, spacedim >::active_cell_iterator, Particle< dim, spacedim > > &particles)
types::particle_index next_free_particle_index
const particle_container::iterator owned_particles_end
particle_iterator end() const
void remove_particle(const particle_iterator &particle)
types::particle_index n_locally_owned_particles() const
void reset_particle_container(particle_container &particles)
types::particle_index n_particles_in_cell(const typename Triangulation< dim, spacedim >::active_cell_iterator &cell) const
void update_ghost_particles()
typename ParticleAccessor< dim, spacedim >::particle_container particle_container
ObserverPointer< const Triangulation< dim, spacedim >, ParticleHandler< dim, spacedim > > triangulation
particle_iterator_range particles_in_cell(const typename Triangulation< dim, spacedim >::active_cell_iterator &cell)
particle_iterator insert_particle(const Particle< dim, spacedim > &particle, const typename Triangulation< dim, spacedim >::active_cell_iterator &cell)
void sort_particles_into_subdomains_and_cells()
void remove_particles(const std::vector< particle_iterator > &particles)
types::particle_index n_global_max_particles_per_cell() const
std::vector< boost::signals2::connection > tria_listeners
particle_container::iterator particle_container_owned_begin() const
void unpack_callback(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const CellStatus status, const boost::iterator_range< std::vector< char >::const_iterator > &data_range)
virtual ~ParticleHandler()
std::vector< typename particle_container::iterator > cells_to_particle_cache
void copy_from(const ParticleHandler< dim, spacedim > &particle_handler)
std::function< const void *(const particle_iterator &, const void *)> load_callback
std::unique_ptr< GridTools::Cache< dim, spacedim > > triangulation_cache
types::particle_index get_next_free_particle_index() const
particle_container particles
IndexSet locally_owned_particle_ids() const
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define Assert(cond, exc)
#define AssertDimension(dim1, dim2)
static ::ExceptionBase & ExcInternalError()
std::vector< index_type > data
boost::signals2::signal< void(const typename Particles::ParticleIterator< dim, spacedim > &particle, const typename Triangulation< dim, spacedim >::active_cell_iterator &cell)> particle_lost