SHMEM_ALLTOALL64(3)
π¦ SHMEM_ALLTOALL64(3) β the Open MPI glow-up β¨
okay bestie, sit down, hydrate, because this is the collective routine where every single PE exchanges a fixed amount of data with literally all the other PEs in the Active set. no cap. the whole girl group is trading fits at once π
it's shmem_alltoall32(3), shmem_alltoall64(3), shmem_alltoalls32(3), and shmem_alltoalls64(3), all living under the Open MPI roof π
β¨ SYNOPSIS β¨
C or C++:
#include <mpp/shmem.h>
void shmem_alltoall32(void *target, const void *source,
size_t nelems, int PE_start, int logPE_stride,
int PE_size, long *pSync)
void shmem_alltoall64(void *target, const void *source,
size_t nelems, int PE_start, int logPE_stride,
int PE_size, long *pSync)
void shmem_alltoalls32(void *target, const void *source,
ptrdiff_t dst, ptrdiff_t sst,
size_t nelems, int PE_start, int logPE_stride,
int PE_size, long *pSync)
void shmem_alltoalls64(void *target, const void *source,
ptrdiff_t dst, ptrdiff_t sst,
size_t nelems, int PE_start, int logPE_stride,
int PE_size, long *pSync)
Fortran:
INCLUDE "mpp/shmem.fh"
INTEGER dst, sst, nelems, PE_root, PE_start, logPE_stride, PE_size
INTEGER pSync(SHMEM_ALLTOALL_SYNC_SIZE)
CALL SHMEM_ALLTOALL32(target, source, nelems,
& PE_start, logPE_stride, PE_size, fIpSync)
CALL SHMEM_ALLTOALL64(target, source, nelems,
& PE_start, logPE_stride, PE_size, pSync)
CALL SHMEM_ALLTOALLS32(target, source, dst, sst, nelems,
& PE_start, logPE_stride, PE_size, pSync)
CALL SHMEM_ALLTOALLS64(target, source, dst, sst, nelems,
& PE_start, logPE_stride, PE_size, pSync)
π DESCRIPTION π
so. the shmem_alltoalls routines are collective routines. every PE in the Active set swaps nelems strided data elements β each one sized 32 bits (that's shmem_alltoalls32) or 64 bits (that's shmem_alltoalls64) β with all the other PEs in the set. everyone eats, nobody is left on read π
the tea on the strides: both strides, dst and sst, must be greater than or equal to 1. and here's the placement rule, write it on your mirror: the sst*jth block sent from PE i to PE j gets placed in the dst*ith block of the dest data object on PE j. iconic choreography, no notes π
and like every OpenSHMEM collective routine ever, these assume that only PEs in the Active set call the routine. if a PE that is not in the Active set calls an OpenSHMEM collective routine, you get undefined behavior. periodt, don't test it π
the values of the arguments dst, sst, nelems, PE_start, logPE_stride, and PE_size must be equal on all PEs in the Active set. the same dest and source data objects, and the same pSync work array, must be passed to all PEs in the Active set. matching outfits only. π
before any PE calls a shmem_alltoalls routine, these conditions must already exist (and yes, synchronization via a barrier or some other method is often needed to make that true):
- the
pSyncarray on all PEs in the Active set is not still in use from a prior call to ashmem_alltoallsroutine - the
destdata object on all PEs in the Active set is ready to accept theshmem_alltoallsdata
upon return from a shmem_alltoalls routine, here's what's true for the local PE:
- its
destsymmetric data object is completely updated and the data has been copied out of thesourcedata object - the values in the
pSyncarray are restored to the original values
the arguments are as follows:
- a symmetric data object with one of the following data types:
shmem_alltoall32: any noncharacter type that has an element size of 32 bits. no Fortran derived types or C/C++ structures are allowed.shmem_alltoall64: any noncharacter type that has an element size of 64 bits. no Fortran derived types or C/C++ structures are allowed.
targetβ a symmetric data object big enough to receive the combined total ofnelemselements from each PE in the Active set. it's giving spacious πsourceβ a symmetric data object holdingnelemselements of data for each PE in the Active set, ordered according to destination PE.dstβ the stride between consecutive elements of thedestdata object. the stride is scaled by the element size. a value of1means contiguous data.dstmust be of typeptrdiff_t. if you're in Fortran, it must be a default integer value.sstβ the stride between consecutive elements of thesourcedata object. the stride is scaled by the element size. a value of1means contiguous data.sstmust be of typeptrdiff_t. if you're in Fortran, it must be a default integer value.nelemsβ the number of elements to exchange for each PE.nelemsmust be of typesize_tfor C/C++. if you're in Fortran, it must be a default integer value.PE_startβ the lowest virtual PE number of the active set of PEs.PE_startmust be of typeinteger. if you're in Fortran, it must be a default integer value.logPE_strideβ the log (base 2) of the stride between consecutive virtual PE numbers in the active set.log_PE_stridemust be of typeinteger. if you're in Fortran, it must be a default integer value.PE_sizeβ the number of PEs in the active set.PE_sizemust be of typeinteger. if you're in Fortran, it must be a default integer value.pSyncβ a symmetric work array. in C/C++,pSyncmust be of typelongand size_SHMEM_ALLTOALL_SYNC_SIZE. in Fortran,pSyncmust be of typeintegerand sizeSHMEM_ALLTOALL_SYNC_SIZE. every element of this array must be initialized with the value_SHMEM_SYNC_VALUE(in C/C++) orSHMEM_SYNC_VALUE(in Fortran) before any of the PEs in the active set entershmem_barrier().
upon return from a shmem_alltoalls routine, the following is true for the local PE: its dest symmetric data object is completely updated and the data has been copied out of the source data object. the values in the pSync array are restored to the original values. (yes bestie, this is the vibe again, it's that important π)
the values of arguments PE_root, PE_start, logPE_stride, and PE_size must be equal on all PEs in the active set. the same target and source data objects and the same pSync work array must be passed to all PEs in the active set.
before any PE calls a alltoall routine, you must make sure the following conditions exist (synchronization via a barrier or some other method is often needed to ensure this):
- the
pSyncarray on all PEs in the active set is not still in use from a prior call to a alltoall routine - the
targetarray on all PEs in the active set is ready to accept the alltoall data
πΈ NOTES πΈ
the terms collective and symmetric are defined in intro_shmem(3). go read her, she's the blueprint π
all SHMEM alltoall routines restore pSync to its original contents. multiple calls to SHMEM routines that use the same pSync array do not require that pSync be reinitialized after the first call. so efficient, so slay β¨
β οΈ you must ensure the that the pSync array is not being updated by any PE in the active set while any of the PEs participates in processing of a SHMEM broadcast routine. be careful to avoid these situations, bestie, this is the part where people get got π
- if the
pSyncarray is initialized at run time, some type of synchronization is needed to ensure that all PEs in the working set have initializedpSyncbefore any of them enter a SHMEM routine called with thepSyncsynchronization array. - a
pSyncarray may be reused on a subsequent SHMEM broadcast routine only if none of the PEs in the active set are still processing a prior SHMEM alltoall routine call that used the samepSyncarray. - in general, this can be ensured only by doing some type of synchronization. however, in the special case of SHMEM routines being called with the same active set, you can allocate two
pSyncarrays and alternate between them on successive calls. two arrays, twin energy, we love a rotation π
π§ EXAMPLES π§
C/C++ example:
#include <shmem.h>
#include <stdio.h>
long pSync[SHMEM_ALLTOALL_SYNC_SIZE];
int main(void)
{
int64_t *source, *dest;
int i, count, pe;
shmem_init();
count = 2;
dest = (int64_t*) shmem_malloc(count * shmem_n_pes() * sizeof(int64_t));
source = (int64_t*) shmem_malloc(count * shmem_n_pes() * sizeof(int64_t));
/* assign source values */
for (pe=0; pe <shmem_n_pes(); pe++){
for (i=0; i<count; i++){
source[(pe*count)+i] = shmem_my_pe() + pe;
dest[(pe*count)+i] = 9999;
}
}
for (i=0; i< SHMEM_ALLTOALLS_SYNC_SIZE; i++) {
pSync[i] = SHMEM_SYNC_VALUE;
}
/* wait for all PEs to initialize pSync */
shmem_barrier_all();
/* alltoalls on all PES */
shmem_alltoalls64(dest, source, 1, 1, count, 0, 0, shmem_n_pes(), pSync);
/* verify results */
for (pe=0; pe<shmem_n_pes(); pe++) {
for (i=0; i<count; i++){
if (dest[(pe*count)+i] != shmem_my_pe() + pe) {
printf("[%d] ERROR: dest[%d]=%ld, should be %d,
shmem_my_pe(),(pe*count)+i,dest[(pe*count)+i],
shmem_n_pes() + pe);
}
}
}
shmem_barrier_all();
shmem_free(dest);
shmem_free(source);
shmem_finalize();
return 0;
}
See also:
π Copyright π
Copyright 2003-2026, The Open MPI Community
and that's the whole page, bestie π«Ά this was SHMEM_ALLTOALL64(3) from Open MPI, dated Sep 10, 2026 β copyright 2003-2026, The Open MPI Community. go forth and alltoall responsibly, ily, mwah πβ¨π
π¦ where this page comes from
- ubuntu stonking: openmpi-doc 5.0.10-6 Β· other
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash on : a language model's rewrite of the original, under its terms