1189 {
1191
1192 const double vol = getMeasure();
1193 auto t_h = getFTensor0FromVec(*
hPtr);
1194 auto t_grad_h = getFTensor1FromMat<SPACE_DIM>(*
gradHPtr);
1195 auto t_g = getFTensor0FromVec(*
gPtr);
1196 auto t_coords = getFTensor1CoordsAtGaussPts();
1197
1198 auto t_base = data.getFTensor0N();
1199 auto t_diff_base = data.getFTensor1DiffN<
SPACE_DIM>();
1200 auto t_w = getFTensor0IntegrationWeight();
1201
1202 for (int gg = 0; gg != nbIntegrationPts; ++gg) {
1203
1204 const double r = t_coords(0);
1206
1207 double set_h;
1209 set_h =
init_h(t_coords(0), t_coords(1), t_coords(2));
1210 else
1211 set_h = t_h;
1212
1213 const double f =
get_f(set_h);
1214
1215 int bb = 0;
1216 for (; bb != nbRows; ++bb) {
1217 locF[bb] += (t_base * alpha) * (t_g - f);
1218 locF[bb] -= (t_diff_base(
i) * (
eta2 * alpha)) * t_grad_h(
i);
1219 ++t_base;
1220 ++t_diff_base;
1221 }
1222
1223 for (; bb < nbRowBaseFunctions; ++bb) {
1224 ++t_base;
1225 ++t_diff_base;
1226 }
1227
1228 ++t_h;
1229 ++t_grad_h;
1230 ++t_g;
1231
1232 ++t_coords;
1233 ++t_w;
1234 }
1235
1237 }
#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()
auto init_h
Initialisation function.
FTensor::Index< 'i', SPACE_DIM > i
constexpr IntegrationType I