MoFEM::Field Struct Reference

Provide data structure for (tensor) field approximation. More...

#include <src/multi_indices/FieldMultiIndices.hpp>

Public Types
using	SequenceDofContainer = multi_index_container< boost::weak_ptr< std::vector< DofEntity > >, indexed_by< sequenced<> >>
typedef std::array< std::array< int, MAX_DOFS_ON_ENTITY >, MBMAXTYPE >	DofsOrderMap

Public Member Functions
	Field (const moab::Interface &moab, const EntityHandle meshset, const boost::shared_ptr< CoordSys > coord_sys_ptr)
	constructor for moab field More...
virtual	~Field ()=default
FieldOrderTable &	getFieldOrderTable ()
	Get the Field Order Table. More...
UId	generateGlobalUniqueIdForTypeLo (const EntityType type, const int owner_proc) const
UId	generateGlobalUniqueIdForTypeHi (const EntityType type, const int owner_proc) const
EntityHandle	getMeshset () const
	Get field meshset. More...
int	getCoordSysDim (const int d=0) const
	Get dimension of general two-point tensor MoFEM::CoordSys::getDim. More...
MoFEMErrorCode	get_E_Base (const double m[]) const
	Get reference base vectors. More...
MoFEMErrorCode	get_E_DualBase (const double m[]) const
	Get reference dual base vectors. More...
MoFEMErrorCode	get_e_Base (const double m[]) const
	Get current dual base vectors. More...
MoFEMErrorCode	get_e_DualBase (const double m[]) const
	Get current dual base vectors. More...
EntityHandle	getCoordSysMeshSet () const
	Returns meshset on which Tags defining coordinate system are stored. More...
std::string	getCoordSysName () const
	Get coordinate system name. More...
boost::string_ref	getCoordSysNameRef () const
	Get coordinate system name. More...
const BitFieldId &	getId () const
	Get unique field id. More...
boost::string_ref	getNameRef () const
	Get string reference to field name. More...
std::string	getName () const
	Get field name. More...
FieldSpace	getSpace () const
	Get field approximation space. More...
FieldApproximationBase	getApproxBase () const
	Get approximation base. More...
FieldCoefficientsNumber	getNbOfCoeffs () const
	Get number of field coefficients. More...
FieldBitNumber	getBitNumber () const
	Get number of set bit in Field ID. Each field has uid, get getBitNumber get number of bit set for given field. Field ID has only one bit set for each field. More...
FieldBitNumber	getBitNumberCalculate () const
	Calculate number of set bit in Field ID. Each field has uid, get getBitNumber get number of bit set for given field. Field ID has only one bit set for each field. More...
SequenceDofContainer &	getDofSequenceContainer () const
	Get reference to sequence data container. More...
std::array< ApproximationOrder, MAX_DOFS_ON_ENTITY > &	getDofOrderMap (const EntityType type) const
	get hash-map relating dof index on entity with its order More...
DofsOrderMap &	getDofOrderMap () const
	get hash-map relating dof index on entity with its order More...
MoFEMErrorCode	rebuildDofsOrderMap () const
const Field *	getFieldRawPtr () const

Static Public Member Functions
static UId	generateGlobalUniqueIdForTypeLo (const char bit_number, const EntityType type, const int owner_proc)
static UId	generateGlobalUniqueIdForTypeHi (const char bit_number, const EntityType type, const int owner_proc)
static FieldBitNumber	getBitNumberCalculate (const BitFieldId &id)
	Calculate number of set bit in Field ID. Each field has uid, get getBitNumber get number of bit set for given field. Field ID has only one bit set for each field. More...

Public Attributes
moab::Interface &	moab
EntityHandle	meshSet
	keeps entities for this meshset More...
boost::shared_ptr< CoordSys >	coordSysPtr
	Pointer to field coordinate system data structure. More...
TagType	th_FieldDataVertsType
Tag	th_FieldDataVerts
	Tag storing field values on vertices in the field. More...
Tag	th_FieldData
	Tag storing field values on entity in the field. More...
Tag	th_AppOrder
	Tag storing approximation order on entity. More...
BitFieldId *	tagId
	tag keeps field id More...
FieldSpace *	tagSpaceData
	tag keeps field space More...
FieldApproximationBase *	tagBaseData
	tag keeps field base More...
FieldCoefficientsNumber *	tagNbCoeffData
const void *	tagName
	tag keeps name of the field More...
int	tagNameSize
	number of bits necessary to keep field name More...
const void *	tagNamePrefixData
	tag keeps name prefix of the field More...
int	tagNamePrefixSize
FieldOrderTable	forderTable
	nb. DOFs table for entities More...
unsigned int	bitNumber

Private Attributes
SequenceDofContainer	sequenceDofContainer
DofsOrderMap	dofOrderMap

Friends
std::ostream &	operator<< (std::ostream &os, const Field &e)

Detailed Description

Provide data structure for (tensor) field approximation.

The Field is intended to provide support for fields, with a strong bias towards supporting first and best the capabilities required for scientific computing applications. Since we work with discrete spaces, data structure has to carry information about type of approximation space, its regularity.

Field data structure storing information about space, approximation base, coordinate systems, etc. It stores additional data needed for book keeping, like tags to data on the mesh.

Each filed has unique ID and name. This data structure is shared between entities on which is spanning and DOFs on those entities.

Definition at line 64 of file FieldMultiIndices.hpp.

Member Typedef Documentation

DofsOrderMap

typedef std::array<std::array<int, MAX_DOFS_ON_ENTITY>, MBMAXTYPE> MoFEM::Field::DofsOrderMap

Definition at line 83 of file FieldMultiIndices.hpp.

SequenceDofContainer

using MoFEM::Field::SequenceDofContainer = multi_index_container< boost::weak_ptr<std::vector<DofEntity> >, indexed_by<sequenced<> >>

Definition at line 76 of file FieldMultiIndices.hpp.

Constructor & Destructor Documentation

Field()

MoFEM::Field::Field	(	const moab::Interface &	moab,
		const EntityHandle	meshset,
		const boost::shared_ptr< CoordSys >	coord_sys_ptr
	)

constructor for moab field

Parameters

meshset

which keeps entities for this field

Definition at line 29 of file FieldMultiIndices.cpp.

    : moab(const_cast<moab::Interface &>(moab)), meshSet(meshset),
      coordSysPtr(coord_sys_ptr), tagId(NULL), tagSpaceData(NULL),
      tagNbCoeffData(NULL), tagName(NULL), tagNameSize(0) {
 
  auto get_tag_data_ptr = [&](const auto name, auto &tag_data) {
    MoFEMFunctionBegin;
    Tag th;
    CHKERR moab.tag_get_handle(name, th);
    CHKERR moab.tag_get_by_ptr(th, &meshset, 1, (const void **)&tag_data);
    MoFEMFunctionReturn(0);
  };
 
  // id
  ierr = get_tag_data_ptr("_FieldId", tagId);
  CHKERRABORT(PETSC_COMM_SELF, ierr);
  // space
  ierr = get_tag_data_ptr("_FieldSpace", tagSpaceData);
  CHKERRABORT(PETSC_COMM_SELF, ierr);
 
  // approx. base
  ierr = get_tag_data_ptr("_FieldBase", tagBaseData);
  CHKERRABORT(PETSC_COMM_SELF, ierr);
 
  auto get_all_tags = [&]() {
    MoFEMFunctionBegin;
    // name
    Tag th_field_name;
    CHKERR moab.tag_get_handle("_FieldName", th_field_name);
    CHKERR moab.tag_get_by_ptr(th_field_name, &meshSet, 1,
                               (const void **)&tagName, &tagNameSize);
    // name prefix
    Tag th_field_name_data_name_prefix;
    CHKERR moab.tag_get_handle("_FieldName_DataNamePrefix",
                               th_field_name_data_name_prefix);
    CHKERR moab.tag_get_by_ptr(th_field_name_data_name_prefix, &meshSet, 1,
                               (const void **)&tagNamePrefixData,
                               &tagNamePrefixSize);
    std::string name_data_prefix((char *)tagNamePrefixData, tagNamePrefixSize);
    // data
    std::string tag_data_name = name_data_prefix + getName();
    CHKERR moab.tag_get_handle(tag_data_name.c_str(), th_FieldData);
    std::string tag_data_name_verts = name_data_prefix + getName() + "V";
    CHKERR moab.tag_get_handle(tag_data_name_verts.c_str(), th_FieldDataVerts);
    CHKERR moab.tag_get_type(th_FieldDataVerts, th_FieldDataVertsType);
    // order
    std::string tag_approximation_order_name = "_App_Order_" + getName();
    CHKERR moab.tag_get_handle(tag_approximation_order_name.c_str(),
                               th_AppOrder);
    // rank
    Tag th_rank;
    std::string Tag_rank_name = "_Field_Rank_" + getName();
    CHKERR moab.tag_get_handle(Tag_rank_name.c_str(), th_rank);
    CHKERR moab.tag_get_by_ptr(th_rank, &meshSet, 1,
                               (const void **)&tagNbCoeffData);
    MoFEMFunctionReturn(0);
  };
 
  ierr = get_all_tags();
  CHKERRABORT(PETSC_COMM_SELF, ierr);
 
  bitNumber = getBitNumberCalculate();
 
  auto reset_entity_order_table = [&]() {
    for (int tt = 0; tt != MBMAXTYPE; ++tt)
      forderTable[tt] = NULL;
  };
 
  auto set_entity_order_table = [&]() {
    switch (*tagBaseData) {
    case AINSWORTH_LEGENDRE_BASE:
    case AINSWORTH_LOBATTO_BASE:
      switch (*tagSpaceData) {
      case H1:
        forderTable[MBVERTEX] = [](int P) -> int { return (P > 0) ? 1 : 0; };
        forderTable[MBEDGE] = [](int P) -> int { return NBEDGE_H1(P); };
        forderTable[MBTRI] = [](int P) -> int { return NBFACETRI_H1(P); };
        forderTable[MBQUAD] = [](int P) -> int { return NBFACEQUAD_H1(P); };
        forderTable[MBTET] = [](int P) -> int { return NBVOLUMETET_H1(P); };
        forderTable[MBPRISM] = [](int P) -> int { return NBVOLUMEPRISM_H1(P); };
        break;
      case HCURL:
        forderTable[MBVERTEX] = [](int P) -> int {
          (void)P;
          return 0;
        };
        forderTable[MBEDGE] = [](int P) -> int {
          return NBEDGE_AINSWORTH_HCURL(P);
        };
        forderTable[MBTRI] = [](int P) -> int {
          return NBFACETRI_AINSWORTH_HCURL(P);
        };
        forderTable[MBTET] = [](int P) -> int {
          return NBVOLUMETET_AINSWORTH_HCURL(P);
        };
        break;
      case HDIV:
        forderTable[MBVERTEX] = [](int P) -> int {
          (void)P;
          return 0;
        };
        forderTable[MBEDGE] = [](int P) -> int {
          (void)P;
          return NBEDGE_HDIV(P);
        };
        forderTable[MBTRI] = [](int P) -> int {
          return NBFACETRI_AINSWORTH_HDIV(P);
        };
        forderTable[MBTET] = [](int P) -> int {
          return NBVOLUMETET_AINSWORTH_HDIV(P);
        };
        break;
      case L2:
        forderTable[MBVERTEX] = [](int P) -> int {
          (void)P;
          return 1;
        };
        forderTable[MBEDGE] = [](int P) -> int { return NBEDGE_L2(P); };
        forderTable[MBTRI] = [](int P) -> int { return NBFACETRI_L2(P); };
        forderTable[MBQUAD] = [](int P) -> int { return NBFACEQUAD_L2(P); };
        forderTable[MBTET] = [](int P) -> int { return NBVOLUMETET_L2(P); };
        break;
      default:
        THROW_MESSAGE("unknown approximation space");
      }
      break;
    case AINSWORTH_BERNSTEIN_BEZIER_BASE:
      switch (*tagSpaceData) {
      case H1:
        forderTable[MBVERTEX] = [](int P) -> int { return (P > 0) ? 1 : 0; };
        forderTable[MBEDGE] = [](int P) -> int { return NBEDGE_H1(P); };
        forderTable[MBTRI] = [](int P) -> int { return NBFACETRI_H1(P); };
        forderTable[MBQUAD] = [](int P) -> int { return NBFACEQUAD_H1(P); };
        forderTable[MBTET] = [](int P) -> int { return NBVOLUMETET_H1(P); };
        forderTable[MBPRISM] = [](int P) -> int { return NBVOLUMEPRISM_H1(P); };
        break;
      case L2:
        forderTable[MBVERTEX] = [](int P) -> int {
          (void)P;
          return 1;
        };
        forderTable[MBEDGE] = [](int P) -> int { return NBEDGE_L2(P); };
        forderTable[MBTRI] = [](int P) -> int { return NBFACETRI_L2(P); };
        forderTable[MBQUAD] = [](int P) -> int { return NBFACEQUAD_L2(P); };
        forderTable[MBTET] = [](int P) -> int { return NBVOLUMETET_L2(P); };
        break;
      default:
        THROW_MESSAGE("unknown approximation space or not yet implemented");
      }
      break;
    case DEMKOWICZ_JACOBI_BASE:
      switch (*tagSpaceData) {
      case H1:
        forderTable[MBVERTEX] = [](int P) -> int { return (P > 0) ? 1 : 0; };
        forderTable[MBEDGE] = [](int P) -> int { return NBEDGE_H1(P); };
        forderTable[MBTRI] = [](int P) -> int { return NBFACETRI_H1(P); };
        forderTable[MBQUAD] = [](int P) -> int { return NBFACEQUAD_H1(P); };
        forderTable[MBTET] = [](int P) -> int { return NBVOLUMETET_H1(P); };
        forderTable[MBPRISM] = [](int P) -> int { return NBVOLUMEPRISM_H1(P); };
        break;
      case HCURL:
        forderTable[MBVERTEX] = [](int P) -> int {
          (void)P;
          return 0;
        };
        forderTable[MBEDGE] = [](int P) -> int {
          return NBEDGE_DEMKOWICZ_HCURL(P);
        };
        forderTable[MBTRI] = [](int P) -> int {
          return NBFACETRI_DEMKOWICZ_HCURL(P);
        };
        forderTable[MBQUAD] = [](int P) -> int {
          return NBFACEQUAD_DEMKOWICZ_HCURL(P);
        };
        forderTable[MBTET] = [](int P) -> int {
          return NBVOLUMETET_DEMKOWICZ_HCURL(P);
        };
        break;
      case HDIV:
        forderTable[MBVERTEX] = [](int P) -> int {
          (void)P;
          return 0;
        };
        forderTable[MBEDGE] = [](int P) -> int {
          (void)P;
          return 0;
        };
        forderTable[MBTRI] = [](int P) -> int {
          return NBFACETRI_DEMKOWICZ_HDIV(P);
        };
        forderTable[MBTET] = [](int P) -> int {
          return NBVOLUMETET_DEMKOWICZ_HDIV(P);
        };
        break;
      case L2:
        forderTable[MBVERTEX] = [](int P) -> int {
          (void)P;
          return 1;
        };
        forderTable[MBEDGE] = [](int P) -> int { return NBEDGE_L2(P); };
        forderTable[MBTRI] = [](int P) -> int { return NBFACETRI_L2(P); };
        forderTable[MBQUAD] = [](int P) -> int { return NBFACEQUAD_L2(P); };
        forderTable[MBTET] = [](int P) -> int { return NBVOLUMETET_L2(P); };
        break;
      default:
        THROW_MESSAGE("unknown approximation space or not yet implemented");
      }
      break;
    case USER_BASE:
      for (int ee = 0; ee < MBMAXTYPE; ee++) {
        forderTable[ee] = [](int P) -> int {
          (void)P;
          return 0;
        };
      }
      break;
    default:
      if (*tagSpaceData != NOFIELD) {
        THROW_MESSAGE("unknown approximation base");
      } else {
        for (EntityType t = MBVERTEX; t < MBMAXTYPE; t++)
          forderTable[t] = [](int P) -> int {
            (void)P;
            return 1;
          };
      }
    }
  };
 
  reset_entity_order_table();
  set_entity_order_table();
  ierr = rebuildDofsOrderMap();
  CHKERRABORT(PETSC_COMM_SELF, ierr);
};

~Field()

virtual MoFEM::Field::~Field ( )

virtualdefault

Member Function Documentation

generateGlobalUniqueIdForTypeHi() [1/2]

static UId MoFEM::Field::generateGlobalUniqueIdForTypeHi ( const char  bit_number, const EntityType  type, const int  owner_proc  )

static

Definition at line 138 of file FieldMultiIndices.hpp.

                                                                   {
    constexpr int ent_shift = 8 * sizeof(EntityHandle);
    return (static_cast<UId>(type) << MB_ID_WIDTH |
            static_cast<UId>(bit_number) << ent_shift |
            static_cast<UId>(owner_proc) << 5 + ent_shift)
           << 9;
  }

generateGlobalUniqueIdForTypeHi() [2/2]

UId MoFEM::Field::generateGlobalUniqueIdForTypeHi	(	const EntityType	type,
		const int	owner_proc
	)		const

Definition at line 148 of file FieldMultiIndices.hpp.

                                                                  {
    return generateGlobalUniqueIdForTypeHi(bitNumber, type, owner_proc);
  }

generateGlobalUniqueIdForTypeLo() [1/2]

static UId MoFEM::Field::generateGlobalUniqueIdForTypeLo ( const char  bit_number, const EntityType  type, const int  owner_proc  )

static

Definition at line 123 of file FieldMultiIndices.hpp.

                                                                   {
    constexpr int ent_shift = 8 * sizeof(EntityHandle);
    return (static_cast<UId>(type) << MB_ID_WIDTH |
            static_cast<UId>(bit_number) << 8 * sizeof(EntityHandle) |
            static_cast<UId>(owner_proc) << 5 + ent_shift)
           << 9;
  }

generateGlobalUniqueIdForTypeLo() [2/2]

UId MoFEM::Field::generateGlobalUniqueIdForTypeLo	(	const EntityType	type,
		const int	owner_proc
	)		const

Definition at line 133 of file FieldMultiIndices.hpp.

                                                                  {
    return generateGlobalUniqueIdForTypeLo(bitNumber, type, owner_proc);
  }

get_E_Base()

MoFEMErrorCode MoFEM::Field::get_E_Base

(

const double

m[]

)

const

Get reference base vectors.

Parameters

Array

where coefficients (covariant) are returned

Returns

Error code

Definition at line 185 of file FieldMultiIndices.hpp.

                                                           {
    MoFEMFunctionBeginHot;
    MoFEMFunctionReturnHot(coordSysPtr->get_E_Base(m));
  }

get_e_Base()

MoFEMErrorCode MoFEM::Field::get_e_Base

(

const double

m[]

)

const

Get current dual base vectors.

Parameters

Array

where coefficients (covariant) are returned

Returns

Error code

Definition at line 205 of file FieldMultiIndices.hpp.

                                                           {
    MoFEMFunctionBeginHot;
    MoFEMFunctionReturnHot(coordSysPtr->get_e_Base(m));
  }

get_E_DualBase()

MoFEMErrorCode MoFEM::Field::get_E_DualBase

(

const double

m[]

)

const

Get reference dual base vectors.

Parameters

Array

where coefficients (contravariant) are returned

Returns

Error code

Definition at line 195 of file FieldMultiIndices.hpp.

                                                               {
    MoFEMFunctionBeginHot;
    MoFEMFunctionReturnHot(coordSysPtr->get_E_DualBase(m));
  }

get_e_DualBase()

MoFEMErrorCode MoFEM::Field::get_e_DualBase

(

const double

m[]

)

const

Get current dual base vectors.

Parameters

Array

where coefficients (covariant) are returned

Returns

Error code

Definition at line 215 of file FieldMultiIndices.hpp.

                                                               {
    MoFEMFunctionBeginHot;
    MoFEMFunctionReturnHot(coordSysPtr->get_e_DualBase(m));
  }

getApproxBase()

FieldApproximationBase MoFEM::Field::getApproxBase

(

)

const

Get approximation base.

Returns

Approximation base

Definition at line 274 of file FieldMultiIndices.hpp.

{ return *tagBaseData; }

getBitNumber()

FieldBitNumber MoFEM::Field::getBitNumber

(

)

const

Get number of set bit in Field ID. Each field has uid, get getBitNumber get number of bit set for given field. Field ID has only one bit set for each field.

Definition at line 299 of file FieldMultiIndices.hpp.

{ return bitNumber; }

getBitNumberCalculate() [1/2]

FieldBitNumber MoFEM::Field::getBitNumberCalculate

(

)

const

Calculate number of set bit in Field ID. Each field has uid, get getBitNumber get number of bit set for given field. Field ID has only one bit set for each field.

Definition at line 320 of file FieldMultiIndices.hpp.

                                                      {
    return getBitNumberCalculate(static_cast<BitFieldId &>(*tagId));
  }

getBitNumberCalculate() [2/2]

static FieldBitNumber MoFEM::Field::getBitNumberCalculate ( const BitFieldId &  id )

static

Calculate number of set bit in Field ID. Each field has uid, get getBitNumber get number of bit set for given field. Field ID has only one bit set for each field.

Definition at line 306 of file FieldMultiIndices.hpp.

                                                                           {
    static_assert(BITFIELDID_SIZE >= 32,
                  "Too many fields allowed, can be more but ...");
    FieldBitNumber b = ffsl(id.to_ulong());
    if (b != 0)
      return b;
    return 0;
  }

getCoordSysDim()

int MoFEM::Field::getCoordSysDim

(

const int

d = 0

)

const

Get dimension of general two-point tensor MoFEM::CoordSys::getDim.

See details here MoFEM::CoordSys::getDim

Definition at line 176 of file FieldMultiIndices.hpp.

                                                   {
    return coordSysPtr->getDim(d);
  }

getCoordSysMeshSet()

EntityHandle MoFEM::Field::getCoordSysMeshSet

(

)

const

Returns meshset on which Tags defining coordinate system are stored.

Returns

Coordinate system EntityHandle

Definition at line 224 of file FieldMultiIndices.hpp.

                                                 {
    return coordSysPtr->getMeshset();
  }

getCoordSysName()

std::string MoFEM::Field::getCoordSysName

(

)

const

Get coordinate system name.

Returns

Coordinate system name

Definition at line 232 of file FieldMultiIndices.hpp.

{ return coordSysPtr->getName(); }

getCoordSysNameRef()

boost::string_ref MoFEM::Field::getCoordSysNameRef

(

)

const

Get coordinate system name.

Returns

Return string_ref with name.

Definition at line 238 of file FieldMultiIndices.hpp.

                                                    {
    return coordSysPtr->getNameRef();
  }

getDofOrderMap() [1/2]

DofsOrderMap& MoFEM::Field::getDofOrderMap

(

)

const

get hash-map relating dof index on entity with its order

Dofs of given field are indexed on entity of the same type, same space, approximation base and number of coefficients, are sorted in the way.

Definition at line 368 of file FieldMultiIndices.hpp.

                                              {
    return const_cast<DofsOrderMap &>(dofOrderMap);
  }

getDofOrderMap() [2/2]

std::array<ApproximationOrder, MAX_DOFS_ON_ENTITY>& MoFEM::Field::getDofOrderMap

(

const EntityType

type

)

const

get hash-map relating dof index on entity with its order

Dofs of given field are indexed on entity of the same type, same space, approximation base and number of coefficients, are sorted in the way.

Definition at line 356 of file FieldMultiIndices.hpp.

                                              {
    return dofOrderMap[type];
  }

getDofSequenceContainer()

SequenceDofContainer& MoFEM::Field::getDofSequenceContainer

(

)

const

Get reference to sequence data container.

In sequence data container data are physically stored. The purpose of this is to allocate DofEntity data in bulk, having only one allocation instead each time entity is inserted. That makes code efficient.

The vector in sequence is destroyed if last entity inside that vector is destroyed. All MoFEM::MoFEMEntities have aliased shared_ptr which points to the vector.

Not all DOFs are starred in this way, currently such cases are considered;

• DOFs on vertices. That is exploited that for H1 space, there is some fixed number of DOFs on each vertex

For other cases, DOFs are stored locally in each MoFEM::MoFEMEntities.

Returns

MoFEM::Field::SequenceDofContainer

Definition at line 343 of file FieldMultiIndices.hpp.

                                                               {
    return sequenceDofContainer;
  }

getFieldOrderTable()

FieldOrderTable& MoFEM::Field::getFieldOrderTable

(

)

Get the Field Order Table.

Returns

FieldOrderTable&

Definition at line 115 of file FieldMultiIndices.hpp.

{ return forderTable; }

getFieldRawPtr()

const Field* MoFEM::Field::getFieldRawPtr

(

)

const

Definition at line 376 of file FieldMultiIndices.hpp.

{ return this; };

getId()

const BitFieldId& MoFEM::Field::getId

(

)

const

Get unique field id.

Returns

Filed ID

Definition at line 246 of file FieldMultiIndices.hpp.

{ return *((BitFieldId *)tagId); }

getMeshset()

EntityHandle MoFEM::Field::getMeshset

(

)

const

Get field meshset.

To meshsets entity are attached Tags which keeps basic information about field. Those information is field name, approximation base, approximation space, id, etc.

In meshset contains entities on which given filed is sparing. Type of entities depended on approximations space.

Returns

EntityHandle

Definition at line 167 of file FieldMultiIndices.hpp.

{ return meshSet; }

getName()

std::string MoFEM::Field::getName

(

)

const

Get field name.

Returns

Field name

Definition at line 260 of file FieldMultiIndices.hpp.

                                   {
    return std::string((char *)tagName, tagNameSize);
  }

getNameRef()

boost::string_ref MoFEM::Field::getNameRef

(

)

const

Get string reference to field name.

Returns

Field name

Definition at line 252 of file FieldMultiIndices.hpp.

                                            {
    return boost::string_ref((char *)tagName, tagNameSize);
  }

getNbOfCoeffs()

FieldCoefficientsNumber MoFEM::Field::getNbOfCoeffs

(

)

const

Get number of field coefficients.

Scalar field has only one coefficient, vector field in 3D has three. In general number determine space needed to keep data on entities. What coefficient means depend on interpretation and associated coordinate system. For example 3 coefficient means could be covariant or contravariant, or mean three temperatures for mixture of solid, air and water, etc.

Definition at line 290 of file FieldMultiIndices.hpp.

                                                       {
    return *tagNbCoeffData;
  };

getSpace()

FieldSpace MoFEM::Field::getSpace

(

)

const

Get field approximation space.

Returns

approximation space

Definition at line 268 of file FieldMultiIndices.hpp.

{ return *tagSpaceData; }

rebuildDofsOrderMap()

MoFEMErrorCode MoFEM::Field::rebuildDofsOrderMap

(

)

const

Definition at line 265 of file FieldMultiIndices.cpp.

                                                {
  MoFEMFunctionBegin;
 
  for (auto t = MBVERTEX; t != MBMAXTYPE; ++t) {
 
    int DD = 0;
    int nb_last_order_dofs = 0;
    const int rank = (*tagNbCoeffData);
 
    if (forderTable[t]) {
 
      for (int oo = 0; oo < MAX_DOFS_ON_ENTITY; ++oo) {
 
        const int nb_order_dofs = forderTable[t](oo);
        const int diff_oo = nb_order_dofs - nb_last_order_dofs;
        if (diff_oo >= 0) {
 
          if ((DD + rank * diff_oo) < MAX_DOFS_ON_ENTITY)
            for (int dd = 0; dd < diff_oo; ++dd)
              for (int rr = 0; rr != rank; ++rr, ++DD)
                dofOrderMap[t][DD] = oo;
          else
            break;
 
          nb_last_order_dofs = nb_order_dofs;
 
        } else {
          break;
        }
      }
    }
 
    std::fill(&dofOrderMap[t][DD], dofOrderMap[t].end(), -1);
  }
 
  MoFEMFunctionReturn(0);
}

Friends And Related Function Documentation

operator<<

std::ostream& operator<< ( std::ostream &  os, const Field &  e  )

friend

Definition at line 303 of file FieldMultiIndices.cpp.

                                                       {
  os << e.getNameRef() << " field_id " << e.getId().to_ulong() << " space "
     << FieldSpaceNames[e.getSpace()] << " approximation base "
     << ApproximationBaseNames[e.getApproxBase()] << " rank "
     << e.getNbOfCoeffs() << " meshset " << e.meshSet;
  return os;
}

Member Data Documentation

bitNumber

unsigned int MoFEM::Field::bitNumber

Field Id is bit set. Each field has only one bit on, bitNumber stores number of set bit

Definition at line 121 of file FieldMultiIndices.hpp.

coordSysPtr

boost::shared_ptr<CoordSys> MoFEM::Field::coordSysPtr

Pointer to field coordinate system data structure.

Definition at line 89 of file FieldMultiIndices.hpp.

dofOrderMap

DofsOrderMap MoFEM::Field::dofOrderMap

mutableprivate

Definition at line 380 of file FieldMultiIndices.hpp.

forderTable

FieldOrderTable MoFEM::Field::forderTable

nb. DOFs table for entities

Definition at line 108 of file FieldMultiIndices.hpp.

meshSet

EntityHandle MoFEM::Field::meshSet

keeps entities for this meshset

Definition at line 87 of file FieldMultiIndices.hpp.

moab

moab::Interface& MoFEM::Field::moab

Definition at line 85 of file FieldMultiIndices.hpp.

sequenceDofContainer

SequenceDofContainer MoFEM::Field::sequenceDofContainer

mutableprivate

Definition at line 379 of file FieldMultiIndices.hpp.

tagBaseData

FieldApproximationBase* MoFEM::Field::tagBaseData

tag keeps field base

Definition at line 98 of file FieldMultiIndices.hpp.

tagId

BitFieldId* MoFEM::Field::tagId

tag keeps field id

Definition at line 96 of file FieldMultiIndices.hpp.

tagName

const void* MoFEM::Field::tagName

tag keeps name of the field

Definition at line 103 of file FieldMultiIndices.hpp.

tagNamePrefixData

const void* MoFEM::Field::tagNamePrefixData

tag keeps name prefix of the field

Definition at line 105 of file FieldMultiIndices.hpp.

tagNamePrefixSize

int MoFEM::Field::tagNamePrefixSize

number of bits necessary to keep field name prefix

Definition at line 106 of file FieldMultiIndices.hpp.

tagNameSize

int MoFEM::Field::tagNameSize

number of bits necessary to keep field name

Definition at line 104 of file FieldMultiIndices.hpp.

tagNbCoeffData

FieldCoefficientsNumber* MoFEM::Field::tagNbCoeffData

tag keeps field rank (dimension, f.e. Temperature field has rank 1, displacements field in 3d has rank 3)

Definition at line 102 of file FieldMultiIndices.hpp.

tagSpaceData

FieldSpace* MoFEM::Field::tagSpaceData

tag keeps field space

Definition at line 97 of file FieldMultiIndices.hpp.

th_AppOrder

Tag MoFEM::Field::th_AppOrder

Tag storing approximation order on entity.

Definition at line 94 of file FieldMultiIndices.hpp.

th_FieldData

Tag MoFEM::Field::th_FieldData

Tag storing field values on entity in the field.

Definition at line 93 of file FieldMultiIndices.hpp.

th_FieldDataVerts

Tag MoFEM::Field::th_FieldDataVerts

Tag storing field values on vertices in the field.

Definition at line 92 of file FieldMultiIndices.hpp.

th_FieldDataVertsType

TagType MoFEM::Field::th_FieldDataVertsType

Definition at line 91 of file FieldMultiIndices.hpp.

---

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

• src/multi_indices/FieldMultiIndices.hpp
• src/multi_indices/impl/FieldMultiIndices.cpp