673 {
675
676 auto cubit_meshset_ptr =
677 m_field.
getInterface<MeshsetsManager>()->getCubitMeshsetPtr(ms_id,
679
680 std::vector<double> block_data;
681 CHKERR cubit_meshset_ptr->getAttributes(block_data);
684 <<
"BLOCKSET is expected to have " <<
FIELD_DIM
685 << " attributes but has size " << block_data.size();
687
689 "Size of attribute in BLOCKSET is too small");
690 }
691 }
692
693 for (
unsigned int ii = 0; ii !=
FIELD_DIM; ++ii) {
694 this->
tForce(ii) = block_data[ii];
695 }
696
698 << "Flux blockset " << cubit_meshset_ptr->getName();
700 << "Number of attributes " << block_data.size();
702 <<
"tForce vector initialised: " << this->
tForce;
703
704 this->entsPtr = boost::make_shared<Range>();
705 CHKERR m_field.
get_moab().get_entities_by_handle(cubit_meshset_ptr->meshset,
706 *(this->entsPtr), true);
707
709}
#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, severity)
Log.
virtual moab::Interface & get_moab()=0
FTensor::Tensor1< double, FIELD_DIM > tForce
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface refernce to pointer of interface.