46 <<
"Your probability density function in the particle generator "
47 "returned a negative value for the following position: "
48 <<
arg1 <<
". Please check your function expression.");
61 template <
int dim,
int spacedim = dim>
86 for (
const auto &cell :
triangulation.active_cell_iterators())
88 if (cell->is_locally_owned())
90 fe_values.reinit(cell);
93 fe_values.get_quadrature_points()[0]);
116 template <
int dim,
int spacedim>
123 Assert(cell->reference_cell().is_hyper_cube() ==
true,
140 for (
unsigned int d = 0;
d < spacedim; ++
d)
150 mapping.transform_real_to_unit_cell(cell, position);
167 for (
unsigned int d = 0;
d < dim; ++
d)
180 template <
int dim,
int spacedim>
192#ifdef DEAL_II_WITH_MPI
193 if (
const auto tria =
207 tria->get_mpi_communicator());
215 for (
const auto &cell :
triangulation.active_cell_iterators() |
236 template <
int dim,
int spacedim>
253 template <
int dim,
int spacedim>
273 template <
int dim,
int spacedim>
285 if (
const auto tria =
318 if (
const auto tria =
334 "The integral of the user prescribed probability "
335 "density function over the domain equals zero, "
336 "deal.II has no way to determine the cell of "
337 "generated particles. Please ensure that the "
338 "provided function is positive in at least a "
339 "part of the domain; also check the syntax of "
390 for (
const auto &cell :
triangulation.active_cell_iterators() |
417 for (
const auto &cell :
triangulation.active_cell_iterators() |
420 for (
unsigned int i = 0;
440 template <
int dim,
int spacedim>
448 const std::vector<std::vector<double>> &properties)
450 const auto &fe = dof_handler.get_fe();
457 const std::map<types::global_dof_index, Point<spacedim>>
474 template <
int dim,
int spacedim>
484 const std::vector<std::vector<double>> &properties)
487 quadrature.get_points();
493 for (
const auto &cell :
triangulation.active_cell_iterators() |
510#include "particles/generators.inst"
unsigned int size() const
Abstract base class for mapping classes.
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
static ::ExceptionBase & ExcNotImplemented()
#define Assert(cond, exc)
#define AssertThrowMPI(error_code)
#define DeclException1(Exception1, type1, outsequence)
static ::ExceptionBase & ExcMessage(std::string arg1)
#define AssertThrow(cond, exc)
@ update_JxW_values
Transformed quadrature weights.
@ update_quadrature_points
Transformed quadrature points.
const unsigned int my_rank
void probabilistic_locations(const Triangulation< dim, spacedim > &triangulation, const Function< spacedim > &probability_density_function, const bool random_cell_selection, const types::particle_index n_particles_to_create, ParticleHandler< dim, spacedim > &particle_handler, const Mapping< dim, spacedim > &mapping=(ReferenceCells::get_hypercube< dim >() .template get_default_linear_mapping< dim, spacedim >()), const unsigned int random_number_seed=5432)
ParticleIterator< dim, spacedim > random_particle_in_cell_insert(const typename Triangulation< dim, spacedim >::active_cell_iterator &cell, const types::particle_index id, std::mt19937 &random_number_generator, ParticleHandler< dim, spacedim > &particle_handler, const Mapping< dim, spacedim > &mapping=(ReferenceCells::get_hypercube< dim >() .template get_default_linear_mapping< dim, spacedim >()))
void quadrature_points(const Triangulation< dim, spacedim > &triangulation, const Quadrature< dim > &quadrature, const std::vector< std::vector< BoundingBox< spacedim > > > &global_bounding_boxes, ParticleHandler< dim, spacedim > &particle_handler, const Mapping< dim, spacedim > &mapping=(ReferenceCells::get_hypercube< dim >() .template get_default_linear_mapping< dim, spacedim >()), const std::vector< std::vector< double > > &properties={})
Particle< dim, spacedim > random_particle_in_cell(const typename Triangulation< dim, spacedim >::active_cell_iterator &cell, const types::particle_index id, std::mt19937 &random_number_generator, const Mapping< dim, spacedim > &mapping=(ReferenceCells::get_hypercube< dim >() .template get_default_linear_mapping< dim, spacedim >()))
void dof_support_points(const DoFHandler< dim, spacedim > &dof_handler, const std::vector< std::vector< BoundingBox< spacedim > > > &global_bounding_boxes, ParticleHandler< dim, spacedim > &particle_handler, const Mapping< dim, spacedim > &mapping=(ReferenceCells::get_hypercube< dim >() .template get_default_linear_mapping< dim, spacedim >()), const ComponentMask &components={}, const std::vector< std::vector< double > > &properties={})
void regular_reference_locations(const Triangulation< dim, spacedim > &triangulation, const std::vector< Point< dim > > &particle_reference_locations, ParticleHandler< dim, spacedim > &particle_handler, const Mapping< dim, spacedim > &mapping=(ReferenceCells::get_hypercube< dim >() .template get_default_linear_mapping< dim, spacedim >()))
SymmetricTensor< 2, dim, Number > d(const Tensor< 2, dim, Number > &F, const Tensor< 2, dim, Number > &dF_dt)
std::pair< T, T > partial_and_total_sum(const T &value, const MPI_Comm comm)
unsigned int this_mpi_process(const MPI_Comm mpi_communicator)
const ::parallel::distributed::Triangulation< dim, spacedim > * triangulation
#define DEAL_II_PARTICLE_INDEX_MPI_TYPE