Skip to content

Instantly share code, notes, and snippets.

@jsquyres
Created February 28, 2015 19:23
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save jsquyres/496e741bd00196ef2d64 to your computer and use it in GitHub Desktop.
Save jsquyres/496e741bd00196ef2d64 to your computer and use it in GitHub Desktop.
Issue #330 allreduce_f.f90 replacement
/*
* Copyright (c) 2004-2005 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) 2011-2012 Cisco Systems, Inc. All rights reserved.
* $COPYRIGHT$
*
* Additional copyrights may follow
*
* $HEADER$
*/
#include "ompi_config.h"
#include "ompi/mpi/fortran/mpif-h/bindings.h"
#include "ompi/mpi/fortran/base/constants.h"
#if OPAL_HAVE_WEAK_SYMBOLS && OMPI_PROFILE_LAYER
#pragma weak PMPI_ALLREDUCE = ompi_allreduce_f
#pragma weak pmpi_allreduce = ompi_allreduce_f
#pragma weak pmpi_allreduce_ = ompi_allreduce_f
#pragma weak pmpi_allreduce__ = ompi_allreduce_f
#pragma weak PMPI_Allreduce_f = ompi_allreduce_f
#pragma weak PMPI_Allreduce_f08 = ompi_allreduce_f
#elif OMPI_PROFILE_LAYER
OMPI_GENERATE_F77_BINDINGS (PMPI_ALLREDUCE,
pmpi_allreduce,
pmpi_allreduce_,
pmpi_allreduce__,
pompi_allreduce_f,
(char *sendbuf, char *recvbuf, MPI_Fint *count, MPI_Fint *datatype, MPI_Fint *op, MPI_Fint *comm, MPI_Fint *ierr),
(sendbuf, recvbuf, count, datatype, op, comm, ierr) )
#endif
#if OPAL_HAVE_WEAK_SYMBOLS
#pragma weak MPI_ALLREDUCE = ompi_allreduce_f
#pragma weak mpi_allreduce = ompi_allreduce_f
#pragma weak mpi_allreduce_ = ompi_allreduce_f
#pragma weak mpi_allreduce__ = ompi_allreduce_f
#pragma weak MPI_Allreduce_f = ompi_allreduce_f
#pragma weak MPI_Allreduce_f08 = ompi_allreduce_f
#endif
#if ! OPAL_HAVE_WEAK_SYMBOLS && ! OMPI_PROFILE_LAYER
OMPI_GENERATE_F77_BINDINGS (MPI_ALLREDUCE,
mpi_allreduce,
mpi_allreduce_,
mpi_allreduce__,
ompi_allreduce_f,
(char *sendbuf, char *recvbuf, MPI_Fint *count, MPI_Fint *datatype, MPI_Fint *op, MPI_Fint *comm, MPI_Fint *ierr),
(sendbuf, recvbuf, count, datatype, op, comm, ierr) )
#endif
#if OMPI_PROFILE_LAYER && ! OPAL_HAVE_WEAK_SYMBOLS
#include "ompi/mpi/fortran/mpif-h/profile/defines.h"
#endif
void ompi_allreduce_f(char *sendbuf, char *recvbuf, MPI_Fint *count,
MPI_Fint *datatype, MPI_Fint *op, MPI_Fint *comm,
MPI_Fint *ierr)
{
int ierr_c;
MPI_Comm c_comm;
MPI_Datatype c_type;
MPI_Op c_op;
c_comm = MPI_Comm_f2c(*comm);
c_type = MPI_Type_f2c(*datatype);
c_op = MPI_Op_f2c(*op);
{
void *tmp = (void*) sendbuf;
#if OPAL_HAVE_WEAK_SYMBOLS
opal_output(0, "Have weak symbols: %p %p %p %p",
(void*) &MPI_FORTRAN_IN_PLACE,
(void*) &mpi_fortran_in_place,
(void*) &mpi_fortran_in_place_,
(void*) &mpi_fortran_in_place__);
#elif OMPI_FORTRAN_CAPS
opal_output(0, "Have CAPS: %p", (void*) &MPI_FORTRAN_IN_PLACE);
#elif OMPI_FORTRAN_PLAIN
opal_output(0, "Have CAPS: %p", (void*) &mpi_fortran_in_place);
#elif OMPI_FORTRAN_SINGLE_UNDERSCORE
opal_output(0, "Have CAPS: %p", (void*) &mpi_fortran_in_place_);
#elif OMPI_FORTRAN_DOUBLE_UNDERSCORE
opal_output(0, "Have CAPS: %p", (void*) &mpi_fortran_in_place__);
#else
#error Impossible!
#endif
opal_output(0, "Fortran MPI_ALLREDUCE: got input buffer %p", tmp);
tmp = OMPI_F2C_IN_PLACE(tmp);
opal_output(0, "Fortran MPI_ALLREDUCE 2: after IN_PLACE %p (MPI_IN_PLACE = %p)", tmp, (void*) MPI_IN_PLACE);
tmp = OMPI_F2C_BOTTOM(tmp);
opal_output(0, "Fortran MPI_ALLREDUCE 3: after IN_PLACE %p (MPI_IN_PLACE = %p)", tmp, (void*) MPI_IN_PLACE);
}
sendbuf = (char *) OMPI_F2C_IN_PLACE(sendbuf);
sendbuf = (char *) OMPI_F2C_BOTTOM(sendbuf);
recvbuf = (char *) OMPI_F2C_BOTTOM(recvbuf);
ierr_c = MPI_Allreduce(sendbuf, recvbuf,
OMPI_FINT_2_INT(*count),
c_type, c_op, c_comm);
if (NULL != ierr) *ierr = OMPI_INT_2_FINT(ierr_c);
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment