c976559877
The log basic bcast was completely broken. The rank 0 gets the hibit set to -1, so it always returned an error.
200 строки
5.9 KiB
C
200 строки
5.9 KiB
C
/*
|
|
* Copyright (c) 2004-2005 The Trustees of Indiana University and Indiana
|
|
* University Research and Technology
|
|
* Corporation. All rights reserved.
|
|
* Copyright (c) 2004-2016 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) 2015 Cisco Systems, Inc. All rights reserved.
|
|
* Copyright (c) 2016 Research Organization for Information Science
|
|
* and Technology (RIST). All rights reserved.
|
|
* $COPYRIGHT$
|
|
*
|
|
* Additional copyrights may follow
|
|
*
|
|
* $HEADER$
|
|
*/
|
|
|
|
#include "ompi_config.h"
|
|
#include "coll_basic.h"
|
|
|
|
#include "mpi.h"
|
|
#include "ompi/constants.h"
|
|
#include "ompi/datatype/ompi_datatype.h"
|
|
#include "ompi/mca/coll/coll.h"
|
|
#include "ompi/mca/coll/base/coll_tags.h"
|
|
#include "coll_basic.h"
|
|
#include "ompi/mca/pml/pml.h"
|
|
#include "opal/util/bit_ops.h"
|
|
|
|
|
|
/*
|
|
* bcast_log_intra
|
|
*
|
|
* Function: - broadcast using O(log(N)) algorithm
|
|
* Accepts: - same arguments as MPI_Bcast()
|
|
* Returns: - MPI_SUCCESS or error code
|
|
*/
|
|
int
|
|
mca_coll_basic_bcast_log_intra(void *buff, int count,
|
|
struct ompi_datatype_t *datatype, int root,
|
|
struct ompi_communicator_t *comm,
|
|
mca_coll_base_module_t *module)
|
|
{
|
|
int i;
|
|
int size;
|
|
int rank;
|
|
int vrank;
|
|
int peer;
|
|
int dim;
|
|
int hibit;
|
|
int mask;
|
|
int err;
|
|
int nreqs;
|
|
ompi_request_t **preq, **reqs;
|
|
|
|
size = ompi_comm_size(comm);
|
|
rank = ompi_comm_rank(comm);
|
|
vrank = (rank + size - root) % size;
|
|
|
|
dim = comm->c_cube_dim;
|
|
hibit = opal_hibit(vrank, dim);
|
|
--dim;
|
|
|
|
/* Receive data from parent in the tree. */
|
|
|
|
if (vrank > 0) {
|
|
assert(hibit >= 0);
|
|
peer = ((vrank & ~(1 << hibit)) + root) % size;
|
|
|
|
err = MCA_PML_CALL(recv(buff, count, datatype, peer,
|
|
MCA_COLL_BASE_TAG_BCAST,
|
|
comm, MPI_STATUS_IGNORE));
|
|
if (MPI_SUCCESS != err) {
|
|
return err;
|
|
}
|
|
}
|
|
|
|
/* Send data to the children. */
|
|
|
|
reqs = coll_base_comm_get_reqs(module->base_data, size);
|
|
if( NULL == reqs ) { return OMPI_ERR_OUT_OF_RESOURCE; }
|
|
|
|
err = MPI_SUCCESS;
|
|
preq = reqs;
|
|
nreqs = 0;
|
|
for (i = hibit + 1, mask = 1 << i; i <= dim; ++i, mask <<= 1) {
|
|
peer = vrank | mask;
|
|
if (peer < size) {
|
|
peer = (peer + root) % size;
|
|
++nreqs;
|
|
|
|
err = MCA_PML_CALL(isend(buff, count, datatype, peer,
|
|
MCA_COLL_BASE_TAG_BCAST,
|
|
MCA_PML_BASE_SEND_STANDARD,
|
|
comm, preq++));
|
|
if (MPI_SUCCESS != err) {
|
|
ompi_coll_base_free_reqs(reqs, nreqs);
|
|
return err;
|
|
}
|
|
}
|
|
}
|
|
|
|
/* Start and wait on all requests. */
|
|
|
|
if (nreqs > 0) {
|
|
|
|
/* Wait for them all. If there's an error, note that we don't
|
|
* care what the error was -- just that there *was* an error.
|
|
* The PML will finish all requests, even if one or more of them
|
|
* fail. i.e., by the end of this call, all the requests are
|
|
* free-able. So free them anyway -- even if there was an
|
|
* error, and return the error after we free everything. */
|
|
|
|
err = ompi_request_wait_all(nreqs, reqs, MPI_STATUSES_IGNORE);
|
|
if( MPI_SUCCESS != err ) {
|
|
ompi_coll_base_free_reqs(reqs, nreqs);
|
|
}
|
|
}
|
|
|
|
/* All done */
|
|
|
|
return err;
|
|
}
|
|
|
|
|
|
/*
|
|
* bcast_lin_inter
|
|
*
|
|
* Function: - broadcast using O(N) algorithm
|
|
* Accepts: - same arguments as MPI_Bcast()
|
|
* Returns: - MPI_SUCCESS or error code
|
|
*/
|
|
int
|
|
mca_coll_basic_bcast_lin_inter(void *buff, int count,
|
|
struct ompi_datatype_t *datatype, int root,
|
|
struct ompi_communicator_t *comm,
|
|
mca_coll_base_module_t *module)
|
|
{
|
|
int i;
|
|
int rsize;
|
|
int err;
|
|
ompi_request_t **reqs = NULL;
|
|
|
|
rsize = ompi_comm_remote_size(comm);
|
|
|
|
if (MPI_PROC_NULL == root) {
|
|
/* do nothing */
|
|
err = OMPI_SUCCESS;
|
|
} else if (MPI_ROOT != root) {
|
|
/* Non-root receive the data. */
|
|
err = MCA_PML_CALL(recv(buff, count, datatype, root,
|
|
MCA_COLL_BASE_TAG_BCAST, comm,
|
|
MPI_STATUS_IGNORE));
|
|
} else {
|
|
reqs = coll_base_comm_get_reqs(module->base_data, rsize);
|
|
if( NULL == reqs ) { return OMPI_ERR_OUT_OF_RESOURCE; }
|
|
|
|
/* root section */
|
|
for (i = 0; i < rsize; i++) {
|
|
err = MCA_PML_CALL(isend(buff, count, datatype, i,
|
|
MCA_COLL_BASE_TAG_BCAST,
|
|
MCA_PML_BASE_SEND_STANDARD,
|
|
comm, &(reqs[i])));
|
|
if (OMPI_SUCCESS != err) {
|
|
ompi_coll_base_free_reqs(reqs, i + 1);
|
|
return err;
|
|
}
|
|
}
|
|
err = ompi_request_wait_all(rsize, reqs, MPI_STATUSES_IGNORE);
|
|
if (OMPI_SUCCESS != err) {
|
|
ompi_coll_base_free_reqs(reqs, rsize);
|
|
}
|
|
}
|
|
|
|
|
|
/* All done */
|
|
return err;
|
|
}
|
|
|
|
|
|
/*
|
|
* bcast_log_inter
|
|
*
|
|
* Function: - broadcast using O(N) algorithm
|
|
* Accepts: - same arguments as MPI_Bcast()
|
|
* Returns: - MPI_SUCCESS or error code
|
|
*/
|
|
int
|
|
mca_coll_basic_bcast_log_inter(void *buff, int count,
|
|
struct ompi_datatype_t *datatype, int root,
|
|
struct ompi_communicator_t *comm,
|
|
mca_coll_base_module_t *module)
|
|
{
|
|
return OMPI_ERR_NOT_IMPLEMENTED;
|
|
}
|