source: XIOS/dev/branch_openmp/extern/src_ep_dev/ep_message.cpp @ 1520

Last change on this file since 1520 was 1520, checked in by yushan, 3 years ago

save dev. TO DO : test with xios

File size: 7.4 KB
Line 
1/*!
2   \file ep_message.cpp
3   \since 2 may 2016
4
5   \brief Definitions of MPI endpoint function: Message_Check
6 */
7
8#include "ep_lib.hpp"
9#include <mpi.h>
10#include "ep_declaration.hpp"
11#include "ep_mpi.hpp"
12
13using namespace std;
14
15extern std::list< ep_lib::MPI_Request* > * EP_PendingRequests;
16#pragma omp threadprivate(EP_PendingRequests)
17
18namespace ep_lib
19{
20  int Request_Check()
21  {
22    if(EP_PendingRequests == 0 ) EP_PendingRequests = new std::list< MPI_Request* >;
23   
24    if(EP_PendingRequests->size() == 0) return 0;
25   
26    MPI_Status status;
27    MPI_Message *message;
28    int probed = false;
29    int recv_count = 0;
30    std::list<MPI_Request* >::iterator it;
31   
32   
33    for(it = EP_PendingRequests->begin(); it!=EP_PendingRequests->end(); it++)
34    { 
35      Message_Check(((*(*it))->comm));
36    }
37
38
39    for(it = EP_PendingRequests->begin(); it!=EP_PendingRequests->end(); )
40    {
41      if((*(*it))->state == 0)
42      {
43        #pragma omp critical (_query0)
44        {
45          MPI_Iprobe((*(*it))->ep_src, (*(*it))->ep_tag, ((*(*it))->comm), &probed, &status);
46          if(probed)
47          {
48            message = new MPI_Message;
49            *message = new ep_message;
50       
51            memcheck("new "<< message <<" : in ep_lib::Request_Check, message = new MPI_Message");
52            memcheck("new "<< *message <<" : in ep_lib::Request_Check, *message = new ep_message");
53         
54         
55            MPI_Improbe((*(*it))->ep_src, (*(*it))->ep_tag, (*(*it))->comm, &probed, message, &status);
56       
57          }
58        }
59     
60       
61        if(probed)
62        {
63          MPI_Get_count(&status, (*(*it))->ep_datatype, &recv_count);
64         
65          MPI_Imrecv((*(*it))->buf, recv_count, (*(*it))->ep_datatype, message, *it);
66          (*(*it))->type = 3;
67          (*(*it))->state = 1;
68
69          memcheck("delete "<< status.mpi_status <<" : in ep_lib::Request_Check, delete status.mpi_status");
70          delete status.mpi_status;         
71
72          memcheck("delete "<< *message <<" : in ep_lib::Request_Check, delete *message");
73          memcheck("delete "<< message <<" : in ep_lib::Request_Check, delete message");
74
75          delete *message;
76          delete message;
77       
78          it++;
79          continue;     
80        }             
81      }
82     
83      if((*(*it))->state == 2)
84      {
85        int ep_rank = ((*(*it))->comm)->ep_comm_ptr->size_rank_info[0].first;
86        memcheck("delete "<< (*(*it)) <<" : in ep_lib::Request_Check, delete (*(*it))");
87       
88       
89        int world_rank;
90        MPI_Comm_rank(MPI_COMM_WORLD, &world_rank);
91        if(world_rank==2) 
92        {
93          printf("ep %d erased one pending request %p\n", world_rank,*(*it));
94        }
95       
96        EP_PendingRequests->erase(it);
97       
98        memcheck("EP_PendingRequests["<<ep_rank<<"]->size() = " << EP_PendingRequests->size());
99        it = EP_PendingRequests->begin();
100        continue;
101      }
102      else it++;
103    }
104  }
105 
106 
107 
108  int Message_Check(MPI_Comm comm)
109  {
110    if(!comm->is_ep) return MPI_SUCCESS;
111
112    if(comm->is_intercomm)
113    {
114      Message_Check_intercomm(comm);
115    }
116   
117    return Message_Check_intracomm(comm);
118
119  }
120 
121 
122  int Message_Check_intracomm(MPI_Comm comm)
123  {
124   
125    int flag = true;
126    ::MPI_Status status;
127    ::MPI_Message message;
128
129    while(flag) // loop until the end of global queue
130    {
131      Debug("Message probing for intracomm\n");
132     
133      #pragma omp critical (_mpi_call)
134      {
135        ::MPI_Iprobe(MPI_ANY_SOURCE, MPI_ANY_TAG, to_mpi_comm(comm->mpi_comm), &flag, &status);
136        if(flag)
137        {
138          Debug("find message in mpi comm \n");
139          ::MPI_Mprobe(status.MPI_SOURCE, status.MPI_TAG, to_mpi_comm(comm->mpi_comm), &message, &status);
140        }
141      }
142
143     
144      if(flag)
145      {
146        MPI_Message msg = new ep_message; 
147        msg->mpi_message = new ::MPI_Message(message);
148
149        memcheck("new "<< msg <<" : in ep_lib::Message_Check, msg = new ep_message");
150        memcheck("new "<< msg->mpi_message <<" : in ep_lib::Message_Check, msg->mpi_message = new ::MPI_Message");
151             
152
153        msg->ep_tag  = bitset<15>(status.MPI_TAG >> 16).to_ulong(); 
154        int src_loc  = bitset<8> (status.MPI_TAG >> 8) .to_ulong(); 
155        int dest_loc = bitset<8> (status.MPI_TAG)           .to_ulong();
156        int src_mpi  = status.MPI_SOURCE;
157             
158        msg->ep_src  = get_ep_rank(comm, src_loc,  src_mpi); 
159
160#ifdef _showinfo
161        printf("status.MPI_TAG = %d, src_loc = %d, dest_loc = %d, ep_tag = %d\n", status.MPI_TAG, src_loc, dest_loc, msg->ep_tag);
162#endif
163
164        msg->mpi_status = new ::MPI_Status(status); 
165        memcheck("new "<< msg->mpi_status <<" : in ep_lib::Message_Check, msg->mpi_status = new ::MPI_Status");
166
167        #pragma omp critical (_query)
168        {
169          #pragma omp flush
170          comm->ep_comm_ptr->comm_list[dest_loc]->ep_comm_ptr->message_queue->push_back(msg); 
171          int dest_mpi = comm->ep_comm_ptr->size_rank_info[2].first;
172          memcheck("message_queue["<<dest_mpi<<","<<dest_loc<<"]->size = "<<comm->ep_comm_ptr->comm_list[dest_loc]->ep_comm_ptr->message_queue->size());
173          #pragma omp flush
174        }
175      }
176    }
177
178    return MPI_SUCCESS;
179  }
180
181
182 
183 
184 
185
186
187  int Message_Check_intercomm(MPI_Comm comm)
188  {
189    if(!comm->ep_comm_ptr->intercomm->mpi_inter_comm) return 0;
190
191    Debug("Message probing for intercomm\n");
192
193    int flag = true;
194    ::MPI_Message message;
195    ::MPI_Status status;
196    int current_ep_rank;
197    MPI_Comm_rank(comm, &current_ep_rank);
198
199    while(flag) // loop until the end of global queue "comm->ep_comm_ptr->intercomm->mpi_inter_comm"
200    {
201      Debug("Message probing for intracomm\n");
202
203      #pragma omp critical (_mpi_call)
204      {
205        ::MPI_Iprobe(MPI_ANY_SOURCE, MPI_ANY_TAG, to_mpi_comm(comm->ep_comm_ptr->intercomm->mpi_inter_comm), &flag, &status);
206        if(flag)
207        {
208          Debug("find message in mpi comm \n");
209          ::MPI_Mprobe(status.MPI_SOURCE, status.MPI_TAG, to_mpi_comm(comm->ep_comm_ptr->intercomm->mpi_inter_comm), &message, &status);
210        }
211      }
212     
213
214      if(flag)
215      {
216
217        MPI_Message msg = new ep_message; 
218        msg->mpi_message = new ::MPI_Message(message);
219
220        memcheck("new "<< msg <<" : in ep_lib::Message_Check, msg = new ep_message");
221        memcheck("new "<< msg->mpi_message <<" : in ep_lib::Message_Check, msg->mpi_message = new ::MPI_Message");
222             
223
224        msg->ep_tag  = bitset<15>(status.MPI_TAG >> 16).to_ulong(); 
225        int src_loc  = bitset<8> (status.MPI_TAG >> 8) .to_ulong(); 
226        int dest_loc = bitset<8> (status.MPI_TAG)           .to_ulong();
227        int src_mpi  = status.MPI_SOURCE;
228             
229        msg->ep_src  = get_ep_rank(comm, src_loc,  src_mpi);   
230#ifdef _showinfo
231        printf("status.MPI_TAG = %d, src_loc = %d, dest_loc = %d, ep_tag = %d\n", status.MPI_TAG, src_loc, dest_loc, msg->ep_tag);
232#endif
233
234        msg->mpi_status = new ::MPI_Status(status); 
235        memcheck("new "<< msg->mpi_status <<" : in ep_lib::Message_Check, msg->mpi_status = new ::MPI_Status");
236
237        #pragma omp critical (_query)
238        {
239          #pragma omp flush
240          comm->ep_comm_ptr->comm_list[dest_loc]->ep_comm_ptr->message_queue->push_back(msg);
241          memcheck("comm->ep_comm_ptr->comm_list["<<dest_loc<<"]->ep_comm_ptr->message_queue->size = "<<comm->ep_comm_ptr->comm_list[dest_loc]->ep_comm_ptr->message_queue->size());
242          #pragma omp flush
243        }
244      }
245    }
246
247    Message_Check_intracomm(comm);
248
249    return MPI_SUCCESS;
250  }
251
252 
253
254}
Note: See TracBrowser for help on using the repository browser.