32 {
35 double *vec_ptr = OP::nF.data().data();
37 int *ind_ptr = data.
getIndices().data().data();
39 }
42 }
45 std::vector<EntityHandle> ents(field_ents.size());
46 std::transform(field_ents.begin(), field_ents.end(), ents.begin(),
47 [](const auto *fe) { return fe->getEnt(); });
48 if (field_ents.empty())
52 auto &moab = OP::getMoab();
55 topo_values.data().data());
56 noalias(topo_values) += OP::nF;
58 OP::nF.data().data());
59 }
61 }
#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.
