1130 {
1132
1133 const double vol = getMeasure();
1134 auto t_h = getFTensor0FromVec(*
hPtr);
1135
1136 auto t_row_diff_base = row_data.getFTensor1DiffN<
SPACE_DIM>();
1137 auto t_w = getFTensor0IntegrationWeight();
1138 auto t_coords = getFTensor1CoordsAtGaussPts();
1139
1140 for (int gg = 0; gg != nbIntegrationPts; gg++) {
1141
1142 const double r = t_coords(0);
1144
1145 double set_h;
1147 set_h =
init_h(t_coords(0), t_coords(1), t_coords(2));
1148 else
1149 set_h = t_h;
1150
1152
1153 int rr = 0;
1154 for (; rr != nbRows; ++rr) {
1155 auto t_col_diff_base = col_data.getFTensor1DiffN<
SPACE_DIM>(gg, 0);
1156
1157 for (int cc = 0; cc != nbCols; ++cc) {
1158 locMat(rr, cc) += (t_row_diff_base(
i) * t_col_diff_base(
i)) *
m;
1159
1160 ++t_col_diff_base;
1161 }
1162
1163 ++t_row_diff_base;
1164 }
1165
1166 for (; rr < nbRowBaseFunctions; ++rr) {
1167 ++t_row_diff_base;
1168 }
1169
1170 ++t_h;
1171 ++t_w;
1172 ++t_coords;
1173 }
1174
1176 }
#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()
FTensor::Index< 'm', SPACE_DIM > m
auto init_h
Initialisation function.
FTensor::Index< 'i', SPACE_DIM > i
constexpr IntegrationType I