11#ifndef __FORCES_AND_SOURCES_CORE__HPP__
12#define __FORCES_AND_SOURCES_CORE__HPP__
14using namespace boost::numeric;
27 typedef boost::function<int(
int order_row,
int order_col,
int order_data)>
31 int order_row,
int order_col,
177 const EntityType type,
178 boost::ptr_vector<EntitiesFieldData::EntData> &data)
const;
189 const EntityType type,
const FieldSpace space,
190 boost::ptr_vector<EntitiesFieldData::EntData> &data)
const;
199 template <EntityType type>
212 template <EntityType type>
223 template <
typename EXTRACTOR>
227 VectorInt &local_nodes_indices, EXTRACTOR &&extractor)
const;
231 const int bit_number)
const;
235 const int bit_number)
const;
237 template <
typename EXTRACTOR>
240 const EntityType type_lo,
241 const EntityType type_hi,
242 EXTRACTOR &&extractor)
const;
246 const EntityType type_lo = MBVERTEX,
247 const EntityType type_hi = MBPOLYHEDRON)
const;
251 const EntityType type_lo = MBVERTEX,
252 const EntityType type_hi = MBPOLYHEDRON)
const;
257 const int bit_number)
const;
261 const int bit_number)
const;
280 const int bit_number)
const;
284 const EntityType type_lo = MBVERTEX,
285 const EntityType type_hi = MBPOLYHEDRON)
const;
295 template <
typename EXTRACTOR>
297 const int bit_number,
298 EXTRACTOR &&extractor)
const;
323 EntityType type,
int side_number,
329 EntityType type,
int side_number,
334 EntityType type,
int side_number,
376 virtual int getRule(
int order_row,
int order_col,
int order_data);
457 const std::array<boost::shared_ptr<EntitiesFieldData>,
LASTSPACE>
464 const std::array<boost::shared_ptr<EntitiesFieldData>,
LASTSPACE>
586 const bool symm =
true);
589 const bool symm =
true);
592 const std::string col_field_name,
const char type,
593 const bool symm =
true);
597 inline boost::shared_ptr<const NumeredEntFiniteElement>
632 const EntityType type);
667 const EntityType type);
691 const EntityType type,
const int side,
709 const EntityType type,
const int side,
791 inline Vec
getTSu()
const;
797 inline Vec
getTSf()
const;
799 inline Mat
getTSA()
const;
801 inline Mat
getTSB()
const;
809 inline double getTSa()
const;
901 std::vector<boost::weak_ptr<NumeredEntFiniteElement>>>;
920 boost::shared_ptr<Range> fe_range =
nullptr,
921 const int verb =
QUIET,
940 const int verb =
QUIET,
955 const int verb =
QUIET,
970 const int verb =
QUIET,
985 boost::shared_ptr<Range> fe_range,
986 const int verb =
QUIET,
1013boost::shared_ptr<const NumeredEntFiniteElement>
1014ForcesAndSourcesCore::UserDataOperator::getNumeredEntFiniteElementPtr()
const {
1018EntityHandle ForcesAndSourcesCore::UserDataOperator::getFEEntityHandle()
const {
1019 return getNumeredEntFiniteElementPtr()->getEnt();
1022int ForcesAndSourcesCore::UserDataOperator::getFEDim()
const {
1026EntityType ForcesAndSourcesCore::UserDataOperator::getFEType()
const {
1030boost::weak_ptr<SideNumber>
1031ForcesAndSourcesCore::UserDataOperator::getSideNumberPtr(
1032 const int side_number,
const EntityType type) {
1033 auto &side_table_by_side_and_type =
1034 ptrFE->numeredEntFiniteElementPtr->getSideNumberTable().get<1>();
1036 side_table_by_side_and_type.find(boost::make_tuple(type, side_number));
1037 if (side_it != side_table_by_side_and_type.end())
1040 return boost::weak_ptr<SideNumber>();
1044ForcesAndSourcesCore::UserDataOperator::getSideEntity(
const int side_number,
1045 const EntityType type) {
1046 if (
auto side_ptr = getSideNumberPtr(side_number, type).lock())
1047 return side_ptr->ent;
1052int ForcesAndSourcesCore::UserDataOperator::getNumberOfNodesOnElement()
const {
1053 return ptrFE->getNumberOfNodes();
1056const FEMethod *ForcesAndSourcesCore::UserDataOperator::getFEMethod()
const {
1060int ForcesAndSourcesCore::UserDataOperator::getOpType()
const {
return opType; }
1062void ForcesAndSourcesCore::UserDataOperator::setOpType(
const OpType type) {
1066void ForcesAndSourcesCore::UserDataOperator::addOpType(
const OpType type) {
1070int ForcesAndSourcesCore::UserDataOperator::getNinTheLoop()
const {
1071 return getFEMethod()->getNinTheLoop();
1074int ForcesAndSourcesCore::UserDataOperator::getLoopSize()
const {
1075 return getFEMethod()->getLoopSize();
1078std::string ForcesAndSourcesCore::UserDataOperator::getFEName()
const {
1079 return getFEMethod()->getFEName();
1083ForcesAndSourcesCore::UserDataOperator::getDataCtx()
const {
1084 return getFEMethod()->data_ctx;
1088ForcesAndSourcesCore::UserDataOperator::getKSPCtx()
const {
1089 return getFEMethod()->ksp_ctx;
1093ForcesAndSourcesCore::UserDataOperator::getSNESCtx()
const {
1094 return getFEMethod()->snes_ctx;
1098ForcesAndSourcesCore::UserDataOperator::getTSCtx()
const {
1099 return getFEMethod()->ts_ctx;
1102Vec ForcesAndSourcesCore::UserDataOperator::getKSPf()
const {
1104 if (getFEMethod()->
ksp_f == PETSC_NULLPTR)
1107 return getFEMethod()->ksp_f;
1110Mat ForcesAndSourcesCore::UserDataOperator::getKSPA()
const {
1112 if (getFEMethod()->
ksp_A == PETSC_NULLPTR)
1115 return getFEMethod()->ksp_A;
1118Mat ForcesAndSourcesCore::UserDataOperator::getKSPB()
const {
1120 if (getFEMethod()->
ksp_B == PETSC_NULLPTR)
1123 return getFEMethod()->ksp_B;
1126Vec ForcesAndSourcesCore::UserDataOperator::getSNESf()
const {
1128 if (getFEMethod()->
snes_f == PETSC_NULLPTR)
1131 return getFEMethod()->snes_f;
1134Vec ForcesAndSourcesCore::UserDataOperator::getSNESx()
const {
1136 if (getFEMethod()->
snes_x == PETSC_NULLPTR)
1139 return getFEMethod()->snes_x;
1142Mat ForcesAndSourcesCore::UserDataOperator::getSNESA()
const {
1144 if (getFEMethod()->
snes_A == PETSC_NULLPTR)
1147 return getFEMethod()->snes_A;
1150Mat ForcesAndSourcesCore::UserDataOperator::getSNESB()
const {
1152 if (getFEMethod()->
snes_B == PETSC_NULLPTR)
1155 return getFEMethod()->snes_B;
1158Vec ForcesAndSourcesCore::UserDataOperator::getTSu()
const {
1160 if (getFEMethod()->
ts_u == PETSC_NULLPTR)
1163 return getFEMethod()->ts_u;
1166Vec ForcesAndSourcesCore::UserDataOperator::getTSu_t()
const {
1168 if (getFEMethod()->
ts_u_t == PETSC_NULLPTR)
1171 return getFEMethod()->ts_u_t;
1174Vec ForcesAndSourcesCore::UserDataOperator::getTSu_tt()
const {
1176 if (getFEMethod()->
ts_u_tt == PETSC_NULLPTR)
1179 return getFEMethod()->ts_u_tt;
1182Vec ForcesAndSourcesCore::UserDataOperator::getTSf()
const {
1184 if (getFEMethod()->
ts_F == PETSC_NULLPTR)
1187 return getFEMethod()->ts_F;
1190Mat ForcesAndSourcesCore::UserDataOperator::getTSA()
const {
1192 if (getFEMethod()->
ts_A == PETSC_NULLPTR)
1195 return getFEMethod()->ts_A;
1198Mat ForcesAndSourcesCore::UserDataOperator::getTSB()
const {
1200 if (getFEMethod()->
ts_B == PETSC_NULLPTR)
1203 return getFEMethod()->ts_B;
1206int ForcesAndSourcesCore::UserDataOperator::getTSstep()
const {
1208 if ((getFEMethod()->
data_ctx & PetscData::PetscData::CtxSetTime).none())
1211 return getFEMethod()->ts_step;
1214double ForcesAndSourcesCore::UserDataOperator::getTStime()
const {
1216 if ((getFEMethod()->
data_ctx & PetscData::PetscData::CtxSetTime).none())
1219 return getFEMethod()->ts_t;
1222double ForcesAndSourcesCore::UserDataOperator::getTStimeStep()
const {
1224 if ((getFEMethod()->
data_ctx & PetscData::PetscData::CtxSetTime).none())
1227 return getFEMethod()->ts_dt;
1230double ForcesAndSourcesCore::UserDataOperator::getTSa()
const {
1237 return getFEMethod()->ts_a;
1240double ForcesAndSourcesCore::UserDataOperator::getTSaa()
const {
1247 return getFEMethod()->ts_aa;
1254auto ForcesAndSourcesCore::UserDataOperator::getFTensor0IntegrationWeight() {
1256 &(getGaussPts()(getGaussPts().size1() - 1, 0)));
1270ForcesAndSourcesCore::UserDataOperator::getSidePtrFE()
const {
1275ForcesAndSourcesCore::UserDataOperator::getRefinePtrFE()
const {
1279MatrixDouble &ForcesAndSourcesCore::UserDataOperator::getCoordsAtGaussPts() {
1283auto ForcesAndSourcesCore::UserDataOperator::getFTensor1CoordsAtGaussPts() {
1285 &getCoordsAtGaussPts()(0, 0), &getCoordsAtGaussPts()(0, 1),
1286 &getCoordsAtGaussPts()(0, 2));
1289double ForcesAndSourcesCore::UserDataOperator::getMeasure()
const {
1293double &ForcesAndSourcesCore::UserDataOperator::getMeasure() {
1303template <
typename E>
1309 const int side_dim, boost::shared_ptr<Range> fe_range,
1311 boost::shared_ptr<AdjCache> adj_cache =
nullptr)
1313 sideFEName(fe_name), sideDim(side_dim), sevLevel(sev),
1314 adjCache(adj_cache), feRange(fe_range) {}
1328 boost::shared_ptr<AdjCache> adj_cache =
nullptr)
1329 :
OpLoopSide(m_field, fe_name, side_dim, nullptr, sev, adj_cache) {}
1333 return loopSide(sideFEName, sideFEPtr.get(), sideDim, 0, feRange,
VERBOSE,
1334 sevLevel, adjCache.get());
1338 return sideFEPtr->getOpPtrVector();
1357template <
typename E>
1363 boost::shared_ptr<Range> fe_range,
1366 rangeFEName(fe_name), sevLevel(sev), feRange(fe_range) {}
1370 return loopRange(rangeFEName, rangeFEPtr.get(), feRange,
VERBOSE, sevLevel);
1374 return rangeFEPtr->getOpPtrVector();
FieldApproximationBase
approximation base
FieldSpace
approximation spaces
@ LASTSPACE
FieldSpace in [ 0, LASTSPACE )
#define THROW_MESSAGE(msg)
Throw MoFEM exception.
multi_index_container< boost::shared_ptr< NumeredDofEntity >, indexed_by< ordered_unique< tag< Unique_mi_tag >, const_mem_fun< NumeredDofEntity::interface_type_DofEntity, UId, &NumeredDofEntity::getLocalUniqueId > >, ordered_non_unique< tag< Part_mi_tag >, member< NumeredDofEntity, unsigned int, &NumeredDofEntity::pArt > >, ordered_non_unique< tag< Idx_mi_tag >, member< NumeredDofEntity, DofIdx, &NumeredDofEntity::dofIdx > >, ordered_non_unique< tag< PetscGlobalIdx_mi_tag >, member< NumeredDofEntity, DofIdx, &NumeredDofEntity::petscGloablDofIdx > >, ordered_non_unique< tag< PetscLocalIdx_mi_tag >, member< NumeredDofEntity, DofIdx, &NumeredDofEntity::petscLocalDofIdx > >, ordered_non_unique< tag< Ent_mi_tag >, const_mem_fun< NumeredDofEntity::interface_type_DofEntity, EntityHandle, &NumeredDofEntity::getEnt > > > > NumeredDofEntity_multiIndex
MultiIndex container keeps NumeredDofEntity.
SeverityLevel
Severity levels.
PetscErrorCode MoFEMErrorCode
MoFEM/PETSc error code.
UBlasVector< int > VectorInt
implementation of Data Operators for Forces and Sources
auto type_from_handle(const EntityHandle h)
get type from entity handle
DEPRECATED typedef ForcesAndSourcesCore ForcesAndSurcesCore
auto dimension_from_handle(const EntityHandle h)
get entity dimension form handle
std::vector< boost::weak_ptr< FieldEntity > > FieldEntity_vector_view
constexpr auto field_name
base operator to do operations at Gauss Pt. level
Deprecated interface functions.
Base face element used to integrate on skeleton.
Data on single entity (This is passed as argument to DataOperator::doWork)
data structure for finite element entity
std::array< boost::ptr_vector< EntData >, MBMAXTYPE > dataOnEntities
structure for User Loop Methods on finite elements
EntityHandle getFEEntityHandle() const
boost::shared_ptr< const NumeredEntFiniteElement > numeredEntFiniteElementPtr
Base face element used to integrate on skeleton.
Base face element used to integrate on skeleton.
void setOpType(const OpType type)
Set operator type.
EntityHandle getFEEntityHandle() const
Return finite element entity handle.
std::map< EntityHandle, std::vector< boost::weak_ptr< NumeredEntFiniteElement > > > AdjCache
MoFEMErrorCode getProblemColIndices(const std::string filed_name, const EntityType type, const int side, VectorInt &indices) const
Get col indices.
int getLoopSize() const
get size of elements in the loop
EntityType getFEType() const
Get dimension of finite element.
const KspMethod::KSPContext getKSPCtx() const
static const char *const OpTypeNames[]
MoFEMErrorCode loopSide(const string &fe_name, ForcesAndSourcesCore *side_fe, const size_t dim, const EntityHandle ent_for_side=0, boost::shared_ptr< Range > fe_range=nullptr, const int verb=QUIET, const LogManager::SeverityLevel sev=Sev::noisy, AdjCache *adj_cache=nullptr)
User calls this function to loop over elements on the side of face. This function calls finite elemen...
MoFEMErrorCode loopRange(const string &fe_name, ForcesAndSourcesCore *range_fe, boost::shared_ptr< Range > fe_range, const int verb=QUIET, const LogManager::SeverityLevel sev=Sev::noisy)
Iterate over range of elements.
int getNinTheLoop() const
get number of finite element in the loop
auto getFTensor1CoordsAtGaussPts()
Get coordinates at integration points assuming linear geometry.
int getNumberOfNodesOnElement() const
Get the number of nodes on finite element.
boost::weak_ptr< SideNumber > getSideNumberPtr(const int side_number, const EntityType type)
Get the side number pointer.
const SnesMethod::SNESContext getSNESCtx() const
MoFEMErrorCode loopParent(const string &fe_name, ForcesAndSourcesCore *parent_fe, const int verb=QUIET, const LogManager::SeverityLevel sev=Sev::noisy)
User calls this function to loop over parent elements. This function calls finite element with its op...
const PetscData::Switches & getDataCtx() const
ForcesAndSourcesCore * ptrFE
EntityHandle getSideEntity(const int side_number, const EntityType type)
Get the side entity.
MatrixDouble & getCoordsAtGaussPts()
Gauss points and weight, matrix (nb. of points x 3)
auto getFTensor0IntegrationWeight()
Get integration weights.
const TSMethod::TSContext getTSCtx() const
double getMeasure() const
get measure of element
std::string getFEName() const
Get name of the element.
ForcesAndSourcesCore * getRefinePtrFE() const
MoFEMErrorCode getProblemRowIndices(const std::string filed_name, const EntityType type, const int side, VectorInt &indices) const
Get row indices.
boost::shared_ptr< const NumeredEntFiniteElement > getNumeredEntFiniteElementPtr() const
Return raw pointer to NumeredEntFiniteElement.
double getTStimeStep() const
void addOpType(const OpType type)
Add operator type.
ForcesAndSourcesCore * getPtrFE() const
int getOpType() const
Get operator types.
OpType
Controls loop over entities on element.
@ OPCOL
operator doWork function is executed on FE columns
@ OPROW
operator doWork function is executed on FE rows
@ OPROWCOL
operator doWork is executed on FE rows &columns
@ OPSPACE
operator do Work is execute on space data
MoFEMErrorCode loopChildren(const string &fe_name, ForcesAndSourcesCore *child_fe, const int verb=QUIET, const LogManager::SeverityLevel sev=Sev::noisy)
User calls this function to loop over parent elements. This function calls finite element with its op...
virtual MoFEMErrorCode setPtrFE(ForcesAndSourcesCore *ptr)
const FEMethod * getFEMethod() const
Return raw pointer to Finite Element Method object.
ForcesAndSourcesCore * getSidePtrFE() const
int getFEDim() const
Get dimension of finite element.
MatrixDouble & getGaussPts()
matrix of integration (Gauss) points for Volume Element
MoFEMErrorCode loopThis(const string &fe_name, ForcesAndSourcesCore *this_fe, const int verb=QUIET, const LogManager::SeverityLevel sev=Sev::noisy)
User calls this function to loop over the same element using a different set of integration points....
structure to get information form mofem into EntitiesFieldData
MoFEMErrorCode loopOverOperators()
Iterate user data operators.
MoFEMErrorCode getProblemNodesRowIndices(const std::string &field_name, VectorInt &nodes_indices) const
EntitiesFieldData & dataNoField
auto & getDataOnElementBySpaceArray()
Get data on entities and space.
MoFEMErrorCode getEntityIndices(EntitiesFieldData &data, const int bit_number, FieldEntity_vector_view &ents_field, const EntityType type_lo, const EntityType type_hi, EXTRACTOR &&extractor) const
MoFEMErrorCode getEntityRowIndices(EntitiesFieldData &data, const int bit_number, const EntityType type_lo=MBVERTEX, const EntityType type_hi=MBPOLYHEDRON) const
MoFEMErrorCode getRowNodesIndices(EntitiesFieldData &data, const int bit_number) const
get row node indices from FENumeredDofEntity_multiIndex
int getMaxRowOrder() const
Get max order of approximation for field in rows.
auto & getUserPolynomialBase()
Get the User Polynomial Base object.
MoFEMErrorCode getEntitySense(EntitiesFieldData &data) const
Get the entity sense (orientation)
EntitiesFieldData & dataHdiv
MoFEMErrorCode getSpacesAndBaseOnEntities(EntitiesFieldData &data) const
Get field approximation space and base on entities.
MoFEMErrorCode getProblemTypeColIndices(const std::string &field_name, EntityType type, int side_number, VectorInt &indices) const
virtual MoFEMErrorCode operator()()
function is run for every finite element
MoFEMErrorCode getFaceNodes(EntitiesFieldData &data) const
Get nodes on faces.
const std::array< boost::shared_ptr< EntitiesFieldData >, LASTSPACE > derivedDataOnElement
Entity data on element entity columns fields.
MoFEMErrorCode getNoFieldRowIndices(EntitiesFieldData &data, const int bit_number) const
get col NoField indices
virtual MoFEMErrorCode setGaussPts(int order_row, int order_col, int order_data)
set user specific integration rule
MoFEMErrorCode getNodesIndices(const int bit_number, FieldEntity_vector_view &ents_field, VectorInt &nodes_indices, VectorInt &local_nodes_indices, EXTRACTOR &&extractor) const
get node indices
MoFEMErrorCode getEntityFieldData(EntitiesFieldData &data, const int bit_number, const EntityType type_lo=MBVERTEX, const EntityType type_hi=MBPOLYHEDRON) const
EntitiesFieldData & dataHcurl
MoFEMErrorCode setRefineFEPtr(const ForcesAndSourcesCore *refine_fe_ptr)
Set the pointer to face element refined.
boost::function< int(int order_row, int order_col, int order_data)> RuleHookFun
auto & getElementPolynomialBase()
Get the Entity Polynomial Base object.
ForcesAndSourcesCore * sidePtrFE
Element to integrate on the sides.
EntitiesFieldData & dataH1
MatrixDouble coordsAtGaussPts
coordinated at gauss points
virtual int getRule(int order_row, int order_col, int order_data)
another variant of getRule
boost::ptr_deque< UserDataOperator > opPtrVector
Vector of finite element users data operators.
MoFEMErrorCode getEntitySense(const EntityType type, boost::ptr_vector< EntitiesFieldData::EntData > &data) const
get sense (orientation) of entity
const std::array< boost::shared_ptr< EntitiesFieldData >, LASTSPACE > dataOnElement
Entity data on element entity rows fields.
auto & getDerivedDataOnElementBySpaceArray()
Get derived data on entities and space.
MoFEMErrorCode getProblemNodesColIndices(const std::string &field_name, VectorInt &nodes_indices) const
MoFEMErrorCode getNoFieldEntityFieldData(EntitiesFieldData &data, const int bit_number, EXTRACTOR &&extractor) const
Get field data on entities where no field is defined.
boost::function< MoFEMErrorCode(ForcesAndSourcesCore *fe_raw_ptr, int order_row, int order_col, int order_data)> GaussHookFun
MoFEMErrorCode calHierarchicalBaseFunctionsOnElement()
Calculate base functions.
MatrixDouble gaussPts
Matrix of integration points.
MoFEMErrorCode calBernsteinBezierBaseFunctionsOnElement()
Calculate Bernstein-Bezier base.
auto & getEntData(const FieldSpace space, const EntityType type, const int side)
Get the entity data.
friend class UserDataOperator
MoFEMErrorCode getEntityDataOrder(EntitiesFieldData &data, const FieldSpace space) const
Get the entity data order for given space.
virtual MoFEMErrorCode preProcess()
function is run at the beginning of loop
MoFEMErrorCode getColNodesIndices(EntitiesFieldData &data, const int bit_number) const
get col node indices from FENumeredDofEntity_multiIndex
MoFEMErrorCode setSideFEPtr(const ForcesAndSourcesCore *side_fe_ptr)
Set the pointer to face element on the side.
int getMaxColOrder() const
Get max order of approximation for field in columns.
EntityType lastEvaluatedElementEntityType
Last evaluated type of element entity.
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
boost::ptr_deque< UserDataOperator > & getOpPtrVector()
Use to push back operator for row operator.
GaussHookFun setRuleHook
Set function to calculate integration rule.
MoFEMErrorCode getBitRefLevelOnData()
MoFEMErrorCode getNodesFieldData(EntitiesFieldData &data, const int bit_number) const
Get data on nodes.
MoFEMErrorCode getEntityColIndices(EntitiesFieldData &data, const int bit_number, const EntityType type_lo=MBVERTEX, const EntityType type_hi=MBPOLYHEDRON) const
boost::shared_ptr< BaseFunction > elementPolynomialBasePtr
Pointer to entity polynomial base.
MoFEMErrorCode getEntityDataOrder(const EntityType type, const FieldSpace space, boost::ptr_vector< EntitiesFieldData::EntData > &data) const
Get the entity data order.
ForcesAndSourcesCore * refinePtrFE
Element to integrate parent or child.
MoFEMErrorCode getNoFieldColIndices(EntitiesFieldData &data, const int bit_number) const
get col NoField indices
int getMaxDataOrder() const
Get max order of approximation for data fields.
boost::shared_ptr< BaseFunction > userPolynomialBasePtr
Pointer to user polynomial base.
MoFEMErrorCode getProblemTypeRowIndices(const std::string &field_name, EntityType type, int side_number, VectorInt &indices) const
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
RuleHookFun getRuleHook
Hook to get rule.
MoFEMErrorCode createDataOnElement(EntityType type)
Create a entity data on element object.
virtual MoFEMErrorCode postProcess()
function is run at the end of loop
EntitiesFieldData & dataL2
KSPContext
pass information about context of KSP/DM for with finite element is computed
Operator for broken loop side.
Copy geometry-related data from one element to other.
Iterate over range of (sub)elements.
OpLoopRange(MoFEM::Interface &m_field, const std::string fe_name, boost::shared_ptr< Range > fe_range, const LogManager::SeverityLevel sev=Sev::noisy)
boost::shared_ptr< E > & getRangeFEPtr()
boost::ptr_deque< UserDataOperator > & getOpPtrVector()
const std::string rangeFEName
const LogManager::SeverityLevel sevLevel
MoFEMErrorCode doWork(int side, EntityType type, EntitiesFieldData::EntData &data)
Operator for linear form, usually to calculate values on right hand side.
boost::shared_ptr< Range > feRange
boost::shared_ptr< E > rangeFEPtr
Element used to execute operators on side of the element.
MoFEMErrorCode doWork(int side, EntityType type, EntitiesFieldData::EntData &data)
Operator for linear form, usually to calculate values on right hand side.
boost::shared_ptr< Range > feRange
boost::shared_ptr< E > & getSideFEPtr()
const std::string sideFEName
boost::shared_ptr< AdjCache > adjCache
boost::shared_ptr< E > sideFEPtr
OpLoopSide(MoFEM::Interface &m_field, const std::string fe_name, const int side_dim, boost::shared_ptr< Range > fe_range, const LogManager::SeverityLevel sev=Sev::noisy, boost::shared_ptr< AdjCache > adj_cache=nullptr)
boost::ptr_deque< UserDataOperator > & getOpPtrVector()
const LogManager::SeverityLevel sevLevel
OpLoopSide(MoFEM::Interface &m_field, const std::string fe_name, const int side_dim, const LogManager::SeverityLevel sev=Sev::noisy, boost::shared_ptr< AdjCache > adj_cache=nullptr)
Construct a new Op Loop Side object.
static constexpr Switches CtxSetA
static constexpr Switches CtxSetX_TT
std::bitset< 8 > Switches
static constexpr Switches CtxSetX_T
static constexpr Switches CtxSetB
Mat & snes_B
preconditioner of jacobian matrix
Mat & snes_A
jacobian matrix
Vec & ts_F
residual vector
Vec & ts_u_tt
second time derivative of state vector
Vec & ts_u_t
time derivative of state vector
Mat & ts_B
Preconditioner for ts_A.
Base volume element used to integrate on skeleton.