995 auto get_body_id = [&](
auto fe_ent) {
997 if (
m.second.find(fe_ent) !=
m.second.end()) {
1007 auto fe_ent = getNumeredEntFiniteElementPtr()->getEnt();
1009 auto body_id = get_body_id(fe_ent);
1013 post_proc_ents, &fe_id);
1015 post_proc_ents, &body_id);
1017 auto nb_gauss_pts = getGaussPts().size2();
1018 auto t_u_h1 = getFTensor1FromMat<3>(*
uH1Ptr);
1019 auto t_u_l2 = getFTensor1FromMat<3>(
commonDataPtr->contactDisp);
1020 auto t_coords = getFTensor1CoordsAtGaussPts();
1023 auto t_x_h1 = getFTensor1FromPtr<3>(&x_h1(0, 0));
1025 auto t_x_l2 = getFTensor1FromPtr<3>(&x_l2(0, 0));
1032 for (
auto gg = 0; gg != nb_gauss_pts; ++gg) {
1033 t_x_h1(
i) = t_coords(
i) + t_u_h1(
i);
1034 t_x_l2(
i) = t_coords(
i) + t_u_l2(
i);
1043 CHKERR moab_post_proc_mesh.set_coords(
1044 &*map_gauss_pts.begin(), map_gauss_pts.size(), &*x_h1.data().begin());
1045 CHKERR moab_post_proc_mesh.tag_set_data(
1046 contactTreePtr->thSmallX, &*map_gauss_pts.begin(), map_gauss_pts.size(),
1047 &*x_h1.data().begin());
1048 CHKERR moab_post_proc_mesh.tag_set_data(
1049 contactTreePtr->thLargeX, &*map_gauss_pts.begin(), map_gauss_pts.size(),
1050 &*coords.data().begin());
1051 CHKERR moab_post_proc_mesh.tag_set_data(
1052 contactTreePtr->thTraction, &*map_gauss_pts.begin(), map_gauss_pts.size(),
1053 &*tractions.data().begin());