1368 {
1370
1371 auto get_body_id = [&](auto fe_ent) {
1373 if (
m.second.find(fe_ent) !=
m.second.end()) {
1375 }
1376 }
1377 return -1;
1378 };
1379
1382
1383 auto fe_ent = getNumeredEntFiniteElementPtr()->getEnt();
1385 auto body_id = get_body_id(fe_ent);
1387
1389 post_proc_ents, &fe_id);
1391 post_proc_ents, &body_id);
1392
1393 auto nb_gauss_pts = getGaussPts().size2();
1394 auto t_u_h1 = getFTensor1FromMat<3>(*
uH1Ptr);
1395 auto t_u_l2 = getFTensor1FromMat<3>(
commonDataPtr->contactDisp);
1396 auto t_coords = getFTensor1CoordsAtGaussPts();
1397
1399 auto t_x_h1 = getFTensor1FromPtr<3>(&x_h1(0, 0));
1401 auto t_x_l2 = getFTensor1FromPtr<3>(&x_l2(0, 0));
1404
1406
1407
1408 for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
1409 t_x_h1(
i) = t_coords(
i) + t_u_h1(
i);
1410 t_x_l2(
i) = t_coords(
i) + t_u_l2(
i);
1411
1412 ++t_coords;
1413 ++t_u_h1;
1414 ++t_u_l2;
1415 ++t_x_h1;
1416 ++t_x_l2;
1417 }
1418
1419 CHKERR moab_post_proc_mesh.set_coords(
1420 &*map_gauss_pts.begin(), map_gauss_pts.size(), &*x_h1.data().begin());
1421 CHKERR moab_post_proc_mesh.tag_set_data(
1422 contactTreePtr->thSmallX, &*map_gauss_pts.begin(), map_gauss_pts.size(),
1423 &*x_h1.data().begin());
1424 CHKERR moab_post_proc_mesh.tag_set_data(
1425 contactTreePtr->thLargeX, &*map_gauss_pts.begin(), map_gauss_pts.size(),
1426 &*coords.data().begin());
1427 CHKERR moab_post_proc_mesh.tag_set_data(
1428 contactTreePtr->thTraction, &*map_gauss_pts.begin(), map_gauss_pts.size(),
1429 &*tractions.data().begin());
1430
1432}
#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()
#define CHKERR
Inline error check.
FTensor::Index< 'i', SPACE_DIM > i
auto id_from_handle(const EntityHandle h)
FTensor::Index< 'm', 3 > m
