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

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

tests in XIOS OK (client, complete, remap, toy)

File size: 5.6 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    if(comm->is_intercomm)
44    {
45      if(src>=0) src = comm->inter_rank_map->at(src);
46    }
47   
48    Debug("MPI_Irecv with EP");
49    int dest_rank;
50    MPI_Comm_rank(comm, &dest_rank);
51   
52    *request = new ep_request;   
53    memcheck("new " << *request <<" : in ep_lib::MPI_Irecv, *request = new ep_request");
54
55    (*request)->mpi_request = new ::MPI_Request;
56    memcheck("new " << (*request)->mpi_request << " : in ep_lib::MPI_Irecv, (*request)->mpi_request = new ::MPI_Request");
57   
58    (*request)->buf = buf;
59    (*request)->comm = comm;
60    (*request)->type = 2;
61    (*request)->probed = false;
62    (*request)->state = 0;
63
64    (*request)->ep_src = src;
65    (*request)->ep_tag = tag;
66    (*request)->ep_datatype = datatype;
67 
68    if(EP_PendingRequests == 0 ) EP_PendingRequests = new std::list< MPI_Request* >;
69
70    EP_PendingRequests->push_back(request); 
71   
72    memcheck("EP_PendingRequests["<<ep_rank<<"]->size() = " << EP_PendingRequests->size());   
73   
74    show_EP_PendingRequests(EP_PendingRequests);
75
76#ifdef _showinfo
77    if(comm->is_intercomm)
78    {
79      int ep_dest_loc  = comm->ep_rank_map->at(dest_rank).first;
80      int ep_src_loc = comm->ep_comm_ptr->intercomm->intercomm_rank_map->at(src).first;
81      int mpi_tag     = tag_combine(tag, ep_src_loc, ep_dest_loc);
82      int mpi_dest    = comm->ep_comm_ptr->intercomm->intercomm_rank_map->at(src).second.first;
83
84      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);
85    }                                                         
86#endif                                     
87
88    return Request_Check();
89  }
90 
91  int MPI_Mrecv(void *buf, int count, MPI_Datatype datatype, MPI_Message *message, MPI_Status *status)
92  {
93    Debug("MPI_Mrecv with MPI/EP");
94
95    status->mpi_status = new ::MPI_Status;
96    memcheck("new " << status->mpi_status << " : in ep_lib::MPI_Mrecv, status->mpi_status = new ::MPI_Status");
97   
98    ::MPI_Mrecv(buf, count, to_mpi_type(datatype), static_cast< ::MPI_Message* >((*message)->mpi_message), to_mpi_status_ptr(*status));
99
100   
101    status->ep_src = (*message)->ep_src;
102    status->ep_datatype = datatype;
103    status->ep_tag = (*message)->ep_tag;
104
105    memcheck("delete " << (*message)->mpi_message << " : in ep_lib::MPI_Mrecv, delete (*message)->mpi_message");
106    delete (*message)->mpi_message;
107    memcheck("delete " << *message << " : in ep_lib::MPI_Imrecv, delete *message");
108    delete *message;
109
110#ifdef _check_sum
111    check_sum_recv(buf, count, datatype, message->ep_src, message->ep_tag);
112#endif
113
114    return Request_Check();
115  }
116
117
118  int MPI_Imrecv(void *buf, int count, MPI_Datatype datatype, MPI_Message *message, MPI_Request *request)
119  {
120    Debug("MPI_Imrecv with MPI/EP");
121
122    (*request)->type = 3;
123    (*request)->ep_datatype = datatype;
124    (*request)->ep_tag = (*message)->ep_tag;
125    (*request)->ep_src = (*message)->ep_src;
126   
127    (*request)->probed = true;
128    (*request)->state = 1;
129       
130    ::MPI_Imrecv(buf, count, to_mpi_type(datatype), to_mpi_message_ptr(*message), to_mpi_request_ptr(*request));               
131   
132    memcheck("delete " << (*message)->mpi_message << " : in ep_lib::MPI_Imrecv, delete (*message)->mpi_message");
133    delete (*message)->mpi_message;
134
135
136#ifdef _check_sum
137    check_sum_recv(buf, count, datatype, message->ep_src, message->ep_tag);
138#endif
139
140   
141    //return Request_Check();
142  }
143
144
145   int MPI_Recv_mpi(void *buf, int count, MPI_Datatype datatype, int src, int tag, MPI_Comm comm, MPI_Status *status)
146  {
147    Debug("MPI_Recv with MPI");
148    status->ep_src = src;
149    status->ep_tag = tag;
150    status->ep_datatype = datatype;
151   
152    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)); 
153  }
154 
155  int MPI_Irecv_mpi(void *buf, int count, MPI_Datatype datatype, int src, int tag, MPI_Comm comm, MPI_Request *request)
156  {
157    Debug("MPI_Irecv with MPI");
158    int dest_rank;
159    MPI_Comm_rank(comm, &dest_rank);
160   
161    *request = new ep_request;
162    memcheck("new "<< *request <<" : in ep_lib::MPI_Irecv, *request = new ep_request");
163
164    (*request)->mpi_request = new ::MPI_Request;
165    memcheck("new "<< (*request)->mpi_request <<" : in ep_lib::MPI_Irecv, (*request)->mpi_request = new ::MPI_Request");
166     
167    (*request)->ep_src = src;
168    (*request)->ep_datatype = datatype;
169    (*request)->type = 2;
170    (*request)->ep_tag = tag;
171   
172    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));
173  }
174}
175
176
177
Note: See TracBrowser for help on using the repository browser.