243 {
244
245 ParallelComm *pcomm =
247
248 MOFEM_LOG(
"EP", Sev::noisy) <<
"get_two_sides_of_crack_surface";
249
250 if (!pcomm->rank()) {
251
252 auto impl = [&](auto &saids) {
254
256
257 auto get_adj = [&](auto e, auto dim) {
260 e, dim, true, adj, moab::Interface::UNION),
261 "get adj");
262 return adj;
263 };
264
265 auto get_conn = [&](auto e) {
268 "get connectivity");
269 return conn;
270 };
271
272 constexpr bool debug =
false;
275 body_ents);
276 auto body_skin =
get_skin(m_field, body_ents);
277 auto body_skin_edges = get_adj(body_skin, 1);
278
279 auto crack_skin =
280 subtract(
get_skin(m_field, crack_faces), body_skin_edges);
281 auto crack_skin_conn = get_conn(crack_skin);
282 auto crack_skin_conn_edges = get_adj(crack_skin_conn, 1);
283 auto crack_edges = get_adj(crack_faces, 1);
284 crack_edges = subtract(crack_edges, crack_skin);
285 auto all_tets = get_adj(crack_edges, 3);
286 crack_edges = subtract(crack_edges, crack_skin_conn_edges);
287 auto crack_conn = get_conn(crack_edges);
288 all_tets.merge(get_adj(crack_conn, 3));
289
294 crack_edges);
295 }
296
297 if (crack_faces.size()) {
298 auto grow = [&](
auto r) {
299 auto crack_faces_conn = get_conn(crack_faces);
301 auto size_r = 0;
302 while (size_r !=
r.size() &&
r.size() > 0) {
304 CHKERR moab.get_connectivity(r,
v,
true);
305 v = subtract(
v, crack_faces_conn);
308 moab::Interface::UNION);
309 r = intersect(r, all_tets);
310 }
312 break;
313 }
314 }
316 };
317
318 Range all_tets_ord = all_tets;
319 while (all_tets.size()) {
320 Range faces = get_adj(unite(saids.first, saids.second), 2);
321 faces = subtract(crack_faces, faces);
322 if (faces.size()) {
324 auto fit = faces.begin();
325 for (; fit != faces.end(); ++fit) {
326 tets = intersect(get_adj(
Range(*fit, *fit), 3), all_tets);
327 if (tets.size() == 2) {
328 break;
329 }
330 }
331 if (tets.empty()) {
332 break;
333 } else {
334 saids.first.insert(tets[0]);
335 saids.first = grow(saids.first);
336 all_tets = subtract(all_tets, saids.first);
337 if (tets.size() == 2) {
338 saids.second.insert(tets[1]);
339 saids.second = grow(saids.second);
340 all_tets = subtract(all_tets, saids.second);
341 }
342 }
343 } else {
344 break;
345 }
346 }
347
348 saids.first = subtract(all_tets_ord, saids.second);
349 saids.second = subtract(all_tets_ord, saids.first);
350 }
351
353 };
354
355 std::pair<Range, Range> saids;
356 if (crack_faces.size())
358 return saids;
359 }
360
361 MOFEM_LOG(
"EP", Sev::noisy) <<
"get_two_sides_of_crack_surface <- done";
362
363 return std::pair<Range, Range>();
364}
#define CHK_THROW_MESSAGE(err, msg)
Check and throw MoFEM exception.
#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 MOFEM_LOG(channel, severity)
Log.
const double v
phase velocity of light in medium (cm/ns)
