Changeset 1087 for XIOS/dev/branch_yushan/extern

Timestamp:

04/12/17 11:55:38 (7 years ago)

Author:

yushan

Message:

save modif from local. context_init/finalize OK

Location:

XIOS/dev/branch_yushan/extern/src_ep_dev

Files:

• ep_fortran.cpp (modified) (3 diffs)
• ep_type.hpp (modified) (1 diff)

XIOS/dev/branch_yushan/extern/src_ep_dev/ep_fortran.cpp

@@ -25 +25 @@
     std::map<std::pair<int, int>, MPI_Comm > ::iterator it;
     
-    it = fc_comm_map.find(std::make_pair(fint, omp_get_thread_num()));
-    if(it == fc_comm_map.end())
+    #pragma omp critical (fc_comm_map)
     {
-      fc_comm_map.insert(std::make_pair( std::make_pair( fint, omp_get_thread_num()) , comm));
-      printf("EP_Comm_c2f : MAP insert: %d, %d, %p\n", fint, omp_get_thread_num(), comm.ep_comm_ptr);
+      it = fc_comm_map.find(std::make_pair(fint, omp_get_thread_num()));
+      if(it == fc_comm_map.end())
+      {
+        fc_comm_map.insert(std::make_pair( std::make_pair( fint, omp_get_thread_num()) , comm));
+        printf("EP_Comm_c2f : MAP insert: %d, %d, %p\n", fint, omp_get_thread_num(), comm.ep_comm_ptr);
+      }
     }
     
@@ -44 +47 @@
     std::map<std::pair<int, int>, MPI_Comm > ::iterator it;
     
+    #pragma omp critical (fc_comm_map)
     it = fc_comm_map.find(std::make_pair(comm, omp_get_thread_num()));
+    
     if(it != fc_comm_map.end())
     {
@@ -52 +57 @@
       return  comm_ptr;
     }
-    else
-    {      
+       
       
-      #ifdef _openmpi
-      ::MPI_Comm base_comm = ::MPI_Comm_f2c(comm);
-      #elif _intelmpi
-      ::MPI_Comm base_comm = (::MPI_Comm)(comm);
-      #endif
+    #ifdef _openmpi
+    ::MPI_Comm base_comm = ::MPI_Comm_f2c(comm);
+    #elif _intelmpi
+    ::MPI_Comm base_comm = (::MPI_Comm)(comm);
+    #endif
 
-      if(base_comm != MPI_COMM_NULL_STD)
+    if(base_comm != MPI_COMM_NULL_STD)
+    {
+      if(omp_get_thread_num() == 0)
       {
-        if(omp_get_thread_num() == 0)
-        {
-          int num_ep = omp_get_num_threads();
-          MPI_Comm *new_comm;
-          MPI_Info info;
-          MPI_Comm_create_endpoints(base_comm, num_ep, info, new_comm);
-          passage = new_comm;
-        }
-        #pragma omp barrier
+        int num_ep = omp_get_num_threads();
+        MPI_Comm *new_comm;
+        MPI_Info info;
+        MPI_Comm_create_endpoints(base_comm, num_ep, info, new_comm);
+        passage = new_comm;
+      }
+      #pragma omp barrier
 
-        MPI_Comm return_comm = passage[omp_get_thread_num()];
-        return return_comm;
+      MPI_Comm return_comm = passage[omp_get_thread_num()];
+      return return_comm;
         
-      }
-      return MPI_COMM_NULL;
-      
     }
+    return MPI_COMM_NULL;
+
   }
-
-  // #ifdef _intelmpi
-
-  // MPI_Fint MPI_Comm_c2f(MPI_Comm comm)
-  // {
-  //   Debug("MPI_Comm_c2f");
-  //   int fint;
-  //   fint = (::MPI_Fint)(comm.mpi_comm);
-    
-  //   std::map<std::pair<int, int>, MPI_Comm  > ::iterator it;
-    
-  //   it = fc_comm_map.find(std::make_pair(fint, omp_get_thread_num()));
-  //   if(it == fc_comm_map.end())
-  //   {
-  //     fc_comm_map.insert(std::make_pair( std::make_pair( fint, omp_get_thread_num()) , comm));
-  //     printf("MAP insert: %d, %d, %p\n", fint, omp_get_thread_num(), &comm);
-  //   }
-    
-  //   MPI_Fint Fint;
-  //   Fint.mpi_fint = fint;
-  //   return Fint;
-    
-  // }
-
-  
-
-
-  // MPI_Comm MPI_Comm_f2c(MPI_Fint comm)
-  // {
-  //   Debug("MPI_Comm_f2c");
-    
-    
-  //   std::map<std::pair<int, int>, MPI_Comm  > ::iterator it;
-    
-  //   it = fc_comm_map.find(std::make_pair(comm.mpi_fint, omp_get_thread_num()));
-  //   if(it != fc_comm_map.end())
-  //   {
-  //     MPI_Comm comm_ptr;
-  //     comm_ptr =  it->second;
-  //     printf("MAP find: %d, %d, %p\n", it->first.first, it->first.second, &comm_ptr);
-  //     return  comm_ptr;
-  //   }
-  //   else
-  //   {      
-  //     MPI_Comm return_comm;
-  //     return_comm.mpi_comm = (::MPI_Comm)(comm.mpi_fint);
-  //     return return_comm;
-  //   }
-  // }
-
-  
-
-  // #elif _openmpi
-  
-  // int MPI_Comm_c2f(MPI_Comm comm)
-  // {
-  //   Debug("MPI_Comm_c2f");
-  //   int fint;
-  //   fint = ::MPI_Comm_c2f(static_cast< ::MPI_Comm>(comm.mpi_comm));
-    
-  //   std::map<std::pair<int, int>, MPI_Comm > ::iterator it;
-    
-  //   it = fc_comm_map.find(std::make_pair(fint, omp_get_thread_num()));
-  //   if(it == fc_comm_map.end())
-  //   {
-  //     fc_comm_map.insert(std::make_pair( std::make_pair( fint, omp_get_thread_num()) , comm));
-  //     printf("MAP insert: %d, %d, %p\n", fint, omp_get_thread_num(), &comm);
-  //   }
-    
-  //   return fint;
-    
-  // }
-
-  // ep_lib::MPI_Comm MPI_Comm_f2c(MPI_Fint comm)
-  // {
-  //   Debug("MPI_Comm_f2c");
-    
-    
-  //   std::map<std::pair<int, int>, MPI_Comm > ::iterator it;
-    
-  //   it = fc_comm_map.find(std::make_pair(comm, omp_get_thread_num()));
-  //   if(it != fc_comm_map.end())
-  //   {
-  //     MPI_Comm comm_ptr;
-  //     comm_ptr =  it->second;
-  //     printf("MAP find: %d, %d, %p\n", it->first.first, it->first.second, &comm_ptr);
-  //     return  comm_ptr;
-  //   }
-  //   else
-  //   {      
-  //     MPI_Comm return_comm;
-  //     return_comm.mpi_comm = (::MPI_Comm)(comm);
-  //     return return_comm;
-  //   }
-  // }
-  // #endif
 
 }

XIOS/dev/branch_yushan/extern/src_ep_dev/ep_type.hpp

@@ -485 +485 @@
 
   static std::map<std::pair<int, int>, MPI_Comm >  fc_comm_map;
+
+  static std::map<std::pair<int, int>, MPI_Comm >  *fc_comm_map_ptr;
+  #pragma omp threadprivate(fc_comm_map_ptr)
             //    <MPI_Fint,thread_num>   EP_Comm