| 1 | /*! |
|---|
| 2 | \file mesh.hpp |
|---|
| 3 | \author Olga Abramkina |
|---|
| 4 | \brief Declaration of class CMesh. |
|---|
| 5 | */ |
|---|
| 6 | |
|---|
| 7 | #ifndef __XIOS_CMesh__ |
|---|
| 8 | #define __XIOS_CMesh__ |
|---|
| 9 | |
|---|
| 10 | #include "array_new.hpp" |
|---|
| 11 | #include "dht_auto_indexing.hpp" |
|---|
| 12 | |
|---|
| 13 | namespace xios { |
|---|
| 14 | |
|---|
| 15 | /// ////////////////////// Déclarations ////////////////////// /// |
|---|
| 16 | |
|---|
| 17 | class CMesh; |
|---|
| 18 | |
|---|
| 19 | ///-------------------------------------------------------------- |
|---|
| 20 | |
|---|
| 21 | /*! |
|---|
| 22 | \class CMesh |
|---|
| 23 | \brief CMesh class. |
|---|
| 24 | */ |
|---|
| 25 | class CMesh |
|---|
| 26 | { |
|---|
| 27 | |
|---|
| 28 | public: |
|---|
| 29 | |
|---|
| 30 | CMesh(void); |
|---|
| 31 | ~CMesh(void); |
|---|
| 32 | |
|---|
| 33 | int nbNodesGlo; |
|---|
| 34 | int nbEdgesGlo; |
|---|
| 35 | |
|---|
| 36 | int node_start; |
|---|
| 37 | int node_count; |
|---|
| 38 | int edge_start; |
|---|
| 39 | int edge_count; |
|---|
| 40 | |
|---|
| 41 | bool nodesAreWritten; |
|---|
| 42 | bool edgesAreWritten; |
|---|
| 43 | bool facesAreWritten; |
|---|
| 44 | |
|---|
| 45 | CArray<double, 1> node_lon; |
|---|
| 46 | CArray<double, 1> node_lat; |
|---|
| 47 | |
|---|
| 48 | CArray<double, 1> edge_lon; |
|---|
| 49 | CArray<double, 1> edge_lat; |
|---|
| 50 | CArray<int, 2> edge_nodes; |
|---|
| 51 | |
|---|
| 52 | CArray<double, 1> face_lon; |
|---|
| 53 | CArray<double, 1> face_lat; |
|---|
| 54 | CArray<int, 2> face_nodes; |
|---|
| 55 | CArray<int, 2> face_edges; |
|---|
| 56 | CArray<int, 2> edge_faces; |
|---|
| 57 | CArray<int, 2> face_faces; |
|---|
| 58 | |
|---|
| 59 | void createMesh(const CArray<double, 1>&, const CArray<double, 1>&, |
|---|
| 60 | const CArray<double, 2>&, const CArray<double, 2>& ); |
|---|
| 61 | |
|---|
| 62 | void createMeshEpsilon(const MPI_Comm&, |
|---|
| 63 | const CArray<double, 1>&, const CArray<double, 1>&, |
|---|
| 64 | const CArray<double, 2>&, const CArray<double, 2>& ); |
|---|
| 65 | |
|---|
| 66 | void getDomainExtended(const int, const MPI_Comm&, |
|---|
| 67 | const CArray<double, 2>&, const CArray<double, 2>&, |
|---|
| 68 | CArray<int, 2>&); |
|---|
| 69 | |
|---|
| 70 | static CMesh* getMesh(StdString, int); |
|---|
| 71 | |
|---|
| 72 | private: |
|---|
| 73 | |
|---|
| 74 | int nbNodes_; |
|---|
| 75 | int nbEdges_; |
|---|
| 76 | int nbFaces_; |
|---|
| 77 | |
|---|
| 78 | static std::map <StdString, CMesh> meshList; |
|---|
| 79 | static std::map <StdString, vector<int> > domainList; |
|---|
| 80 | CClientClientDHTSizet* pNodeGlobalIndex; // pointer to a map <nodeHash, nodeIdxGlo> |
|---|
| 81 | CClientClientDHTSizet* pEdgeGlobalIndex; // pointer to a map <edgeHash, edgeIdxGlo> |
|---|
| 82 | void getNghbFacesNodeType(const MPI_Comm&, const CArray<double, 2>&, const CArray<double, 2>&, CArray<int, 2>&); |
|---|
| 83 | void getNghbFacesEdgeType(const MPI_Comm&, const CArray<double, 2>&, const CArray<double, 2>&, CArray<int, 2>&); |
|---|
| 84 | |
|---|
| 85 | vector<size_t> createHashes (const double, const double); |
|---|
| 86 | |
|---|
| 87 | size_t nodeIndex (double, double); // redundant in parallel version with epsilon precision |
|---|
| 88 | boost::unordered_map <size_t, size_t> hashed_map_nodes; // redundant in parallel version with epsilon precision |
|---|
| 89 | boost::unordered_map <pair<double,double>, int> map_nodes; // redundant in parallel version with epsilon precision |
|---|
| 90 | boost::unordered_map <pair<int,int>, int> map_edges; // redundant in parallel version with epsilon precision |
|---|
| 91 | |
|---|
| 92 | }; |
|---|
| 93 | |
|---|
| 94 | } // namespace xios |
|---|
| 95 | |
|---|
| 96 | #endif //__XIOS_CMesh__ |
|---|
