/* -*- Mode: C; c-basic-offset:4 ; indent-tabs-mode:nil ; -*- */
/*
*
* (C) 2001 by Argonne National Laboratory.
* See COPYRIGHT in top-level directory.
*/
#include "mpiimpl.h"
/* -- Begin Profiling Symbol Block for routine MPI_Op_create */
#if defined(HAVE_PRAGMA_WEAK)
#pragma weak MPI_Op_create = PMPI_Op_create
#elif defined(HAVE_PRAGMA_HP_SEC_DEF)
#pragma _HP_SECONDARY_DEF PMPI_Op_create MPI_Op_create
#elif defined(HAVE_PRAGMA_CRI_DUP)
#pragma _CRI duplicate MPI_Op_create as PMPI_Op_create
#elif defined(HAVE_WEAK_ATTRIBUTE)
int MPI_Op_create(MPI_User_function * user_fn, int commute, MPI_Op * op)
__attribute__ ((weak, alias("PMPI_Op_create")));
#endif
/* -- End Profiling Symbol Block */
/* Define MPICH_MPI_FROM_PMPI if weak symbols are not supported to build
the MPI routines */
#ifndef MPICH_MPI_FROM_PMPI
#undef MPI_Op_create
#define MPI_Op_create PMPI_Op_create
#ifndef MPIR_OP_PREALLOC
#define MPIR_OP_PREALLOC 16
#endif
/* Preallocated op objects */
MPIR_Op MPIR_Op_builtin[MPIR_OP_N_BUILTIN] = { {0} };
MPIR_Op MPIR_Op_direct[MPIR_OP_PREALLOC] = { {0} };
MPIR_Object_alloc_t MPIR_Op_mem = { 0, 0, 0, 0, MPIR_OP,
sizeof(MPIR_Op),
MPIR_Op_direct,
MPIR_OP_PREALLOC,
};
#ifdef HAVE_CXX_BINDING
void MPII_Op_set_cxx(MPI_Op op, void (*opcall) (void))
{
MPIR_Op *op_ptr;
MPIR_Op_get_ptr(op, op_ptr);
op_ptr->language = MPIR_LANG__CXX;
MPIR_Process.cxx_call_op_fn = (void (*)(const void *, void *, int,
MPI_Datatype, MPI_User_function *)) opcall;
}
#endif
#if defined(HAVE_FORTRAN_BINDING) && !defined(HAVE_FINT_IS_INT)
/* Normally, the C and Fortran versions are the same, by design in the
MPI Standard. However, if MPI_Fint and int are not the same size (e.g.,
MPI_Fint was made 8 bytes but int is 4 bytes), then the C and Fortran
versions must be distinquished. */
void MPII_Op_set_fc(MPI_Op op)
{
MPIR_Op *op_ptr;
MPIR_Op_get_ptr(op, op_ptr);
op_ptr->language = MPIR_LANG__FORTRAN;
}
#endif
#undef FUNCNAME
#define FUNCNAME MPIR_Op_create_impl
#undef FCNAME
#define FCNAME MPL_QUOTE(FUNCNAME)
int MPIR_Op_create_impl(MPI_User_function * user_fn, int commute, MPI_Op * op)
{
MPIR_Op *op_ptr;
int mpi_errno = MPI_SUCCESS;
op_ptr = (MPIR_Op *) MPIR_Handle_obj_alloc(&MPIR_Op_mem);
/* --BEGIN ERROR HANDLING-- */
if (!op_ptr) {
mpi_errno =
MPIR_Err_create_code(MPI_SUCCESS, MPIR_ERR_RECOVERABLE, FCNAME, __LINE__, MPI_ERR_OTHER,
"**nomem", "**nomem %s", "MPI_Op");
goto fn_fail;
}
/* --END ERROR HANDLING-- */
op_ptr->language = MPIR_LANG__C;
op_ptr->kind = commute ? MPIR_OP_KIND__USER : MPIR_OP_KIND__USER_NONCOMMUTE;
op_ptr->function.c_function = (void (*)(const void *, void *,
const int *, const MPI_Datatype *)) user_fn;
MPIR_Object_set_ref(op_ptr, 1);
MPIR_OBJ_PUBLISH_HANDLE(*op, op_ptr->handle);
#ifdef MPID_Op_commit_hook
MPID_Op_commit_hook(op_ptr);
#endif
fn_exit:
return mpi_errno;
fn_fail:
goto fn_exit;
}
#endif /* MPICH_MPI_FROM_PMPI */
#undef FUNCNAME
#define FUNCNAME MPI_Op_create
#undef FCNAME
#define FCNAME MPL_QUOTE(FUNCNAME)
/*@
MPI_Op_create - Creates a user-defined combination function handle
Input Parameters:
+ user_fn - user defined function (function)
- commute - true if commutative; false otherwise. (logical)
Output Parameters:
. op - operation (handle)
Notes on the user function:
The calling list for the user function type is
.vb
typedef void (MPI_User_function) (void * a,
void * b, int * len, MPI_Datatype *);
.ve
where the operation is 'b[i] = a[i] op b[i]', for 'i=0,...,len-1'. A pointer
to the datatype given to the MPI collective computation routine (i.e.,
'MPI_Reduce', 'MPI_Allreduce', 'MPI_Scan', or 'MPI_Reduce_scatter') is also
passed to the user-specified routine.
.N ThreadSafe
.N Fortran
.N collops
.N Errors
.N MPI_SUCCESS
.seealso: MPI_Op_free
@*/
int MPI_Op_create(MPI_User_function * user_fn, int commute, MPI_Op * op)
{
int mpi_errno = MPI_SUCCESS;
MPIR_FUNC_TERSE_STATE_DECL(MPID_STATE_MPI_OP_CREATE);
MPIR_ERRTEST_INITIALIZED_ORDIE();
MPID_THREAD_CS_ENTER(GLOBAL, MPIR_THREAD_GLOBAL_ALLFUNC_MUTEX);
MPIR_FUNC_TERSE_ENTER(MPID_STATE_MPI_OP_CREATE);
/* ... body of routine ... */
mpi_errno = MPIR_Op_create_impl(user_fn, commute, op);
if (mpi_errno)
goto fn_fail;
/* ... end of body of routine ... */
fn_exit:
MPIR_FUNC_TERSE_EXIT(MPID_STATE_MPI_OP_CREATE);
MPID_THREAD_CS_EXIT(GLOBAL, MPIR_THREAD_GLOBAL_ALLFUNC_MUTEX);
return mpi_errno;
fn_fail:
/* --BEGIN ERROR HANDLING-- */
#ifdef HAVE_ERROR_CHECKING
{
mpi_errno =
MPIR_Err_create_code(mpi_errno, MPIR_ERR_RECOVERABLE, FCNAME, __LINE__, MPI_ERR_OTHER,
"**mpi_op_create", "**mpi_op_create %p %d %p", user_fn, commute,
op);
}
#endif
mpi_errno = MPIR_Err_return_comm(0, FCNAME, mpi_errno);
goto fn_exit;
/* --END ERROR HANDLING-- */
}