22 {
24
25 auto cubit_meshset_ptr =
26 m_field.
getInterface<MeshsetsManager>()->getCubitMeshsetPtr(ms_id,
28 std::vector<double> block_data;
29 CHKERR cubit_meshset_ptr->getAttributes(block_data);
30 if (block_data.size() != 2) {
32 "Expected that block has two attribute");
33 }
34 normal_stiffness = block_data[0];
35 tangent_stiffness = block_data[1];
36
39 << "Normal stiffness " << normal_stiffness;
41 << "Tangent stiffness " << tangent_stiffness;
42
43 ents = boost::make_shared<Range>();
45 m_field.
get_moab().get_entities_by_handle(cubit_meshset_ptr->meshset,
46 *(ents), true);
47
51
53 }
#define MOFEM_LOG_SEVERITY_SYNC(comm, severity)
Synchronise "SYNC" on curtain severity level.
#define MOFEM_TAG_AND_LOG(channel, severity, tag)
Tag and log in channel.
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
@ MOFEM_DATA_INCONSISTENCY
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define CHKERR
Inline error check.
#define MOFEM_LOG_CHANNEL(channel)
Set and reset channel.
virtual moab::Interface & get_moab()=0
virtual MPI_Comm & get_comm() const =0
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface refernce to pointer of interface.