1375 {
1377
1378 auto get_body_id = [&](auto fe_ent) {
1380 if (
m.second.find(fe_ent) !=
m.second.end()) {
1382 }
1383 }
1384 return -1;
1385 };
1386
1389
1390 auto fe_ent = getNumeredEntFiniteElementPtr()->getEnt();
1392 auto body_id = get_body_id(fe_ent);
1394
1396 post_proc_ents, &fe_id);
1398 post_proc_ents, &body_id);
1399
1400 auto nb_gauss_pts = getGaussPts().size2();
1401 auto t_u_h1 = getFTensor1FromMat<3>(*
uH1Ptr);
1402 auto t_u_l2 = getFTensor1FromMat<3>(
commonDataPtr->contactDisp);
1403 auto t_coords = getFTensor1CoordsAtGaussPts();
1404
1406 auto t_x_h1 = getFTensor1FromPtr<3>(&x_h1(0, 0));
1408 auto t_x_l2 = getFTensor1FromPtr<3>(&x_l2(0, 0));
1411
1413
1414
1415 for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
1416 t_x_h1(
i) = t_coords(
i) + t_u_h1(
i);
1417 t_x_l2(
i) = t_coords(
i) + t_u_l2(
i);
1418
1419 ++t_coords;
1420 ++t_u_h1;
1421 ++t_u_l2;
1422 ++t_x_h1;
1423 ++t_x_l2;
1424 }
1425
1426 CHKERR moab_post_proc_mesh.set_coords(
1427 &*map_gauss_pts.begin(), map_gauss_pts.size(), &*x_h1.data().begin());
1428 CHKERR moab_post_proc_mesh.tag_set_data(
1429 contactTreePtr->thSmallX, &*map_gauss_pts.begin(), map_gauss_pts.size(),
1430 &*x_h1.data().begin());
1431 CHKERR moab_post_proc_mesh.tag_set_data(
1432 contactTreePtr->thLargeX, &*map_gauss_pts.begin(), map_gauss_pts.size(),
1433 &*coords.data().begin());
1434 CHKERR moab_post_proc_mesh.tag_set_data(
1435 contactTreePtr->thTraction, &*map_gauss_pts.begin(), map_gauss_pts.size(),
1436 &*tractions.data().begin());
1437
1439}
#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
