45 {
48 double *vec_ptr = OP::nF.data().data();
50 int *ind_ptr = data.
getIndices().data().data();
52 }
55 }
58 std::vector<EntityHandle> ents(field_ents.size());
59 std::transform(field_ents.begin(), field_ents.end(), ents.begin(),
60 [](const auto *fe) { return fe->getEnt(); });
61 if (field_ents.empty())
65 auto &moab = OP::getMoab();
68 topo_values.data().data());
69 noalias(topo_values) += OP::nF;
71 OP::nF.data().data());
72 }
74 }
#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 MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define CHKERR
Inline error check.
UBlasVector< double > VectorDouble
auto type_from_handle(const EntityHandle h)
get type from entity handle
MoFEMErrorCode VecSetValues(Vec V, const EntitiesFieldData::EntData &data, const double *ptr, InsertMode iora)
Assemble PETSc vector.
const VectorFieldEntities & getFieldEntities() const
Get field entities (const version)
const VectorInt & getIndices() const
Get global indices of degrees of freedom on entity.
