4660 {
4662
4663 if (OP::entsPtr) {
4664 if (OP::entsPtr->find(this->getFEEntityHandle()) == OP::entsPtr->end())
4666 }
4667
4668#ifndef NDEBUG
4669 if (!brokenBaseSideData) {
4671 }
4672#endif
4673
4674 auto do_work_rhs = [this](int row_side, EntityType row_type,
4677
4679 if (!OP::nbRows)
4681
4682 OP::nbIntegrationPts = OP::getGaussPts().size2();
4683
4684 OP::nbRowBaseFunctions = OP::getNbOfBaseFunctions(row_data);
4685
4686 OP::locF.resize(OP::nbRows, false);
4687 OP::locF.clear();
4688
4689 CHKERR this->iNtegrate(row_data);
4690
4691 CHKERR this->aSsemble(row_data);
4693 };
4694
4695 switch (OP::opType) {
4696 case OP::OPSPACE:
4697 for (auto &bd : *brokenBaseSideData) {
4698 this->sourceVec =
4699 boost::shared_ptr<MatrixDouble>(brokenBaseSideData, &bd.getFlux());
4700 CHKERR do_work_rhs(bd.getSide(), bd.getType(), bd.getData());
4701 this->sourceVec.reset();
4702 }
4703 break;
4704 default:
4706 (std::string("wrong op type ") +
4707 OpBaseDerivativesBase::OpTypeNames[OP::opType])
4708 .c_str());
4709 }
4710
4712}
#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.
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.