1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
|
/*
* Copyright (c) 2004-2007 The Trustees of Indiana University and Indiana
* University Research and Technology
* Corporation. All rights reserved.
* Copyright (c) 2004-2020 The University of Tennessee and The University
* of Tennessee Research Foundation. All rights
* reserved.
* Copyright (c) 2004-2005 High Performance Computing Center Stuttgart,
* University of Stuttgart. All rights reserved.
* Copyright (c) 2004-2005 The Regents of the University of California.
* All rights reserved.
* Copyright (c) 2006 University of Houston. All rights reserved.
* Copyright (c) 2006-2012 Cisco Systems, Inc. All rights reserved.
* Copyright (c) 2015 Research Organization for Information Science
* and Technology (RIST). All rights reserved.
* $COPYRIGHT$
*
* Additional copyrights may follow
*
* $HEADER$
*/
#include "ompi_config.h"
#include <stdio.h>
#include "ompi/mpi/c/bindings.h"
#include "ompi/runtime/params.h"
#include "ompi/group/group.h"
#include "ompi/errhandler/errhandler.h"
#include "ompi/communicator/communicator.h"
#if OMPI_BUILD_MPI_PROFILING
#if OPAL_HAVE_WEAK_SYMBOLS
#pragma weak MPI_Group_range_incl = PMPI_Group_range_incl
#endif
#define MPI_Group_range_incl PMPI_Group_range_incl
#endif
static const char FUNC_NAME[] = "MPI_Group_range_incl";
int MPI_Group_range_incl(MPI_Group group, int n_triplets, int ranges[][3],
MPI_Group *new_group)
{
int err, i,indx;
int group_size;
int * elements_int_list;
/* can't act on NULL group */
if( MPI_PARAM_CHECK ) {
OMPI_ERR_INIT_FINALIZE(FUNC_NAME);
if ( (MPI_GROUP_NULL == group) || (NULL == group) ||
(NULL == new_group) ) {
return OMPI_ERRHANDLER_NOHANDLE_INVOKE(MPI_ERR_GROUP,
FUNC_NAME);
}
group_size = ompi_group_size ( group);
elements_int_list = (int *) malloc(sizeof(int) * (group_size+1));
if (NULL == elements_int_list) {
return OMPI_ERRHANDLER_NOHANDLE_INVOKE(MPI_ERR_OTHER, FUNC_NAME);
}
for (i = 0; i < group_size; i++) {
elements_int_list[i] = -1;
}
for ( i=0; i < n_triplets; i++) {
if ((0 > ranges[i][0]) || (ranges[i][0] > group_size)) {
goto error_rank;
}
if ((0 > ranges[i][1]) || (ranges[i][1] > group_size)) {
goto error_rank;
}
if (ranges[i][2] == 0) {
goto error_rank;
}
if ((ranges[i][0] < ranges[i][1])) {
if (ranges[i][2] < 0) {
goto error_rank;
}
/* positive stride */
for (indx = ranges[i][0]; indx <= ranges[i][1]; indx += ranges[i][2]) {
/* make sure rank has not already been selected */
if (elements_int_list[indx] != -1) {
goto error_rank;
}
elements_int_list[indx] = i;
}
} else if (ranges[i][0] > ranges[i][1]) {
if (ranges[i][2] > 0) {
goto error_rank;
}
/* negative stride */
for (indx = ranges[i][0]; indx >= ranges[i][1]; indx += ranges[i][2]) {
/* make sure rank has not already been selected */
if (elements_int_list[indx] != -1) {
goto error_rank;
}
elements_int_list[indx] = i;
}
} else {
/* first_rank == last_rank */
indx = ranges[i][0];
if (elements_int_list[indx] != -1) {
goto error_rank;
}
elements_int_list[indx] = i;
}
}
free ( elements_int_list);
}
err = ompi_group_range_incl ( group, n_triplets, ranges, new_group );
OMPI_ERRHANDLER_NOHANDLE_RETURN(err, err, FUNC_NAME );
error_rank:
free(elements_int_list);
return OMPI_ERRHANDLER_NOHANDLE_INVOKE(MPI_ERR_RANK, FUNC_NAME);
}
|
