/* -*- 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-- */

}