980 {
982
983#ifndef NDEBUG
986 "Topological data pointer is null");
987#endif
988
990 if (!nb_dofs)
992
993 const int nb_integration_pts = data.
getN().size1();
994
995 auto t_w = getFTensor0IntegrationWeight();
996 const int nb_base_functions = data.
getN().size2();
998 auto t_obj_du_gamma =
999 topoData->getFTensorObjDDisplacement(nb_integration_pts);
1000
1002
1003 auto get_ftensor1 = [](
auto &
v) {
1005 &
v[0], &
v[1], &
v[2]);
1006 };
1007
1008 for (int gg = 0; gg != nb_integration_pts; ++gg) {
1009 const double a = t_w * getMeasure();
1010 auto t_nf = get_ftensor1(
nF);
1011
1012 int bb = 0;
1013 for (; bb != nb_dofs /
SPACE_DIM; ++bb) {
1014 t_nf(
i) +=
a * t_row_base_fun * t_obj_du_gamma(
i);
1015 ++t_nf;
1016 ++t_row_base_fun;
1017 }
1018 for (; bb != nb_base_functions; ++bb)
1019 ++t_row_base_fun;
1020
1021 ++t_w;
1022 ++t_obj_du_gamma;
1023 }
1024
1026}
#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