114 {
117 double *vec_ptr = OP::locF.data().data();
119 int *ind_ptr = data.
getIndices().data().data();
121 }
124 std::vector<EntityHandle> ents(field_ents.size());
125 std::transform(field_ents.begin(), field_ents.end(), ents.begin(),
126 [](const auto *fe) { return fe->getEnt(); });
127 if (field_ents.empty())
131 auto &moab = getMoab();
134 topo_values.data().data());
135 topo_values += OP::locF;
137 OP::locF.data().data());
138 }
140 }
#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.
