Changeset 715 for XIOS/trunk/src/node/domain.cpp

Timestamp:

10/05/15 14:43:48 (9 years ago)

Author:

mhnguyen

Message:

Temporary modification of auto-generate domain. This should be changed to a better solution

+) Divide a rectangular domain into parts being equal to number of procs

Test
+) On local with 2 procs, mode server and attached.
+) test_remap passes and result is correct

File:

• XIOS/trunk/src/node/domain.cpp (modified) (6 diffs)

XIOS/trunk/src/node/domain.cpp

@@ -35 +35 @@
       , isClientAfterTransformationChecked(false), hasLonLat(false)
       , lonvalue_client(), latvalue_client(), bounds_lon_client(), bounds_lat_client()
-      , srcObject_(0)
+      , srcObject_(0), isRedistributed_(false)
    { /* Ne rien faire de plus */ }
 
@@ -45 +45 @@
       , isClientAfterTransformationChecked(false), hasLonLat(false)
       , lonvalue_client(), latvalue_client(), bounds_lon_client(), bounds_lat_client()
-      , srcObject_(0)
+      , srcObject_(0), isRedistributed_(false)
    { /* Ne rien faire de plus */ }
 
@@ -171 +171 @@
    void CDomain::redistribute(int nbLocalDomain)
    {
+     if (this->isRedistributed_) return;
      if (type_attr::rectilinear == type)
      {
@@ -270 +271 @@
         // Now fill other attributes
         fillInRectilinearLonLat();
-//        fillInRectilinearBoundLonLat();
+        this->isRedistributed_ = true;
+        info <<"now, we are here " << std::endl;
+        info << "domain " << this->getId() << " ni " << ni.getValue() << " nj " << nj.getValue() << std::endl;
      }
    }
@@ -327 +330 @@
        }
    }
+
+   /*!
+     Temporary function to verify whether a rectilinear domain is created automatically.
+   The domain is distributed into number of parts which are equal to number of clients (intracomm)
+   */
+   void CDomain::checkGenerate()
+   {
+     TransMapTypes trans = this->getAllTransformations();
+     TransMapTypes::const_iterator it = trans.begin(), ite = trans.end();
+     int transOrder = 0;
+     for (; it != ite; ++it, ++transOrder)
+     {
+       ETranformationType transType = it->first;
+       if ((TRANS_GENERATE_RECTILINEAR_DOMAIN == transType) && (0 == transOrder))
+       {
+         CContext* context = CContext::getCurrent();
+         CContextClient* client = context->client;
+         int nbClient;
+         MPI_Comm_size(client->intraComm,&nbClient);
+         this->redistribute(nbClient);
+         break;
+       }
+     }
+   }
+
 
    void CDomain::checkDomain(void)
@@ -942 +970 @@
      CContext* context=CContext::getCurrent();
 
+      this->checkGenerate();
       this->checkDomain();
       this->checkBounds();