MoFEM::ForcesAndSourcesCore Struct Reference

structure to get information form mofem into DataForcesAndSourcesCore More...

#include <src/finite_elements/ForcesAndSourcesCore.hpp>

Inheritance diagram for MoFEM::ForcesAndSourcesCore:

Collaboration diagram for MoFEM::ForcesAndSourcesCore:

Classes
struct	UserDataOperator
	Data operator to do calculations at integration points.Is inherited and implemented by user to do calculations. It can be used in many different ways but typically is used to integrate matrices (f.e. stiffness matrix) and the right hand vector (f.e. force vector). More...

Public Types
typedef boost::function< int(int order_row, int order_col, int order_data)>	RuleHookFun
Public Types inherited from MoFEM::KspMethod
enum	KSPContext { CTX_SETFUNCTION, CTX_OPERATORS, CTX_KSPNONE }
	pass information about context of KSP/DM for with finite element is computed More...
Public Types inherited from MoFEM::SnesMethod
enum	SNESContext { CTX_SNESSETFUNCTION, CTX_SNESSETJACOBIAN, CTX_SNESNONE }
Public Types inherited from MoFEM::TSMethod
enum	TSContext {   CTX_TSSETRHSFUNCTION, CTX_TSSETRHSJACOBIAN, CTX_TSSETIFUNCTION, CTX_TSSETIJACOBIAN,   CTX_TSTSMONITORSET, CTX_TSNONE }

Public Member Functions
	ForcesAndSourcesCore (Interface &m_field)
virtual	~ForcesAndSourcesCore ()
MoFEMErrorCode	getNumberOfNodes (int &num_nodes) const
int	getMaxDataOrder () const
	Get max order of approximation for data fields. More...
int	getMaxRowOrder () const
	Get max order of approximation for field in rows. More...
int	getMaxColOrder () const
	Get max order of approximation for field in columns. More...
MoFEMErrorCode	getSense (EntityType type, boost::ptr_vector< DataForcesAndSourcesCore::EntData > &data) const
	get sense (orientation) of entity More...
MoFEMErrorCode	getDataOrder (const EntityType type, const FieldSpace space, boost::ptr_vector< DataForcesAndSourcesCore::EntData > &data) const
	get maximal approximation order of approximation on the entity More...
MoFEMErrorCode	getDataOrderSpaceAndBase (const std::string &field_name, const EntityType type, boost::ptr_vector< DataForcesAndSourcesCore::EntData > &data) const
	get maximal approximation order on entity More...
MoFEMErrorCode	getEdgesSense (DataForcesAndSourcesCore &data) const
MoFEMErrorCode	getTrisSense (DataForcesAndSourcesCore &data) const
MoFEMErrorCode	getQuadSense (DataForcesAndSourcesCore &data) const
MoFEMErrorCode	getEdgesDataOrder (DataForcesAndSourcesCore &data, const FieldSpace space) const
MoFEMErrorCode	getTrisDataOrder (DataForcesAndSourcesCore &data, const FieldSpace space) const
MoFEMErrorCode	getQuadDataOrder (DataForcesAndSourcesCore &data, const FieldSpace space) const
MoFEMErrorCode	getTetDataOrder (DataForcesAndSourcesCore &data, const FieldSpace space) const
MoFEMErrorCode	getPrismDataOrder (DataForcesAndSourcesCore &data, const FieldSpace space) const
MoFEMErrorCode	getEdgesDataOrderSpaceAndBase (DataForcesAndSourcesCore &data, const std::string &field_name) const
MoFEMErrorCode	getTrisDataOrderSpaceAndBase (DataForcesAndSourcesCore &data, const std::string &field_name) const
MoFEMErrorCode	getQuadDataOrderSpaceAndBase (DataForcesAndSourcesCore &data, const std::string &field_name) const
MoFEMErrorCode	getTetDataOrderSpaceAndBase (DataForcesAndSourcesCore &data, const std::string &field_name) const
MoFEMErrorCode	getPrismDataOrderSpaceAndBase (DataForcesAndSourcesCore &data, const std::string &field_name) const
MoFEMErrorCode	getNodesIndices (const boost::string_ref field_name, FENumeredDofEntity_multiIndex &dofs, VectorInt &nodes_indices, VectorInt &local_nodes_indices) const
	get node indices More...
MoFEMErrorCode	getTypeIndices (const boost::string_ref field_name, FENumeredDofEntity_multiIndex &dofs, EntityType type, int side_number, VectorInt &indices, VectorInt &local_indices) const
	get indices by type (generic function) More...
MoFEMErrorCode	getTypeIndices (const boost::string_ref field_name, FENumeredDofEntity_multiIndex &dofs, EntityType type, boost::ptr_vector< DataForcesAndSourcesCore::EntData > &data) const
	get indices by type (generic function) More...
MoFEMErrorCode	getRowNodesIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get row node indices from FENumeredDofEntity_multiIndex More...
MoFEMErrorCode	getColNodesIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get col node indices from FENumeredDofEntity_multiIndex More...
MoFEMErrorCode	getEdgesRowIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get Edges row indices from FENumeredDofEntity_multiIndex More...
MoFEMErrorCode	getEdgesColIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get Edges col indices from FENumeredDofEntity_multiIndex More...
MoFEMErrorCode	getTrisRowIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get Tris row indices from FENumeredDofEntity_multiIndex More...
MoFEMErrorCode	getTrisColIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get Tris col indices from FENumeredDofEntity_multiIndex More...
MoFEMErrorCode	getTetsRowIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get Tets row indices from FENumeredDofEntity_multiIndex More...
MoFEMErrorCode	getTetsColIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get Tets col indices from FENumeredDofEntity_multiIndex More...
MoFEMErrorCode	getQuadRowIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get Quad row indices from FENumeredDofEntity_multiIndex More...
MoFEMErrorCode	getQuadColIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get Quad col indices from FENumeredDofEntity_multiIndex More...
MoFEMErrorCode	getPrismRowIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get Prism row indices from FENumeredDofEntity_multiIndex More...
MoFEMErrorCode	getPrismColIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get Prism col indices from FENumeredDofEntity_multiIndex More...
MoFEMErrorCode	getNoFieldIndices (const std::string &field_name, FENumeredDofEntity_multiIndex &dofs, VectorInt &nodes_indices) const
	get NoField indices More...
MoFEMErrorCode	getNoFieldRowIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get col NoField indices More...
MoFEMErrorCode	getNoFieldColIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get col NoField indices More...
MoFEMErrorCode	getNodesFieldData (const boost::string_ref field_name, FEDofEntity_multiIndex &dofs, VectorDouble &nodes_data, VectorDofs &nodes_dofs, FieldSpace &space, FieldApproximationBase &base) const
	Get field data on nodes. More...
MoFEMErrorCode	getTypeFieldData (const boost::string_ref field_name, FEDofEntity_multiIndex &dofs, EntityType type, int side_number, VectorDouble &ent_field_data, VectorDofs &ent_field_dofs) const
	Get field data on entities. More...
MoFEMErrorCode	getTypeFieldData (const boost::string_ref field_name, FEDofEntity_multiIndex &dofs, EntityType type, boost::ptr_vector< DataForcesAndSourcesCore::EntData > &data) const
MoFEMErrorCode	getNoFieldFieldData (const boost::string_ref field_name, FEDofEntity_multiIndex &dofs, VectorDouble &ent_field_data, VectorDofs &ent_field_dofs) const
MoFEMErrorCode	getNoFieldFieldData (DataForcesAndSourcesCore &data, const boost::string_ref field_name) const
MoFEMErrorCode	getNodesFieldData (DataForcesAndSourcesCore &data, const std::string &field_name) const
	Get data on nodes. More...
MoFEMErrorCode	getEdgesFieldData (DataForcesAndSourcesCore &data, const std::string &field_name) const
MoFEMErrorCode	getTrisFieldData (DataForcesAndSourcesCore &data, const std::string &field_name) const
MoFEMErrorCode	getQuadFieldData (DataForcesAndSourcesCore &data, const std::string &field_name) const
MoFEMErrorCode	getTetsFieldData (DataForcesAndSourcesCore &data, const std::string &field_name) const
MoFEMErrorCode	getPrismFieldData (DataForcesAndSourcesCore &data, const std::string &field_name) const
MoFEMErrorCode	getFaceTriNodes (DataForcesAndSourcesCore &data) const
	Get nodes on triangles. More...
MoFEMErrorCode	getSpacesAndBaseOnEntities (DataForcesAndSourcesCore &data) const
	Get field approximation space and base on entities. More...
MoFEMErrorCode	getProblemNodesIndices (const std::string &field_name, const NumeredDofEntity_multiIndex &dofs, VectorInt &nodes_indices) const
	get indices of nodal indices which are declared for problem but not this particular element More...
MoFEMErrorCode	getProblemTypeIndices (const std::string &field_name, const NumeredDofEntity_multiIndex &dofs, EntityType type, int side_number, VectorInt &indices) const
	get indices by type (generic function) which are declared for problem but not this particular element More...
MoFEMErrorCode	getProblemNodesRowIndices (const std::string &field_name, VectorInt &nodes_indices) const
MoFEMErrorCode	getProblemTypeRowIndices (const std::string &field_name, EntityType type, int side_number, VectorInt &indices) const
MoFEMErrorCode	getProblemNodesColIndices (const std::string &field_name, VectorInt &nodes_indices) const
MoFEMErrorCode	getProblemTypeColIndices (const std::string &field_name, EntityType type, int side_number, VectorInt &indices) const
virtual int	getRule (int order)
	set integration rule for finite element More...
virtual int	getRule (int order_row, int order_col, int order_data)
	another variant of getRule More...
virtual MoFEMErrorCode	setGaussPts (int order)
	It will be removed in the future use other variant. More...
virtual MoFEMErrorCode	setGaussPts (int order_row, int order_col, int order_data)
	set user specific integration rule More...
boost::ptr_vector< UserDataOperator > &	getOpPtrVector ()
	Use to push back operator for row operator. More...
virtual MoFEMErrorCode	preProcess ()
	function is run at the beginning of loop More...
virtual MoFEMErrorCode	operator() ()
	function is run for every finite element More...
virtual MoFEMErrorCode	postProcess ()
	function is run at the end of loop More...
Public Member Functions inherited from MoFEM::FEMethod
MoFEMErrorCode	query_interface (const MOFEMuuid &uuid, UnknownInterface **iface) const
	FEMethod ()
template<class MULTIINDEX >
MULTIINDEX::iterator	get_begin (const MULTIINDEX &index, const std::string &field_name, const EntityType type, const int side_number) const
template<class MULTIINDEX >
MULTIINDEX::iterator	get_end (const MULTIINDEX &index, const std::string &field_name, const EntityType type, const int side_number) const
template<class MULTIINDEX >
MULTIINDEX::iterator	get_begin (const MULTIINDEX &index, const std::string &field_name, const EntityType type) const
template<class MULTIINDEX >
MULTIINDEX::iterator	get_end (const MULTIINDEX &index, const std::string &field_name, const EntityType type) const
template<class MULTIINDEX >
MULTIINDEX::iterator	get_begin (const MULTIINDEX &index, const std::string &field_name) const
template<class MULTIINDEX >
MULTIINDEX::iterator	get_end (const MULTIINDEX &index, const std::string &field_name) const
template<class MULTIINDEX >
MULTIINDEX::iterator	get_begin (const MULTIINDEX &index, const EntityHandle ent) const
template<class MULTIINDEX >
MULTIINDEX::iterator	get_end (const MULTIINDEX &index, const EntityHandle ent) const
template<class MULTIINDEX >
MULTIINDEX::iterator	get_begin (const MULTIINDEX &index, const std::string &field_name, const EntityHandle ent) const
template<class MULTIINDEX >
MULTIINDEX::iterator	get_end (const MULTIINDEX &index, const std::string &field_name, const EntityHandle ent) const
Public Member Functions inherited from MoFEM::BasicMethod
	BasicMethod ()
virtual	~BasicMethod ()
int	getNinTheLoop () const
	get number of evaluated element in the loop More...
int	getLoopSize () const
	get loop size More...
MoFEMErrorCode	copyBasicMethod (const BasicMethod &basic)
Public Member Functions inherited from MoFEM::KspMethod
	KspMethod ()
virtual	~KspMethod ()
MoFEMErrorCode	setKspCtx (const KSPContext &ctx)
	set operator type More...
MoFEMErrorCode	setKsp (KSP ksp)
	set solver More...
MoFEMErrorCode	copyKsp (const KspMethod &ksp)
	copy data form another method More...
Public Member Functions inherited from MoFEM::UnknownInterface
template<class IFACE >
MoFEMErrorCode	registerInterface (const MOFEMuuid &uuid, bool error_if_registration_failed=true)
	Register interface. More...
template<class IFACE , bool VERIFY = false>
MoFEMErrorCode	getInterface (const MOFEMuuid &uuid, IFACE *&iface) const
	Get interface by uuid and return reference to pointer of interface. More...
template<class IFACE >
MoFEMErrorCode	getInterface (IFACE *&iface) const
	Get interface refernce to pointer of interface. More...
template<class IFACE >
MoFEMErrorCode	getInterface (IFACE **const iface) const
	Get interface pointer to pointer of interface. More...
template<class IFACE , typename boost::enable_if< boost::is_pointer< IFACE >, int >::type = 0>
IFACE	getInterface () const
	Get interface pointer to pointer of interface. More...
template<class IFACE , typename boost::enable_if< boost::is_reference< IFACE >, int >::type = 0>
IFACE	getInterface () const
	Get reference to interface. More...
template<class IFACE >
IFACE *	getInterface () const
	Function returning pointer to interface. More...
virtual	~UnknownInterface ()
virtual MoFEMErrorCode	getLibVersion (Version &version) const
	Get library version. More...
virtual const MoFEMErrorCode	getFileVersion (moab::Interface &moab, Version &version) const
	Get database major version. More...
virtual MoFEMErrorCode	getInterfaceVersion (Version &version) const
	Get database major version. More...
template<>
MoFEMErrorCode	getInterface (const MOFEMuuid &uuid, UnknownInterface *&iface) const
Public Member Functions inherited from MoFEM::SnesMethod
	SnesMethod ()
virtual	~SnesMethod ()
MoFEMErrorCode	setSnesCtx (const SNESContext &ctx)
	Set SNES context. More...
MoFEMErrorCode	setSnes (SNES snes)
	Set SNES instance. More...
MoFEMErrorCode	copySnes (const SnesMethod &snes)
	Copy snes data. More...
Public Member Functions inherited from MoFEM::TSMethod
	TSMethod ()
virtual	~TSMethod ()
MoFEMErrorCode	setTsCtx (const TSContext &ctx)
	Set Ts context. More...
MoFEMErrorCode	copyTs (const TSMethod &ts)
	Copy TS solver data. More...
MoFEMErrorCode	setTs (TS _ts)
	Set TS solver. More...

Public Attributes
Interface &	mField
RuleHookFun	getRuleHook
	Hook to get rule. More...
boost::ptr_vector< UserDataOperator >	opPtrVector
Public Attributes inherited from MoFEM::FEMethod
std::string	feName
boost::shared_ptr< const NumeredEntFiniteElement >	numeredEntFiniteElementPtr
boost::shared_ptr< const FENumeredDofEntity_multiIndex >	rowPtr
boost::shared_ptr< const FENumeredDofEntity_multiIndex >	colPtr
boost::shared_ptr< const FEDofEntity_multiIndex >	dataPtr
Public Attributes inherited from MoFEM::BasicMethod
int	nInTheLoop
int	loopSize
int	rAnk
int	sIze
const RefEntity_multiIndex *	refinedEntitiesPtr
const RefElement_multiIndex *	refinedFiniteElementsPtr
const Problem *	problemPtr
const Field_multiIndex *	fieldsPtr
const FieldEntity_multiIndex *	entitiesPtr
const DofEntity_multiIndex *	dofsPtr
const FiniteElement_multiIndex *	finiteElementsPtr
const EntFiniteElement_multiIndex *	finiteElementsEntitiesPtr
const FieldEntityEntFiniteElementAdjacencyMap_multiIndex *	adjacenciesPtr
boost::function< MoFEMErrorCode()>	preProcessHook
boost::function< MoFEMErrorCode()>	postProcessHook
boost::function< MoFEMErrorCode()>	operatorHook
Public Attributes inherited from MoFEM::KspMethod
KSPContext	ksp_ctx
	Context. More...
KSP	ksp
	KSP solver. More...
Vec	ksp_f
	the right hand side vector More...
Mat	ksp_A
	matrix More...
Mat	ksp_B
	preconditioner matrix More...
Public Attributes inherited from MoFEM::SnesMethod
SNESContext	snes_ctx
SNES	snes
Vec	snes_x
Vec	snes_f
Mat	snes_A
Mat	snes_B
Public Attributes inherited from MoFEM::TSMethod
TSContext	ts_ctx
TS	ts
Vec	ts_u
Vec	ts_u_t
Vec	ts_F
Mat	ts_A
Mat	ts_B
PetscInt	ts_step
PetscReal	ts_a
PetscReal	ts_t

Additional Inherited Members
Protected Member Functions inherited from MoFEM::UnknownInterface
boost::typeindex::type_index	getClassIdx (const MOFEMuuid &uid) const
	Get type name for interface Id. More...
MOFEMuuid	getUId (const boost::typeindex::type_index &class_idx) const
	Get interface Id for class name. More...

Detailed Description

structure to get information form mofem into DataForcesAndSourcesCore

Examples:

hello_world.cpp.

Definition at line 36 of file ForcesAndSourcesCore.hpp.

Member Typedef Documentation

RuleHookFun

typedef boost::function<int (int order_row,int order_col,int order_data)> MoFEM::ForcesAndSourcesCore::RuleHookFun

Definition at line 311 of file ForcesAndSourcesCore.hpp.

Constructor & Destructor Documentation

ForcesAndSourcesCore()

MoFEM::ForcesAndSourcesCore::ForcesAndSourcesCore

(

Interface &

m_field

)

Definition at line 39 of file ForcesAndSourcesCore.hpp.

                                          :
  mField(m_field),
  getRuleHook(0) {
  }

~ForcesAndSourcesCore()

virtual MoFEM::ForcesAndSourcesCore::~ForcesAndSourcesCore ( )

virtual

Definition at line 43 of file ForcesAndSourcesCore.hpp.

                                  {

  }

Member Function Documentation

getColNodesIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getColNodesIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get col node indices from FENumeredDofEntity_multiIndex

Definition at line 505 of file ForcesAndSourcesCore.cpp.

                                                                            {

  MoFEMFunctionBeginHot;
  ierr = getNodesIndices(field_name,
                         const_cast<FENumeredDofEntity_multiIndex &>(
                             numeredEntFiniteElementPtr->getColsDofs()),
                         data.dataOnEntities[MBVERTEX][0].getIndices(),
                         data.dataOnEntities[MBVERTEX][0].getLocalIndices());
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getDataOrder()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getDataOrder	(	const EntityType	type,
		const FieldSpace	space,
		boost::ptr_vector< DataForcesAndSourcesCore::EntData > &	data
	)		const

get maximal approximation order of approximation on the entity

Definition at line 194 of file ForcesAndSourcesCore.cpp.

                                                                  {
  MoFEMFunctionBeginHot;
  try {
    SideNumber_multiIndex &side_table = const_cast<SideNumber_multiIndex &>(
        numeredEntFiniteElementPtr->getSideNumberTable());
    if (data.size() < side_table.get<2>().count(type)) {
      // prims has 9 edges, some of edges for "flat" prism are not active
      SETERRQ2(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
               "data inconsistency %d < %d", data.size(),
               side_table.get<2>().count(type));
    }
    for (unsigned int side = 0; side < data.size(); side++) {
      data[side].getDataOrder() = 0;
    }
    FEDofEntity_multiIndex::index<Composite_EntType_and_Space_mi_tag>::type
        &data_dofs = const_cast<FEDofEntity_multiIndex::index<
            Composite_EntType_and_Space_mi_tag>::type &>(
            numeredEntFiniteElementPtr->getDataDofs()
                .get<Composite_EntType_and_Space_mi_tag>());
    FEDofEntity_multiIndex::index<
        Composite_EntType_and_Space_mi_tag>::type::iterator dit,
        hi_dit;
    dit = data_dofs.lower_bound(boost::make_tuple(type, space));
    hi_dit = data_dofs.upper_bound(boost::make_tuple(type, space));
    for (; dit != hi_dit; dit++) {
      ApproximationOrder ent_order = (*dit)->getMaxOrder();
      int side_number = (*dit)->sideNumberPtr->side_number;
      if (side_number < 0) {
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "data inconsistency");
      }
      data[side_number].getDataOrder() =
          data[side_number].getDataOrder() > ent_order
              ? data[side_number].getDataOrder()
              : ent_order;
      if ((*dit)->sideNumberPtr->brother_side_number != -1) {
        if (data.size() <
            (unsigned int)(*dit)->sideNumberPtr->brother_side_number) {
          SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                  "data inconsistency");
        }
        data[(*dit)->sideNumberPtr->brother_side_number].getDataOrder() =
            data[side_number].getDataOrder();
      }
    }
  } catch (std::exception &ex) {
    std::ostringstream ss;
    ss << "thorw in method: " << ex.what() << " at line " << __LINE__
       << " in file " << __FILE__;
    SETERRQ(PETSC_COMM_SELF, MOFEM_STD_EXCEPTION_THROW, ss.str().c_str());
  }
  MoFEMFunctionReturnHot(0);
}

getDataOrderSpaceAndBase()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getDataOrderSpaceAndBase	(	const std::string &	field_name,
		const EntityType	type,
		boost::ptr_vector< DataForcesAndSourcesCore::EntData > &	data
	)		const

get maximal approximation order on entity

Definition at line 300 of file ForcesAndSourcesCore.cpp.

                                                                  {
  MoFEMFunctionBeginHot;

  SideNumber_multiIndex &side_table = const_cast<SideNumber_multiIndex &>(
      numeredEntFiniteElementPtr->getSideNumberTable());
  if (data.size() < side_table.get<2>().count(type)) {
    // prims has 9 edges, some of edges for "flat" prism are not active
    SETERRQ2(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "data inconsistency %d < %d", data.size(),
             side_table.get<2>().count(type));
  }

  for (unsigned int side = 0; side != data.size(); ++side) {
    data[side].getDataOrder() = 0;
    data[side].getBase() = NOBASE;
    data[side].getSpace() = NOSPACE;
  }

  // get dofs by name, type and side
  const FEDofEntity_multiIndex::index<
      Composite_Name_Type_And_Side_Number_mi_tag>::type &data_dofs =
      numeredEntFiniteElementPtr->getDataDofs()
          .get<Composite_Name_Type_And_Side_Number_mi_tag>();
  FEDofEntity_multiIndex::index<
      Composite_Name_Type_And_Side_Number_mi_tag>::type::iterator dit,
      hi_dit;
  dit = data_dofs.lower_bound(boost::make_tuple(field_name, type, 0));
  if (dit == data_dofs.end()) {
    MoFEMFunctionReturnHot(0);
  }
  hi_dit =
      data_dofs.upper_bound(boost::make_tuple(field_name, type, data.size()));

  for (; dit != hi_dit;) {

    // std::cerr << ApproximationBaseNames[dit->getApproxBase()] << std::endl;

    int side_number = (*dit)->sideNumberPtr->side_number;
    if (data[side_number].getDataOrder()) {
      const int nb_dofs_on_ent = (*dit)->getNbDofsOnEnt();
      for (int i = 0; i != nb_dofs_on_ent; ++i, ++dit) {
      }
      continue;
    }

    ApproximationOrder ent_order = (*dit)->getMaxOrder();
    if (PetscUnlikely(side_number < 0)) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Side number is not set to dof entity");
    }
    data[side_number].getBase() = (*dit)->getApproxBase();
    data[side_number].getSpace() = (*dit)->getSpace();
    data[side_number].getDataOrder() =
        data[side_number].getDataOrder() > ent_order
            ? data[side_number].getDataOrder()
            : ent_order;
    if ((*dit)->sideNumberPtr->brother_side_number != -1) {
      if (data.size() <
          (unsigned int)(*dit)->sideNumberPtr->brother_side_number) {
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "data inconsistency");
      }
      data[(*dit)->sideNumberPtr->brother_side_number].getBase() =
          data[side_number].getBase();
      data[(*dit)->sideNumberPtr->brother_side_number].getSpace() =
          data[side_number].getSpace();
      data[(*dit)->sideNumberPtr->brother_side_number].getDataOrder() =
          data[side_number].getDataOrder();
    }

    const int nb_dofs_on_ent = (*dit)->getNbDofsOnEnt();
    for (int i = 0; i != nb_dofs_on_ent; ++i, ++dit) {
    }
  }

  MoFEMFunctionReturnHot(0);
}

getEdgesColIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getEdgesColIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get Edges col indices from FENumeredDofEntity_multiIndex

Definition at line 631 of file ForcesAndSourcesCore.cpp.

                                                                            {

  MoFEMFunctionBeginHot;
  ierr = getTypeIndices(field_name,
                        const_cast<FENumeredDofEntity_multiIndex &>(
                            numeredEntFiniteElementPtr->getColsDofs()),
                        MBEDGE, data.dataOnEntities[MBEDGE]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getEdgesDataOrder()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getEdgesDataOrder	(	DataForcesAndSourcesCore &	data,
		const FieldSpace	space
	)		const

Definition at line 251 of file ForcesAndSourcesCore.cpp.

                                                                      {

  MoFEMFunctionBeginHot;
  ierr = getDataOrder(MBEDGE, space, data.dataOnEntities[MBEDGE]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getEdgesDataOrderSpaceAndBase()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getEdgesDataOrderSpaceAndBase	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

Definition at line 380 of file ForcesAndSourcesCore.cpp.

                                                                       {

  MoFEMFunctionBeginHot;
  ierr =
      getDataOrderSpaceAndBase(field_name, MBEDGE, data.dataOnEntities[MBEDGE]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getEdgesFieldData()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getEdgesFieldData	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

Definition at line 1167 of file ForcesAndSourcesCore.cpp.

                                                                           {

  MoFEMFunctionBeginHot;
  ierr = getTypeFieldData(field_name,
                          const_cast<FEDofEntity_multiIndex &>(
                              numeredEntFiniteElementPtr->getDataDofs()),
                          MBEDGE, data.dataOnEntities[MBEDGE]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getEdgesRowIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getEdgesRowIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get Edges row indices from FENumeredDofEntity_multiIndex

Definition at line 618 of file ForcesAndSourcesCore.cpp.

                                                                            {

  MoFEMFunctionBeginHot;
  ierr = getTypeIndices(field_name,
                        const_cast<FENumeredDofEntity_multiIndex &>(
                            numeredEntFiniteElementPtr->getRowsDofs()),
                        MBEDGE, data.dataOnEntities[MBEDGE]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getEdgesSense()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getEdgesSense

(

DataForcesAndSourcesCore &

data

)

const

get edge sense )orientation) in respect finite element entity

Parameters

data	structure where results are stored

Returns

error code

Definition at line 139 of file ForcesAndSourcesCore.cpp.

                                                                        {

  MoFEMFunctionBeginHot;
  ierr = getSense(MBEDGE, data.dataOnEntities[MBEDGE]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getFaceTriNodes()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getFaceTriNodes

(

DataForcesAndSourcesCore &

data

)

const

Get nodes on triangles.

Definition at line 1244 of file ForcesAndSourcesCore.cpp.

                                                                          {
  //
  MoFEMFunctionBeginHot;
  // PetscAttachDebugger();
  data.facesNodes.resize(4, 3, false);
  SideNumber_multiIndex &side_table = const_cast<SideNumber_multiIndex &>(
      numeredEntFiniteElementPtr->getSideNumberTable());
  SideNumber_multiIndex::nth_index<1>::type::iterator siit =
      side_table.get<1>().lower_bound(boost::make_tuple(MBTRI, 0));
  SideNumber_multiIndex::nth_index<1>::type::iterator hi_siit =
      side_table.get<1>().upper_bound(boost::make_tuple(MBTRI, 4));
  if (distance(siit, hi_siit) != 4) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "Should be 4 triangles on tet, side_table not initialized");
  }
  const int cannonical_face_sense_p1[4][3] = {
      {0, 1, 3},
      {1, 2, 3},
      {0, 3, 2} /**/,
      {0, 2, 1} /**/}; //secon index is offset (positive sense)
  const int cannonical_face_sense_m1[4][3] = {
      {0, 3, 1},
      {1, 3, 2},
      {0, 2, 3},
      {0, 1, 2}}; //second index is offset (negative sense
  for (; siit != hi_siit; siit++) {
    const boost::shared_ptr<SideNumber> side = *siit;
    int face_conn[3] = {-1, -1, -1};
    if (side->offset == 0) {
      face_conn[0] = side->sense == 1
                         ? cannonical_face_sense_p1[(int)side->side_number][0]
                         : cannonical_face_sense_m1[(int)side->side_number][0];
      face_conn[1] = side->sense == 1
                         ? cannonical_face_sense_p1[(int)side->side_number][1]
                         : cannonical_face_sense_m1[(int)side->side_number][1];
      face_conn[2] = side->sense == 1
                         ? cannonical_face_sense_p1[(int)side->side_number][2]
                         : cannonical_face_sense_m1[(int)side->side_number][2];
    }
    if (side->offset == 1) {
      face_conn[0] =
          side->sense == 1
              ? cannonical_face_sense_p1[(int)side->side_number][1]
              : cannonical_face_sense_m1[(int)side->side_number][2] /**/;
      face_conn[1] = side->sense == 1
                         ? cannonical_face_sense_p1[(int)side->side_number][2]
                         : cannonical_face_sense_m1[(int)side->side_number][0];
      face_conn[2] = side->sense == 1
                         ? cannonical_face_sense_p1[(int)side->side_number][0]
                         : cannonical_face_sense_m1[(int)side->side_number][1];
    }
    if (side->offset == 2) {
      face_conn[0] =
          side->sense == 1
              ? cannonical_face_sense_p1[(int)side->side_number][2]
              : cannonical_face_sense_m1[(int)side->side_number][1] /**/;
      face_conn[1] = side->sense == 1
                         ? cannonical_face_sense_p1[(int)side->side_number][0]
                         : cannonical_face_sense_m1[(int)side->side_number][2];
      face_conn[2] = side->sense == 1
                         ? cannonical_face_sense_p1[(int)side->side_number][1]
                         : cannonical_face_sense_m1[(int)side->side_number][0];
    }
    for (int nn = 0; nn < 3; nn++)
      data.facesNodes(side->side_number, nn) = face_conn[nn];
    {
      const EntityHandle *conn_tet;
      int num_nodes_tet;
      EntityHandle ent = numeredEntFiniteElementPtr->getEnt();
      rval = mField.get_moab().get_connectivity(ent, conn_tet, num_nodes_tet,
                                                true);
      CHKERRQ_MOAB(rval);
      if (num_nodes_tet != 4)
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "data inconsistency");
      int num_nodes_face;
      const EntityHandle *conn_face;
      rval = mField.get_moab().get_connectivity(side->ent, conn_face,
                                                num_nodes_face, true);
      CHKERRQ_MOAB(rval);
      if (num_nodes_face != 3)
        SETERRQ(PETSC_COMM_SELF, 1, "data inconsistency");
      if (conn_face[0] != conn_tet[data.facesNodes(side->side_number, 0)])
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "data inconsistency");
      if (conn_face[1] != conn_tet[data.facesNodes(side->side_number, 1)])
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "data inconsistency");
      if (conn_face[2] != conn_tet[data.facesNodes(side->side_number, 2)])
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "data inconsistency");
    }
  }
  MoFEMFunctionReturnHot(0);
}

getMaxColOrder()

int MoFEM::ForcesAndSourcesCore::getMaxColOrder

(

)

const

Get max order of approximation for field in columns.

Definition at line 190 of file ForcesAndSourcesCore.cpp.

                                               {
  return getMaxOrder(numeredEntFiniteElementPtr->getColsDofs());
}

getMaxDataOrder()

int MoFEM::ForcesAndSourcesCore::getMaxDataOrder

(

)

const

Get max order of approximation for data fields.

Method getMaxDataOrder () return maximal order on entities, for all data on the element. So for example if finite element is triangle, and triangle base function have order 4 and on edges base function have order 2, this function return 4.

If finite element has for example 2 or more approximated fields, for example Pressure (order 3) and displacement field (order 5), this function returns 5.

Definition at line 182 of file ForcesAndSourcesCore.cpp.

                                                {
  return getMaxOrder(numeredEntFiniteElementPtr->getDataDofs());
}

getMaxRowOrder()

int MoFEM::ForcesAndSourcesCore::getMaxRowOrder

(

)

const

Get max order of approximation for field in rows.

Definition at line 186 of file ForcesAndSourcesCore.cpp.

                                               {
  return getMaxOrder(numeredEntFiniteElementPtr->getRowsDofs());
}

getNodesFieldData() [1/2]

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNodesFieldData	(	const boost::string_ref	field_name,
		FEDofEntity_multiIndex &	dofs,
		VectorDouble &	nodes_data,
		VectorDofs &	nodes_dofs,
		FieldSpace &	space,
		FieldApproximationBase &	base
	)		const

Get field data on nodes.

Parameters

field_name	Name of field
dofs	Dofs (element) multi index
nodes_data	Returned DOFs values
nodes_dofs	Vector of pointers to DOFs data structure
space	Get space on nodes (Only H! is valid)
base	Get base on nodes

Returns

Error code

Definition at line 936 of file ForcesAndSourcesCore.cpp.

                                        {

  MoFEMFunctionBeginHot;
  try {
    FEDofEntityByNameAndType::iterator dit, hi_dit, it;
    dit = dofs.get<Composite_Name_And_Type_mi_tag>().lower_bound(
        boost::make_tuple(field_name, MBVERTEX));
    hi_dit = dofs.get<Composite_Name_And_Type_mi_tag>().upper_bound(
        boost::make_tuple(field_name, MBVERTEX));

    int num_nodes;
    ierr = getNumberOfNodes(num_nodes);
    CHKERRG(ierr);
    int max_nb_dofs = 0;
    if (dit != hi_dit) {
      max_nb_dofs = (*dit)->getNbOfCoeffs() * num_nodes;
    }

    if (distance(dit, hi_dit) != max_nb_dofs) {
      nodes_data.resize(max_nb_dofs, false);
      nodes_data.clear();
      nodes_dofs.resize(max_nb_dofs, false);
    } else {
      int size = distance(dit, hi_dit);
      nodes_data.resize(size, false);
      nodes_dofs.resize(size, false);
    }

    if (dit != hi_dit) {
      space = (*dit)->getSpace();
      base = (*dit)->getApproxBase();
    }

    for (; dit != hi_dit; dit++) {
      FieldData val = (*dit)->getFieldData();
      int side_number = (*dit)->sideNumberPtr->side_number;
      if (side_number == -1) {
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "data inconsistency");
      }
      int pos =
          side_number * (*dit)->getNbOfCoeffs() + (*dit)->getDofCoeffIdx();
      nodes_data[pos] = val;
      nodes_dofs[pos] = *dit;
      int brother_side_number = (*dit)->sideNumberPtr->brother_side_number;
      if (brother_side_number != -1) {
        if (nodes_data.size() <
            (unsigned int)(brother_side_number * (*dit)->getNbOfCoeffs() +
                           (*dit)->getNbOfCoeffs())) {
          nodes_data.resize(brother_side_number * (*dit)->getNbOfCoeffs() +
                            (*dit)->getNbOfCoeffs());
          nodes_dofs.resize(brother_side_number * (*dit)->getNbOfCoeffs() +
                            (*dit)->getNbOfCoeffs());
        }
        int brother_pos = brother_side_number * (*dit)->getNbOfCoeffs() +
                          (*dit)->getDofCoeffIdx();
        nodes_data[brother_pos] = val;
        nodes_dofs[brother_pos] = *dit;
      }
    }

  } catch (std::exception &ex) {
    std::ostringstream ss;
    ss << "thorw in method: " << ex.what() << " at line " << __LINE__
       << " in file " << __FILE__;
    SETERRQ(PETSC_COMM_SELF, MOFEM_STD_EXCEPTION_THROW, ss.str().c_str());
  }
  MoFEMFunctionReturnHot(0);
}

getNodesFieldData() [2/2]

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNodesFieldData	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

Get data on nodes.

Parameters

data	Data structure
field_name	Field name

Returns

Error code

Definition at line 1010 of file ForcesAndSourcesCore.cpp.

                                                                           {
  MoFEMFunctionBeginHot;

  ierr = getNodesFieldData(field_name,
                           const_cast<FEDofEntity_multiIndex &>(
                               numeredEntFiniteElementPtr->getDataDofs()),
                           data.dataOnEntities[MBVERTEX][0].getFieldData(),
                           data.dataOnEntities[MBVERTEX][0].getFieldDofs(),
                           data.dataOnEntities[MBVERTEX][0].getSpace(),
                           data.dataOnEntities[MBVERTEX][0].getBase());
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getNodesIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNodesIndices	(	const boost::string_ref	field_name,
		FENumeredDofEntity_multiIndex &	dofs,
		VectorInt &	nodes_indices,
		VectorInt &	local_nodes_indices
	)		const

get node indices

Definition at line 432 of file ForcesAndSourcesCore.cpp.

                                                                    {

  MoFEMFunctionBeginHot;

  FENumeredDofEntityByNameAndType::iterator dit, hi_dit, it;
  dit = dofs.get<Composite_Name_And_Type_mi_tag>().lower_bound(
      boost::make_tuple(field_name, MBVERTEX));
  hi_dit = dofs.get<Composite_Name_And_Type_mi_tag>().upper_bound(
      boost::make_tuple(field_name, MBVERTEX));

  int num_nodes;
  ierr = getNumberOfNodes(num_nodes);
  CHKERRG(ierr);
  int max_nb_dofs = 0;
  if (dit != hi_dit) {
    max_nb_dofs = (*dit)->getNbOfCoeffs() * num_nodes;
  }

  if (distance(dit, hi_dit) != max_nb_dofs) {
    nodes_indices.resize(max_nb_dofs, false);
    local_nodes_indices.resize(max_nb_dofs, false);
    for (int dd = 0; dd < max_nb_dofs; dd++) {
      nodes_indices[dd] = -1;
      local_nodes_indices[dd] = -1;
    }
  } else {
    nodes_indices.resize(distance(dit, hi_dit), false);
    local_nodes_indices.resize(distance(dit, hi_dit), false);
  }

  for (; dit != hi_dit; dit++) {
    const int idx = (*dit)->getPetscGlobalDofIdx();
    const int local_idx = (*dit)->getPetscLocalDofIdx();
    const int side_number = (*dit)->sideNumberPtr->side_number;
    const int pos =
        side_number * (*dit)->getNbOfCoeffs() + (*dit)->getDofCoeffIdx();
    nodes_indices[pos] = idx;
    local_nodes_indices[pos] = local_idx;
    const int brother_side_number = (*dit)->sideNumberPtr->brother_side_number;
    if (brother_side_number != -1) {
      if (nodes_indices.size() <
          (unsigned int)(brother_side_number * (*dit)->getNbOfCoeffs() +
                         (*dit)->getNbOfCoeffs())) {
        nodes_indices.resize(brother_side_number * (*dit)->getNbOfCoeffs() +
                             (*dit)->getNbOfCoeffs());
      }
      const int elem_idx = brother_side_number * (*dit)->getNbOfCoeffs() +
                           (*dit)->getDofCoeffIdx();
      nodes_indices[elem_idx] = idx;
      local_nodes_indices[elem_idx] = local_idx;
    }
  }

  MoFEMFunctionReturnHot(0);
}

getNoFieldColIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNoFieldColIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get col NoField indices

Definition at line 788 of file ForcesAndSourcesCore.cpp.

                                                                       {

  MoFEMFunctionBeginHot;
  if (data.dataOnEntities[MBENTITYSET].size() == 0) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
  }
  ierr = getNoFieldIndices(field_name,
                           const_cast<FENumeredDofEntity_multiIndex &>(
                               numeredEntFiniteElementPtr->getColsDofs()),
                           data.dataOnEntities[MBENTITYSET][0].getIndices());
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getNoFieldFieldData() [1/2]

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNoFieldFieldData	(	const boost::string_ref	field_name,
		FEDofEntity_multiIndex &	dofs,
		VectorDouble &	ent_field_data,
		VectorDofs &	ent_field_dofs
	)		const

Definition at line 1130 of file ForcesAndSourcesCore.cpp.

                                                                    {
  //
  MoFEMFunctionBeginHot;
  FEDofEntityByFieldName::iterator dit, hi_dit;
  dit = dofs.get<FieldName_mi_tag>().lower_bound(field_name);
  hi_dit = dofs.get<FieldName_mi_tag>().upper_bound(field_name);
  int size = distance(dit, hi_dit);
  ent_field_data.resize(size, false);
  ent_field_dofs.resize(size, false);
  for (; dit != hi_dit; dit++) {
    int idx = (*dit)->getDofCoeffIdx();
    ent_field_data[idx] = (*dit)->getFieldData();
    ent_field_dofs[idx] = *dit;
  }
  MoFEMFunctionReturnHot(0);
}

getNoFieldFieldData() [2/2]

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNoFieldFieldData	(	DataForcesAndSourcesCore &	data,
		const boost::string_ref	field_name
	)		const

Definition at line 1149 of file ForcesAndSourcesCore.cpp.

                                                                            {

  MoFEMFunctionBeginHot;
  if (data.dataOnEntities[MBENTITYSET].size() == 0) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
  }
  ierr =
      getNoFieldFieldData(field_name,
                          const_cast<FEDofEntity_multiIndex &>(
                              numeredEntFiniteElementPtr->getDataDofs()),
                          data.dataOnEntities[MBENTITYSET][0].getFieldData(),
                          data.dataOnEntities[MBENTITYSET][0].getFieldDofs());
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getNoFieldIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNoFieldIndices	(	const std::string &	field_name,
		FENumeredDofEntity_multiIndex &	dofs,
		VectorInt &	nodes_indices
	)		const

get NoField indices

Definition at line 756 of file ForcesAndSourcesCore.cpp.

                                                                  {
  //
  MoFEMFunctionBeginHot;
  FENumeredDofEntityByFieldName::iterator dit, hi_dit;
  dit = dofs.get<FieldName_mi_tag>().lower_bound(field_name);
  hi_dit = dofs.get<FieldName_mi_tag>().upper_bound(field_name);
  indices.resize(distance(dit, hi_dit));
  for (; dit != hi_dit; dit++) {
    int idx = (*dit)->getPetscGlobalDofIdx();
    indices[(*dit)->getDofCoeffIdx()] = idx;
  }
  MoFEMFunctionReturnHot(0);
}

getNoFieldRowIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNoFieldRowIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get col NoField indices

Definition at line 772 of file ForcesAndSourcesCore.cpp.

                                                                       {

  MoFEMFunctionBeginHot;
  // EntityType fe_type = numeredEntFiniteElementPtr->getEntType();
  if (data.dataOnEntities[MBENTITYSET].size() == 0) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
  }
  ierr = getNoFieldIndices(field_name,
                           const_cast<FENumeredDofEntity_multiIndex &>(
                               numeredEntFiniteElementPtr->getRowsDofs()),
                           data.dataOnEntities[MBENTITYSET][0].getIndices());
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getNumberOfNodes()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNumberOfNodes

(

int &

num_nodes

)

const

Definition at line 69 of file ForcesAndSourcesCore.cpp.

                                                                          {
  MoFEMFunctionBeginHot;

  EntityHandle ent = numeredEntFiniteElementPtr->getEnt();
  switch (mField.get_moab().type_from_handle(ent)) {
  case MBVERTEX:
    num_nodes = 1;
    break;
  case MBEDGE:
    num_nodes = 2;
    break;
  case MBTRI:
    num_nodes = 3;
    break;
  case MBQUAD:
    num_nodes = 4;
    break;
  case MBTET:
    num_nodes = 4;
    break;
  case MBPRISM:
    num_nodes = 6;
    break;
  default:
    SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED, "not implemented");
  }

  MoFEMFunctionReturnHot(0);
}

getOpPtrVector()

boost::ptr_vector<UserDataOperator>& MoFEM::ForcesAndSourcesCore::getOpPtrVector

(

)

Use to push back operator for row operator.

It can be used to calculate nodal forces or other quantities on the mesh.

Examples:

forces_and_sources_testing_edge_element.cpp, and MagneticElement.hpp.

Definition at line 620 of file ForcesAndSourcesCore.hpp.

{ return opPtrVector; }

getPrismColIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getPrismColIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get Prism col indices from FENumeredDofEntity_multiIndex

Definition at line 743 of file ForcesAndSourcesCore.cpp.

                                                                            {

  MoFEMFunctionBeginHot;
  ierr = getTypeIndices(field_name,
                        const_cast<FENumeredDofEntity_multiIndex &>(
                            numeredEntFiniteElementPtr->getColsDofs()),
                        MBPRISM, data.dataOnEntities[MBPRISM]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getPrismDataOrder()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getPrismDataOrder	(	DataForcesAndSourcesCore &	data,
		const FieldSpace	space
	)		const

Definition at line 291 of file ForcesAndSourcesCore.cpp.

                                                                      {

  MoFEMFunctionBeginHot;
  ierr = getDataOrder(MBPRISM, space, data.dataOnEntities[MBPRISM]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getPrismDataOrderSpaceAndBase()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getPrismDataOrderSpaceAndBase	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

Definition at line 420 of file ForcesAndSourcesCore.cpp.

                                                                       {

  MoFEMFunctionBeginHot;
  ierr = getDataOrderSpaceAndBase(field_name, MBPRISM,
                                  data.dataOnEntities[MBPRISM]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getPrismFieldData()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getPrismFieldData	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

Definition at line 1224 of file ForcesAndSourcesCore.cpp.

                                                                           {

  MoFEMFunctionBeginHot;
  if (data.dataOnEntities[MBPRISM].size() == 0) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
  }
  ierr = getTypeFieldData(field_name,
                          const_cast<FEDofEntity_multiIndex &>(
                              numeredEntFiniteElementPtr->getDataDofs()),
                          MBPRISM, 0,
                          data.dataOnEntities[MBPRISM][0].getFieldData(),
                          data.dataOnEntities[MBPRISM][0].getFieldDofs());
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getPrismRowIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getPrismRowIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get Prism row indices from FENumeredDofEntity_multiIndex

Definition at line 730 of file ForcesAndSourcesCore.cpp.

                                                                            {

  MoFEMFunctionBeginHot;
  ierr = getTypeIndices(field_name,
                        const_cast<FENumeredDofEntity_multiIndex &>(
                            numeredEntFiniteElementPtr->getRowsDofs()),
                        MBPRISM, data.dataOnEntities[MBPRISM]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getProblemNodesColIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getProblemNodesColIndices	(	const std::string &	field_name,
		VectorInt &	nodes_indices
	)		const

Definition at line 912 of file ForcesAndSourcesCore.cpp.

                                                                 {
  MoFEMFunctionBeginHot;

  ierr = getProblemNodesIndices(field_name, *(problemPtr->numeredDofsCols),
                                nodes_indices);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getProblemNodesIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getProblemNodesIndices	(	const std::string &	field_name,
		const NumeredDofEntity_multiIndex &	dofs,
		VectorInt &	nodes_indices
	)		const

get indices of nodal indices which are declared for problem but not this particular element

Definition at line 805 of file ForcesAndSourcesCore.cpp.

                                    {
  MoFEMFunctionBeginHot;

  nodes_indices.resize(0);

  SideNumber_multiIndex &side_table = const_cast<SideNumber_multiIndex &>(
      numeredEntFiniteElementPtr->getSideNumberTable());
  SideNumber_multiIndex::nth_index<1>::type::iterator siit =
      side_table.get<1>().lower_bound(boost::make_tuple(MBVERTEX, 0));
  SideNumber_multiIndex::nth_index<1>::type::iterator hi_siit =
      side_table.get<1>().lower_bound(boost::make_tuple(MBVERTEX, 10000));

  int nn = 0;
  for (; siit != hi_siit; siit++, nn++) {

    if (siit->get()->side_number == -1)
      continue;

    const EntityHandle ent = siit->get()->ent;
    NumeredDofEntity_multiIndex::index<
        Composite_Name_And_Ent_And_EntDofIdx_mi_tag>::type::iterator dit,
        hi_dit;
    dit = dofs.get<Composite_Name_And_Ent_And_EntDofIdx_mi_tag>().lower_bound(
        boost::make_tuple(field_name, ent, 0));
    hi_dit =
        dofs.get<Composite_Name_And_Ent_And_EntDofIdx_mi_tag>().upper_bound(
            boost::make_tuple(field_name, ent, 10000)); /// very large number

    if (dit != hi_dit) {

      if (!nn) {
        nodes_indices.resize((*dit)->getNbOfCoeffs() * distance(siit, hi_siit));
      }
      for (; dit != hi_dit; dit++) {
        nodes_indices[siit->get()->side_number * (*dit)->getNbOfCoeffs() +
                      (*dit)->getDofCoeffIdx()] =
            (*dit)->getPetscGlobalDofIdx();
      }
    }
  }

  MoFEMFunctionReturnHot(0);
}

getProblemNodesRowIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getProblemNodesRowIndices	(	const std::string &	field_name,
		VectorInt &	nodes_indices
	)		const

Definition at line 890 of file ForcesAndSourcesCore.cpp.

                                                                 {
  MoFEMFunctionBeginHot;

  ierr = getProblemNodesIndices(field_name, *(problemPtr->numeredDofsRows),
                                nodes_indices);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getProblemTypeColIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getProblemTypeColIndices	(	const std::string &	field_name,
		EntityType	type,
		int	side_number,
		VectorInt &	indices
	)		const

Definition at line 923 of file ForcesAndSourcesCore.cpp.

                                                                         {
  MoFEMFunctionBeginHot;

  ierr = getProblemTypeIndices(field_name, *(problemPtr->numeredDofsCols), type,
                               side_number, indices);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getProblemTypeIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getProblemTypeIndices	(	const std::string &	field_name,
		const NumeredDofEntity_multiIndex &	dofs,
		EntityType	type,
		int	side_number,
		VectorInt &	indices
	)		const

get indices by type (generic function) which are declared for problem but not this particular element

Definition at line 851 of file ForcesAndSourcesCore.cpp.

                                                                {
  MoFEMFunctionBeginHot;

  indices.resize(0);

  SideNumber_multiIndex &side_table = const_cast<SideNumber_multiIndex &>(
      numeredEntFiniteElementPtr->getSideNumberTable());
  SideNumber_multiIndex::nth_index<1>::type::iterator siit =
      side_table.get<1>().lower_bound(boost::make_tuple(type, side_number));
  SideNumber_multiIndex::nth_index<1>::type::iterator hi_siit =
      side_table.get<1>().upper_bound(boost::make_tuple(type, side_number));

  for (; siit != hi_siit; siit++) {

    if (siit->get()->side_number == -1)
      continue;

    const EntityHandle ent = siit->get()->ent;
    NumeredDofEntity_multiIndex::index<
        Composite_Name_And_Ent_And_EntDofIdx_mi_tag>::type::iterator dit,
        hi_dit;
    dit = dofs.get<Composite_Name_And_Ent_And_EntDofIdx_mi_tag>().lower_bound(
        boost::make_tuple(field_name, ent, 0));
    hi_dit =
        dofs.get<Composite_Name_And_Ent_And_EntDofIdx_mi_tag>().upper_bound(
            boost::make_tuple(field_name, ent, 10000)); /// very large number

    indices.resize(distance(dit, hi_dit));
    for (; dit != hi_dit; dit++) {

      indices[(*dit)->getEntDofIdx()] = (*dit)->getPetscGlobalDofIdx();
    }
  }

  MoFEMFunctionReturnHot(0);
}

getProblemTypeRowIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getProblemTypeRowIndices	(	const std::string &	field_name,
		EntityType	type,
		int	side_number,
		VectorInt &	indices
	)		const

Definition at line 901 of file ForcesAndSourcesCore.cpp.

                                                                         {
  MoFEMFunctionBeginHot;

  ierr = getProblemTypeIndices(field_name, *(problemPtr->numeredDofsRows), type,
                               side_number, indices);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getQuadColIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getQuadColIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get Quad col indices from FENumeredDofEntity_multiIndex

Definition at line 717 of file ForcesAndSourcesCore.cpp.

                                                                           {

  MoFEMFunctionBeginHot;
  ierr = getTypeIndices(field_name,
                        const_cast<FENumeredDofEntity_multiIndex &>(
                            numeredEntFiniteElementPtr->getColsDofs()),
                        MBQUAD, data.dataOnEntities[MBQUAD]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getQuadDataOrder()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getQuadDataOrder	(	DataForcesAndSourcesCore &	data,
		const FieldSpace	space
	)		const

Definition at line 271 of file ForcesAndSourcesCore.cpp.

                                                                     {

  MoFEMFunctionBeginHot;
  ierr = getDataOrder(MBQUAD, space, data.dataOnEntities[MBQUAD]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getQuadDataOrderSpaceAndBase()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getQuadDataOrderSpaceAndBase	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

Definition at line 400 of file ForcesAndSourcesCore.cpp.

                                                                       {

  MoFEMFunctionBeginHot;
  ierr =
      getDataOrderSpaceAndBase(field_name, MBQUAD, data.dataOnEntities[MBQUAD]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getQuadFieldData()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getQuadFieldData	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

Definition at line 1193 of file ForcesAndSourcesCore.cpp.

                                                                          {

  MoFEMFunctionBeginHot;
  ierr = getTypeFieldData(field_name,
                          const_cast<FEDofEntity_multiIndex &>(
                              numeredEntFiniteElementPtr->getDataDofs()),
                          MBQUAD, data.dataOnEntities[MBQUAD]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getQuadRowIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getQuadRowIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get Quad row indices from FENumeredDofEntity_multiIndex

Definition at line 704 of file ForcesAndSourcesCore.cpp.

                                                                           {

  MoFEMFunctionBeginHot;
  ierr = getTypeIndices(field_name,
                        const_cast<FENumeredDofEntity_multiIndex &>(
                            numeredEntFiniteElementPtr->getRowsDofs()),
                        MBQUAD, data.dataOnEntities[MBQUAD]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getQuadSense()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getQuadSense

(

DataForcesAndSourcesCore &

data

)

const

Definition at line 157 of file ForcesAndSourcesCore.cpp.

                                                                       {

  MoFEMFunctionBeginHot;
  ierr = getSense(MBQUAD, data.dataOnEntities[MBQUAD]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getRowNodesIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getRowNodesIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get row node indices from FENumeredDofEntity_multiIndex

Definition at line 491 of file ForcesAndSourcesCore.cpp.

                                                                            {

  MoFEMFunctionBeginHot;
  ierr = getNodesIndices(field_name,
                         const_cast<FENumeredDofEntity_multiIndex &>(
                             numeredEntFiniteElementPtr->getRowsDofs()),
                         data.dataOnEntities[MBVERTEX][0].getIndices(),
                         data.dataOnEntities[MBVERTEX][0].getLocalIndices());
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getRule() [1/2]

virtual int MoFEM::ForcesAndSourcesCore::getRule ( int  order )

virtual

set integration rule for finite element

This function is overloaded by the user. The integration rule is set such that specific operator implemented by the user is integrated accurately. For example if user implement bilinear operator

\[ b(u,v) = \int_\mathcal{T} \frac{\partial u_i}{\partial x_j}\frac{\partial v_i}{\partial x_j} \textrm{d}\mathcal{T} \]

then if \(u\) and \(v\) are polynomial of given order, then exact integral would be

int getRule(int order) { return 2*(order-1); };

The integration points and weights are set appropriately for given entity type and integration rule from quad.c

Method ForcesAndSourcesCore::getRule takes at argument takes maximal polynomial order set on the element on all fields defined on the element. If a user likes to have more control, another variant of this function can be called which distinguishing between field orders on rows, columns and data, the i.e. first argument of a bilinear form, the second argument of bilinear form and field coefficents on the element.

Note

If user set rule to -1 or any other negative integer, then method ForcesAndSourcesCore::setGaussPts is called. In that method user can implement own (specific) integration method.

Bug:

this function should be const

Reimplemented in PostProcFaceOnRefinedMesh, BoneRemodeling::DensityMapFe, MoFEM::VolumeElementForcesAndSourcesCoreOnSide, GelModule::Gel::GelFE, SmallStrainPlasticity::MyVolumeFE, PostProcTemplateVolumeOnRefinedMesh< MoFEM::VolumeElementForcesAndSourcesCore >, KelvinVoigtDamper::DamperFE, NitscheMethod::MyVolumeFE, ThermalElement::MyTriFE, MixTransport::MixTransportElement::MyTriFE, MagneticElement::TriFE, NonlinearElasticElement::MyVolumeFE, NeummanForcesSurface::MyTriangleFE, MoFEM::NormElement::MyTriFE, NeummanForcesSurfaceComplexForLazy::MyTriangleSpatialFE, ThermalElement::MyVolumeFE, CellEngineering::FaceElement, MixTransport::MixTransportElement::MyVolumeFE, MoFEM::NormElement::MyVolumeFE, OptimalMassTransport::TriFE, MagneticElement::VolumeFE, BoneRemodeling::Remodeling::Fe, AnalyticalDirichletBC::ApproxField::MyTriFE, AnalyticalDirichletHelmholtzBC::ApproxField::MyTriFE, AnalyticalDirihletBC::ApproxField::MyTriFE, FluidPressure::MyTriangleFE, NitscheMethod::MyFace, OptimalMassTransport::VolumeFE, AnalyticalDirihletBC::ApproxField::MyTriFE, CohesiveElement::CohesiveInterfaceElement::MyPrism, EdgeForce::MyFE, BodyForceConstantField::MyVolumeFE, ThermalStressElement::MyVolumeFE, TestFE, and CellEngineering::DispMapFe.

Definition at line 309 of file ForcesAndSourcesCore.hpp.

{ return 2*order; }

getRule() [2/2]

virtual int MoFEM::ForcesAndSourcesCore::getRule ( int  order_row, int  order_col, int  order_data  )

virtual

another variant of getRule

Parameters

order_row	order of base function on row
order_col	order of base function on columns
order_data	order of base function approximating data

Returns

integration rule

Bug:

this function should be const

Definition at line 330 of file ForcesAndSourcesCore.hpp.

                                                                  {
    return getRuleHook? getRuleHook(order_row,order_col,order_data) : getRule(order_data);
  }

getSense()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getSense	(	EntityType	type,
		boost::ptr_vector< DataForcesAndSourcesCore::EntData > &	data
	)		const

get sense (orientation) of entity

Parameters

type	type of entity
data	entity data

Returns

error code

Definition at line 101 of file ForcesAndSourcesCore.cpp.

                                                                  {
  MoFEMFunctionBeginHot;
  try {
    SideNumber_multiIndex &side_table = const_cast<SideNumber_multiIndex &>(
        numeredEntFiniteElementPtr->getSideNumberTable());
    if (data.size() < side_table.get<2>().count(type)) {
      // prims has 9 edges, some of edges for "flat" prism are not active
      SETERRQ2(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
               "data inconsistency %u < %u", data.size(),
               side_table.get<2>().count(type));
    }
    SideNumber_multiIndex::nth_index<2>::type::iterator siit =
        side_table.get<2>().lower_bound(type);
    SideNumber_multiIndex::nth_index<2>::type::iterator hi_siit =
        side_table.get<2>().upper_bound(type);
    for (; siit != hi_siit; siit++) {
      data[siit->get()->side_number].getSense() = siit->get()->sense;
      if (siit->get()->brother_side_number != -1) {
        if (data.size() < (unsigned)siit->get()->brother_side_number) {
          SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                  "data inconsistency");
        }
        data[siit->get()->brother_side_number].getSense() = siit->get()->sense;
      }
    }
  } catch (std::exception &ex) {
    std::ostringstream ss;
    ss << "thorw in method: " << ex.what() << " at line " << __LINE__
       << " in file " << __FILE__;
    SETERRQ(PETSC_COMM_SELF, MOFEM_STD_EXCEPTION_THROW, ss.str().c_str());
  }

  MoFEMFunctionReturnHot(0);
}

getSpacesAndBaseOnEntities()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getSpacesAndBaseOnEntities

(

DataForcesAndSourcesCore &

data

)

const

Get field approximation space and base on entities.

Definition at line 1342 of file ForcesAndSourcesCore.cpp.

                                          {
  MoFEMFunctionBeginHot;
  if (nInTheLoop == 0) {
    data.sPace.reset();
    data.bAse.reset();
    for (EntityType t = MBVERTEX; t != MBMAXTYPE; ++t) {
      data.spacesOnEntities[t].reset();
      data.basesOnEntities[t].reset();
    }
    for (int s = 0; s != LASTSPACE; ++s) {
      data.basesOnSpaces[s].reset();
    }
  }
  for (_IT_GET_FEDATA_DOFS_FOR_LOOP_(this, dof)) {
    if (dof->get()->getEntDofIdx() != 0)
      continue;
    const EntityType type = dof->get()->getEntType();
    const FieldSpace space = dof->get()->getSpace();
    const FieldApproximationBase approx = dof->get()->getApproxBase();
    data.sPace.set(space);
    data.bAse.set(approx);
    data.spacesOnEntities[type].set(space);
    data.basesOnEntities[type].set(approx);
    data.basesOnSpaces[space].set(approx);
  }
  MoFEMFunctionReturnHot(0);
}

getTetDataOrder()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getTetDataOrder	(	DataForcesAndSourcesCore &	data,
		const FieldSpace	space
	)		const

Definition at line 281 of file ForcesAndSourcesCore.cpp.

                                                                    {

  MoFEMFunctionBeginHot;
  ierr = getDataOrder(MBTET, space, data.dataOnEntities[MBTET]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getTetDataOrderSpaceAndBase()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getTetDataOrderSpaceAndBase	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

Definition at line 410 of file ForcesAndSourcesCore.cpp.

                                                                       {

  MoFEMFunctionBeginHot;
  ierr =
      getDataOrderSpaceAndBase(field_name, MBTET, data.dataOnEntities[MBTET]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getTetsColIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getTetsColIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get Tets col indices from FENumeredDofEntity_multiIndex

Definition at line 687 of file ForcesAndSourcesCore.cpp.

                                                                           {

  MoFEMFunctionBeginHot;
  if (data.dataOnEntities[MBTET].size() == 0) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
  }
  ierr = getTypeIndices(field_name,
                        const_cast<FENumeredDofEntity_multiIndex &>(
                            numeredEntFiniteElementPtr->getColsDofs()),
                        MBTET, 0, data.dataOnEntities[MBTET][0].getIndices(),
                        data.dataOnEntities[MBTET][0].getLocalIndices());
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getTetsFieldData()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getTetsFieldData	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

Definition at line 1206 of file ForcesAndSourcesCore.cpp.

                                                                          {

  MoFEMFunctionBeginHot;
  if (data.dataOnEntities[MBTET].size() == 0) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
  }
  ierr =
      getTypeFieldData(field_name,
                       const_cast<FEDofEntity_multiIndex &>(
                           numeredEntFiniteElementPtr->getDataDofs()),
                       MBTET, 0, data.dataOnEntities[MBTET][0].getFieldData(),
                       data.dataOnEntities[MBTET][0].getFieldDofs());
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getTetsRowIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getTetsRowIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get Tets row indices from FENumeredDofEntity_multiIndex

Definition at line 670 of file ForcesAndSourcesCore.cpp.

                                                                           {

  MoFEMFunctionBeginHot;
  if (data.dataOnEntities[MBTET].size() == 0) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
  }
  ierr = getTypeIndices(field_name,
                        const_cast<FENumeredDofEntity_multiIndex &>(
                            numeredEntFiniteElementPtr->getRowsDofs()),
                        MBTET, 0, data.dataOnEntities[MBTET][0].getIndices(),
                        data.dataOnEntities[MBTET][0].getLocalIndices());
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getTrisColIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getTrisColIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get Tris col indices from FENumeredDofEntity_multiIndex

Definition at line 657 of file ForcesAndSourcesCore.cpp.

                                                                           {

  MoFEMFunctionBeginHot;
  ierr = getTypeIndices(field_name,
                        const_cast<FENumeredDofEntity_multiIndex &>(
                            numeredEntFiniteElementPtr->getColsDofs()),
                        MBTRI, data.dataOnEntities[MBTRI]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getTrisDataOrder()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getTrisDataOrder	(	DataForcesAndSourcesCore &	data,
		const FieldSpace	space
	)		const

Definition at line 261 of file ForcesAndSourcesCore.cpp.

                                                                     {

  MoFEMFunctionBeginHot;
  ierr = getDataOrder(MBTRI, space, data.dataOnEntities[MBTRI]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getTrisDataOrderSpaceAndBase()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getTrisDataOrderSpaceAndBase	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

Definition at line 390 of file ForcesAndSourcesCore.cpp.

                                                                       {

  MoFEMFunctionBeginHot;
  ierr =
      getDataOrderSpaceAndBase(field_name, MBTRI, data.dataOnEntities[MBTRI]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getTrisFieldData()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getTrisFieldData	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

Definition at line 1180 of file ForcesAndSourcesCore.cpp.

                                                                          {

  MoFEMFunctionBeginHot;
  ierr = getTypeFieldData(field_name,
                          const_cast<FEDofEntity_multiIndex &>(
                              numeredEntFiniteElementPtr->getDataDofs()),
                          MBTRI, data.dataOnEntities[MBTRI]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getTrisRowIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getTrisRowIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get Tris row indices from FENumeredDofEntity_multiIndex

Definition at line 644 of file ForcesAndSourcesCore.cpp.

                                                                           {

  MoFEMFunctionBeginHot;
  ierr = getTypeIndices(field_name,
                        const_cast<FENumeredDofEntity_multiIndex &>(
                            numeredEntFiniteElementPtr->getRowsDofs()),
                        MBTRI, data.dataOnEntities[MBTRI]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getTrisSense()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getTrisSense

(

DataForcesAndSourcesCore &

data

)

const

get triangle sense (orientation)

Parameters

data	structure where results are stored

Returns

error code

Definition at line 148 of file ForcesAndSourcesCore.cpp.

                                                                       {

  MoFEMFunctionBeginHot;
  ierr = getSense(MBTRI, data.dataOnEntities[MBTRI]);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getTypeFieldData() [1/2]

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getTypeFieldData	(	const boost::string_ref	field_name,
		FEDofEntity_multiIndex &	dofs,
		EntityType	type,
		int	side_number,
		VectorDouble &	ent_field_data,
		VectorDofs &	ent_field_dofs
	)		const

Get field data on entities.

Parameters

field_name	Field name
dofs	Dofs (element) multi index
type	Entity type
side_number	Side number (Local number of entity on element in canonical order)
ent_field_data	Vector of DOFs values on entities
ent_field_dofs	Vector of pointers to DOFs data structure

Returns

Error code

Definition at line 1025 of file ForcesAndSourcesCore.cpp.

                                      {
  //
  MoFEMFunctionBeginHot;
  FEDofEntity_multiIndex::index<
      Composite_Name_Type_And_Side_Number_mi_tag>::type::iterator dit,
      hi_dit;
  dit = dofs.get<Composite_Name_Type_And_Side_Number_mi_tag>().lower_bound(
      boost::make_tuple(field_name, type, side_number));
  if (dit == dofs.get<Composite_Name_Type_And_Side_Number_mi_tag>().end()) {
    ent_field_data.resize(0, false);
    ent_field_dofs.resize(0, false);
    MoFEMFunctionReturnHot(0);
  }
  hi_dit = dofs.get<Composite_Name_Type_And_Side_Number_mi_tag>().upper_bound(
      boost::make_tuple(field_name, type, side_number));
  if (dit != hi_dit) {
    ent_field_data.resize((*dit)->getNbDofsOnEnt(), false);
    ent_field_dofs.resize((*dit)->getNbDofsOnEnt(), false);
    for (; dit != hi_dit; dit++) {
      const FieldData val = (*dit)->getFieldData();
      const int idx = (*dit)->getEntDofIdx();
      ent_field_data[idx] = val;
      ent_field_dofs[idx] = *dit;
    }
  } else {
    ent_field_data.resize(0, false);
    ent_field_dofs.resize(0, false);
  }
  MoFEMFunctionReturnHot(0);
}

getTypeFieldData() [2/2]

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getTypeFieldData	(	const boost::string_ref	field_name,
		FEDofEntity_multiIndex &	dofs,
		EntityType	type,
		boost::ptr_vector< DataForcesAndSourcesCore::EntData > &	data
	)		const

Definition at line 1059 of file ForcesAndSourcesCore.cpp.

                                                                  {

  MoFEMFunctionBeginHot;
  SideNumber_multiIndex &side_table = const_cast<SideNumber_multiIndex &>(
      numeredEntFiniteElementPtr->getSideNumberTable());
  // if(data.size() < side_table.get<2>().count(type)) {
  //   SETERRQ(PETSC_COMM_SELF,MOFEM_DATA_INCONSISTENCY,"data inconsistency");
  // }
  SideNumber_multiIndex::nth_index<2>::type::iterator siit, hi_siit;
  siit = side_table.get<2>().lower_bound(type);
  hi_siit = side_table.get<2>().upper_bound(type);
  const boost::tuple<const boost::string_ref, EntityType> tuple =
      boost::make_tuple(field_name, type);
  FEDofEntityByNameAndType::iterator dit, hi_dit;
  dit = dofs.get<Composite_Name_And_Type_mi_tag>().lower_bound(tuple);
  if (dit == dofs.get<Composite_Name_And_Type_mi_tag>().end()) {
    for (SideNumber_multiIndex::nth_index<2>::type::iterator siiit = siit;
         siiit != hi_siit; siiit++) {
      const int side_number = siiit->get()->side_number;
      data[side_number].getFieldData().resize(0, false);
      data[side_number].getFieldDofs().resize(0, false);
    }
    MoFEMFunctionReturnHot(0);
  }
  for (SideNumber_multiIndex::nth_index<2>::type::iterator siiit = siit;
       siiit != hi_siit; siiit++) {
    const int side_number = siiit->get()->side_number;
    data[side_number].semaphore = false;
  }
  hi_dit = dofs.get<Composite_Name_And_Type_mi_tag>().upper_bound(tuple);
  for (; dit != hi_dit; dit++) {
    const int side = dit->get()->sideNumberPtr->side_number;
    const int nb_dofs_on_ent = (*dit)->getNbDofsOnEnt();
    VectorDouble &ent_field_data = data[side].getFieldData();
    VectorDofs &ent_field_dofs = data[side].getFieldDofs();
    if (!data[side].semaphore) {
      data[side].semaphore = true;
      ent_field_data.resize(nb_dofs_on_ent, false);
      ent_field_dofs.resize(nb_dofs_on_ent, false);
    }
    if (!nb_dofs_on_ent) {
      continue;
    }
    const int idx = dit->get()->getEntDofIdx();
    ent_field_data[idx] = dit->get()->getFieldData();
    ent_field_dofs[idx] = *dit;
  }
  for (; siit != hi_siit; siit++) {
    const int side_number = siit->get()->side_number;
    if (!data[side_number].semaphore) {
      data[side_number].getFieldData().resize(0, false);
      data[side_number].getFieldDofs().resize(0, false);
    }
    if (siit->get()->brother_side_number != -1) {
      if (data.size() < (unsigned int)siit->get()->brother_side_number) {
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "data inconsistency");
      }
      ierr = getTypeFieldData(
          field_name, dofs, type, side_number,
          data[siit->get()->brother_side_number].getFieldData(),
          data[siit->get()->brother_side_number].getFieldDofs());
      CHKERRG(ierr);
    }
  }
  MoFEMFunctionReturnHot(0);
}

getTypeIndices() [1/2]

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getTypeIndices	(	const boost::string_ref	field_name,
		FENumeredDofEntity_multiIndex &	dofs,
		EntityType	type,
		int	side_number,
		VectorInt &	indices,
		VectorInt &	local_indices
	)		const

get indices by type (generic function)

Definition at line 518 of file ForcesAndSourcesCore.cpp.

                                    {
  //
  MoFEMFunctionBeginHot;
  FENumeredDofEntity_multiIndex::index<
      Composite_Name_Type_And_Side_Number_mi_tag>::type::iterator dit,
      hi_dit;
  dit = dofs.get<Composite_Name_Type_And_Side_Number_mi_tag>().lower_bound(
      boost::make_tuple(field_name, type, side_number));
  if (dit == dofs.get<Composite_Name_Type_And_Side_Number_mi_tag>().end()) {
    indices.resize(0);
    local_indices.resize(0);
    MoFEMFunctionReturnHot(0);
  }
  hi_dit = dofs.get<Composite_Name_Type_And_Side_Number_mi_tag>().upper_bound(
      boost::make_tuple(field_name, type, side_number));
  if (dit != hi_dit) {
    indices.resize((*dit)->getNbDofsOnEnt(), false);
    local_indices.resize((*dit)->getNbDofsOnEnt(), false);
    for (; dit != hi_dit; dit++) {
      int idx = (*dit)->getPetscGlobalDofIdx();
      int elem_idx = (*dit)->getEntDofIdx();
      indices[elem_idx] = idx;
      int local_idx = (*dit)->getPetscLocalDofIdx();
      local_indices[elem_idx] = local_idx;
    }
  } else {
    indices.resize(0);
    local_indices.resize(0);
  }
  MoFEMFunctionReturnHot(0);
}

getTypeIndices() [2/2]

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getTypeIndices	(	const boost::string_ref	field_name,
		FENumeredDofEntity_multiIndex &	dofs,
		EntityType	type,
		boost::ptr_vector< DataForcesAndSourcesCore::EntData > &	data
	)		const

get indices by type (generic function)

Definition at line 553 of file ForcesAndSourcesCore.cpp.

                                                                  {

  MoFEMFunctionBeginHot;
  SideNumber_multiIndex &side_table = const_cast<SideNumber_multiIndex &>(
      numeredEntFiniteElementPtr->getSideNumberTable());
  SideNumber_multiIndex::nth_index<2>::type::iterator siit =
      side_table.get<2>().lower_bound(type);
  SideNumber_multiIndex::nth_index<2>::type::iterator hi_siit =
      side_table.get<2>().upper_bound(type);
  const boost::tuple<const boost::string_ref, EntityType> tuple =
      boost::make_tuple(field_name, type);
  FENumeredDofEntityByNameAndType::iterator dit, hi_dit;
  dit = dofs.get<Composite_Name_And_Type_mi_tag>().lower_bound(tuple);
  if (dit == dofs.get<Composite_Name_And_Type_mi_tag>().end()) {
    for (SideNumber_multiIndex::nth_index<2>::type::iterator siiit = siit;
         siiit != hi_siit; siiit++) {
      data[siiit->get()->side_number].getIndices().resize(0, false);
      data[siiit->get()->side_number].getLocalIndices().resize(0, false);
    }
    MoFEMFunctionReturnHot(0);
  }
  for (SideNumber_multiIndex::nth_index<2>::type::iterator siiit = siit;
       siiit != hi_siit; siiit++) {
    const int side_number = siiit->get()->side_number;
    data[side_number].semaphore = false;
  }
  hi_dit = dofs.get<Composite_Name_And_Type_mi_tag>().upper_bound(tuple);
  for (; dit != hi_dit; dit++) {
    const int side = dit->get()->sideNumberPtr->side_number;
    const int nb_dofs_on_ent = dit->get()->getNbDofsOnEnt();
    VectorInt &indices = data[side].getIndices();
    VectorInt &local_indices = data[side].getLocalIndices();
    if (!data[side].semaphore) {
      data[side].semaphore = true;
      indices.resize(nb_dofs_on_ent, false);
      local_indices.resize(nb_dofs_on_ent, false);
    }
    if (!nb_dofs_on_ent) {
      continue;
    }
    const int idx = dit->get()->getEntDofIdx();
    indices[idx] = dit->get()->getPetscGlobalDofIdx();
    local_indices[idx] = dit->get()->getPetscLocalDofIdx();
  }
  for (; siit != hi_siit; siit++) {
    const int side_number = siit->get()->side_number;
    if (!data[side_number].semaphore) {
      data[side_number].getIndices().resize(0, false);
      data[side_number].getLocalIndices().resize(0, false);
    }
    if (siit->get()->brother_side_number != -1) {
      ierr = getTypeIndices(
          field_name, dofs, type, side_number,
          data[siit->get()->brother_side_number].getIndices(),
          data[siit->get()->brother_side_number].getLocalIndices());
      CHKERRG(ierr);
    }
  }
  MoFEMFunctionReturnHot(0);
}

operator()()

virtual MoFEMErrorCode MoFEM::ForcesAndSourcesCore::operator() ( )

virtual

function is run for every finite element

It is used to calculate element local matrices and assembly. It can be used for post-processing.

Reimplemented from MoFEM::BasicMethod.

Reimplemented in MoFEM::FaceElementForcesAndSourcesCore, MoFEM::VolumeElementForcesAndSourcesCore, MoFEM::FatPrismElementForcesAndSourcesCore, MoFEM::FlatPrismElementForcesAndSourcesCore, MoFEM::EdgeElementForcesAndSourcesCore, NeummanForcesSurfaceComplexForLazy::MyTriangleSpatialFE, and MoFEM::VertexElementForcesAndSourcesCore.

Definition at line 629 of file ForcesAndSourcesCore.hpp.

                                      {
    MoFEMFunctionBeginHot;
    if(operatorHook) {
      ierr = operatorHook(); CHKERRG(ierr);
    }
    MoFEMFunctionReturnHot(0);
  }

postProcess()

virtual MoFEMErrorCode MoFEM::ForcesAndSourcesCore::postProcess ( )

virtual

function is run at the end of loop

It is used to assembly matrices and vectors, calculating global variables, f.e. total internal energy, ect.

Iterating over dofs: Example1 iterating over dofs in row by name of the field for(IT_GET_FEROW_BY_NAME_DOFS_FOR_LOOP(this,"DISPLACEMENT",it)) { ... }

Reimplemented from MoFEM::BasicMethod.

Reimplemented in PostProcFaceOnRefinedMesh, PostProcFatPrismOnRefinedMesh, SolidShellModule::SolidShellPrismElement::PostProcFatPrismOnTriangleOnRefinedMesh, PostProcTemplateVolumeOnRefinedMesh< MoFEM::VolumeElementForcesAndSourcesCore >, GelModule::Gel::GelFE, MoFEM::FatPrismElementForcesAndSourcesCore, KelvinVoigtDamper::DamperFE, SolidShellModule::SolidShellPrismElement::SolidShellError, SolidShellModule::SolidShellPrismElement::SolidShell, Smoother::MyVolumeFE, and NonlinearElasticElement::MyVolumeFE.

Definition at line 636 of file ForcesAndSourcesCore.hpp.

                                       {
    MoFEMFunctionBeginHot;
    if(postProcessHook) {
      ierr = postProcessHook(); CHKERRG(ierr);
    }
    MoFEMFunctionReturnHot(0);
  }

preProcess()

virtual MoFEMErrorCode MoFEM::ForcesAndSourcesCore::preProcess ( )

virtual

function is run at the beginning of loop

It is used to zeroing matrices and vectors, calculation of shape functions on reference element, preprocessing boundary conditions, etc.

Reimplemented from MoFEM::BasicMethod.

Reimplemented in PostProcFaceOnRefinedMesh, PostProcFatPrismOnRefinedMesh, PostProcTemplateVolumeOnRefinedMesh< MoFEM::VolumeElementForcesAndSourcesCore >, GelModule::Gel::GelFE, MoFEM::FatPrismElementForcesAndSourcesCore, KelvinVoigtDamper::DamperFE, SolidShellModule::SolidShellPrismElement::SolidShellError, SolidShellModule::SolidShellPrismElement::SolidShell, NeummanForcesSurfaceComplexForLazy::MyTriangleSpatialFE, Smoother::MyVolumeFE, NonlinearElasticElement::MyVolumeFE, and FluidPressure::MyTriangleFE.

Definition at line 622 of file ForcesAndSourcesCore.hpp.

                                      {
    MoFEMFunctionBeginHot;
    if(preProcessHook) {
      ierr = preProcessHook(); CHKERRG(ierr);
    }
    MoFEMFunctionReturnHot(0);
  }

setGaussPts() [1/2]

virtual MoFEMErrorCode MoFEM::ForcesAndSourcesCore::setGaussPts ( int  order )

virtual

It will be removed in the future use other variant.

Reimplemented in PostProcFaceOnRefinedMesh, PostProcTemplateVolumeOnRefinedMesh< MoFEM::VolumeElementForcesAndSourcesCore >, MoFEM::VolumeElementForcesAndSourcesCoreOnSide, and NitscheMethod::MyVolumeFE.

Definition at line 336 of file ForcesAndSourcesCore.hpp.

                                                {
    MoFEMFunctionBeginHot;
    SETERRQ(PETSC_COMM_SELF,MOFEM_NOT_IMPLEMENTED,"sorry, not implemented");
    MoFEMFunctionReturnHot(0);
  }

setGaussPts() [2/2]

virtual MoFEMErrorCode MoFEM::ForcesAndSourcesCore::setGaussPts ( int  order_row, int  order_col, int  order_data  )

virtual

set user specific integration rule

This function allows for user defined integration rule. The key is to called matrix gaussPts, which is used by other MoFEM procedures. Matrix has number of rows equal to problem dimension plus one, where last index is used to store weight values. Number of columns is equal to number of integration points.

Note

This function is called if method ForcesAndSourcesCore::getRule is returning integer -1 or any other negative integer.

User sets

MatrixDouble gaussPts;

where

gaussPts.resize(dim+1,nb_gauss_pts);

number rows represents local coordinates of integration points in reference element, where last index in row is for integration weight.

Definition at line 364 of file ForcesAndSourcesCore.hpp.

                                                                                 {

    MoFEMFunctionBeginHot;
    ierr = setGaussPts(order_data); CHKERRG(ierr);
    MoFEMFunctionReturnHot(0);
  }

Member Data Documentation

getRuleHook

RuleHookFun MoFEM::ForcesAndSourcesCore::getRuleHook

Hook to get rule.

Todo:

check preferred format how works with gcc and clang, see http://www.boost.org/doc/libs/1_64_0/doc/html/function/tutorial.html#idp247873024

Definition at line 319 of file ForcesAndSourcesCore.hpp.

mField

Interface& MoFEM::ForcesAndSourcesCore::mField

Definition at line 38 of file ForcesAndSourcesCore.hpp.

opPtrVector

boost::ptr_vector<UserDataOperator> MoFEM::ForcesAndSourcesCore::opPtrVector

Definition at line 613 of file ForcesAndSourcesCore.hpp.

---

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

• src/finite_elements/ForcesAndSourcesCore.hpp
• src/finite_elements/impl/ForcesAndSourcesCore.cpp