source: XIOS/dev/dev_trunk_omp/extern/src_ep_dev/ep_recv.cpp @ 1646

Last change on this file since 1646 was 1646, checked in by yushan, 5 years ago

branch merged with trunk @1645. arch file (ep&mpi) added for ADA

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