-/* Copyright (c) 2009-2021. The SimGrid Team. All rights reserved. */
+/* Copyright (c) 2009-2022. The SimGrid Team. All rights reserved. */
/* This program is free software; you can redistribute it and/or modify it
* under the terms of the license (GNU LGPL) which comes with this package. */
((b).value) *= ((a).value); \
((b).index) *= ((a).index); \
}
-#define LAND_OP(a, b) (b) = (a) && (b)
-#define LOR_OP(a, b) (b) = (a) || (b)
-#define LXOR_OP(a, b) (b) = bool(a) != bool(b)
+#define LAND_OP(a, b) (b) = static_cast<std::remove_reference_t<decltype(b)>>((a) && (b))
+#define LOR_OP(a, b) (b) = static_cast<std::remove_reference_t<decltype(b)>>((a) || (b))
+#define LXOR_OP(a, b) (b) = static_cast<std::remove_reference_t<decltype(b)>>(bool(a) != bool(b))
#define BAND_OP(a, b) (b) &= (a)
#define BOR_OP(a, b) (b) |= (a)
#define BXOR_OP(a, b) (b) ^= (a)
} \
}
+#define APPLY_BEGIN_OP_LOOP() \
+ MPI_Datatype datatype_base = *datatype; \
+ while (datatype_base->duplicated_datatype() != MPI_DATATYPE_NULL) \
+ datatype_base = datatype_base->duplicated_datatype();
+
#define APPLY_OP_LOOP(dtype, type, op) \
- if (*datatype == (dtype)) { \
+ if (datatype_base == (dtype)) { \
APPLY_FUNC(a, b, length, type, op) \
} else
static void max_func(void *a, void *b, int *length, MPI_Datatype * datatype)
{
+ APPLY_BEGIN_OP_LOOP()
APPLY_BASIC_OP_LOOP(MAX_OP)
APPLY_FLOAT_OP_LOOP(MAX_OP)
APPLY_END_OP_LOOP(MAX_OP)
static void min_func(void *a, void *b, int *length, MPI_Datatype * datatype)
{
+ APPLY_BEGIN_OP_LOOP()
APPLY_BASIC_OP_LOOP(MIN_OP)
APPLY_FLOAT_OP_LOOP(MIN_OP)
APPLY_END_OP_LOOP(MIN_OP)
static void sum_func(void *a, void *b, int *length, MPI_Datatype * datatype)
{
+ APPLY_BEGIN_OP_LOOP()
APPLY_BASIC_OP_LOOP(SUM_OP)
APPLY_FLOAT_OP_LOOP(SUM_OP)
APPLY_COMPLEX_OP_LOOP(SUM_OP)
static void prod_func(void *a, void *b, int *length, MPI_Datatype * datatype)
{
+ APPLY_BEGIN_OP_LOOP()
APPLY_BASIC_OP_LOOP(PROD_OP)
APPLY_FLOAT_OP_LOOP(PROD_OP)
APPLY_COMPLEX_OP_LOOP(PROD_OP)
static void land_func(void *a, void *b, int *length, MPI_Datatype * datatype)
{
+ APPLY_BEGIN_OP_LOOP()
APPLY_BASIC_OP_LOOP(LAND_OP)
APPLY_FLOAT_OP_LOOP(LAND_OP)
APPLY_BOOL_OP_LOOP(LAND_OP)
static void lor_func(void *a, void *b, int *length, MPI_Datatype * datatype)
{
+ APPLY_BEGIN_OP_LOOP()
APPLY_BASIC_OP_LOOP(LOR_OP)
APPLY_FLOAT_OP_LOOP(LOR_OP)
APPLY_BOOL_OP_LOOP(LOR_OP)
static void lxor_func(void *a, void *b, int *length, MPI_Datatype * datatype)
{
+ APPLY_BEGIN_OP_LOOP()
APPLY_BASIC_OP_LOOP(LXOR_OP)
APPLY_FLOAT_OP_LOOP(LXOR_OP)
APPLY_BOOL_OP_LOOP(LXOR_OP)
static void band_func(void *a, void *b, int *length, MPI_Datatype * datatype)
{
+ APPLY_BEGIN_OP_LOOP()
APPLY_BASIC_OP_LOOP(BAND_OP)
APPLY_BOOL_OP_LOOP(BAND_OP)
APPLY_BYTE_OP_LOOP(BAND_OP)
static void bor_func(void *a, void *b, int *length, MPI_Datatype * datatype)
{
+ APPLY_BEGIN_OP_LOOP()
APPLY_BASIC_OP_LOOP(BOR_OP)
APPLY_BOOL_OP_LOOP(BOR_OP)
APPLY_BYTE_OP_LOOP(BOR_OP)
static void bxor_func(void *a, void *b, int *length, MPI_Datatype * datatype)
{
+ APPLY_BEGIN_OP_LOOP()
APPLY_BASIC_OP_LOOP(BXOR_OP)
APPLY_BOOL_OP_LOOP(BXOR_OP)
APPLY_BYTE_OP_LOOP(BXOR_OP)
static void minloc_func(void *a, void *b, int *length, MPI_Datatype * datatype)
{
+ APPLY_BEGIN_OP_LOOP()
APPLY_PAIR_OP_LOOP(MINLOC_OP)
APPLY_END_OP_LOOP(MINLOC_OP)
}
static void maxloc_func(void *a, void *b, int *length, MPI_Datatype * datatype)
{
+ APPLY_BEGIN_OP_LOOP()
APPLY_PAIR_OP_LOOP(MAXLOC_OP)
APPLY_END_OP_LOOP(MAXLOC_OP)
}
/* obviously a no-op */
}
-
-#define CREATE_MPI_OP(name, func, types) \
- SMPI_Op _XBT_CONCAT(smpi_MPI_, name)(&(func) /* func */, true, true, types);
+#define CREATE_MPI_OP(name, func, types) \
+ SMPI_Op _XBT_CONCAT(smpi_MPI_, name)(&(func), true, true, types, _XBT_STRINGIFY(MPI_##name));
#define MAX_TYPES DT_FLAG_C_INTEGER|DT_FLAG_F_INTEGER|DT_FLAG_FP|DT_FLAG_MULTILANG
#define LAND_TYPES DT_FLAG_C_INTEGER|DT_FLAG_FP|DT_FLAG_LOGICAL|DT_FLAG_MULTILANG
CREATE_MPI_OP(REPLACE, replace_func, 0)
CREATE_MPI_OP(NO_OP, no_func, 0)
-namespace simgrid{
-namespace smpi{
+namespace simgrid::smpi {
void Op::apply(const void* invec, void* inoutvec, const int* len, MPI_Datatype datatype) const
{
- if (smpi_cfg_privatization() == SmpiPrivStrategies::MMAP) {
- // we need to switch as the called function may silently touch global variables
- XBT_DEBUG("Applying operation, switch to the right data frame ");
- smpi_switch_data_segment(simgrid::s4u::Actor::self());
- }
+ // we need to switch as the called function may silently touch global variables
+ smpi_switch_data_segment(simgrid::s4u::Actor::self());
if (not smpi_process()->replaying() && *len > 0) {
+ XBT_DEBUG("Applying operation of length %d from %p and from/to %p", *len, invec, inoutvec);
if (not is_fortran_op_)
this->func_(const_cast<void*>(invec), inoutvec, const_cast<int*>(len), &datatype);
else{
- XBT_DEBUG("Applying operation of length %d from %p and from/to %p", *len, invec, inoutvec);
int tmp = datatype->c2f();
/* Unfortunately, the C and Fortran version of the MPI standard do not agree on the type here,
thus the reinterpret_cast. */
}
}
-}
-}
+} // namespace simgrid::smpi