source: XIOS/dev/branch_openmp/extern/ep_dev/ep_recv.cpp @ 1503

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

rank_map is passed from vector to map, in order to have more flexibility in comm_split

File size: 4.7 KB
Line 
1/*!
2   \file ep_recv.cpp
3   \since 2 may 2016
4
5   \brief Definitions of MPI receive functions: MPI_Recv, MPI_Mrecv, MPI_Irecv, MPI_Imrecv
6 */
7
8
9#include "ep_lib.hpp"
10#include <mpi.h>
11#include "ep_declaration.hpp"
12#include "ep_mpi.hpp"
13
14using namespace std;
15
16extern std::list< ep_lib::MPI_Request* > * EP_PendingRequests;
17#pragma omp threadprivate(EP_PendingRequests)
18
19namespace ep_lib
20{
21
22  int MPI_Recv(void *buf, int count, MPI_Datatype datatype, int src, int tag, MPI_Comm comm, MPI_Status *status)
23  {
24
25    if(!comm->is_ep) return MPI_Recv_mpi(buf, count, datatype, src, tag, comm, status);
26   
27    Debug("MPI_Recv with EP");
28
29    MPI_Request request;
30    MPI_Irecv(buf, count, datatype, src, tag, comm, &request);
31    MPI_Wait(&request, status);
32
33    return 0;
34  }
35
36 
37
38
39  int MPI_Irecv(void *buf, int count, MPI_Datatype datatype, int src, int tag, MPI_Comm comm, MPI_Request *request)
40  {
41    if(!comm->is_ep) return MPI_Irecv_mpi(buf, count, datatype, src, tag, comm, request);
42   
43    Debug("MPI_Irecv with EP");
44    int dest_rank;
45    MPI_Comm_rank(comm, &dest_rank);
46   
47    *request = new ep_request;   
48    memcheck("new " << *request <<" : in ep_lib::MPI_Irecv, *request = new ep_request");
49
50    (*request)->mpi_request = new ::MPI_Request;
51    memcheck("new " << (*request)->mpi_request << " : in ep_lib::MPI_Irecv, (*request)->mpi_request = new ::MPI_Request");
52   
53    (*request)->buf = buf;
54    (*request)->comm = comm;
55    (*request)->type = 2;
56    (*request)->state = 0;
57   
58
59    (*request)->ep_src = src;
60    (*request)->ep_tag = tag;
61    (*request)->ep_datatype = datatype;
62 
63    if(EP_PendingRequests == 0 ) EP_PendingRequests = new std::list< MPI_Request* >;
64
65    EP_PendingRequests->push_back(request); 
66    int ep_rank = comm->ep_comm_ptr->size_rank_info[0].first;
67    memcheck("EP_PendingRequests["<<ep_rank<<"]->size() = " << EP_PendingRequests->size());                                                             
68                                         
69    return Request_Check();
70  }
71 
72  int MPI_Mrecv(void *buf, int count, MPI_Datatype datatype, MPI_Message *message, MPI_Status *status)
73  {
74    Debug("MPI_Mrecv with MPI/EP");
75
76    status->mpi_status = new ::MPI_Status;
77    memcheck("new " << status->mpi_status << " : in ep_lib::MPI_Mrecv, status->mpi_status = new ::MPI_Status");
78   
79    ::MPI_Mrecv(buf, count, to_mpi_type(datatype), static_cast< ::MPI_Message* >((*message)->mpi_message), to_mpi_status_ptr(*status));
80
81   
82    status->ep_src = (*message)->ep_src;
83    status->ep_datatype = datatype;
84    status->ep_tag = (*message)->ep_tag;
85
86    memcheck("delete " << (*message)->mpi_message << " : in ep_lib::MPI_Mrecv, delete (*message)->mpi_message");
87    delete (*message)->mpi_message;
88
89#ifdef _check_sum
90    check_sum_recv(buf, count, datatype, message->ep_src, message->ep_tag);
91#endif
92
93    return Request_Check();
94  }
95
96
97  int MPI_Imrecv(void *buf, int count, MPI_Datatype datatype, MPI_Message *message, MPI_Request *request)
98  {
99    Debug("MPI_Imrecv with MPI/EP");
100
101    (*request)->type = 3;
102    (*request)->ep_datatype = datatype;
103    (*request)->ep_tag = (*message)->ep_tag;
104    (*request)->ep_src = (*message)->ep_src;
105       
106    ::MPI_Imrecv(buf, count, to_mpi_type(datatype), to_mpi_message_ptr(*message), to_mpi_request_ptr(*request));               
107   
108    memcheck("delete " << (*message)->mpi_message << " : in ep_lib::MPI_Imrecv, delete (*message)->mpi_message");
109    delete (*message)->mpi_message;
110
111#ifdef _check_sum
112    check_sum_recv(buf, count, datatype, message->ep_src, message->ep_tag);
113#endif
114
115   
116    return Request_Check();
117  }
118
119
120   int MPI_Recv_mpi(void *buf, int count, MPI_Datatype datatype, int src, int tag, MPI_Comm comm, MPI_Status *status)
121  {
122    Debug("MPI_Recv with MPI");
123    status->ep_src = src;
124    status->ep_tag = tag;
125    status->ep_datatype = datatype;
126   
127    return ::MPI_Recv(buf, count, to_mpi_type(datatype), src<0? MPI_ANY_SOURCE : src, tag<0? MPI_ANY_TAG: tag, to_mpi_comm(comm->mpi_comm), to_mpi_status_ptr(*status)); 
128  }
129 
130  int MPI_Irecv_mpi(void *buf, int count, MPI_Datatype datatype, int src, int tag, MPI_Comm comm, MPI_Request *request)
131  {
132    Debug("MPI_Irecv with MPI");
133    int dest_rank;
134    MPI_Comm_rank(comm, &dest_rank);
135   
136    *request = new ep_request;
137    memcheck("new "<< *request <<" : in ep_lib::MPI_Irecv, *request = new ep_request");
138
139    (*request)->mpi_request = new ::MPI_Request;
140    memcheck("new "<< (*request)->mpi_request <<" : in ep_lib::MPI_Irecv, (*request)->mpi_request = new ::MPI_Request");
141     
142    (*request)->ep_src = src;
143    (*request)->ep_datatype = datatype;
144    (*request)->type = 2;
145    (*request)->ep_tag = tag;
146   
147    return ::MPI_Irecv(buf, count, to_mpi_type(datatype), src<0? MPI_ANY_SOURCE : src, tag<0? MPI_ANY_TAG: tag, to_mpi_comm(comm->mpi_comm), to_mpi_request_ptr(*request));
148  }
149}
150
151
152
Note: See TracBrowser for help on using the repository browser.