240 {
242
243
244 if (type == MBVERTEX) {
245 int nb_gauss_pts = data.getN().size1();
247 for (int gg = 0; gg != nb_gauss_pts; gg++) {
249 gg);
251 inner_prod(grad_at_gauss_pt, grad_at_gauss_pt);
252 }
254 for (int gg = 0; gg != nb_gauss_pts; gg++) {
256 }
258
260 for (int gg = 0; gg != nb_gauss_pts; gg++) {
262 gg);
263 for (int rr = 0; rr != 2; rr++) {
266 }
267 }
268 } else {
269
271 for (int gg = 0; gg != nb_gauss_pts; gg++) {
273 gg);
274 for (int rr = 0; rr != 2; rr++) {
277 }
278 }
279 }
280 }
281
282
283
284
285
286
287
288
290 }
#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()
MatrixDouble gradU
nb_gauss_pts x 2 gradients at integration pts
VectorDouble normGradU2
size of nb_gauss_pts, norm of gradient
VectorDouble aN
size of nb_gauss_pts,
MatrixDouble aNpow3byGradU
nb_gauss_pts x 2,
MatrixDouble aNbyGradU
nb_gauss_pts x 2,