3508 {
3510
3511#ifndef NDEBUG
3512 auto fe_mi_ptr = getFEMethod()->numeredEntFiniteElementPtr;
3513 auto pstatus = fe_mi_ptr->getPStatus();
3514 if (pstatus) {
3515 auto owner = fe_mi_ptr->getOwnerProc();
3517 << "OpFaceMaterialForce: owner proc is not 0, owner proc: " << owner
3518 << " " << getPtrFE()->mField.get_comm_rank();
3519 }
3520#endif
3521
3523
3526 double face_pressure = 0.;
3527 auto t_T = getFTensor1FromMat<SPACE_DIM>(
3529 auto t_p =
3531 auto t_w = getFTensor0IntegrationWeight();
3532 for (auto gg = 0; gg != getGaussPts().size2(); ++gg) {
3533 t_face_T(
I) += t_w * t_T(
I);
3534 face_pressure += t_w * t_p;
3535 ++t_w;
3536 ++t_T;
3537 ++t_p;
3538 }
3539 t_face_T(
I) *= getMeasure();
3540 face_pressure *= getMeasure();
3541
3542 auto get_tag = [&](auto name, auto dim) {
3543 auto &moab = getPtrFE()->mField.get_moab();
3545 double def_val[] = {0., 0., 0.};
3546 CHK_MOAB_THROW(moab.tag_get_handle(name, dim, MB_TYPE_DOUBLE, tag,
3547 MB_TAG_CREAT | MB_TAG_SPARSE, def_val),
3548 "create tag");
3549 return tag;
3550 };
3551
3552 auto set_tag = [&](auto &&tag, auto ptr) {
3553 auto &moab = getPtrFE()->mField.get_moab();
3554 auto face = getPtrFE()->getFEEntityHandle();
3555 CHK_MOAB_THROW(moab.tag_set_data(tag, &face, 1, ptr),
"set tag");
3556 };
3557
3558 set_tag(get_tag("MaterialForce", 3), &t_face_T(0));
3559 set_tag(get_tag("FacePressure", 1), &face_pressure);
3560
3562}
#define FTENSOR_INDEX(DIM, I)
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
#define CHK_MOAB_THROW(err, msg)
Check error code of MoAB function and throw MoFEM exception.
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define MOFEM_LOG(channel, severity)
Log.
static auto getFTensor0FromVec(ublas::vector< T, A > &data)
Get tensor rank 0 (scalar) form data vector.
constexpr IntegrationType I