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

Last change on this file since 931 was 931, checked in by oabramkina, 6 years ago

Mesh connectivity:

Functions for determining local and global face-to-face connectivity added. Both types of neighboring cells (those that share a node or those that share an edge) are considered.

File size: 3.5 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      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      CClientClientDHTSizet* pNodeGlobalIndex;                    // pointer to a map <nodeHash, nodeIdxGlo>
85      CClientClientDHTSizet* pEdgeGlobalIndex;                    // pointer to a map <edgeHash, edgeIdxGlo>
86      void getGloNghbFacesNodeType(const MPI_Comm&, const CArray<int, 1>&, const CArray<double, 2>&, const CArray<double, 2>&, CArray<int, 2>&);
87      void getGloNghbFacesEdgeType(const MPI_Comm&, const CArray<int, 1>&, const CArray<double, 2>&, const CArray<double, 2>&, CArray<int, 2>&);
88      void getLocNghbFacesNodeType(const CArray<int, 1>&, const CArray<double, 2>&, const CArray<double, 2>&, CArray<int, 2>&, CArray<int, 1>&);
89      void getLocNghbFacesEdgeType(const CArray<int, 1>&, const CArray<double, 2>&, const CArray<double, 2>&, CArray<int, 2>&, CArray<int, 1>&);
90
91      vector<size_t> createHashes (const double, const double);
92
93      size_t nodeIndex (double, double);                           // redundant in parallel version with epsilon precision
94      boost::unordered_map <size_t, size_t> hashed_map_nodes;      // redundant in parallel version with epsilon precision
95      boost::unordered_map <pair<double,double>, int> map_nodes;   // redundant in parallel version with epsilon precision
96      boost::unordered_map <pair<int,int>, int> map_edges;         // redundant in parallel version with epsilon precision
97
98  }; 
99
100} // namespace xios
101
102#endif //__XIOS_CMesh__
Note: See TracBrowser for help on using the repository browser.