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-2005 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 Cisco Systems, Inc. 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 OPAL_HAVE_WEAK_SYMBOLS && OMPI_PROFILING_DEFINES
#pragma weak MPI_Group_range_excl = PMPI_Group_range_excl
#endif
#if OMPI_PROFILING_DEFINES
#include "ompi/mpi/c/profile/defines.h"
#endif
static const char FUNC_NAME[] = "MPI_Group_range_excl";
int MPI_Group_range_excl(MPI_Group group, int n_triplets, int ranges[][3],
MPI_Group *new_group)
{
int err, i, group_size, index;
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_INVOKE(MPI_COMM_WORLD, 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_INVOKE(MPI_COMM_WORLD, 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 (index = ranges[i][0]; index <= ranges[i][1]; index += ranges[i][2]) {
/* make sure rank has not already been selected */
if (elements_int_list[index] != -1) {
goto error_rank;
}
elements_int_list[index] = i;
}
} else if (ranges[i][0] > ranges[i][1]) {
if (ranges[i][2] > 0) {
goto error_rank;
}
/* negative stride */
for (index = ranges[i][0]; index >= ranges[i][1]; index += ranges[i][2]) {
/* make sure rank has not already been selected */
if (elements_int_list[index] != -1) {
goto error_rank;
}
elements_int_list[index] = i;
}
} else {
/* first_rank == last_rank */
index = ranges[i][0];
if (elements_int_list[index] != -1) {
goto error_rank;
}
elements_int_list[index] = i;
}
}
free (elements_int_list);
}
OPAL_CR_ENTER_LIBRARY();
err = ompi_group_range_excl(group,n_triplets,ranges,new_group);
OMPI_ERRHANDLER_RETURN(err, MPI_COMM_WORLD,err,FUNC_NAME);
error_rank:
free(elements_int_list);
return OMPI_ERRHANDLER_INVOKE(MPI_COMM_WORLD, MPI_ERR_RANK, FUNC_NAME);
}
|
