source: XIOS/dev/branch_yushan/src/node/mesh.hpp @ 1109

Last change on this file since 1109 was 1109, checked in by yushan, 7 years ago

test_omp OK

File size: 3.7 KB
Line 
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
13namespace 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 ep_lib::MPI_Comm&,
63                             const CArray<double, 1>&, const CArray<double, 1>&,
64                             const CArray<double, 2>&, const CArray<double, 2>& );
65
66      void getGlobalNghbFaces(const int, const ep_lib::MPI_Comm&, const CArray<int, 1>&,
67                              const CArray<double, 2>&, const CArray<double, 2>&,
68                              CArray<int, 2>&);
69
70      void getLocalNghbFaces(const int, const CArray<int, 1>&,
71                             const CArray<double, 2>&, const CArray<double, 2>&,
72                             CArray<int, 2>&, CArray<int, 1>&);
73           
74      static CMesh* getMesh(StdString, int);
75
76    private:
77
78      int nbNodes_;
79      int nbEdges_;
80      int nbFaces_;
81
82      static std::map <StdString, CMesh> meshList;
83      static std::map <StdString, vector<int> > domainList;
84
85      static std::map <StdString, CMesh> *meshList_ptr;
86      static std::map <StdString, vector<int> > *domainList_ptr;
87      #pragma omp threadprivate(meshList_ptr, domainList_ptr)
88   
89      CClientClientDHTSizet* pNodeGlobalIndex;                    // pointer to a map <nodeHash, nodeIdxGlo>
90      CClientClientDHTSizet* pEdgeGlobalIndex;                    // pointer to a map <edgeHash, edgeIdxGlo>
91      void getGloNghbFacesNodeType(const ep_lib::MPI_Comm&, const CArray<int, 1>&, const CArray<double, 2>&, const CArray<double, 2>&, CArray<int, 2>&);
92      void getGloNghbFacesEdgeType(const ep_lib::MPI_Comm&, const CArray<int, 1>&, const CArray<double, 2>&, const CArray<double, 2>&, CArray<int, 2>&);
93      void getLocNghbFacesNodeType(const CArray<int, 1>&, const CArray<double, 2>&, const CArray<double, 2>&, CArray<int, 2>&, CArray<int, 1>&);
94      void getLocNghbFacesEdgeType(const CArray<int, 1>&, const CArray<double, 2>&, const CArray<double, 2>&, CArray<int, 2>&, CArray<int, 1>&);
95
96      vector<size_t> createHashes (const double, const double);
97
98      size_t nodeIndex (double, double);                           // redundant in parallel version with epsilon precision
99      boost::unordered_map <size_t, size_t> hashed_map_nodes;      // redundant in parallel version with epsilon precision
100      boost::unordered_map <pair<double,double>, int> map_nodes;   // redundant in parallel version with epsilon precision
101      boost::unordered_map <pair<int,int>, int> map_edges;         // redundant in parallel version with epsilon precision
102
103  }; 
104
105} // namespace xios
106
107#endif //__XIOS_CMesh__
Note: See TracBrowser for help on using the repository browser.