22 #ifndef _aspect_initial_temperature_interface_h 23 #define _aspect_initial_temperature_interface_h 29 #include <deal.II/base/point.h> 30 #include <deal.II/base/parameter_handler.h> 32 #include <boost/core/demangle.hpp> 38 template <
int dim>
class SimulatorAccess;
46 namespace InitialTemperature
63 virtual ~Interface() =
default;
78 double initial_temperature (
const Point<dim> &position)
const = 0;
99 parse_parameters (ParameterHandler &prm);
117 ~Manager ()
override;
133 parse_parameters (ParameterHandler &prm);
141 initial_temperature (
const Point<dim> &position)
const;
162 register_initial_temperature (
const std::string &name,
163 const std::string &description,
164 void (*declare_parameters_function) (ParameterHandler &),
165 std::unique_ptr<Interface<dim>> (*factory_function) ());
172 const std::vector<std::string> &
173 get_active_initial_temperature_names ()
const;
179 const std::list<std::unique_ptr<Interface<dim>>> &
180 get_active_initial_temperature_conditions ()
const;
191 template <
typename InitialTemperatureType,
192 typename =
typename std::enable_if_t<std::is_base_of<Interface<dim>,InitialTemperatureType>::value>>
194 has_matching_initial_temperature_model ()
const;
207 template <
typename InitialTemperatureType,
208 typename =
typename std::enable_if_t<std::is_base_of<Interface<dim>,InitialTemperatureType>::value>>
209 const InitialTemperatureType &
210 get_matching_initial_temperature_model ()
const;
231 <<
"Could not find entry <" 233 <<
"> among the names of registered initial temperature objects.");
259 template <
typename InitialTemperatureType,
typename>
264 for (
const auto &p : initial_temperature_objects)
265 if (Plugins::plugin_type_matches<InitialTemperatureType>(*p))
272 template <
typename InitialTemperatureType,
typename>
274 const InitialTemperatureType &
277 AssertThrow(has_matching_initial_temperature_model<InitialTemperatureType> (),
278 ExcMessage(
"You asked InitialTemperature::Manager::get_initial_temperature_model() for a " 279 "initial temperature model of type <" + boost::core::demangle(
typeid(InitialTemperatureType).name()) +
"> " 280 "that could not be found in the current model. Activate this " 281 "initial temperature model in the input file."));
283 typename std::list<std::unique_ptr<Interface<dim>>>::const_iterator initial_temperature_model;
284 for (
typename std::list<std::unique_ptr<
Interface<dim>>>::const_iterator
285 p = initial_temperature_objects.begin();
286 p != initial_temperature_objects.end(); ++p)
287 if (Plugins::plugin_type_matches<InitialTemperatureType>(*(*p)))
288 return Plugins::get_plugin_as_type<InitialTemperatureType>(*(*p));
291 return Plugins::get_plugin_as_type<InitialTemperatureType>(*(*initial_temperature_model));
314 #define ASPECT_REGISTER_INITIAL_TEMPERATURE_MODEL(classname,name,description) \ 315 template class classname<2>; \ 316 template class classname<3>; \ 317 namespace ASPECT_REGISTER_INITIAL_TEMPERATURE_MODEL_ ## classname \ 319 aspect::internal::Plugins::RegisterHelper<aspect::InitialTemperature::Interface<2>,classname<2>> \ 320 dummy_ ## classname ## _2d (&aspect::InitialTemperature::Manager<2>::register_initial_temperature, \ 321 name, description); \ 322 aspect::internal::Plugins::RegisterHelper<aspect::InitialTemperature::Interface<3>,classname<3>> \ 323 dummy_ ## classname ## _3d (&aspect::InitialTemperature::Manager<3>::register_initial_temperature, \ 324 name, description); \ bool has_matching_initial_temperature_model() const
std::string get_valid_model_names_pattern()
void write_plugin_graph(std::ostream &output_stream)
const InitialTemperatureType & get_matching_initial_temperature_model() const
std::vector< std::string > model_names
std::vector< aspect::Utilities::Operator > model_operators
std::list< std::unique_ptr< Interface< dim > > > initial_temperature_objects
void declare_parameters(ParameterHandler &prm)
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.")