Detailed Description

Definition at line 60 of file meshPartitioner.H.

Public Member Functions
	meshPartition (int pidx, const std::vector< int > &glbElmPartedIdx, const std::vector< int > &glbElmConn, MeshType_t inMshType)

	~meshPartition ()=default

std::vector< double >	getCrds (const std::vector< double > &crds) const

std::vector< int >	getConns () const

std::vector< double >	getElmSlns (const std::vector< double > &slns) const

std::vector< double >	getElmSlnsVec (const std::vector< double > &slns, int nComp) const

std::map< int, int >	getPartToGlobNodeMap () const

std::map< int, int >	getPartToGlobElmMap () const

Public Attributes
int	nNde

int	nElm

Private Attributes
int	pIdx

int	nNdeElm

std::vector< int >	globNdeIdx

std::vector< int >	globElmIdx

std::vector< int >	partElmConn

std::map< int, int >	ndeIdxGlobToPart

std::map< int, int >	ndeIdxPartToGlob

std::map< int, int >	elmIdxGlobToPart

std::map< int, int >	elmIdxPartToGlob

MeshType_t	mshType

Constructor & Destructor Documentation

◆ meshPartition()

meshPartition::meshPartition	(	int	pidx,
		const std::vector< int > &	glbElmPartedIdx,
		const std::vector< int > &	glbElmConn,
		MeshType_t	inMshType
	)

Definition at line 42 of file meshPartitioner.C.

References elmIdxGlobToPart, elmIdxPartToGlob, globElmIdx, globNdeIdx, MESH_TETRA_4, MESH_TRI_3, mshType, ndeIdxGlobToPart, ndeIdxPartToGlob, nElm, nNde, nNdeElm, partElmConn, and pIdx.

Referenced by meshPartitioner::buildPartitions().

 {
   mshType = inMshType;
   pIdx = pidx;
   nNde = 0;
   nElm = 0;
   if (mshType == MeshType_t::MESH_TETRA_4)
     nNdeElm = 4;
   else if (mshType == MeshType_t::MESH_TRI_3)
     nNdeElm = 3;
   // finding partition elements
   int glbElmIdx = 1;
   for (int ie : glbElmPartedIdx)
   {
     if (ie == pIdx)
     {
       nElm++;
       elmIdxGlobToPart[glbElmIdx] = nElm;
       elmIdxPartToGlob[nElm] = glbElmIdx;
       globElmIdx.push_back(glbElmIdx);
       //std::cout << "My element = " << glbElmIdx << std::endl;
       // attaching nodes of this element to the current partition
       for (int iNde = 0; iNde < nNdeElm; iNde++)
       {
         int glbNdeIdx = glbElmConn[(glbElmIdx - 1) * nNdeElm + iNde];
         auto it = ndeIdxGlobToPart.find(glbNdeIdx);
         int partNdeIdx;
         if (it == ndeIdxGlobToPart.end())
         {
           nNde++;
           partNdeIdx = nNde;
           ndeIdxGlobToPart[glbNdeIdx] = nNde;
           ndeIdxPartToGlob[nNde] = glbNdeIdx;
           globNdeIdx.push_back(glbNdeIdx);
           //std::cout << "partition " << pIdx
           //          << " global node id " << glbNdeIdx
           //          << " is local " << nNde
           //          << std::endl;
         }
         else
         {
           partNdeIdx = it->second;
         }
         partElmConn.push_back(partNdeIdx);
       }
     }
     glbElmIdx++;
   }
   //std::cout << "Min connectivity passed to partitioner = "
   //          << *std::min_element(glbElmConn.begin(), glbElmConn.end())
   //          << std::endl;
 
   // Debug information
   //std::cout << "I am partition " << pidx << std::endl
   //          << " cycled on " << glbElmIdx - 1 << " global elements "
   //          << std::endl
   //          << " my min global element idx is "
   //          << *std::min_element(globElmIdx.begin(), globElmIdx.end())
   //          << std::endl
   //          << " my max global element idx is "
   //          << *std::max_element(globElmIdx.begin(), globElmIdx.end())
   //          << std::endl
   //          << " my nNde = " << nNde << " my nElm = " << nElm << std::endl;
 }

◆ ~meshPartition()

meshPartition::~meshPartition ( )

default

Member Function Documentation

◆ getConns()

std::vector<int> meshPartition::getConns ( ) const

inline

Definition at line 71 of file meshPartitioner.H.

71 { return partElmConn; }

meshPartition::partElmConn

std::vector< int > partElmConn

Definition: meshPartitioner.H:87

◆ getCrds()

std::vector< double > meshPartition::getCrds ( const std::vector< double > & crds ) const

Definition at line 111 of file meshPartitioner.C.

References globNdeIdx.

 {
   std::vector<double> x;
   for (int in : globNdeIdx)
     x.push_back(crds[in - 1]);
   return x;
 }

◆ getElmSlns()

std::vector< double > meshPartition::getElmSlns ( const std::vector< double > & slns ) const

Definition at line 120 of file meshPartitioner.C.

References globElmIdx.

 {
   std::vector<double> x;
   for (int ie : globElmIdx)
     x.push_back(slns[ie - 1]);
   return x;
 }

◆ getElmSlnsVec()

std::vector< double > meshPartition::getElmSlnsVec	(	const std::vector< double > &	slns,
		int	nComp
	)		const

Definition at line 130 of file meshPartitioner.C.

References globElmIdx.

 {
   std::vector<double> x;
   for (int ie : globElmIdx)
     for (int iComp = 0; iComp < nComp; iComp++)
       x.push_back(slns[(ie - 1) * nComp + iComp]);
   return x;
 }