Flow123d
jenkins-Flow123d-linux-release-multijob-198
|
Class template representing a field with values dependent on: point, element, and region. More...
#include <field.hh>
Classes | |
class | FactoryBase |
struct | SharedData |
Public Types | |
typedef FieldAlgorithmBase < spacedim, Value > | FieldBaseType |
typedef std::shared_ptr < FieldBaseType > | FieldBasePtr |
typedef FieldAlgorithmBase < spacedim, Value >::Point | Point |
Public Member Functions | |
Field () | |
Field (const string &name, bool bc=false) | |
Field (const Field &other) | |
Field & | operator= (const Field &other) |
IT::AbstractRecord & | get_input_type () override |
IT::Record & | get_multifield_input_type () override |
auto | disable_where (const Field< spacedim, typename FieldValue< spacedim >::Enum > &control_field, const vector< FieldEnum > &value_list) -> Field & |
void | set_mesh (const Mesh &mesh) override |
bool | is_constant (Region reg) override |
void | set_field (const RegionSet &domain, FieldBasePtr field, double time=0.0) |
void | set_field (const RegionSet &domain, const Input::AbstractRecord &a_rec, double time=0.0) |
bool | set_time (const TimeStep &time) override |
void | copy_from (const FieldCommon &other) override |
void | output (OutputTime *stream) override |
FieldResult | field_result (ElementAccessor< spacedim > &elm) const |
virtual Value::return_type const & | value (const Point &p, const ElementAccessor< spacedim > &elm) const |
virtual void | value_list (const std::vector< Point > &point_list, const ElementAccessor< spacedim > &elm, std::vector< typename Value::return_type > &value_list) const |
void | add_factory (std::shared_ptr< FactoryBase > factory) |
Public Member Functions inherited from FieldCommon | |
TYPEDEF_ERR_INFO (EI_Time, double) | |
TYPEDEF_ERR_INFO (EI_Field, std::string) | |
DECLARE_INPUT_EXCEPTION (ExcNonascendingTime,<< "Non-ascending time: "<< EI_Time::val<< " for field "<< EI_Field::qval<< ".\n") | |
DECLARE_INPUT_EXCEPTION (ExcMissingDomain,<< "Missing domain specification (region, r_id, or r_set) in the field descriptor:") | |
DECLARE_EXCEPTION (ExcFieldMeshDifference,<< "Two copies of the field "<< EI_Field::qval<< "call set_mesh with different arguments.\n") | |
FieldCommon & | name (const string &name) |
FieldCommon & | description (const string &description) |
FieldCommon & | input_default (const string &input_default) |
FieldCommon & | units (const UnitSI &units) |
Set basic units of the field. More... | |
FieldCommon & | input_selection (const Input::Type::Selection *element_selection) |
FieldCommon & | output_type (OutputTime::DiscreteSpace rt) |
FieldCommon & | flags (FieldFlag::Flags::Mask mask) |
FieldCommon & | flags_add (FieldFlag::Flags::Mask mask) |
void | set_components (const std::vector< string > &names) |
void | set_input_list (const Input::Array &list) |
void | set_limit_side (LimitSide side) |
const std::string & | input_name () const |
const std::string & | name () const |
const std::string | description () const |
const std::string & | input_default () const |
const UnitSI & | units () const |
OutputTime::DiscreteSpace | output_type () const |
bool | is_bc () const |
unsigned int | n_comp () const |
const Mesh * | mesh () const |
LimitSide | limit_side () const |
FieldFlag::Flags & | flags () |
double | time () const |
void | mark_input_times (TimeMark::Type mark_type) |
bool | changed () const |
void | set_component_index (unsigned int idx) |
bool | is_multifield () const |
virtual | ~FieldCommon () |
Static Public Attributes | |
static constexpr bool | is_enum_valued = boost::is_same<typename Value::element_type, FieldEnum>::value |
static const unsigned int | space_dim = spacedim |
Protected Types | |
typedef pair< double, FieldBasePtr > | HistoryPoint |
Pair: time, pointer to FieldBase instance. More... | |
typedef boost::circular_buffer < HistoryPoint > | RegionHistory |
Nearest history of one region. More... | |
typedef Field< spacedim, typename FieldValue< spacedim > ::Enum > | ControlField |
Protected Types inherited from FieldCommon | |
enum | TimeStatus { TimeStatus::changed, TimeStatus::constant, TimeStatus::unknown } |
Protected Member Functions | |
IT::AbstractRecord | make_input_tree () |
-------— end helper function template More... | |
void | update_history (const TimeStep &time) |
void | check_initialized_region_fields_ () |
Protected Member Functions inherited from FieldCommon | |
FieldCommon () | |
FieldCommon (const FieldCommon &other) | |
void | set_history_changed () |
Protected Attributes | |
std::shared_ptr< SharedData > | data_ |
std::shared_ptr< ControlField > | no_check_control_field_ |
std::vector< FieldBasePtr > | region_fields_ |
std::vector< std::shared_ptr < FactoryBase > > | factories_ |
Protected Attributes inherited from FieldCommon | |
std::string | name_ |
std::shared_ptr< SharedData > | shared_ |
LimitSide | limit_side_ |
TimeStatus | set_time_result_ |
Status of history . More... | |
double | last_time_ = -numeric_limits<double>::infinity() |
OutputTime::DiscreteSpace | type_of_output_data_ = OutputTime::ELEM_DATA |
unsigned int | component_index_ |
bool | multifield_ |
FieldFlag::Flags | flags_ = FieldFlag::declare_input & FieldFlag::equation_input & FieldFlag::allow_output |
Field flags. Default setting is "an equation input field, that can read from user input, and can be written to output". More... | |
Friends | |
template<int dim, class Val > | |
class | MultiField |
Additional Inherited Members | |
Static Public Member Functions inherited from FieldCommon | |
static IT::Record | field_descriptor_record (const string &record_name) |
Static Protected Attributes inherited from FieldCommon | |
static const unsigned int | history_length_limit_ =3 |
Class template representing a field with values dependent on: point, element, and region.
By "field" we mean a mapping of a a pair (Point, Time) to a Value
, where Point is from spacedim
dimensional ambient space, Time is real number (set by set_time
method), and Value
type representing range of the field, which can be: real scalar, integer scalar (a discrete value), real vector of fixed (compile time) size, real vector of runtime size, or a matrix of fixed dimensions. Extensions to vectors or matrices of integers, or to variable tensors are possible. For vector and matrix values we use classes provided by Armadillo library for linear algebra. The Value
template parameter should FieldValue<> template, usual choices are: FieldValue<spacedim>::Scalar, FieldValue<spacedim>::Integer, FieldValue<spacedim>::Enum, FieldValue<spacedim>::VectorFixed, FieldValue<spacedim>::TensorFixed deprecated choices: FieldValue<spacedim>::Vector, FieldValue<spacedim>::VectorEnum.
This class assign particular fields (instances of descendants of FiledBase) to the regions. It keeps a table of pointers to fields for every possible bulk region index (very same functionality, but for boundary regions is provided by BCField
class). This class has interface very similar to FiledBase, however key methods value
, and value_list
are not virtual in this class by contrast these methods are inlined to minimize overhead for simplest fields like FieldConstant.
TODO: currently it works only for spacedim==3 since we have only mesh in 3D ambient space.
|
protected |
If this pointer is set, turn off check of initialization in the set_time
method on the regions where the method get_constant_enum_value
of the control field returns value from no_check_values_
. This field is private copy, its set_time method is called from the set_Time method of actual object.
typedef std::shared_ptr< FieldBaseType > Field< spacedim, Value >::FieldBasePtr |
typedef FieldAlgorithmBase<spacedim, Value> Field< spacedim, Value >::FieldBaseType |
|
protected |
typedef FieldAlgorithmBase<spacedim, Value>::Point Field< spacedim, Value >::Point |
|
protected |
Default constructor.
Definition at line 23 of file field.impl.hh.
Field< spacedim, Value >::Field | ( | const string & | name, |
bool | bc = false |
||
) |
Definition at line 36 of file field.impl.hh.
Field< spacedim, Value >::Field | ( | const Field< spacedim, Value > & | other | ) |
Copy constructor. Keeps shared history, declaration data, mesh.
Definition at line 52 of file field.impl.hh.
void Field< spacedim, Value >::add_factory | ( | std::shared_ptr< FactoryBase > | factory | ) |
Add a new factory for creating Field algorithms on individual regions. The last factory is tried first, the last one is always the default implementation Field<...>::FactoryBase.
The Field<...> object keeps a list of such factories. When the instance of a new field algorithm has to be created from the input field descriptor, we pass through the list of factories backward and let factories to create the field algorithm instance from the actual input field descriptor. The first instance (non-null pointer) is used.
Definition at line 469 of file field.impl.hh.
|
protected |
Check that whole field list (region_fields_
) is set, possibly use default values for unset regions.
Definition at line 412 of file field.impl.hh.
|
overridevirtual |
Check that other has same type and assign from it.
Implements FieldCommon.
Definition at line 312 of file field.impl.hh.
auto Field< spacedim, Value >::disable_where | ( | const Field< spacedim, typename FieldValue< spacedim >::Enum > & | control_field, |
const vector< FieldEnum > & | value_list | ||
) | -> Field & |
By this method you can allow that the field need not to be set on regions (and times) where the given control_field
is FieldConstant and has value in given value_list
. We check this in the set_time method. Through this mechanism we can switch of e.g. boundary data fields according to the type of the boundary condition.
Definition at line 162 of file field.impl.hh.
|
inline |
Returns true, if field is currently set to a time in which it is discontinuous. Special field values spatially constant. Could allow optimization of tensor multiplication and tensor or vector addition. field_result_ should be set in constructor and in set_time method of particular Field implementation. We return value result_none
, if the field is not initialized on the region of the given element accessor elm
.
Definition at line 334 of file field.impl.hh.
|
overridevirtual |
Returns reference to input type of particular field instance, this is static member input_type
of the corresponding FieldBase class (i.e. with same template parameters). However, for fields returning "Enum" we have to create whole unique Input::Type hierarchy using following method meka_input_tree
. every instance since every such field use different Selection for initialization, even if all returns just unsigned int.
Implements FieldCommon.
Definition at line 101 of file field.impl.hh.
|
overridevirtual |
Returns input type for MultiField instance. TODO: temporary solution, see multifield_
Implements FieldCommon.
Definition at line 121 of file field.impl.hh.
|
overridevirtual |
Direct read access to the table of Field pointers on regions. Implementation of FieldCommonBase::is_constant()
.
Implements FieldCommon.
Definition at line 202 of file field.impl.hh.
|
protected |
-------— end helper function template
For fields returning "enum", i.e. with Value
== FieldEnum, the input type (meaning whole input_Type tree of the field) depends on the Input::Type::Selection object that represents particular C enum type. Therefore, we have to create whole tree for the selection that was set through FieldBaseCommon::set_selection()
method.
Definition at line 151 of file field.impl.hh.
Field< spacedim, Value > & Field< spacedim, Value >::operator= | ( | const Field< spacedim, Value > & | other | ) |
Assignment operator. Same properties as copy constructor.
Question: do we really need this, isn't copy constructor enough? Answer: It is necessary in (usual) case when Field instance is created as the class member but is filled later by assignment possibly from other class.
Definition at line 69 of file field.impl.hh.
|
overridevirtual |
Implementation of FieldCommonBase::output().
Implements FieldCommon.
Definition at line 324 of file field.impl.hh.
void Field< spacedim, Value >::set_field | ( | const RegionSet & | domain, |
FieldBasePtr | field, | ||
double | time = 0.0 |
||
) |
Assigns given field
to all regions in given region set domain
. Field is added to the history with given time and possibly used in the next call of the set_time method. Caller is responsible for correct construction of given field.
Use this method only if necessary.
Default time simplify setting steady fields.
Definition at line 210 of file field.impl.hh.
void Field< spacedim, Value >::set_field | ( | const RegionSet & | domain, |
const Input::AbstractRecord & | a_rec, | ||
double | time = 0.0 |
||
) |
Same as before but the field is first created using FieldBase::function_factory(), from given abstract record accessor a_rec
.
Definition at line 237 of file field.impl.hh.
|
overridevirtual |
Set mesh pointer and resize region arrays.
Implements abstract method.
Reimplemented from FieldCommon.
Definition at line 172 of file field.impl.hh.
|
overridevirtual |
Check that whole field list is set, possibly use default values for unset regions and call set_time for every field in the field list.
Returns true if the field has been changed.
Implements FieldCommon.
Definition at line 249 of file field.impl.hh.
|
protected |
Read input into regions_history_
possibly pop some old values from the history queue to keep its size less then history_length_limit_
.
Definition at line 343 of file field.impl.hh.
|
inlinevirtual |
Returns one value in one given point on
an element given by ElementAccessor elm
. It returns reference to he actual value in order to avoid temporaries for vector and tensor values.
|
inlinevirtual |
Returns std::vector of scalar values in several points at once. The base class implements trivial implementation using the value(,,)
method. This is not optimal as it involves lot of virtual calls, but this overhead can be negligible for more complex fields as Python of Formula.
Definition at line 326 of file field.hh.
|
protected |
|
protected |
|
protected |
|
protected |
|
static |