3821 {
3823
3824 if (OP::entsPtr) {
3825 if (OP::entsPtr->find(this->getFEEntityHandle()) == OP::entsPtr->end())
3827 }
3828
3829#ifndef NDEBUG
3830 if (!brokenBaseSideData) {
3832 }
3833#endif
3834
3835 auto do_work_rhs = [this](int row_side, EntityType row_type,
3838
3840 if (!OP::nbRows)
3842
3843 OP::nbIntegrationPts = OP::getGaussPts().size2();
3844
3845 OP::nbRowBaseFunctions = OP::getNbOfBaseFunctions(row_data);
3846
3847 OP::locF.resize(OP::nbRows, false);
3848 OP::locF.clear();
3849
3851
3852 CHKERR this->aSsemble(row_data);
3854 };
3855
3856 switch (OP::opType) {
3857 case OP::OPSPACE:
3858 for (auto &bd : *brokenBaseSideData) {
3859 this->sourceVec =
3860 boost::shared_ptr<MatrixDouble>(brokenBaseSideData, &bd.getFlux());
3861 CHKERR do_work_rhs(bd.getSide(), bd.getType(), bd.getData());
3862 this->sourceVec.reset();
3863 }
3864 break;
3865 default:
3867 (std::string("wrong op type ") +
3868 OpBaseDerivativesBase::OpTypeNames[OP::opType])
3869 .c_str());
3870 }
3871
3873}
#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 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 CHKERR
Inline error check.
MoFEMErrorCode iNtegrate(EntitiesFieldData::EntData &row_data)
Data on single entity (This is passed as argument to DataOperator::doWork)
const VectorInt & getIndices() const
Get global indices of degrees of freedom on entity.