976 {
978
979#ifndef NDEBUG
982 "Topological data pointer is null");
983#endif
984
986 if (!nb_dofs)
988
989 const int nb_integration_pts = data.
getN().size1();
990
991 auto t_w = getFTensor0IntegrationWeight();
992 const int nb_base_functions = data.
getN().size2();
994 auto t_obj_du_gamma =
995 topoData->getFTensorObjDDisplacement(nb_integration_pts);
996
998
999 auto get_ftensor1 = [](
auto &
v) {
1001 &
v[0], &
v[1], &
v[2]);
1002 };
1003
1004 for (int gg = 0; gg != nb_integration_pts; ++gg) {
1005 const double a = t_w * getMeasure();
1006 auto t_nf = get_ftensor1(
nF);
1007
1008 int bb = 0;
1009 for (; bb != nb_dofs /
SPACE_DIM; ++bb) {
1010 t_nf(
i) +=
a * t_row_base_fun * t_obj_du_gamma(
i);
1011 ++t_nf;
1012 ++t_row_base_fun;
1013 }
1014 for (; bb != nb_base_functions; ++bb)
1015 ++t_row_base_fun;
1016
1017 ++t_w;
1018 ++t_obj_du_gamma;
1019 }
1020
1022}
#define FTENSOR_INDEX(DIM, I)
#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 ...
@ MOFEM_DATA_INCONSISTENCY
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
FTensor::Index< 'i', SPACE_DIM > i
const double v
phase velocity of light in medium (cm/ns)
boost::shared_ptr< TopologicalData > topoData
FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > getFTensor0N(const FieldApproximationBase base)
Get base function as Tensor0.
MatrixDouble & getN(const FieldApproximationBase base)
get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb....
const VectorInt & getIndices() const
Get global indices of degrees of freedom on entity.
VectorDouble nF
local right hand side vector