#include "users_modules/eshelbian_plasticity/src/EshelbianContact.hpp"

Inheritance diagram for EshelbianPlasticity::ContactTree:

Collaboration diagram for EshelbianPlasticity::ContactTree:

Classes
struct	FaceData

Public Types
using	Base = PostProcBrokenMeshInMoabBase<FaceElementForcesAndSourcesCore>

using	MapFaceData = std::map<EntityHandle, std::vector<FaceData>>

Public Member Functions
	ContactTree (MoFEM::Interface &m_field, boost::shared_ptr< moab::Core > core_mesh_ptr, int max_order, std::map< int, Range > &&body_map)

MoFEMErrorCode	preProcess ()

MoFEMErrorCode	postProcess ()

boost::shared_ptr< OrientedBoxTreeTool > &	getTreeSurfPtr ()

auto	getRootSetSurf ()

int	getMaxLevel () const

auto	findFaceDataVecPtr (EntityHandle fe_ent)

auto	getFaceDataPtr (std::vector< FaceData >::iterator &it, int gg, std::vector< FaceData > *vec_ptr)

Protected Member Functions
MoFEMErrorCode	buildTree (Range &ents)

Protected Attributes
boost::shared_ptr< OrientedBoxTreeTool >	treeSurfPtr

EntityHandle	rootSetSurf

Tag	thEleId

Tag	thBodyId

Tag	thSmallX

Tag	thLargeX

Tag	thTraction

MapFaceData	shadowDataMap

std::map< int, Range >	bodyMap

const int	maxOrder

Friends
struct	OpMoveNode

struct	OpTreeSearch

Detailed Description

Definition at line 34 of file EshelbianContact.hpp.

Member Typedef Documentation

◆ Base

using EshelbianPlasticity::ContactTree::Base = PostProcBrokenMeshInMoabBase<FaceElementForcesAndSourcesCore>

Definition at line 37 of file EshelbianContact.hpp.

◆ MapFaceData

using EshelbianPlasticity::ContactTree::MapFaceData = std::map<EntityHandle, std::vector<FaceData>>

Definition at line 74 of file EshelbianContact.hpp.

Constructor & Destructor Documentation

◆ ContactTree()

EshelbianPlasticity::ContactTree::ContactTree	(	MoFEM::Interface &	m_field,
		boost::shared_ptr< moab::Core >	core_mesh_ptr,
		int	max_order,
		std::map< int, Range > &&	body_map )

< So far only linear geometry is implemented for contact

Definition at line 961 of file EshelbianContact.cpp.

    : Base(m_field, core_mesh_ptr, "contact"), maxOrder(max_order),
      bodyMap(body_map) {
 
  auto ref_ele_ptr = boost::make_shared<PostProcGenerateRefMesh<MBTRI>>();
  ref_ele_ptr->hoNodes =
      PETSC_FALSE; ///< So far only linear geometry is implemented for contact
 
  CHK_THROW_MESSAGE(ref_ele_ptr->getOptions(optionsPrefix), "getOptions");
  CHK_THROW_MESSAGE(ref_ele_ptr->generateReferenceElementMesh(),
                    "Error when generating reference element");
 
  MOFEM_LOG("EP", Sev::inform) << "Contact hoNodes " << ref_ele_ptr->hoNodes
      ? "true"
      : "false";
  MOFEM_LOG("EP", Sev::inform)
      << "Contact maxOrder " << ref_ele_ptr->defMaxLevel;
 
  refElementsMap[MBTRI] = ref_ele_ptr;
 
  int def_ele_id = -1;
  CHKERR getPostProcMesh().tag_get_handle("ELE_ID", 1, MB_TYPE_INTEGER, thEleId,
                                          MB_TAG_CREAT | MB_TAG_DENSE,
                                          &def_ele_id);
  CHKERR getPostProcMesh().tag_get_handle("BODY_ID", 1, MB_TYPE_INTEGER,
                                          thBodyId, MB_TAG_CREAT | MB_TAG_DENSE,
                                          &def_ele_id);
 
  // std::vector<int> def_ids(3 * nb_bases, 0);
  // CHKERR getPostProcMesh().tag_get_handle("IDS", 3 * nb_bases,
  // MB_TYPE_INTEGER,
  //                                         thIds, MB_TAG_CREAT | MB_TAG_DENSE,
  //                                         &*def_ids.begin());
  // std::vector<double> def_coeffs(3 * nb_bases, 0);
  // CHKERR getPostProcMesh().tag_get_handle("COEFF", 3 * nb_bases,
  // MB_TYPE_DOUBLE,
  //                                         thCoeff, MB_TAG_CREAT |
  //                                         MB_TAG_DENSE,
  //                                         &*def_coeffs.begin());
  // std::vector<double> def_basses(nb_bases, 0);
  // CHKERR getPostProcMesh().tag_get_handle("BASES", nb_bases, MB_TYPE_DOUBLE,
  //                                         thBases, MB_TAG_CREAT |
  //                                         MB_TAG_DENSE,
  //                                         &*def_basses.begin());
 
  std::array<double, 3> def_small_x{0., 0., 0.};
  CHKERR getPostProcMesh().tag_get_handle("x", 3, MB_TYPE_DOUBLE, thSmallX,
                                          MB_TAG_CREAT | MB_TAG_DENSE,
                                          def_small_x.data());
  CHKERR getPostProcMesh().tag_get_handle("X", 3, MB_TYPE_DOUBLE, thLargeX,
                                          MB_TAG_CREAT | MB_TAG_DENSE,
                                          def_small_x.data());
 
  std::array<double, 3> def_tractions{0., 0., 0.};
  CHKERR getPostProcMesh().tag_get_handle(
      "TRACTION", 3, MB_TYPE_DOUBLE, thTraction, MB_TAG_CREAT | MB_TAG_DENSE,
      &*def_tractions.begin());
}

Member Function Documentation

◆ buildTree()

MoFEMErrorCode EshelbianPlasticity::ContactTree::buildTree ( Range & ents )

protected

Definition at line 1071 of file EshelbianContact.cpp.

                                                 {
  MoFEMFunctionBegin;
  // MOFEM_LOG("SELF", Sev::inform) << ents << endl;
 
  treeSurfPtr = boost::shared_ptr<OrientedBoxTreeTool>(
      new OrientedBoxTreeTool(&getPostProcMesh(), "ROOTSETSURF", true));
  CHKERR treeSurfPtr->build(ents, rootSetSurf);
 
  // CHKERR treeSurfPtr->stats(rootSetSurf, std::cerr);
 
  MoFEMFunctionReturn(0);
}

◆ findFaceDataVecPtr()

auto EshelbianPlasticity::ContactTree::findFaceDataVecPtr ( EntityHandle fe_ent )

inline

Definition at line 76 of file EshelbianContact.hpp.

                                                      {
    auto &map_face_data = shadowDataMap;
    auto it = map_face_data.find(fe_ent);
    if (it == map_face_data.end()) {
      return (std::vector<FaceData> *)nullptr;
    }
    return &(it->second);
  }

◆ getFaceDataPtr()

auto EshelbianPlasticity::ContactTree::getFaceDataPtr	(	std::vector< FaceData >::iterator &	it,
		int	gg,
		std::vector< FaceData > *	vec_ptr )

inline

Definition at line 85 of file EshelbianContact.hpp.

                                                           {
    FaceData *face_data_ptr = nullptr;
    if (it != vec_ptr->end()) {
      if (it->gaussPtNb == gg) {
        face_data_ptr = &(*it);
        ++it;
      }
    }
    return face_data_ptr;
  }

◆ getMaxLevel()

int EshelbianPlasticity::ContactTree::getMaxLevel ( ) const

inline

Definition at line 49 of file EshelbianContact.hpp.

49{ return refElementsMap.at(MBTRI)->defMaxLevel; }

◆ getRootSetSurf()

auto EshelbianPlasticity::ContactTree::getRootSetSurf ( )

inline

Definition at line 48 of file EshelbianContact.hpp.

48{ return rootSetSurf; }

◆ getTreeSurfPtr()

boost::shared_ptr< OrientedBoxTreeTool > & EshelbianPlasticity::ContactTree::getTreeSurfPtr ( )

inline

Definition at line 45 of file EshelbianContact.hpp.

                                                                {
    return treeSurfPtr;
  }

◆ postProcess()

MoFEMErrorCode EshelbianPlasticity::ContactTree::postProcess ( )

Definition at line 1028 of file EshelbianContact.cpp.

                                        {
  MoFEMFunctionBegin;
 
  CHKERR Base::postProcess();
  shadowDataMap.clear();
 
  PetscBarrier(nullptr);
 
  ParallelComm *pcomm_post_proc_mesh =
      ParallelComm::get_pcomm(&getPostProcMesh(), MYPCOMM_INDEX);
  if (pcomm_post_proc_mesh == nullptr)
    SETERRQ(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY, "PComm not allocated");
 
  // MOFEM_LOG("EP", Sev::verbose) << "ContactTree broadcast entities";
 
  auto brodacts = [&](auto &brodacts_ents) {
    MoFEMFunctionBegin;
    Range recived_ents;
    for (auto rr = 0; rr != mField.get_comm_size(); ++rr) {
      pcomm_post_proc_mesh->broadcast_entities(
          rr, rr == mField.get_comm_rank() ? brodacts_ents : recived_ents,
          false, true);
    }
    MoFEMFunctionReturn(0);
  };
 
  Range brodacts_ents;
  CHKERR getPostProcMesh().get_entities_by_dimension(0, 2, brodacts_ents, true);
  CHKERR brodacts(brodacts_ents);
 
  // MOFEM_LOG("EP", Sev::verbose) << "ContactTree build tree";
 
  Range ents;
  CHKERR getPostProcMesh().get_entities_by_dimension(0, 2, ents, true);
  CHKERR buildTree(ents);
 
#ifndef NDEBUG
  CHKERR writeFile("debug_tree.h5m");
#endif // NDEBUG
 
  MoFEMFunctionReturn(0);
}

◆ preProcess()

MoFEMErrorCode EshelbianPlasticity::ContactTree::preProcess ( )

Definition at line 1022 of file EshelbianContact.cpp.

                                       {
  MoFEMFunctionBegin;
  CHKERR Base::preProcess();
  MoFEMFunctionReturn(0);
}

Friends And Related Symbol Documentation

◆ OpMoveNode

friend struct OpMoveNode

friend

Definition at line 51 of file EshelbianContact.hpp.

◆ OpTreeSearch

friend struct OpTreeSearch

friend

Definition at line 52 of file EshelbianContact.hpp.

Member Data Documentation

◆ bodyMap

std::map<int, Range> EshelbianPlasticity::ContactTree::bodyMap

protected

Definition at line 113 of file EshelbianContact.hpp.

◆ maxOrder

const int EshelbianPlasticity::ContactTree::maxOrder

protected

Definition at line 115 of file EshelbianContact.hpp.

◆ rootSetSurf

EntityHandle EshelbianPlasticity::ContactTree::rootSetSurf

protected

Definition at line 100 of file EshelbianContact.hpp.

◆ shadowDataMap

MapFaceData EshelbianPlasticity::ContactTree::shadowDataMap

protected

Definition at line 111 of file EshelbianContact.hpp.

◆ thBodyId

Tag EshelbianPlasticity::ContactTree::thBodyId

protected

Definition at line 106 of file EshelbianContact.hpp.

◆ thEleId

Tag EshelbianPlasticity::ContactTree::thEleId

protected

Definition at line 102 of file EshelbianContact.hpp.

◆ thLargeX

Tag EshelbianPlasticity::ContactTree::thLargeX

protected

Definition at line 108 of file EshelbianContact.hpp.

◆ thSmallX

Tag EshelbianPlasticity::ContactTree::thSmallX

protected

Definition at line 107 of file EshelbianContact.hpp.

◆ thTraction

Tag EshelbianPlasticity::ContactTree::thTraction

protected

Definition at line 109 of file EshelbianContact.hpp.

◆ treeSurfPtr

boost::shared_ptr<OrientedBoxTreeTool> EshelbianPlasticity::ContactTree::treeSurfPtr

protected

Definition at line 99 of file EshelbianContact.hpp.

The documentation for this struct was generated from the following files:

users_modules/eshelbian_plasticity/src/EshelbianContact.hpp
users_modules/eshelbian_plasticity/src/impl/EshelbianContact.cpp

Classes

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Friends

Detailed Description

Member Typedef Documentation

◆ Base

◆ MapFaceData

Constructor & Destructor Documentation

◆ ContactTree()

Member Function Documentation

◆ buildTree()

◆ findFaceDataVecPtr()

◆ getFaceDataPtr()

◆ getMaxLevel()

◆ getRootSetSurf()

◆ getTreeSurfPtr()

◆ postProcess()

◆ preProcess()

Friends And Related Symbol Documentation

◆ OpMoveNode

◆ OpTreeSearch

Member Data Documentation

◆ bodyMap

◆ maxOrder

◆ rootSetSurf

◆ shadowDataMap

◆ thBodyId

◆ thEleId

◆ thLargeX

◆ thSmallX

◆ thTraction

◆ treeSurfPtr