101 {
104 double *vec_ptr = OP::locF.data().data();
106 int *ind_ptr = data.
getIndices().data().data();
108 }
111 std::vector<EntityHandle> ents(field_ents.size());
112 std::transform(field_ents.begin(), field_ents.end(), ents.begin(),
113 [](const auto *fe) { return fe->getEnt(); });
114 if (field_ents.empty())
118 auto &moab = getMoab();
121 topo_values.data().data());
122 topo_values += OP::locF;
124 OP::locF.data().data());
125 }
127 }
#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.
