source: XIOS/trunk/src/node/mesh.hpp @ 2250

Last change on this file since 2250 was 1639, checked in by yushan, 5 years ago

revert erroneous commit on trunk

File size: 3.1 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 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 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      typedef std::pair<double, double> pairDouble;
79      typedef std::pair<int, int> pairInt;
80
81      int nbNodes_;
82      int nbEdges_;
83      int nbFaces_;
84
85      static std::map <StdString, CMesh> meshList;
86      static std::map <StdString, vector<int> > domainList;
87      CClientClientDHTSizet* pNodeGlobalIndex;                    // pointer to a map <nodeHash, nodeIdxGlo>
88      CClientClientDHTSizet* pEdgeGlobalIndex;                    // pointer to a map <edgeHash, edgeIdxGlo>
89      void getGloNghbFacesNodeType(const MPI_Comm&, const CArray<int, 1>&, const CArray<double, 2>&, const CArray<double, 2>&, CArray<int, 2>&);
90      void getGloNghbFacesEdgeType(const MPI_Comm&, const CArray<int, 1>&, const CArray<double, 2>&, const CArray<double, 2>&, CArray<int, 2>&);
91      void getLocNghbFacesNodeType(const CArray<int, 1>&, const CArray<double, 2>&, const CArray<double, 2>&, CArray<int, 2>&, CArray<int, 1>&);
92      void getLocNghbFacesEdgeType(const CArray<int, 1>&, const CArray<double, 2>&, const CArray<double, 2>&, CArray<int, 2>&, CArray<int, 1>&);
93
94      vector<size_t> createHashes (const double, const double);
95
96  }; 
97
98} // namespace xios
99
100#endif //__XIOS_CMesh__
Note: See TracBrowser for help on using the repository browser.