source: XIOS/trunk/src/node/domain.hpp @ 809

Last change on this file since 809 was 809, checked in by ymipsl, 7 years ago

Fix problems for interpolation onto regular domain.

YM

  • Property copyright set to
    Software name : XIOS (Xml I/O Server)
    http://forge.ipsl.jussieu.fr/ioserver
    Creation date : January 2009
    Licence : CeCCIL version2
    see license file in root directory : Licence_CeCILL_V2-en.txt
    or http://www.cecill.info/licences/Licence_CeCILL_V2-en.html
    Holder : CEA/LSCE (Laboratoire des Sciences du CLimat et de l'Environnement)
    CNRS/IPSL (Institut Pierre Simon Laplace)
    Project Manager : Yann Meurdesoif
    yann.meurdesoif@cea.fr
File size: 7.4 KB
Line 
1#ifndef __XIOS_CDomain__
2#define __XIOS_CDomain__
3
4/// XIOS headers ///
5#include "xios_spl.hpp"
6#include "group_factory.hpp"
7
8#include "declare_group.hpp"
9#include "declare_ref_func.hpp"
10#include "event_client.hpp"
11#include "event_server.hpp"
12#include "buffer_in.hpp"
13#include "array_new.hpp"
14#include "attribute_array.hpp"
15#include "attribute_enum.hpp"
16#include "transformation.hpp"
17#include "transformation_enum.hpp"
18
19namespace xios {
20
21   /// ////////////////////// Déclarations ////////////////////// ///
22
23   class CDomainGroup;
24   class CDomainAttributes;
25   class CDomain;
26
27   ///--------------------------------------------------------------
28
29   // Declare/Define CDomainAttribute
30   BEGIN_DECLARE_ATTRIBUTE_MAP(CDomain)
31#  include "domain_attribute.conf"
32#  include "domain_attribute_private.conf"
33   END_DECLARE_ATTRIBUTE_MAP(CDomain)
34
35   ///--------------------------------------------------------------
36
37   class CDomain
38      : public CObjectTemplate<CDomain>
39      , public CDomainAttributes
40   {
41         enum EEventId
42         {
43           EVENT_ID_SERVER_ATTRIBUT, EVENT_ID_INDEX, EVENT_ID_LON, EVENT_ID_LAT, EVENT_ID_AREA
44         } ;
45
46         /// typedef ///
47         typedef CObjectTemplate<CDomain>   SuperClass;
48         typedef CDomainAttributes SuperClassAttribute;
49
50      public:
51
52         typedef CDomainAttributes RelAttributes;
53         typedef CDomainGroup      RelGroup;
54         typedef CTransformation<CDomain>::TransformationMapTypes TransMapTypes;
55
56         /// Constructeurs ///
57         CDomain(void);
58         explicit CDomain(const StdString & id);
59         CDomain(const CDomain & domain);       // Not implemented yet.
60         CDomain(const CDomain * const domain); // Not implemented yet.
61
62         static CDomain* createDomain();
63
64         virtual void parse(xml::CXMLNode & node);
65
66         /// Vérifications ///
67         void checkAttributes(void);
68
69         void checkAttributesOnClient();
70         void checkAttributesOnClientAfterTransformation();
71
72         void checkEligibilityForCompressedOutput(void);
73
74         void sendCheckedAttributes();
75
76         bool hasTransformation();
77         void solveInheritanceTransformation();
78         TransMapTypes getAllTransformations();
79         void redistribute(int nbLocalDomain);
80
81      public:
82         const std::set<StdString> & getRelFiles(void) const;
83         bool IsWritten(const StdString & filename) const;
84         bool isWrittenCompressed(const StdString& filename) const;
85
86         const std::vector<int>& getIndexesToWrite(void) const;
87         int getNumberWrittenIndexes() const;
88         int getTotalNumberWrittenIndexes() const;
89         int getOffsetWrittenIndexes() const;
90
91         std::map<int, StdSize> getAttributesBufferSize();
92
93         bool isEmpty(void) const;
94         bool isDistributed(void) const;
95         bool isCompressible(void) const;
96
97         int global_zoom_ibegin, global_zoom_ni;
98         int global_zoom_jbegin, global_zoom_nj;
99
100         int ni_srv,ibegin_srv,iend_srv ;
101         int zoom_ni_srv,zoom_ibegin_srv,zoom_iend_srv ;
102
103         int nj_srv,jbegin_srv,jend_srv ;
104         int zoom_nj_srv,zoom_jbegin_srv,zoom_jend_srv ;
105
106         CArray<double, 1> lonvalue_srv, latvalue_srv ;
107         CArray<double, 2> bounds_lon_srv, bounds_lat_srv ;
108         CArray<double, 1> lonvalue_client, latvalue_client;
109         CArray<double, 2> bounds_lon_client, bounds_lat_client;
110         CArray<double, 1> area_srv;
111
112        vector<int> connectedServer ; // list of connected server
113        vector<int> nbSenders ; // for each communication with a server, number of communicating client
114        vector<int> nbDataSrv ; // size of data to send to each server
115        vector< vector<int> > i_indSrv ; // for each server, i global index to send
116        vector< vector<int> > j_indSrv ; // for each server, j global index to send
117
118      public:
119         /// Mutateur ///
120         void addRelFile(const StdString & filename);
121         void addRelFileCompressed(const StdString& filename);
122         void completeLonLatClient(void);
123         void sendServerAttribut(void) ;
124         void sendLonLatArea(void);
125         void computeConnectedServer(void) ;
126         
127         void AllgatherRectilinearLonLat(CArray<double,1>& lon, CArray<double,1>& lat, 
128                                         CArray<double,1>& lon_g, CArray<double,1>& lat_g);
129         
130         void fillInRectilinearBoundLonLat(CArray<double,1>& lon, CArray<double,1>& lat, 
131                                           CArray<double,2>& boundsLon, CArray<double,2>& boundsLat);
132         void fillInRectilinearLonLat();
133
134         static bool dispatchEvent(CEventServer& event);
135         static void recvServerAttribut(CEventServer& event);
136         static void recvIndex(CEventServer& event);
137         static void recvLon(CEventServer& event);
138         static void recvLat(CEventServer& event);
139         static void recvArea(CEventServer& event);
140         void recvServerAttribut(CBufferIn& buffer);
141         void recvIndex(int rank, CBufferIn& buffer);
142         void recvLon(int rank, CBufferIn& buffer);
143         void recvLat(int rank, CBufferIn& buffer);
144         void recvArea(int rank, CBufferIn& buffer);
145
146         /// Destructeur ///
147         virtual ~CDomain(void);
148
149         /// Accesseurs statiques ///
150         static StdString GetName(void);
151         static StdString GetDefName(void);
152
153         static ENodeType GetType(void);
154         const std::map<int, vector<size_t> >& getIndexServer() const;
155         CArray<int, 2> local_mask;
156         bool isCurvilinear ;
157         bool hasBounds ;
158         bool hasArea;
159         bool hasLonLat;
160      private:
161         void checkDomain(void);
162         void checkLocalIDomain(void);
163         void checkLocalJDomain(void);
164
165         void checkMask(void);
166         void checkDomainData(void);
167         void checkCompression(void);
168
169         void checkBounds(void);
170         void checkArea(void);
171         void checkLonLat();
172
173         void checkTransformations();
174         void setTransformations(const TransMapTypes&);
175         void computeNGlobDomain();
176
177         void sendIndex();
178         void sendArea();
179         void sendLonLat();
180
181       private:
182         bool isChecked;
183         std::set<StdString> relFiles, relFilesCompressed;
184         bool isClientChecked; // Verify whether all attributes of domain on the client side are good
185         bool isClientAfterTransformationChecked;
186         std::map<int, CArray<int,1> > indiSrv, indjSrv;
187         std::map<int,int> nbConnectedClients_; // Mapping of number of communicating client to a server
188         std::map<int, vector<size_t> > indSrv_; // Global index of each client sent to server
189         std::map<int, vector<int> > indWrittenSrv_; // Global written index of each client sent to server
190         std::vector<int> indexesToWrite;
191         int numberWrittenIndexes_, totalNumberWrittenIndexes_, offsetWrittenIndexes_;
192         std::vector<int> connectedServerRank_;
193         bool isDistributed_;
194         //! True if and only if the data defined on the domain can be outputted in a compressed way
195         bool isCompressible_;
196         bool isRedistributed_;
197         TransMapTypes transformationMap_;
198         std::vector<int> nGlobDomain_;
199         bool isUnstructed_;
200
201         DECLARE_REF_FUNC(Domain,domain)
202
203   }; // class CDomain
204
205   ///--------------------------------------------------------------
206
207   // Declare/Define CDomainGroup and CDomainDefinition
208   DECLARE_GROUP(CDomain);
209
210   ///--------------------------------------------------------------
211
212} // namespace xios
213
214#endif // __XIOS_CDomain__
Note: See TracBrowser for help on using the repository browser.