68 {
70
80 mesh_nodals_positions);
81 }
82
83
85
87
88 if (
bit->getName().compare(0, 14,
"FLUID_PRESSURE") == 0) {
89
90
91 std::vector<double> attributes;
93 if (attributes.size() < 7) {
94 SETERRQ1(PETSC_COMM_SELF, 1,
95 "not enough block attributes to deffine fluid pressure "
96 "element, attributes.size() = %d ",
97 attributes.size());
98 }
108
111 true);
115
116 Range tets_skin_tris;
117 CHKERR skin.find_skin(0, tets,
false, tets_skin_tris);
119 std::ostringstream ss;
122
124 setOfFluids[
bit->getMeshsetId()].tRis, MBTRI,
"FLUID_PRESSURE_FE");
125 }
126 }
127
129}
#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.
virtual MoFEMErrorCode add_finite_element(const std::string &fe_name, enum MoFEMTypes bh=MF_EXCL, int verb=DEFAULT_VERBOSITY)=0
add finite element
virtual MoFEMErrorCode modify_finite_element_add_field_col(const std::string &fe_name, const std::string name_row)=0
set field col which finite element use
virtual MoFEMErrorCode add_ents_to_finite_element_by_type(const EntityHandle entities, const EntityType type, const std::string &name, const bool recursive=true)=0
add entities to finite element
virtual MoFEMErrorCode modify_finite_element_add_field_data(const std::string &fe_name, const std::string name_filed)=0
set finite element field data
virtual MoFEMErrorCode modify_finite_element_add_field_row(const std::string &fe_name, const std::string name_row)=0
set field row which finite element use
virtual bool check_field(const std::string &name) const =0
check if field is in database
#define _IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(MESHSET_MANAGER, CUBITBCTYPE, IT)
Iterator that loops over a specific Cubit MeshSet having a particular BC meshset in a moFEM field.
constexpr auto field_name
std::map< MeshSetId, FluidData > setOfFluids
virtual moab::Interface & get_moab()=0
virtual MPI_Comm & get_comm() const =0