Changeset 1146 for XIOS/dev/branch_yushan_merged/extern/remap/src

Timestamp:

05/24/17 16:59:40 (7 years ago)

Author:

yushan

Message:

save modif

Location:

XIOS/dev/branch_yushan_merged/extern/remap/src

Files:

• clipper.cpp (modified) (2 diffs)
• mapper.cpp (modified) (4 diffs)
• timerRemap.cpp (modified) (2 diffs)
• timerRemap.hpp (modified) (1 diff)

XIOS/dev/branch_yushan_merged/extern/remap/src/clipper.cpp

@@ -1047 +1047 @@
 
   //create a new edge array ...
-  TEdge *edges = new TEdge [highI +1];
+  TEdge *edges = new TEdge[highI+1];
 
   bool IsFlat = true;
@@ -4274 +4274 @@
 {
   //The equation of a line in general form (Ax + By + C = 0)
-  //given 2 points (x¹,y¹) & (x²,y²) is ...
-  //(y¹ - y²)x + (x² - x¹)y + (y² - y¹)x¹ - (x² - x¹)y¹ = 0
-  //A = (y¹ - y²); B = (x² - x¹); C = (y² - y¹)x¹ - (x² - x¹)y¹
-  //perpendicular distance of point (x³,y³) = (Ax³ + By³ + C)/Sqrt(A² + B²)
+  //given 2 points (x,y) & (x,y) is ...
+  //(y - y)x + (x - x)y + (y - y)x - (x - x)y = 0
+  //A = (y - y); B = (x - x); C = (y - y)x - (x - x)y
+  //perpendicular distance of point (x,y) = (Ax + By + C)/Sqrt(A + B)
   //see http://en.wikipedia.org/wiki/Perpendicular_distance
   double A = double(ln1.Y - ln2.Y);

XIOS/dev/branch_yushan_merged/extern/remap/src/mapper.cpp

@@ -548 +548 @@
     }
 
+    MPI_Waitall(nbSendRequest, sendRequest, status);
     MPI_Waitall(nbRecvRequest, recvRequest, status);
-    MPI_Waitall(nbSendRequest, sendRequest, status);
 
     for (int rank = 0; rank < mpiSize; rank++)
@@ -620 +620 @@
     }
 
+    MPI_Waitall(nbSendRequest, sendRequest, status);
     MPI_Waitall(nbRecvRequest, recvRequest, status);
-    MPI_Waitall(nbSendRequest, sendRequest, status);
-
+ 
     int nbNeighbourNodes = 0;
     for (int rank = 0; rank < mpiSize; rank++)
@@ -803 +803 @@
         }
     }
-
+ 
+    MPI_Waitall(nbSendRequest, sendRequest, status);
     MPI_Waitall(nbRecvRequest, recvRequest, status);
-    MPI_Waitall(nbSendRequest, sendRequest, status);
+    
     char **sendBuffer2 = new char*[mpiSize];
     char **recvBuffer2 = new char*[mpiSize];
@@ -886 +887 @@
         }
     }
-
+    
+    MPI_Waitall(nbSendRequest, sendRequest, status);
     MPI_Waitall(nbRecvRequest, recvRequest, status);
-    MPI_Waitall(nbSendRequest, sendRequest, status);
+   
 
     delete [] sendRequest;

XIOS/dev/branch_yushan_merged/extern/remap/src/timerRemap.cpp

@@ -9 +9 @@
 using namespace std;
 
-map<string,CTimer*> CTimer::allTimer;
+map<string,CTimer*> *CTimer::allTimer = 0;
 
 CTimer::CTimer(const string& name_) : name(name_)
@@ -55 +55 @@
 CTimer& CTimer::get(const string name)
 {
+        if(allTimer == 0) allTimer = new map<string,CTimer*>;
         map<string,CTimer*>::iterator it;
-        it=allTimer.find(name);
-        if (it==allTimer.end()) it=allTimer.insert(pair<string,CTimer*>(name,new CTimer(name))).first;
+        it=(*allTimer).find(name);
+        if (it==(*allTimer).end()) it=(*allTimer).insert(pair<string,CTimer*>(name,new CTimer(name))).first;
         return *(it->second);
 }

XIOS/dev/branch_yushan_merged/extern/remap/src/timerRemap.hpp

@@ -26 +26 @@
     double getCumulatedTime(void);
     void print(void);
-    static map<string,CTimer*> allTimer;
+    //static map<string,CTimer*> allTimer;
+    static map<string,CTimer*> *allTimer;
+    #pragma omp threadprivate(allTimer)
+    
     static double getTime(void);
     static CTimer& get(string name);