source: XIOS/dev/branch_yushan_merged/extern/src_ep_dev/ep_test.cpp @ 1149

Last change on this file since 1149 was 1149, checked in by yushan, 7 years ago

bug fixed in mpi_comm_split. Key needs to be specifify.

File size: 2.8 KB
Line 
1/*!
2   \file ep_test.cpp
3   \since 2 may 2016
4
5   \brief Definitions of MPI test function: MPI_Test, MPI_Testsome, MPI_Testany, MPI_Testall
6 */
7
8#include "ep_lib.hpp"
9#include <mpi.h>
10#include "ep_declaration.hpp"
11
12using namespace std;
13
14
15namespace ep_lib {
16
17  int MPI_Test(MPI_Request *request, int *flag, MPI_Status *status)
18  {
19    Debug("MPI_Test with EP");
20
21    *flag = false;
22
23    if(request->type == 1)      // isend
24    {
25      ::MPI_Request mpi_request = static_cast< ::MPI_Request >(request->mpi_request);
26      ::MPI_Status mpi_status;
27      ::MPI_Test(&mpi_request, flag, &mpi_status);
28     
29      if(*flag) 
30      {
31        status->mpi_status = new ::MPI_Status(mpi_status);
32        status->ep_src = request->ep_src;
33        status->ep_tag = request->ep_tag;
34        status->ep_datatype = request->ep_datatype;
35      }
36
37      return 0;
38    }
39
40    if(request->type == 2)   // irecv message not probed
41    {
42      Message_Check(request->comm);
43      #pragma omp flush
44      MPI_Message message;
45      MPI_Improbe(request->ep_src, request->ep_tag, request->comm, flag, &message, status);
46      if(*flag)
47      {
48        int count;
49        MPI_Get_count(status, request->ep_datatype, &count);
50        MPI_Imrecv(request->buf, count, request->ep_datatype, &message, request);
51        MPI_Test(request, flag, status);
52      }
53      return 0;
54    }
55
56    if(request->type == 3)  // imrecv
57    {
58      ::MPI_Request *mpi_request = static_cast< ::MPI_Request* >(&(request->mpi_request));
59      ::MPI_Status mpi_status;
60     
61      ::MPI_Errhandler_set(MPI_COMM_WORLD_STD, MPI_ERRORS_RETURN);
62      int error_code = ::MPI_Test(mpi_request, flag, &mpi_status);
63      if (error_code != MPI_SUCCESS) {
64     
65         char error_string[BUFSIZ];
66         int length_of_error_string, error_class;
67     
68         ::MPI_Error_class(error_code, &error_class);
69         ::MPI_Error_string(error_class, error_string, &length_of_error_string);
70         printf("%s\n", error_string);
71      }
72     
73      if(*flag)
74      {
75        status->mpi_status = new ::MPI_Status(mpi_status);
76        status->ep_src = request->ep_src;
77        status->ep_tag = request->ep_tag;
78        status->ep_datatype = request->ep_datatype;
79        //int count;
80        //MPI_Get_count(status, request->ep_datatype, &count);
81        //check_sum_recv(request->buf, count, request->ep_datatype, request->ep_src, request->ep_tag, request->comm, 2);
82      }
83
84      status->ep_src = request->ep_src;
85      status->ep_tag = request->ep_tag;
86      status->ep_datatype = request->ep_datatype;
87
88     
89
90      return 0;
91    }
92               
93  }
94
95
96  int MPI_Testall(int count, MPI_Request *array_of_requests, int *flag, MPI_Status *array_of_statuses)
97  {
98    Debug("MPI_Testall with EP");
99    *flag = true;
100    int i=0;
101    while(*flag && i<count )
102    {
103      MPI_Test(&array_of_requests[i], flag, &array_of_statuses[i]);
104      i++;
105    }
106  }
107
108
109}
Note: See TracBrowser for help on using the repository browser.