332 auto refined_ents_ptr = m_field.get_ref_ents();
333 auto refined_finite_elements_ptr = m_field.get_ref_finite_elements();
334 ReadUtilIface *read_util;
337 CHKERR m_field.get_moab().query_interface(read_util);
339 Range ents_to_set_bit;
344 auto get_parent_ents_view = [&](
const auto type_child,
345 const auto type_parent) {
347 refined_ents_ptr->get<Composite_EntType_and_ParentEntType_mi_tag>();
349 ref_ents.equal_range(boost::make_tuple(type_child, type_parent));
353 using I = decltype(range.first);
355 boost::function<
bool(
I it)> tester = [&](
I it) {
356 return ((*it)->getBitRefLevel() &
bit).any();
360 ref_parent_ents_view.insert(
f, s);
366 return ref_parent_ents_view;
369 auto ref_parent_ents_view = get_parent_ents_view(MBVERTEX, MBEDGE);
371 auto tets = _tets.subset_by_type(MBTET);
373 auto get_ele_parent_view = [&]() {
374 auto &ref_ents = refined_finite_elements_ptr->get<Ent_mi_tag>();
376 for (
auto p = tets.pair_begin(); p != tets.pair_end(); ++p) {
377 auto tmi = ref_ents.lower_bound(p->first);
378 auto tme = ref_ents.upper_bound(p->second);
379 ref_ele_parent_view.insert(tmi, tme);
381 return ref_ele_parent_view;
384 auto ref_ele_parent_view = get_ele_parent_view();
386 std::vector<EntityHandle> parent_tets_refinded;
387 std::vector<EntityHandle> vertices_of_created_tets;
388 std::vector<BitRefEdges> parent_edges_bit_vec;
389 std::vector<int> nb_new_tets_vec;
390 std::vector<int> sub_type_vec;
392 parent_tets_refinded.reserve(tets.size());
393 vertices_of_created_tets.reserve(4 * tets.size());
394 parent_edges_bit_vec.reserve(tets.size());
395 nb_new_tets_vec.reserve(tets.size());
396 sub_type_vec.reserve(tets.size());
399 for (
auto tit : tets) {
404 CHKERR moab.get_connectivity(tit, conn, num_nodes,
true);
409 int split_edges[] = {-1, -1, -1, -1, -1, -1};
411 CHKERR set_edge_bits(moab, ref_parent_ents_view, tit, parent_edges_bit,
412 edge_new_nodes, split_edges);
414 if (parent_edges_bit.count()) {
418 std::copy(&conn[0], &conn[4], &_conn_[0]);
422 std::fill(&new_tets_conns[0], &new_tets_conns[8 * 4],
no_handle);
424 int sub_type = -1, nb_new_tets = 0;
425 switch (parent_edges_bit.count()) {
428 tet_type_1(_conn_, split_edges[0], edge_new_nodes[split_edges[0]],
434 tet_type_2(_conn_, split_edges, edge_new_nodes, new_tets_conns);
435 if (sub_type & (4 | 8 | 16)) {
438 }
else if (sub_type == 1) {
446 tet_type_3(_conn_, split_edges, edge_new_nodes, new_tets_conns);
450 }
else if (sub_type <= 7) {
457 tet_type_4(_conn_, split_edges, edge_new_nodes, new_tets_conns);
461 }
else if (sub_type <= 7) {
467 sub_type =
tet_type_5(moab, _conn_, edge_new_nodes, new_tets_conns);
472 tet_type_6(moab, _conn_, edge_new_nodes, new_tets_conns);
480 auto &ref_ele_by_parent_and_ref_edges =
482 .get<Composite_ParentEnt_And_BitsOfRefinedEdges_mi_tag>();
483 auto it_by_ref_edges = ref_ele_by_parent_and_ref_edges.equal_range(
484 boost::make_tuple(tit, parent_edges_bit.to_ulong()));
486 std::vector<EntityHandle> ents_to_set_bit_vec;
487 if (std::distance(it_by_ref_edges.first, it_by_ref_edges.second) ==
488 static_cast<size_t>(nb_new_tets)) {
489 ents_to_set_bit_vec.reserve(nb_new_tets);
490 for (
int tt = 0; it_by_ref_edges.first != it_by_ref_edges.second;
491 it_by_ref_edges.first++, tt++) {
492 auto tet = (*it_by_ref_edges.first)->getEnt();
493 ents_to_set_bit_vec.emplace_back(tet);
497 auto ref_tet_it = refined_ents_ptr->find(tet);
498 if (ref_tet_it == refined_ents_ptr->end()) {
500 "Tetrahedron should be in database");
504 ents_to_set_bit.insert_list(ents_to_set_bit_vec.begin(),
505 ents_to_set_bit_vec.end());
511 parent_tets_refinded.emplace_back(tit);
512 for (
int tt = 0; tt != nb_new_tets; ++tt) {
513 for (
auto nn : {0, 1, 2, 3})
514 vertices_of_created_tets.emplace_back(new_tets_conns[4 * tt + nn]);
516 parent_edges_bit_vec.emplace_back(parent_edges_bit);
517 nb_new_tets_vec.emplace_back(nb_new_tets);
518 sub_type_vec.emplace_back(sub_type);
523 RefEntity_multiIndex::iterator tit_miit;
524 tit_miit = refined_ents_ptr->find(tit);
525 if (tit_miit == refined_ents_ptr->end())
527 "can not find this tet");
530 ents_to_set_bit.insert(tit);
537 const int nb_new_tets = vertices_of_created_tets.size() / 4;
538 read_util->get_element_connect(nb_new_tets, 4, MBTET, 0, start_e, conn_moab);
539 std::copy(vertices_of_created_tets.begin(), vertices_of_created_tets.end(),
541 CHKERR read_util->update_adjacencies(start_e, nb_new_tets, 4, conn_moab);
542 ents_to_set_bit.insert(start_e, start_e + nb_new_tets - 1);
546 for (
auto d : {1, 2}) {
547 CHKERR moab.get_adjacencies(ents_to_set_bit,
d,
true, ents,
548 moab::Interface::UNION);
551 SetParent set_parent;
554 for (
int idx = 0; idx != parent_tets_refinded.size(); ++idx) {
557 const BitRefEdges &parent_edges_bit = parent_edges_bit_vec[idx];
558 const int nb_new_tets = nb_new_tets_vec[idx];
560 std::array<EntityHandle, 8> ref_tets;
561 for (
int tt = 0; tt != nb_new_tets; ++tt, ++start_e)
562 ref_tets[tt] = start_e;
567 nb_new_tets, &parent_edges_bit);
573 map_ref_nodes_by_edges;
576 CHKERR moab.get_adjacencies(&tit, 1, 1,
false, tit_edges);
577 for (
auto edge : tit_edges) {
578 auto eit = ref_parent_ents_view.find(edge);
579 if (eit != ref_parent_ents_view.end()) {
580 map_ref_nodes_by_edges[(*eit)->getParentEnt()] = eit->get()->getEnt();
587 CHKERR moab.get_adjacencies(&tit, 1, 2,
false, tit_faces);
588 Range edges_nodes[6], faces_nodes[4];
592 Range::iterator eit = tit_edges.begin();
593 for (
int ee = 0; eit != tit_edges.end(); eit++, ee++) {
594 CHKERR moab.get_connectivity(&*eit, 1, edges_nodes[ee],
true);
595 auto map_miit = map_ref_nodes_by_edges.find(*eit);
596 if (map_miit != map_ref_nodes_by_edges.end()) {
597 edges_nodes[ee].insert(map_miit->second);
603 Range::iterator fit = tit_faces.begin();
604 for (
int ff = 0; fit != tit_faces.end(); fit++, ff++) {
605 CHKERR moab.get_connectivity(&*fit, 1, faces_nodes[ff],
true);
608 CHKERR moab.get_adjacencies(&*fit, 1, 1,
false, fit_edges);
609 for (Range::iterator eit2 = fit_edges.begin(); eit2 != fit_edges.end();
611 auto map_miit = map_ref_nodes_by_edges.find(*eit2);
612 if (map_miit != map_ref_nodes_by_edges.end()) {
613 faces_nodes[ff].insert(map_miit->second);
620 CHKERR moab.get_connectivity(&tit, 1, tet_nodes,
true);
621 for (std::map<EntityHandle, EntityHandle>::iterator map_miit =
622 map_ref_nodes_by_edges.begin();
623 map_miit != map_ref_nodes_by_edges.end(); map_miit++) {
624 tet_nodes.insert(map_miit->second);
628 CHKERR moab.get_adjacencies(ref_tets.data(), nb_new_tets, 1,
false,
629 ref_edges, moab::Interface::UNION);
631 for (Range::iterator reit = ref_edges.begin(); reit != ref_edges.end();
633 Range ref_edges_nodes;
634 CHKERR moab.get_connectivity(&*reit, 1, ref_edges_nodes,
true);
638 for (; ee < 6; ee++) {
642 if (intersect(edges_nodes[ee], ref_edges_nodes).size() == 2) {
644 CHKERR set_parent(*reit, edge, refined_ents_ptr,
cOre);
652 for (; ff < 4; ff++) {
655 if (intersect(faces_nodes[ff], ref_edges_nodes).size() == 2) {
657 CHKERR set_parent(*reit, face, refined_ents_ptr,
cOre);
666 if (intersect(tet_nodes, ref_edges_nodes).size() == 2) {
667 CHKERR set_parent(*reit, tit, refined_ents_ptr,
cOre);
674 "impossible refined edge");
678 CHKERR moab.get_adjacencies(ref_tets.data(), nb_new_tets, 2,
false,
679 ref_faces, moab::Interface::UNION);
680 Tag th_interface_side;
681 const int def_side[] = {0};
682 CHKERR moab.tag_get_handle(
"INTERFACE_SIDE", 1, MB_TYPE_INTEGER,
684 MB_TAG_CREAT | MB_TAG_SPARSE, def_side);
686 for (
auto rfit : ref_faces) {
687 Range ref_faces_nodes;
688 CHKERR moab.get_connectivity(&rfit, 1, ref_faces_nodes,
true);
691 for (; ff < 4; ff++) {
693 if (intersect(faces_nodes[ff], ref_faces_nodes).size() == 3) {
695 CHKERR set_parent(rfit, face, refined_ents_ptr,
cOre);
698 CHKERR moab.tag_get_data(th_interface_side, &face, 1, &side);
699 CHKERR moab.tag_set_data(th_interface_side, &rfit, 1, &side);
707 if (intersect(tet_nodes, ref_faces_nodes).size() == 3) {
708 CHKERR set_parent(rfit, tit, refined_ents_ptr,
cOre);
712 "impossible refined face");
717 CHKERR set_parent(refined_ents_ptr);
718 CHKERR m_field.getInterface<BitRefManager>()->setBitRefLevel(
719 ents_to_set_bit,
bit,
true, verb, &ents);
Public Member Functions inherited from MoFEM::UnknownInterface