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

Last change on this file since 929 was 929, checked in by oabramkina, 5 years ago

Mesh connectivity:

Fixing a bug for UGRID in case meshes containing cells with different number of vertexes.

Function for determining face neighbors of a local domain added.

File size: 2.9 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 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__
Note: See TracBrowser for help on using the repository browser.