4674 {
4676
4679
4683 auto &m_field = this->getPtrFE()->mField;
4684 CHKERR m_field.get_moab().get_adjacencies(&prism_ent, 1, 2,
false, tri_ents,
4685 moab::Interface::UNION);
4686 tri_ents = tri_ents.subset_by_type(MBTRI);
4689 }
4690
4691 const int nb_gauss_pts = data.
getN().size1();
4692
4694 const double area_s =
4696
4699
4700 auto t_lagrange_slave =
4703 auto t_w = getFTensor0IntegrationWeightSlave();
4704
4705 for (int gg = 0; gg != nb_gauss_pts; gg++) {
4706 const double val_s = t_w * area_s;
4708 bool gap_below_tolerance =
4711 if (gap_below_tolerance) {
4713 } else {
4717 }
4721 }
4722 }
4723
4724 ++t_lagrange_slave;
4725 ++t_gap_gp;
4726 ++t_w;
4727 }
4728
4729 constexpr std::array<int, 2> indices = {
4732 };
4734 indices.data(), &
vecR[0], ADD_VALUES);
4735
4737}
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define CHKERR
Inline error check.
static auto getFTensor0FromVec(ublas::vector< T, A > &data)
Get tensor rank 0 (scalar) form data vector.
MoFEMErrorCode VecSetValues(Vec V, const EntitiesFieldData::EntData &data, const double *ptr, InsertMode iora)
Assemble PETSc vector.
MatrixDouble & getN(const FieldApproximationBase base)
get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb....
const VectorDouble & getFieldData() const
get dofs values
const VectorInt & getIndices() const
Get global indices of dofs on entity.