X-Git-Url: http://info.iut-bm.univ-fcomte.fr/pub/gitweb/simgrid.git/blobdiff_plain/c57842fe897f53b46f0be3da87f7c996674be7d6..5ca17dc381bb7fba4f2689143e2080f7549ee0ab:/teshsuite/smpi/mpich3-test/rma/strided_getacc_indexed.c diff --git a/teshsuite/smpi/mpich3-test/rma/strided_getacc_indexed.c b/teshsuite/smpi/mpich3-test/rma/strided_getacc_indexed.c index e3293a1eff..88a08adbe7 100644 --- a/teshsuite/smpi/mpich3-test/rma/strided_getacc_indexed.c +++ b/teshsuite/smpi/mpich3-test/rma/strided_getacc_indexed.c @@ -6,7 +6,7 @@ /* One-Sided MPI 2-D Strided Accumulate Test * - * Author: James Dinan + * Author: James Dinan * Date : December, 2010 * * This code performs N strided put operations followed by get operations into @@ -28,7 +28,8 @@ #define SUB_YDIM 2 #define ITERATIONS 10 -int main(int argc, char **argv) { +int main(int argc, char **argv) +{ int i, j, rank, nranks, peer, bufsize, errors; double *win_buf, *src_buf, *dst_buf; MPI_Win buf_win; @@ -43,47 +44,47 @@ int main(int argc, char **argv) { MPI_Alloc_mem(bufsize, MPI_INFO_NULL, &src_buf); MPI_Alloc_mem(bufsize, MPI_INFO_NULL, &dst_buf); - for (i = 0; i < XDIM*YDIM; i++) { + for (i = 0; i < XDIM * YDIM; i++) { *(win_buf + i) = -1.0; - *(src_buf + i) = 1.0 + rank; + *(src_buf + i) = 1.0 + rank; } MPI_Win_create(win_buf, bufsize, 1, MPI_INFO_NULL, MPI_COMM_WORLD, &buf_win); - peer = (rank+1) % nranks; + peer = (rank + 1) % nranks; /* Perform ITERATIONS strided accumulate operations */ for (i = 0; i < ITERATIONS; i++) { - int idx_rem[SUB_YDIM]; - int blk_len[SUB_YDIM]; - MPI_Datatype src_type, dst_type; - - for (j = 0; j < SUB_YDIM; j++) { - idx_rem[j] = j*XDIM; - blk_len[j] = SUB_XDIM; - } - - MPI_Type_indexed(SUB_YDIM, blk_len, idx_rem, MPI_DOUBLE, &src_type); - MPI_Type_indexed(SUB_YDIM, blk_len, idx_rem, MPI_DOUBLE, &dst_type); - - MPI_Type_commit(&src_type); - MPI_Type_commit(&dst_type); - - /* PUT */ - MPI_Win_lock(MPI_LOCK_EXCLUSIVE, peer, 0, buf_win); - MPI_Get_accumulate(src_buf, 1, src_type, dst_buf, 1, src_type, peer, 0, - 1, dst_type, MPI_REPLACE, buf_win); - MPI_Win_unlock(peer, buf_win); - - /* GET */ - MPI_Win_lock(MPI_LOCK_EXCLUSIVE, peer, 0, buf_win); - MPI_Get_accumulate(src_buf, 1, src_type, dst_buf, 1, src_type, peer, 0, - 1, dst_type, MPI_NO_OP, buf_win); - MPI_Win_unlock(peer, buf_win); - - MPI_Type_free(&src_type); - MPI_Type_free(&dst_type); + int idx_rem[SUB_YDIM]; + int blk_len[SUB_YDIM]; + MPI_Datatype src_type, dst_type; + + for (j = 0; j < SUB_YDIM; j++) { + idx_rem[j] = j * XDIM; + blk_len[j] = SUB_XDIM; + } + + MPI_Type_indexed(SUB_YDIM, blk_len, idx_rem, MPI_DOUBLE, &src_type); + MPI_Type_indexed(SUB_YDIM, blk_len, idx_rem, MPI_DOUBLE, &dst_type); + + MPI_Type_commit(&src_type); + MPI_Type_commit(&dst_type); + + /* PUT */ + MPI_Win_lock(MPI_LOCK_EXCLUSIVE, peer, 0, buf_win); + MPI_Get_accumulate(src_buf, 1, src_type, dst_buf, 1, src_type, peer, 0, + 1, dst_type, MPI_REPLACE, buf_win); + MPI_Win_unlock(peer, buf_win); + + /* GET */ + MPI_Win_lock(MPI_LOCK_EXCLUSIVE, peer, 0, buf_win); + MPI_Get_accumulate(src_buf, 1, src_type, dst_buf, 1, src_type, peer, 0, + 1, dst_type, MPI_NO_OP, buf_win); + MPI_Win_unlock(peer, buf_win); + + MPI_Type_free(&src_type); + MPI_Type_free(&dst_type); } MPI_Barrier(MPI_COMM_WORLD); @@ -93,40 +94,40 @@ int main(int argc, char **argv) { MPI_Win_lock(MPI_LOCK_EXCLUSIVE, rank, 0, buf_win); errors = 0; for (i = 0; i < SUB_XDIM; i++) { - for (j = 0; j < SUB_YDIM; j++) { - const double actual = *(win_buf + i + j*XDIM); - const double expected = (1.0 + ((rank+nranks-1)%nranks)); - if (fabs(actual - expected) > 1.0e-10) { - SQUELCH( printf("%d: Data validation failed at [%d, %d] expected=%f actual=%f\n", - rank, j, i, expected, actual); ); - errors++; - fflush(stdout); + for (j = 0; j < SUB_YDIM; j++) { + const double actual = *(win_buf + i + j * XDIM); + const double expected = (1.0 + ((rank + nranks - 1) % nranks)); + if (fabs(actual - expected) > 1.0e-10) { + SQUELCH(printf("%d: Data validation failed at [%d, %d] expected=%f actual=%f\n", + rank, j, i, expected, actual);); + errors++; + fflush(stdout); + } } - } } for (i = SUB_XDIM; i < XDIM; i++) { - for (j = 0; j < SUB_YDIM; j++) { - const double actual = *(win_buf + i + j*XDIM); - const double expected = -1.0; - if (fabs(actual - expected) > 1.0e-10) { - SQUELCH( printf("%d: Data validation failed at [%d, %d] expected=%f actual=%f\n", - rank, j, i, expected, actual); ); - errors++; - fflush(stdout); + for (j = 0; j < SUB_YDIM; j++) { + const double actual = *(win_buf + i + j * XDIM); + const double expected = -1.0; + if (fabs(actual - expected) > 1.0e-10) { + SQUELCH(printf("%d: Data validation failed at [%d, %d] expected=%f actual=%f\n", + rank, j, i, expected, actual);); + errors++; + fflush(stdout); + } } - } } for (i = 0; i < XDIM; i++) { - for (j = SUB_YDIM; j < YDIM; j++) { - const double actual = *(win_buf + i + j*XDIM); - const double expected = -1.0; - if (fabs(actual - expected) > 1.0e-10) { - SQUELCH( printf("%d: Data validation failed at [%d, %d] expected=%f actual=%f\n", - rank, j, i, expected, actual); ); - errors++; - fflush(stdout); + for (j = SUB_YDIM; j < YDIM; j++) { + const double actual = *(win_buf + i + j * XDIM); + const double expected = -1.0; + if (fabs(actual - expected) > 1.0e-10) { + SQUELCH(printf("%d: Data validation failed at [%d, %d] expected=%f actual=%f\n", + rank, j, i, expected, actual);); + errors++; + fflush(stdout); + } } - } } MPI_Win_unlock(rank, buf_win); @@ -135,7 +136,7 @@ int main(int argc, char **argv) { MPI_Free_mem(src_buf); MPI_Free_mem(dst_buf); - MTest_Finalize( errors ); + MTest_Finalize(errors); MPI_Finalize(); - return MTestReturnValue( errors ); + return MTestReturnValue(errors); }