22 #ifndef _aspect_boundary_composition_interface_h 23 #define _aspect_boundary_composition_interface_h 30 #include <deal.II/base/parameter_handler.h> 31 #include <deal.II/distributed/tria.h> 33 #include <boost/core/demangle.hpp> 45 namespace BoundaryComposition
62 virtual ~Interface() =
default;
102 boundary_composition (
const types::boundary_id boundary_indicator,
103 const Point<dim> &position,
104 const unsigned int compositional_field)
const = 0;
124 parse_parameters (ParameterHandler &prm);
140 ~Manager ()
override;
169 parse_parameters (ParameterHandler &prm);
177 boundary_composition (
const types::boundary_id boundary_indicator,
178 const Point<dim> &position,
179 const unsigned int compositional_field)
const;
200 register_boundary_composition (
const std::string &name,
201 const std::string &description,
202 void (*declare_parameters_function) (ParameterHandler &),
203 std::unique_ptr<Interface<dim>> (*factory_function) ());
210 const std::vector<std::string> &
211 get_active_boundary_composition_names ()
const;
217 const std::vector<std::unique_ptr<Interface<dim>>> &
218 get_active_boundary_composition_conditions ()
const;
229 template <
typename BoundaryCompositionType,
230 typename =
typename std::enable_if_t<std::is_base_of<Interface<dim>,BoundaryCompositionType>::value>>
232 has_matching_boundary_composition_model ()
const;
245 template <
typename BoundaryCompositionType,
246 typename =
typename std::enable_if_t<std::is_base_of<Interface<dim>,BoundaryCompositionType>::value>>
247 const BoundaryCompositionType &
248 get_matching_boundary_composition_model ()
const;
254 const std::set<types::boundary_id> &
255 get_fixed_composition_boundary_indicators()
const;
262 allows_fixed_composition_on_outflow_boundaries()
const;
283 <<
"Could not find entry <" 285 <<
"> among the names of registered boundary composition objects.");
323 template <
typename BoundaryCompositionType,
typename>
328 for (
const auto &p : boundary_composition_objects)
329 if (Plugins::plugin_type_matches<BoundaryCompositionType>(*p))
338 template <
typename BoundaryCompositionType,
typename>
340 const BoundaryCompositionType &
343 AssertThrow(has_matching_boundary_composition_model<BoundaryCompositionType> (),
344 ExcMessage(
"You asked BoundaryComposition::Manager::get_boundary_composition_model() for a " 345 "boundary composition model of type <" + boost::core::demangle(
typeid(BoundaryCompositionType).name()) +
"> " 346 "that could not be found in the current model. Activate this " 347 "boundary composition model in the input file."));
349 typename std::vector<std::unique_ptr<Interface<dim>>>::const_iterator boundary_composition_model;
350 for (
typename std::vector<std::unique_ptr<
Interface<dim>>>::const_iterator
351 p = boundary_composition_objects.begin();
352 p != boundary_composition_objects.end(); ++p)
353 if (Plugins::plugin_type_matches<BoundaryCompositionType>(*(*p)))
354 return Plugins::get_plugin_as_type<BoundaryCompositionType>(*(*p));
357 return Plugins::get_plugin_as_type<BoundaryCompositionType>(*(*boundary_composition_model));
381 #define ASPECT_REGISTER_BOUNDARY_COMPOSITION_MODEL(classname, name, description) \ 382 template class classname<2>; \ 383 template class classname<3>; \ 384 namespace ASPECT_REGISTER_BOUNDARY_COMPOSITION_MODEL_ ## classname \ 386 aspect::internal::Plugins::RegisterHelper<aspect::BoundaryComposition::Interface<2>,classname<2>> \ 387 dummy_ ## classname ## _2d (&aspect::BoundaryComposition::Manager<2>::register_boundary_composition, \ 388 name, description); \ 389 aspect::internal::Plugins::RegisterHelper<aspect::BoundaryComposition::Interface<3>,classname<3>> \ 390 dummy_ ## classname ## _3d (&aspect::BoundaryComposition::Manager<3>::register_boundary_composition, \ 391 name, description); \ std::vector< aspect::Utilities::Operator > model_operators
bool has_matching_boundary_composition_model() const
std::vector< std::unique_ptr< Interface< dim > > > boundary_composition_objects
void write_plugin_graph(std::ostream &output_stream)
const BoundaryCompositionType & get_matching_boundary_composition_model() const
bool allow_fixed_composition_on_outflow_boundaries
std::vector< std::string > model_names
void declare_parameters(ParameterHandler &prm)
std::string get_valid_model_names_pattern()
std::set< types::boundary_id > fixed_composition_boundary_indicators
DeclException1(ProbabilityFunctionNegative, Point< dim >,<< "Your probability density function in the particle generator " "returned a negative probability density for the following position: "<< arg1<< ". Please check your function expression.")