PMesh.H

Go to the documentation of this file.

/*
                  Supported Cell Types

        4                                              6   
       /|\             5       8---------7            /|\ 
      / | \         /  |       | \       |\          / | \
     /  |  \       /  /        |  \      | \        /  |  \
    /   |   \     /   |        |   5--------\6     /   3   \
   1----|----3   4----|-3      4---|-----3   |    4---------5
    \   |   /     \  /   \      \  |      \  |    |  /   \  |
     \  |  /       \ |    \      \ |       \ |    | /     \ |
      \ | /         \|     \      \|        \|    |/       \|
        2            1------2      1---------2    1---------2
 
*/
#ifndef _PMESH_H_
#define _PMESH_H_
#include "Mesh.H"
#include "COMM.H"

namespace Mesh {  
  struct BorderData {
    std::vector<Mesh::IndexType> SendAp; 
    std::vector<Mesh::IndexType> SendAi;
    std::vector<Mesh::IndexType> RecvAp;
    std::vector<Mesh::IndexType> RecvAi; 
    std::vector<Mesh::IndexType> border_dof_to_global;
    std::vector<double> SendBuffer;
    std::vector<double> RecvBuffer;
  };
  struct Border {
    Mesh::IndexType rpart; 
    Mesh::IndexVec nsend;  
    Mesh::IndexVec nrecv;  
    Mesh::IndexType recvsize; 
    Mesh::IndexType sendsize; 
    Mesh::IndexType remote_dofcount; 
    BorderData data; 
  };

  struct PartInfo {
    Mesh::IndexType npart;   
    Mesh::IndexType part;    
    Mesh::IndexType nelem;   
    Mesh::IndexType nnodes;  
    Mesh::IndexType nborder; 
    Mesh::IndexType nshared; 
    Mesh::IndexType nown;    
    Mesh::IndexType nlocal;  
    Mesh::IndexType doffset; 
  };
  
  class Partition {
  private:
    IRAD::Comm::CommunicatorObject *_communicator;
  public:
    Partition() : _communicator(NULL), _mesh(),_info(),_borders() {};
    MeshUtilityObject _mesh;
    //    NodalCoordinates    _nc;
    //    Connectivity        _ec;
    PartInfo            _info;
    std::vector<Border> _borders;
    int Read(const std::string &meshname,int id);
    int Read(const std::string &MeshName,IRAD::Comm::CommunicatorObject &comm,
             bool allow_n2m,std::ostream &ErrOut);
    //    IRAD::Comm::CommunicatorObject *CommunicatorPtr(){ return(_communicator); };
    //    std::ostream &Report(std::ostream &Ostr);
    MeshUtilityObject &Mesh() { return (_mesh); };
    const MeshUtilityObject &Mesh() const { return (_mesh); };
    PartInfo &Info() { return (_info); };
    const PartInfo &Info() const { return(_info); };
    std::vector<Border> &Borders() { return(_borders); };
    const std::vector<Border> &Borders() const { return(_borders); };
    IRAD::Comm::CommunicatorObject *CommunicatorPtr(){ return(_communicator); };
     //    std::ostream &Report(std::ostream &Ostr);
    int GetBorderElements(std::vector<Mesh::IndexType> &be) const;
  };
  
}
#endif