Flow123d  release_2.2.0-33-g759111d
Classes | Public Member Functions | Static Public Member Functions | Protected Types | Protected Member Functions | Protected Attributes | Static Protected Attributes | Friends | List of all members
FieldCommon Class Referenceabstract

Common abstract parent of all Field<...> classes. More...

#include <field_common.hh>

Inheritance diagram for FieldCommon:
Inheritance graph
[legend]
Collaboration diagram for FieldCommon:
Collaboration graph
[legend]

Classes

struct  MessageData
 Store data of one initialization message. More...
 
struct  SharedData
 

Public Member Functions

 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 or r_id) in the field descriptor:")
 
 DECLARE_EXCEPTION (ExcFieldMeshDifference,<< "Two copies of the field "<< EI_Field::qval<< "call set_mesh with different arguments.\n")
 
FieldCommonname (const string &name)
 
FieldCommondescription (const string &description)
 
FieldCommoninput_default (const string &input_default)
 
FieldCommonunits (const UnitSI &units)
 Set basic units of the field. More...
 
FieldCommonset_limits (double min, double max=std::numeric_limits< double >::max())
 
FieldCommoninput_selection (Input::Type::Selection element_selection)
 
FieldCommonoutput_type (OutputTime::DiscreteSpace rt)
 
FieldCommonflags (FieldFlag::Flags::Mask mask)
 
FieldCommonflags_add (FieldFlag::Flags::Mask mask)
 
void set_components (const std::vector< string > &names)
 
virtual void set_mesh (const Mesh &mesh)
 
virtual void set_input_list (const Input::Array &list)=0
 
const std::string & input_name () const
 
const std::string & name () const
 
const std::string description () const
 
const std::string & input_default () const
 
const UnitSIunits () const
 
std::pair< double, double > limits () const
 
OutputTime::DiscreteSpace output_type () const
 
bool is_bc () const
 
unsigned int n_comp () const
 
const Meshmesh () const
 
FieldFlag::Flagsflags ()
 
double time () const
 
bool is_jump_time ()
 
unsigned int input_list_size () const
 
virtual bool is_constant (Region reg)=0
 
virtual FieldResult field_result (RegionSet region_set) const =0
 Indicates special field states. More...
 
virtual std::string get_value_attribute () const =0
 
bool changed () const
 
virtual IT::Instance get_input_type ()=0
 
virtual IT::Array get_multifield_input_type ()=0
 
void mark_input_times (const TimeGovernor &tg)
 
virtual bool set_time (const TimeStep &time, LimitSide limit_side)=0
 
virtual void copy_from (const FieldCommon &other)=0
 
virtual void field_output (std::shared_ptr< OutputTime > stream)=0
 
virtual void observe_output (std::shared_ptr< Observe > observe)=0
 
void set_component_index (unsigned int idx)
 
bool is_multifield () const
 
virtual ~FieldCommon ()
 

Static Public Member Functions

static IT::Record field_descriptor_record (const string &record_name)
 
static const std::string field_descriptor_record_description (const string &record_name)
 
static bool print_message_table (ostream &stream, std::string equation_name)
 

Protected Types

enum  TimeStatus { TimeStatus::changed, TimeStatus::constant, TimeStatus::unknown }
 

Protected Member Functions

 FieldCommon ()
 
 FieldCommon (const FieldCommon &other)
 
void set_history_changed ()
 

Protected Attributes

std::string name_
 
std::shared_ptr< SharedDatashared_
 
TimeStatus set_time_result_
 
double last_time_ = -numeric_limits<double>::infinity()
 
LimitSide last_limit_side_ = LimitSide::left
 
bool is_jump_time_
 
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...
 

Static Protected Attributes

static const unsigned int history_length_limit_ =3
 
static std::vector< MessageDatamessages_data_ = std::vector<FieldCommon::MessageData>()
 Vector of data of initialization messages. More...
 

Friends

std::ostream & operator<< (std::ostream &stream, const FieldCommon &field)
 

Detailed Description

Common abstract parent of all Field<...> classes.

We need common ancestor in order to keep a list of all fields in one EqData object and allow collective operations like set_time or init_from_input.

Definition at line 60 of file field_common.hh.

Member Enumeration Documentation

enum FieldCommon::TimeStatus
strongprotected

Result of last set time method

Enumerator
changed 
constant 
unknown 

Definition at line 562 of file field_common.hh.

Constructor & Destructor Documentation

FieldCommon::~FieldCommon ( )
virtual

Virtual destructor.

Definition at line 114 of file field_common.cc.

FieldCommon::FieldCommon ( )
protected

Private default constructor. Should be used only through Field<...>

Definition at line 26 of file field_common.cc.

FieldCommon::FieldCommon ( const FieldCommon other)
protected

Private copy constructor. Should be used only through Field<...>

Definition at line 41 of file field_common.cc.

Member Function Documentation

bool FieldCommon::changed ( ) const
inline

Returns true if set_time_result_ is not TimeStatus::constant. Returns the same value as last set_time method.

Definition at line 323 of file field_common.hh.

Here is the caller graph for this function:

virtual void FieldCommon::copy_from ( const FieldCommon other)
pure virtual

Check that other is instance of the same Field<..> class and perform assignment. Polymorphic copy.

The copy is performed only if *this have set flag 'input_copy'. If *this have set also the flag 'decare_input' the copy is performed only if the input_list is empty.

Implemented in Field< spacedim, Value >, Field< 3, FieldValue< 3 >::Scalar >, Field< 3, FieldValue< 3 >::VectorFixed >, Field< spacedim, FieldValue< 3 >::Scalar >, Field< spacedim, FieldValue< 3 >::Enum >, Field< 3, FieldValue< 3 >::Integer >, Field< 3, FieldValue< 3 >::TensorFixed >, MultiField< spacedim, Value >, MultiField< 3, FieldValue< 3 >::Scalar >, MultiField< spacedim, FieldValue< 3 >::Enum >, MultiField< spacedim, FieldValue< 3 >::Scalar >, MultiField< 3, FieldValue< 3 >::TensorFixed >, and MultiField< 3, FieldValue< 3 >::Enum >.

Here is the caller graph for this function:

FieldCommon::DECLARE_EXCEPTION ( ExcFieldMeshDifference  ,
<< "Two copies of the field "<< EI_Field::qval<< "call set_mesh with different arguments.\n"   
)
FieldCommon::DECLARE_INPUT_EXCEPTION ( ExcNonascendingTime  ,
<< "Non-ascending time: "<< EI_Time::val<< " for field "<< EI_Field::qval<< ".\n"   
)
FieldCommon::DECLARE_INPUT_EXCEPTION ( ExcMissingDomain  ,
<< "Missing domain specification (region or r_id) in the field descriptor:"   
)
FieldCommon& FieldCommon::description ( const string &  description)
inline

Set description of the field, used for description of corresponding key in documentation.

Definition at line 104 of file field_common.hh.

Here is the caller graph for this function:

const std::string FieldCommon::description ( ) const
inline

Definition at line 223 of file field_common.hh.

IT::Record FieldCommon::field_descriptor_record ( const string &  record_name)
static

Common part of the field descriptor. To get finished record one has to add keys for individual fields. This is done automatically using FieldSet::get_input_type().

Definition at line 55 of file field_common.cc.

Here is the caller graph for this function:

const std::string FieldCommon::field_descriptor_record_description ( const string &  record_name)
static

Create description of field descriptor record.

Definition at line 68 of file field_common.cc.

Here is the caller graph for this function:

virtual void FieldCommon::field_output ( std::shared_ptr< OutputTime stream)
pure virtual
virtual FieldResult FieldCommon::field_result ( RegionSet  region_set) const
pure virtual

Indicates special field states.

Extension of the previous method. Return possible values from the enum FieldResult, see description there. The initial state is field_none, if the field is correctly set on all regions of the region_set given as parameter we return state field_other or even more particular result.

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. Other possible results are: result_zeros, result_eye, result_ones, result_constant, result_other see FieldResult for explanation.

Multifield return most particular value that holds for all its subfields.

Implemented in Field< spacedim, Value >, Field< 3, FieldValue< 3 >::Scalar >, Field< 3, FieldValue< 3 >::VectorFixed >, Field< spacedim, FieldValue< 3 >::Scalar >, Field< spacedim, FieldValue< 3 >::Enum >, Field< 3, FieldValue< 3 >::Integer >, Field< 3, FieldValue< 3 >::TensorFixed >, MultiField< spacedim, Value >, MultiField< 3, FieldValue< 3 >::Scalar >, MultiField< spacedim, FieldValue< 3 >::Enum >, MultiField< spacedim, FieldValue< 3 >::Scalar >, MultiField< 3, FieldValue< 3 >::TensorFixed >, and MultiField< 3, FieldValue< 3 >::Enum >.

Here is the caller graph for this function:

FieldCommon& FieldCommon::flags ( FieldFlag::Flags::Mask  mask)
inline

Set given mask to the field flags, ignoring default setting. Default setting is declare_input & equation_input & allow_output.

Definition at line 168 of file field_common.hh.

Here is the caller graph for this function:

FieldFlag::Flags& FieldCommon::flags ( )
inline

Definition at line 252 of file field_common.hh.

Here is the caller graph for this function:

FieldCommon& FieldCommon::flags_add ( FieldFlag::Flags::Mask  mask)
inline

Add given mask to the field flags.

Definition at line 174 of file field_common.hh.

Here is the caller graph for this function:

virtual IT::Instance FieldCommon::get_input_type ( )
pure virtual
virtual IT::Array FieldCommon::get_multifield_input_type ( )
pure virtual
virtual std::string FieldCommon::get_value_attribute ( ) const
pure virtual
FieldCommon& FieldCommon::input_default ( const string &  input_default)
inline

Set default value for the field's key from which the default constant valued field will be constructed.

During the first call of the set_time method, we check that the field is defined on all regions. On regions where it is not set yet, we use given dflt string to get particular instance of FieldBase<> (see check_initialized_region_fields_). The default string is interpreted in the same way as if it appears in the input file as the value of the field. In particular it can be whole record with TYPE of the field etc. Most common choice is however mere constant.

Definition at line 116 of file field_common.hh.

Here is the caller graph for this function:

const std::string& FieldCommon::input_default ( ) const
inline

Definition at line 226 of file field_common.hh.

Here is the caller graph for this function:

unsigned int FieldCommon::input_list_size ( ) const
inline

Returns number of field descriptors containing the field.

Definition at line 273 of file field_common.hh.

const std::string& FieldCommon::input_name ( ) const
inline

Getters.

Definition at line 217 of file field_common.hh.

Here is the caller graph for this function:

FieldCommon& FieldCommon::input_selection ( Input::Type::Selection  element_selection)
inline

For the fields returning "Enum", we have to pass the Input::Type::Selection object to the field implementations.

We must save raw pointer since selection may not be yet initialized (during static initialization phase).

Definition at line 148 of file field_common.hh.

Here is the caller graph for this function:

bool FieldCommon::is_bc ( ) const
inline

Definition at line 243 of file field_common.hh.

Here is the caller graph for this function:

virtual bool FieldCommon::is_constant ( Region  reg)
pure virtual

If the field on given region reg exists and is of type FieldConstant<...> the method method returns true otherwise it returns false. Then one can call ElementAccessor<spacedim>(mesh(), reg ) to construct an ElementAccessor elm pointing to "virtual" element on which Field::value returns constant value. Unlike the Field<>::field_result method, this one provides no value, so it have common header (arguments, return type) and could be part of FieldCommon and FieldSet which is useful in some applications.

TODO:Current implementation use virtual functions and can be prohibitively slow if called for every element. If this becomes necessary it is possible to incorporate such test into set_time method and in this method just return precomputed result.

Implemented in Field< spacedim, Value >, Field< 3, FieldValue< 3 >::Scalar >, Field< 3, FieldValue< 3 >::VectorFixed >, Field< spacedim, FieldValue< 3 >::Scalar >, Field< spacedim, FieldValue< 3 >::Enum >, Field< 3, FieldValue< 3 >::Integer >, Field< 3, FieldValue< 3 >::TensorFixed >, MultiField< spacedim, Value >, MultiField< 3, FieldValue< 3 >::Scalar >, MultiField< spacedim, FieldValue< 3 >::Enum >, MultiField< spacedim, FieldValue< 3 >::Scalar >, MultiField< 3, FieldValue< 3 >::TensorFixed >, and MultiField< 3, FieldValue< 3 >::Enum >.

Here is the caller graph for this function:

bool FieldCommon::is_jump_time ( )
inline

Returns true if the field change algorithm for the current time set through the set_time method. This happen for all times in the field descriptors on the input of this particular field.

Definition at line 266 of file field_common.hh.

Here is the caller graph for this function:

bool FieldCommon::is_multifield ( ) const
inline

Return multifield_ flag. TODO: temporary solution

Definition at line 422 of file field_common.hh.

Here is the caller graph for this function:

std::pair<double, double> FieldCommon::limits ( ) const
inline

Definition at line 235 of file field_common.hh.

Here is the caller graph for this function:

void FieldCommon::mark_input_times ( const TimeGovernor tg)

Pass through the input array input_list_, collect all times where the field could change and put appropriate time marks into global TimeMarks object. Introduced time marks have both given mark_type and type_input() type.

Further development:

  • we have to distinguish "jump" times and "smooth" times

Definition at line 77 of file field_common.cc.

Here is the caller graph for this function:

const Mesh* FieldCommon::mesh ( ) const
inline

Definition at line 249 of file field_common.hh.

Here is the caller graph for this function:

unsigned int FieldCommon::n_comp ( ) const
inline

Definition at line 246 of file field_common.hh.

Here is the caller graph for this function:

FieldCommon& FieldCommon::name ( const string &  name)
inline

Set name of the field. In fact there are two attributes set by this method.

The first is name used to identify the field as part of a FieldSet or MultiField objects. This name is permanent and can be set only by this method. Can be accessed by name() method. This name is also used at output.

The second is input_name_ that determines appropriate key name in the input field descriptor. This name is also set by this method, but is stored in the internal shared space which is overwritten during call of copy_from method or assignment operator. Can be accessed by input_name() mathod.

Definition at line 97 of file field_common.hh.

Here is the caller graph for this function:

const std::string& FieldCommon::name ( ) const
inline

Definition at line 220 of file field_common.hh.

Here is the caller graph for this function:

virtual void FieldCommon::observe_output ( std::shared_ptr< Observe observe)
pure virtual
FieldCommon& FieldCommon::output_type ( OutputTime::DiscreteSpace  rt)
inline

Output discrete space used in the output() method. Can be different for different field copies. one can choose between: data constant on elements, linear data given in nodes, and discontinuous linear data.

If not set explicitly by this method, the default value is OutputTime::ELEM_DATA

Definition at line 161 of file field_common.hh.

Here is the caller graph for this function:

OutputTime::DiscreteSpace FieldCommon::output_type ( ) const
inline

Definition at line 240 of file field_common.hh.

Here is the caller graph for this function:

bool FieldCommon::print_message_table ( ostream &  stream,
std::string  equation_name 
)
static

Print stored messages to table.

Return true if messages_data_ vector is nonempty and clear its.

Definition at line 91 of file field_common.cc.

Here is the caller graph for this function:

void FieldCommon::set_component_index ( unsigned int  idx)
inline

Sets component_index_

Definition at line 413 of file field_common.hh.

Here is the caller graph for this function:

void FieldCommon::set_components ( const std::vector< string > &  names)
inline

Set vector of component names. Set number of components for run-time sized vectors. This is used latter when we construct objects derived from FieldBase<...>.

n_comp_ is constant zero for fixed values, this zero is set by Field<...> constructors

Definition at line 184 of file field_common.hh.

Here is the caller graph for this function:

void FieldCommon::set_history_changed ( )
inlineprotected

Invalidate last time in order to force set_time method update region_fields_.

Definition at line 457 of file field_common.hh.

Here is the caller graph for this function:

virtual void FieldCommon::set_input_list ( const Input::Array list)
pure virtual

Set the data list from which field will read its input. It is list of "field descriptors". When reading from the input list we consider only field descriptors containing key of named by the field name. These field descriptors has to have times forming ascending sequence.

The list is used by set_time method to set field on individual regions to actual FieldBase descendants.

Implemented in Field< spacedim, Value >, Field< 3, FieldValue< 3 >::Scalar >, Field< 3, FieldValue< 3 >::VectorFixed >, Field< spacedim, FieldValue< 3 >::Scalar >, Field< spacedim, FieldValue< 3 >::Enum >, Field< 3, FieldValue< 3 >::Integer >, Field< 3, FieldValue< 3 >::TensorFixed >, MultiField< spacedim, Value >, MultiField< 3, FieldValue< 3 >::Scalar >, MultiField< spacedim, FieldValue< 3 >::Enum >, MultiField< spacedim, FieldValue< 3 >::Scalar >, MultiField< 3, FieldValue< 3 >::TensorFixed >, and MultiField< 3, FieldValue< 3 >::Enum >.

Here is the caller graph for this function:

FieldCommon& FieldCommon::set_limits ( double  min,
double  max = std::numeric_limits<double>::max() 
)
inline

Set limits of value of the field.

Definition at line 135 of file field_common.hh.

Here is the caller graph for this function:

virtual void FieldCommon::set_mesh ( const Mesh mesh)
inlinevirtual
virtual bool FieldCommon::set_time ( const TimeStep time,
LimitSide  limit_side 
)
pure virtual

Abstract method to update field to the new time level. Implemented by in class template Field<...>.

Return true if the value of the field was changed on some region. The returned value is also stored in changed_during_set_time data member.

Default values helps when creating steady field. Note that default TimeGovernor constructor set time to 0.0.

Different field copies can be set to different times.

TODO: update following: Set side of limit when calling set_time with jump time, i.e. time where the field change implementation on some region. Wee assume that implementations prescribe only smooth fields. This method invalidate result of changed() so it should be called just before set_time. Can be different for different field copies.

Implemented in Field< spacedim, Value >, Field< 3, FieldValue< 3 >::Scalar >, Field< 3, FieldValue< 3 >::VectorFixed >, Field< spacedim, FieldValue< 3 >::Scalar >, Field< spacedim, FieldValue< 3 >::Enum >, Field< 3, FieldValue< 3 >::Integer >, Field< 3, FieldValue< 3 >::TensorFixed >, MultiField< spacedim, Value >, MultiField< 3, FieldValue< 3 >::Scalar >, MultiField< spacedim, FieldValue< 3 >::Enum >, MultiField< spacedim, FieldValue< 3 >::Scalar >, MultiField< 3, FieldValue< 3 >::TensorFixed >, and MultiField< 3, FieldValue< 3 >::Enum >.

Here is the caller graph for this function:

double FieldCommon::time ( ) const
inline

Returns time set by last call of set_time method. Can be different for different field copies.

Definition at line 259 of file field_common.hh.

Here is the caller graph for this function:

FieldCommon::TYPEDEF_ERR_INFO ( EI_Time  ,
double   
)
FieldCommon::TYPEDEF_ERR_INFO ( EI_Field  ,
std::string   
)
FieldCommon& FieldCommon::units ( const UnitSI units)
inline

Set basic units of the field.

Currently, we use it only during output and we represents units just by a string.

TODO: Particular class for representing and conversion of various units would be more appropriate. This can allow specification of the units on the inptu, automatic conversion and the same on the output. Possibly this allow using Boost::Units library, however, it seems to introduce lot of boilerplate code. But can increase correctness of the calculations.

Definition at line 129 of file field_common.hh.

Here is the caller graph for this function:

const UnitSI& FieldCommon::units ( ) const
inline

Definition at line 229 of file field_common.hh.

Here is the caller graph for this function:

Friends And Related Function Documentation

std::ostream& operator<< ( std::ostream &  stream,
const FieldCommon field 
)
friend

Stream output operator

Definition at line 618 of file field_common.hh.

Member Data Documentation

unsigned int FieldCommon::component_index_
protected

Specify if the field is part of a MultiField and which component it is

Definition at line 596 of file field_common.hh.

Field flags. Default setting is "an equation input field, that can read from user input, and can be written to output".

Definition at line 610 of file field_common.hh.

const unsigned int FieldCommon::history_length_limit_ =3
staticprotected

Maximum number of FieldBase objects we store per one region.

Definition at line 607 of file field_common.hh.

bool FieldCommon::is_jump_time_
protected

Set to true by the set_time method the field algorithm change on any region. Accessible through the is_jump_time method.

Definition at line 586 of file field_common.hh.

LimitSide FieldCommon::last_limit_side_ = LimitSide::left
protected

Definition at line 580 of file field_common.hh.

double FieldCommon::last_time_ = -numeric_limits<double>::infinity()
protected

Last set time. Can be different for different field copies. Store also time limit, since the field may be discontinuous.

Definition at line 579 of file field_common.hh.

std::vector< FieldCommon::MessageData > FieldCommon::messages_data_ = std::vector<FieldCommon::MessageData>()
staticprotected

Vector of data of initialization messages.

Definition at line 613 of file field_common.hh.

bool FieldCommon::multifield_
protected

Flag determining if object is Multifield or Field. TODO: temporary solution, goal is to make these two classes to behave similarly

Definition at line 602 of file field_common.hh.

std::string FieldCommon::name_
protected

Name that identifies the field in the field_set. By default this is same as shared_->input_name_.

Definition at line 552 of file field_common.hh.

TimeStatus FieldCommon::set_time_result_
protected

Status of history.

Definition at line 573 of file field_common.hh.

std::shared_ptr<SharedData> FieldCommon::shared_
protected

Data shared among copies of the same input field.

Definition at line 557 of file field_common.hh.

OutputTime::DiscreteSpace FieldCommon::type_of_output_data_ = OutputTime::ELEM_DATA
protected

Output data type used in the output() method. Can be different for different field copies.

Definition at line 591 of file field_common.hh.


The documentation for this class was generated from the following files: