MoFEM::ForcesAndSourcesCore::UserDataOperator Struct Reference

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...

#include <src/finite_elements/ForcesAndSourcesCore.hpp>

Inheritance diagram for MoFEM::ForcesAndSourcesCore::UserDataOperator:

Collaboration diagram for MoFEM::ForcesAndSourcesCore::UserDataOperator:

Public Types
enum	OpType { OPROW = 1 << 0, OPCOL = 1 << 1, OPROWCOL = 1 << 2, OPLAST = 1 << 3 }
	Controls loop over entities on element. More...

Public Member Functions
	UserDataOperator (const FieldSpace space, const char type=OPLAST, const bool symm=true)
	UserDataOperator (const std::string &field_name, const char type, const bool symm=true)
	UserDataOperator (const std::string &row_field_name, const std::string &col_field_name, const char type, const bool symm=true)
boost::shared_ptr< const NumeredEntFiniteElement >	getNumeredEntFiniteElementPtr () const
	Return raw pointer to NumeredEntFiniteElement. More...
EntityHandle	getFEEntityHandle () const
	Return finite element entity handle. More...
boost::weak_ptr< SideNumber >	getSideNumberPtr (const int side_number, const EntityType type)
	Get the side number pointer. More...
EntityHandle	getSideEntity (const int side_number, const EntityType type)
	Get the side entity. More...
int	getNumberOfNodesOnElement ()
	Get the number of nodes on finite element. More...
MoFEMErrorCode	getProblemRowIndices (const std::string filed_name, const EntityType type, const int side, VectorInt &indices) const
	Get row indices. More...
MoFEMErrorCode	getProblemColIndices (const std::string filed_name, const EntityType type, const int side, VectorInt &indices) const
	Get col indices. More...
const FEMethod *	getFEMethod () const
	Return raw pointer to Finite Element Method object. More...
int	getOpType () const
	Get operator types. More...
void	setOpType (const OpType type)
	Set operator type. More...
void	addOpType (const OpType type)
	Add operator type. More...
int	getNinTheLoop () const
	get number of finite element in the loop More...
int	getLoopSize () const
	get size of elements in the loop More...
const std::string &	getFEName () const
	Get name of the element. More...
Accessing KSP
const PetscData::Switches &	getDataCtx () const
const KspMethod::KSPContext	getKSPCtx () const
const SnesMethod::SNESContext	getSNESCtx () const
const TSMethod::TSContext	getTSCtx () const
Vec	getKSPf () const
Mat	getKSPA () const
Mat	getKSPB () const
Accessing SNES
Vec	getSNESf () const
Vec	getSNESx () const
Mat	getSNESA () const
Mat	getSNESB () const
DEPRECATED Vec	getSnesF () const
DEPRECATED Vec	getSnesX () const
DEPRECATED Mat	getSnesA () const
DEPRECATED Mat	getSnesB () const
Accessing TS
Vec	getTSu () const
Vec	getTSu_t () const
Vec	getTSu_tt () const
Vec	getTSf () const
Mat	getTSA () const
Mat	getTSB () const
int	getTSstep () const
double	getTStime () const
double	getTSa () const
Base funtions and integration points
MatrixDouble &	getGaussPts ()
	matrix of integration (Gauss) points for Volume Element More...
auto	getFTensor0IntegrationWeight ()
	Get integration weights. More...
Deprecated (do not use)
DEPRECATED MoFEMErrorCode	getPorblemRowIndices (const std::string filed_name, const EntityType type, const int side, VectorInt &indices) const
Public Member Functions inherited from MoFEM::DataOperator
	DataOperator (const bool symm=true)
virtual	~DataOperator ()=default
virtual MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)
	Operator for bi-linear form, usually to calculate values on left hand side. More...
virtual MoFEMErrorCode	opLhs (DataForcesAndSourcesCore &row_data, DataForcesAndSourcesCore &col_data)
virtual MoFEMErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)
	Operator for linear form, usually to calculate values on right hand side. More...
virtual MoFEMErrorCode	opRhs (DataForcesAndSourcesCore &data, const bool error_if_no_base=false)
bool	getSymm () const
	Get if operator uses symmetry of DOFs or not. More...
void	setSymm ()
	set if operator is executed taking in account symmetry More...
void	unSetSymm ()
	unset if operator is executed for non symmetric problem More...

Public Attributes
char	opType
std::string	rowFieldName
std::string	colFieldName
FieldSpace	sPace
Public Attributes inherited from MoFEM::DataOperator
bool	sYmm
	If true assume that matrix is symmetric structure. More...
std::array< bool, MBMAXTYPE >	doEntities
	If true operator is executed for entity. More...
bool &	doVertices
	\deprectaed If false skip vertices More...
bool &	doEdges
	\deprectaed If false skip edges More...
bool &	doQuads
	\deprectaed More...
bool &	doTris
	\deprectaed More...
bool &	doTets
	\deprectaed More...
bool &	doPrisms
	\deprectaed More...

Protected Member Functions
virtual MoFEMErrorCode	setPtrFE (ForcesAndSourcesCore *ptr)
ForcesAndSourcesCore *	getPtrFE () const
ForcesAndSourcesCore *	getSidePtrFE () const

Protected Attributes
ForcesAndSourcesCore *	ptrFE

Private Member Functions
MoFEMErrorCode	loopSide (const string &fe_name, ForcesAndSourcesCore *side_fe, const size_t dim)
	User call this function to loop over elements on the side of face. This function calls finite element with is operator to do calculations. More...

Friends
class	ForcesAndSourcesCore
class	EdgeElementForcesAndSourcesCoreBase
class	FaceElementForcesAndSourcesCoreBase
class	ContactPrismElementForcesAndSourcesCore

Detailed Description

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).

Note: It is assumed that operator is executed for symmetric problem. That means that is executed for not repeating entities on finite element. For example on triangle we have nodes, 3 edges and one face. Because of symmetry calculations are for: nodes-nodes, nodes-edge0, nodes-edge_1, nodes-edge_2, nodes-face, edges_1-edges_1, edges_1-edge_1, edge_1-edge_2, edge_1-edge_1, edge_1-edge_2, edge_2-edge_2, edge_1-face, edge_1-face, edges_3-face, face - face

In case of non symmetric problem in addition calculations of lower off diagonal terms. F.e. edge_1-edge_0, esges_3-edge_0, edge_3-edge_1,

In that case class variable UserDataOperator::sYmm = false;

NoteL: By default sYmm is set for symmetric problems

Examples

continuity_check_on_skeleton_with_simple_2d.cpp, forces_and_sources_testing_users_base.cpp, hcurl_curl_operator.cpp, hdiv_divergence_operator.cpp, hello_world.cpp, mesh_smoothing.cpp, prism_elements_from_surface.cpp, prism_polynomial_approximation.cpp, quad_polynomial_approximation.cpp, simple_elasticity.cpp, and simple_interface.cpp.

Definition at line 84 of file ForcesAndSourcesCore.hpp.

Member Enumeration Documentation

OpType

enum MoFEM::ForcesAndSourcesCore::UserDataOperator::OpType

Controls loop over entities on element.

OPROW is used if row vector is assembled OPCOL is usually used if column vector is assembled OPROWCOL is usually used for assemble matrices.

For typical problem like Bubnov-Galerkin OPROW and OPCOL are the same. In more general case for example for non-square matrices columns and rows could have different numeration and/or different set of fields.

Enumerator
OPROW
OPCOL
OPROWCOL
OPLAST

Definition at line 98 of file ForcesAndSourcesCore.hpp.

                {
      OPROW = 1 << 0,
      OPCOL = 1 << 1,
      OPROWCOL = 1 << 2,
      OPLAST = 1 << 3
    };

Constructor & Destructor Documentation

UserDataOperator() [1/3]

MoFEM::ForcesAndSourcesCore::UserDataOperator::UserDataOperator	(	const FieldSpace	space,
		const char	type = OPLAST,
		const bool	symm = true
	)

This Constructor is used typically when some modification base shape functions on some approx. space is applied. Operator is run for all data on space.

User has no access to field data from this operator.

Definition at line 1375 of file ForcesAndSourcesCore.cpp.

    : DataOperator(symm), opType(type), sPace(space), ptrFE(nullptr) {}

UserDataOperator() [2/3]

MoFEM::ForcesAndSourcesCore::UserDataOperator::UserDataOperator	(	const std::string &	field_name,
		const char	type,
		const bool	symm = true
	)

Definition at line 1380 of file ForcesAndSourcesCore.cpp.

    : DataOperator(symm), opType(type), rowFieldName(field_name),
      colFieldName(field_name), sPace(LASTSPACE), ptrFE(nullptr) {}

UserDataOperator() [3/3]

MoFEM::ForcesAndSourcesCore::UserDataOperator::UserDataOperator	(	const std::string &	row_field_name,
		const std::string &	col_field_name,
		const char	type,
		const bool	symm = true
	)

Definition at line 1385 of file ForcesAndSourcesCore.cpp.

    : DataOperator(symm), opType(type), rowFieldName(row_field_name),
      colFieldName(col_field_name), sPace(LASTSPACE), ptrFE(nullptr) {}

Member Function Documentation

addOpType()

void MoFEM::ForcesAndSourcesCore::UserDataOperator::addOpType

(

const OpType

type

)

Add operator type.

Definition at line 918 of file ForcesAndSourcesCore.hpp.

                                                                      {
  opType |= type;
}

getDataCtx()

const PetscData::Switches & MoFEM::ForcesAndSourcesCore::UserDataOperator::getDataCtx

(

)

const

Definition at line 935 of file ForcesAndSourcesCore.hpp.

                                                       {
  return getFEMethod()->data_ctx;
}

getFEEntityHandle()

EntityHandle MoFEM::ForcesAndSourcesCore::UserDataOperator::getFEEntityHandle

(

)

const

Return finite element entity handle.

Returns

Finite element entity handle

Examples

UnsaturatedFlow.hpp.

Definition at line 876 of file ForcesAndSourcesCore.hpp.

                                                                           {
  return getNumeredEntFiniteElementPtr()->getEnt();
}

getFEMethod()

const FEMethod * MoFEM::ForcesAndSourcesCore::UserDataOperator::getFEMethod

(

)

const

Return raw pointer to Finite Element Method object.

Examples

PoissonOperators.hpp, and UnsaturatedFlow.hpp.

Definition at line 908 of file ForcesAndSourcesCore.hpp.

                                                                        {
  return ptrFE;
}

getFEName()

const std::string & MoFEM::ForcesAndSourcesCore::UserDataOperator::getFEName

(

)

const

Get name of the element.

Definition at line 930 of file ForcesAndSourcesCore.hpp.

                                                                       {
  return getFEMethod()->feName;
}

getFTensor0IntegrationWeight()

auto MoFEM::ForcesAndSourcesCore::UserDataOperator::getFTensor0IntegrationWeight

(

)

Get integration weights.

auto t_w = getFTensor0IntegrationWeight();
for(int gg = 0; gg!=getGaussPts.size2(); ++gg) {
 // integrate something
 ++t_w;
}

Returns

FTensor::Tensor0<FTensor::PackPtr<double *, 1>>

Examples

PoissonOperators.hpp, prism_polynomial_approximation.cpp, quad_polynomial_approximation.cpp, and UnsaturatedFlow.hpp.

Definition at line 1022 of file ForcesAndSourcesCore.hpp.

                                                                        {
  return FTensor::Tensor0<FTensor::PackPtr<double *, 1>>(
      &(getGaussPts()(getGaussPts().size1() - 1, 0)));
}

getGaussPts()

MatrixDouble & MoFEM::ForcesAndSourcesCore::UserDataOperator::getGaussPts

(

)

matrix of integration (Gauss) points for Volume Element

For triangle: columns 0,1 are x,y coordinates respectively and column 2 is a weight value for example getGaussPts()(1,13) returns y coordinate of 13th Gauss point on particular volume element

For tetrahedron: columns 0,1,2 are x,y,z coordinates respectively and column 3 is a weight value for example getGaussPts()(1,13) returns y coordinate of 13th Gauss point on particular volume element

Examples

MagneticElement.hpp, PoissonOperators.hpp, and UnsaturatedFlow.hpp.

Definition at line 1018 of file ForcesAndSourcesCore.hpp.

                                                                {
  return static_cast<ForcesAndSourcesCore *>(ptrFE)->gaussPts;
}

getKSPA()

Mat MoFEM::ForcesAndSourcesCore::UserDataOperator::getKSPA

(

)

const

Definition at line 958 of file ForcesAndSourcesCore.hpp.

                                                        {
  return getFEMethod()->ksp_A;
}

getKSPB()

Mat MoFEM::ForcesAndSourcesCore::UserDataOperator::getKSPB

(

)

const

Definition at line 962 of file ForcesAndSourcesCore.hpp.

                                                        {
  return getFEMethod()->ksp_B;
}

getKSPCtx()

const KspMethod::KSPContext MoFEM::ForcesAndSourcesCore::UserDataOperator::getKSPCtx

(

)

const

Definition at line 940 of file ForcesAndSourcesCore.hpp.

                                                      {
  return getFEMethod()->ksp_ctx;
}

getKSPf()

Vec MoFEM::ForcesAndSourcesCore::UserDataOperator::getKSPf

(

)

const

Definition at line 954 of file ForcesAndSourcesCore.hpp.

                                                        {
  return getFEMethod()->ksp_f;
}

getLoopSize()

int MoFEM::ForcesAndSourcesCore::UserDataOperator::getLoopSize

(

)

const

get size of elements in the loop

Returns

loop size

Definition at line 926 of file ForcesAndSourcesCore.hpp.

                                                            {
  return getFEMethod()->getLoopSize();
}

getNinTheLoop()

int MoFEM::ForcesAndSourcesCore::UserDataOperator::getNinTheLoop

(

)

const

get number of finite element in the loop

Returns

number of finite element

Definition at line 922 of file ForcesAndSourcesCore.hpp.

                                                              {
  return getFEMethod()->getNinTheLoop();
}

getNumberOfNodesOnElement()

int MoFEM::ForcesAndSourcesCore::UserDataOperator::getNumberOfNodesOnElement

(

)

Get the number of nodes on finite element.

Returns

int

Definition at line 902 of file ForcesAndSourcesCore.hpp.

                                                                    {
  int num_nodes;
  CHKERR ptrFE->getNumberOfNodes(num_nodes);
  return num_nodes;
}

getNumeredEntFiniteElementPtr()

boost::shared_ptr< const NumeredEntFiniteElement > MoFEM::ForcesAndSourcesCore::UserDataOperator::getNumeredEntFiniteElementPtr

(

)

const

Return raw pointer to NumeredEntFiniteElement.

Examples

ElasticityMixedFormulation.hpp, and UnsaturatedFlow.hpp.

Definition at line 872 of file ForcesAndSourcesCore.hpp.

                                                                          {
  return ptrFE->numeredEntFiniteElementPtr;
};

getOpType()

int MoFEM::ForcesAndSourcesCore::UserDataOperator::getOpType

(

)

const

Get operator types.

Returns

Return operator type

Definition at line 912 of file ForcesAndSourcesCore.hpp.

{ return opType; }

getPorblemRowIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::UserDataOperator::getPorblemRowIndices	(	const std::string	filed_name,
		const EntityType	type,
		const int	side,
		VectorInt &	indices
	)		const

Definition at line 1027 of file ForcesAndSourcesCore.hpp.

                              {
  return getProblemRowIndices(filed_name, type, side, indices);
}

getProblemColIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::UserDataOperator::getProblemColIndices	(	const std::string	filed_name,
		const EntityType	type,
		const int	side,
		VectorInt &	indices
	)		const

Get col indices.

Field could be or not declared for this element but is declared for problem

Parameters

field_name
type	entity type
side	side number, any number if type is MBVERTEX

Returns

indices

NOTE: Using those indices to assemble matrix will result in error if new non-zero values need to be created.

Definition at line 1280 of file ForcesAndSourcesCore.cpp.

                              {
  MoFEMFunctionBegin;
  if (ptrFE == NULL)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");

  switch (type) {
  case MBVERTEX:
    CHKERR ptrFE->getProblemNodesColIndices(field_name, indices);
    break;
  default:
    CHKERR ptrFE->getProblemTypeColIndices(field_name, type, side, indices);
  }
  MoFEMFunctionReturn(0);
}

getProblemRowIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::UserDataOperator::getProblemRowIndices	(	const std::string	filed_name,
		const EntityType	type,
		const int	side,
		VectorInt &	indices
	)		const

Get row indices.

Field could be or not declared for this element but is declared for problem

Parameters

field_name
type	entity type
side	side number, any number if type is MBVERTEX

Returns

indices

NOTE: Using those indices to assemble matrix will result in error if new non-zero values need to be created.

Definition at line 1263 of file ForcesAndSourcesCore.cpp.

                              {
  MoFEMFunctionBegin;
  if (ptrFE == NULL)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");

  switch (type) {
  case MBVERTEX:
    CHKERR ptrFE->getProblemNodesRowIndices(field_name, indices);
    break;
  default:
    CHKERR ptrFE->getProblemTypeRowIndices(field_name, type, side, indices);
  }
  MoFEMFunctionReturn(0);
}

getPtrFE()

ForcesAndSourcesCore * MoFEM::ForcesAndSourcesCore::UserDataOperator::getPtrFE ( ) const

protected

Definition at line 1033 of file ForcesAndSourcesCore.hpp.

                                                                           {
  return ptrFE;
}

getSideEntity()

EntityHandle MoFEM::ForcesAndSourcesCore::UserDataOperator::getSideEntity	(	const int	side_number,
		const EntityType	type
	)

Get the side entity.

Note

For vertex is expection. Side basses in argument of function doWork is zero. For other entity types side can be used as argument of this function.

MoFEMErrorCode doWork(int side, EntityType type,
                    DataForcesAndSourcesCore::EntData &data) {
 MoFEMFunctionBegin;

 if (type == MBVERTEX) {
   for (int n = 0; n != number_of_nodes; ++n)
     EntityHandle ent = getSideEntity(n, type);

     // Do somthing

 } else {
   EntityHandle ent = getSideEntity(side, type);

   // Do somthing

 }

 MoFEMFunctionReturn(0);
}

Parameters

side_number
type

Definition at line 894 of file ForcesAndSourcesCore.hpp.

                                                                             {
  if (auto side_ptr = getSideNumberPtr(side_number, type).lock())
    return side_ptr->ent;
  else
    return 0;
}

getSideNumberPtr()

boost::weak_ptr< SideNumber > MoFEM::ForcesAndSourcesCore::UserDataOperator::getSideNumberPtr	(	const int	side_number,
		const EntityType	type
	)

Get the side number pointer.

Note

For vertex is expection. Side basses in argument of function doWork is zero. For other entity types side can be used as argument of this function.

Parameters

side_number
type

Returns

boost::weak_ptr<SideNumber>

Definition at line 881 of file ForcesAndSourcesCore.hpp.

                                                  {
  auto &side_table_by_side_and_type =
      ptrFE->numeredEntFiniteElementPtr->getSideNumberTable().get<1>();
  auto side_it =
      side_table_by_side_and_type.find(boost::make_tuple(type, side_number));
  if (side_it != side_table_by_side_and_type.end())
    return *side_it;
  else
    return boost::weak_ptr<SideNumber>();
}

getSidePtrFE()

ForcesAndSourcesCore * MoFEM::ForcesAndSourcesCore::UserDataOperator::getSidePtrFE ( ) const

protected

Definition at line 1038 of file ForcesAndSourcesCore.hpp.

                                                         {
  return ptrFE->sidePtrFE;
}

getSNESA()

Mat MoFEM::ForcesAndSourcesCore::UserDataOperator::getSNESA

(

)

const

Definition at line 974 of file ForcesAndSourcesCore.hpp.

                                                         {
  return getFEMethod()->snes_A;
}

getSnesA()

DEPRECATED Mat MoFEM::ForcesAndSourcesCore::UserDataOperator::getSnesA

(

)

const

Deprecated:

Use getSNESA intead

Definition at line 310 of file ForcesAndSourcesCore.hpp.

{ return getSNESA(); }

getSNESB()

Mat MoFEM::ForcesAndSourcesCore::UserDataOperator::getSNESB

(

)

const

Definition at line 978 of file ForcesAndSourcesCore.hpp.

                                                         {
  return getFEMethod()->snes_B;
}

getSnesB()

DEPRECATED Mat MoFEM::ForcesAndSourcesCore::UserDataOperator::getSnesB

(

)

const

Deprecated:

Use getSNESB intead

Definition at line 313 of file ForcesAndSourcesCore.hpp.

{ return getSNESB(); }

getSNESCtx()

const SnesMethod::SNESContext MoFEM::ForcesAndSourcesCore::UserDataOperator::getSNESCtx

(

)

const

Definition at line 945 of file ForcesAndSourcesCore.hpp.

                                                       {
  return getFEMethod()->snes_ctx;
}

getSNESf()

Vec MoFEM::ForcesAndSourcesCore::UserDataOperator::getSNESf

(

)

const

Definition at line 966 of file ForcesAndSourcesCore.hpp.

                                                         {
  return getFEMethod()->snes_f;
}

getSnesF()

DEPRECATED Vec MoFEM::ForcesAndSourcesCore::UserDataOperator::getSnesF

(

)

const

Deprecated:

Use getSNESF intead

Definition at line 304 of file ForcesAndSourcesCore.hpp.

{ return getSNESf(); }

getSNESx()

Vec MoFEM::ForcesAndSourcesCore::UserDataOperator::getSNESx

(

)

const

Definition at line 970 of file ForcesAndSourcesCore.hpp.

                                                         {
  return getFEMethod()->snes_x;
}

getSnesX()

DEPRECATED Vec MoFEM::ForcesAndSourcesCore::UserDataOperator::getSnesX

(

)

const

Deprecated:

Use getSNESX intead

Definition at line 307 of file ForcesAndSourcesCore.hpp.

{ return getSNESx(); }

getTSA()

Mat MoFEM::ForcesAndSourcesCore::UserDataOperator::getTSA

(

)

const

Definition at line 998 of file ForcesAndSourcesCore.hpp.

                                                       {
  return getFEMethod()->ts_A;
}

getTSa()

double MoFEM::ForcesAndSourcesCore::UserDataOperator::getTSa

(

)

const

Definition at line 1014 of file ForcesAndSourcesCore.hpp.

                                                          {
  return getFEMethod()->ts_a;
}

getTSB()

Mat MoFEM::ForcesAndSourcesCore::UserDataOperator::getTSB

(

)

const

Definition at line 1002 of file ForcesAndSourcesCore.hpp.

                                                       {
  return getFEMethod()->ts_B;
}

getTSCtx()

const TSMethod::TSContext MoFEM::ForcesAndSourcesCore::UserDataOperator::getTSCtx

(

)

const

Definition at line 950 of file ForcesAndSourcesCore.hpp.

                                                     {
  return getFEMethod()->ts_ctx;
}

getTSf()

Vec MoFEM::ForcesAndSourcesCore::UserDataOperator::getTSf

(

)

const

Definition at line 994 of file ForcesAndSourcesCore.hpp.

                                                       {
  return getFEMethod()->ts_F;
}

getTSstep()

int MoFEM::ForcesAndSourcesCore::UserDataOperator::getTSstep

(

)

const

Definition at line 1006 of file ForcesAndSourcesCore.hpp.

                                                          {
  return getFEMethod()->ts_step;
}

getTStime()

double MoFEM::ForcesAndSourcesCore::UserDataOperator::getTStime

(

)

const

Definition at line 1010 of file ForcesAndSourcesCore.hpp.

                                                             {
  return getFEMethod()->ts_t;
}

getTSu()

Vec MoFEM::ForcesAndSourcesCore::UserDataOperator::getTSu

(

)

const

Definition at line 982 of file ForcesAndSourcesCore.hpp.

                                                       {
  return getFEMethod()->ts_u;
}

getTSu_t()

Vec MoFEM::ForcesAndSourcesCore::UserDataOperator::getTSu_t

(

)

const

Definition at line 986 of file ForcesAndSourcesCore.hpp.

                                                         {
  return getFEMethod()->ts_u_t;
}

getTSu_tt()

Vec MoFEM::ForcesAndSourcesCore::UserDataOperator::getTSu_tt

(

)

const

Definition at line 990 of file ForcesAndSourcesCore.hpp.

                                                          {
  return getFEMethod()->ts_u_tt;
}

loopSide()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::UserDataOperator::loopSide ( const string &  fe_name, ForcesAndSourcesCore *  side_fe, const size_t  dim  )

private

User call this function to loop over elements on the side of face. This function calls finite element with is operator to do calculations.

Parameters

fe_name
side_fe
dim

Returns

MoFEMErrorCode

Definition at line 1305 of file ForcesAndSourcesCore.cpp.

                                                                        {
  MoFEMFunctionBegin;
  const EntityHandle ent = getNumeredEntFiniteElementPtr()->getEnt();
  const Problem *problem_ptr = getFEMethod()->problemPtr;
  Range adjacent_ents;
  CHKERR ptrFE->mField.getInterface<BitRefManager>()->getAdjacenciesAny(
      ent, side_dim, adjacent_ents);
  typedef NumeredEntFiniteElement_multiIndex::index<
      Composite_Name_And_Ent_mi_tag>::type FEByComposite;
  FEByComposite &numered_fe =
      problem_ptr->numeredFiniteElements->get<Composite_Name_And_Ent_mi_tag>();

  side_fe->feName = fe_name;

  CHKERR side_fe->setSideFEPtr(ptrFE);
  CHKERR side_fe->copyBasicMethod(*getFEMethod());
  CHKERR side_fe->copyKsp(*getFEMethod());
  CHKERR side_fe->copySnes(*getFEMethod());
  CHKERR side_fe->copyTs(*getFEMethod());

  CHKERR side_fe->preProcess();

  int nn = 0;
  side_fe->loopSize = adjacent_ents.size();
  for (Range::iterator vit = adjacent_ents.begin(); vit != adjacent_ents.end();
       vit++) {
    FEByComposite::iterator miit =
        numered_fe.find(boost::make_tuple(fe_name, *vit));
    if (miit != numered_fe.end()) {
      side_fe->nInTheLoop = nn++;
      side_fe->numeredEntFiniteElementPtr = *miit;
      side_fe->dataFieldEntsPtr = (*miit)->sPtr->data_field_ents_view;
      side_fe->rowFieldEntsPtr = (*miit)->sPtr->row_field_ents_view;
      side_fe->colFieldEntsPtr = (*miit)->sPtr->col_field_ents_view;
      side_fe->dataPtr = (*miit)->sPtr->data_dofs;
      side_fe->rowPtr = (*miit)->rows_dofs;
      side_fe->colPtr = (*miit)->cols_dofs;
      CHKERR (*side_fe)();
    }
  }

  CHKERR side_fe->postProcess();

  MoFEMFunctionReturn(0);
}

setOpType()

void MoFEM::ForcesAndSourcesCore::UserDataOperator::setOpType

(

const OpType

type

)

Set operator type.

Parameters

Operator

type

Definition at line 914 of file ForcesAndSourcesCore.hpp.

                                                                      {
  opType = type;
}

setPtrFE()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::UserDataOperator::setPtrFE ( ForcesAndSourcesCore *  ptr )

protectedvirtual

Reimplemented in MoFEM::FaceElementForcesAndSourcesCoreBase::UserDataOperator, MoFEM::VolumeElementForcesAndSourcesCoreBase::UserDataOperator, MoFEM::FatPrismElementForcesAndSourcesCore::UserDataOperator, MoFEM::VolumeElementForcesAndSourcesCoreOnSideBase::UserDataOperator, MoFEM::FlatPrismElementForcesAndSourcesCore::UserDataOperator, MoFEM::EdgeElementForcesAndSourcesCoreBase::UserDataOperator, MoFEM::FaceElementForcesAndSourcesCoreOnSideBase::UserDataOperator, and MoFEM::VertexElementForcesAndSourcesCore::UserDataOperator.

Definition at line 1417 of file ForcesAndSourcesCore.cpp.

                                                                        {
  MoFEMFunctionBeginHot;
  if (!(ptrFE = dynamic_cast<ForcesAndSourcesCore *>(ptr)))
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "User operator and finite element do not work together");
  MoFEMFunctionReturnHot(0);
}

Friends And Related Function Documentation

ContactPrismElementForcesAndSourcesCore

friend class ContactPrismElementForcesAndSourcesCore

friend

Definition at line 412 of file ForcesAndSourcesCore.hpp.

EdgeElementForcesAndSourcesCoreBase

friend class EdgeElementForcesAndSourcesCoreBase

friend

Definition at line 410 of file ForcesAndSourcesCore.hpp.

FaceElementForcesAndSourcesCoreBase

friend class FaceElementForcesAndSourcesCoreBase

friend

Definition at line 411 of file ForcesAndSourcesCore.hpp.

ForcesAndSourcesCore

friend class ForcesAndSourcesCore

friend

Definition at line 409 of file ForcesAndSourcesCore.hpp.

Member Data Documentation

colFieldName

std::string MoFEM::ForcesAndSourcesCore::UserDataOperator::colFieldName

Definition at line 107 of file ForcesAndSourcesCore.hpp.

opType

char MoFEM::ForcesAndSourcesCore::UserDataOperator::opType

Definition at line 105 of file ForcesAndSourcesCore.hpp.

ptrFE

ForcesAndSourcesCore* MoFEM::ForcesAndSourcesCore::UserDataOperator::ptrFE

protected

Definition at line 387 of file ForcesAndSourcesCore.hpp.

rowFieldName

std::string MoFEM::ForcesAndSourcesCore::UserDataOperator::rowFieldName

Definition at line 106 of file ForcesAndSourcesCore.hpp.

sPace

FieldSpace MoFEM::ForcesAndSourcesCore::UserDataOperator::sPace

Definition at line 108 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