source: XIOS/dev/branch_openmp/extern/src_ep_dev/ep_recv.cpp @ 1220

Last change on this file since 1220 was 1220, checked in by yushan, 4 years ago

test_remap_omp tested on ADA except two fields

File size: 3.4 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
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 
21
22  int MPI_Recv(void *buf, int count, MPI_Datatype datatype, int src, int tag, MPI_Comm comm, MPI_Status *status)
23  {
24    int dest_rank;
25    MPI_Comm_rank(comm, &dest_rank);
26
27    if(!comm.is_ep)
28    {
29      ::MPI_Comm mpi_comm = static_cast< ::MPI_Comm >(comm.mpi_comm);
30      ::MPI_Status mpi_status;
31      ::MPI_Recv(buf, count, static_cast< ::MPI_Datatype >(datatype), src<0? MPI_ANY_SOURCE : src, tag<0? MPI_ANY_TAG: tag, mpi_comm, &mpi_status);
32
33      status->ep_src = src;
34      status->ep_tag = tag;
35      status->ep_datatype = datatype;
36
37      return 0; 
38    }
39
40    Message_Check(comm);
41
42    MPI_Request request;
43    MPI_Irecv(buf, count, datatype, src, tag, comm, &request);
44    MPI_Wait(&request, status);
45
46    return 0;
47  }
48
49
50
51
52  int MPI_Irecv(void *buf, int count, MPI_Datatype datatype, int src, int tag, MPI_Comm comm, MPI_Request *request)
53  {
54
55    Debug("MPI_Irecv with EP");
56    int dest_rank;
57    MPI_Comm_rank(comm, &dest_rank);
58    int world_rank;
59    MPI_Comm_rank(MPI_COMM_WORLD_STD, &world_rank);
60
61    if(!comm.is_ep)
62    {
63      ::MPI_Request mpi_request;
64      ::MPI_Comm mpi_comm = static_cast< ::MPI_Comm > (comm.mpi_comm);
65      ::MPI_Irecv(buf, count, static_cast< ::MPI_Datatype> (datatype), src<0? MPI_ANY_SOURCE : src, tag<0? MPI_ANY_TAG: tag, mpi_comm, &mpi_request);
66
67      request->mpi_request = mpi_request;
68      request->ep_src = src;
69      request->ep_datatype = datatype;
70      request->ep_tag = tag;
71    }
72
73    request->mpi_request = MPI_REQUEST_NULL_STD;
74    request->buf = buf;
75    request->comm = comm;
76    request->type = 2;
77
78    request->ep_src = src;
79    request->ep_tag = tag;
80    request->ep_datatype = datatype;
81
82
83
84    /* With Improbe*/
85    Message_Check(comm);
86
87    if(EP_PendingRequests == 0 ) 
88    {
89      EP_PendingRequests = new std::list< MPI_Request* >;
90    }
91
92
93    EP_PendingRequests->push_back(request);
94
95    Request_Check();
96
97    return 0;
98  }
99
100  int MPI_Imrecv(void *buf, int count, MPI_Datatype datatype, MPI_Message *message, MPI_Request *request)
101  {
102    Debug("MPI_Imrecv");
103
104    request->type = 3;
105
106    ::MPI_Request mpi_request;
107    ::MPI_Message mpi_message = static_cast< ::MPI_Message >(message->mpi_message);
108               
109    ::MPI_Imrecv(buf, count, static_cast< ::MPI_Datatype>(datatype), &mpi_message, &mpi_request);
110
111    request->mpi_request = mpi_request;
112    request->ep_datatype = datatype;
113    request->ep_tag = message->ep_tag;
114    request->ep_src = message->ep_src;
115
116    return 0;
117  }
118
119
120  int MPI_Mrecv(void *buf, int count, MPI_Datatype datatype, MPI_Message *message, MPI_Status *status)
121  {
122    Debug("EP Mrecv called\n");
123
124    ::MPI_Status mpi_status;
125    ::MPI_Message mpi_message = static_cast< ::MPI_Message >(message->mpi_message);
126   
127    ::MPI_Mrecv(buf, count, static_cast< ::MPI_Datatype>(datatype), &mpi_message, &mpi_status);
128
129    status->mpi_status = new ::MPI_Status(mpi_status);
130    status->ep_src = message->ep_src;
131    status->ep_datatype = datatype;
132    status->ep_tag = message->ep_tag;
133
134    //check_sum_recv(buf, count, datatype, message->ep_src, message->ep_tag);
135
136    return 0;
137  }
138
139}
140
141
142
Note: See TracBrowser for help on using the repository browser.