X-Git-Url: http://info.iut-bm.univ-fcomte.fr/pub/gitweb/simgrid.git/blobdiff_plain/fc0a76c21f0bceb0314aee5d62f785ea8009495e..5ca17dc381bb7fba4f2689143e2080f7549ee0ab:/teshsuite/smpi/mpich3-test/rma/strided_acc_subarray.c diff --git a/teshsuite/smpi/mpich3-test/rma/strided_acc_subarray.c b/teshsuite/smpi/mpich3-test/rma/strided_acc_subarray.c index c8f850c6a7..b1c667c739 100644 --- a/teshsuite/smpi/mpich3-test/rma/strided_acc_subarray.c +++ b/teshsuite/smpi/mpich3-test/rma/strided_acc_subarray.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 accumulates into a 2d patch of a shared array. The @@ -22,13 +22,14 @@ #include "mpitest.h" #include "squelch.h" -#define XDIM 1024 +#define XDIM 1024 #define YDIM 1024 #define SUB_XDIM 512 #define SUB_YDIM 512 #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; MPI_Win buf_win; @@ -42,44 +43,44 @@ int main(int argc, char **argv) { MPI_Alloc_mem(bufsize, MPI_INFO_NULL, &win_buf); MPI_Alloc_mem(bufsize, MPI_INFO_NULL, &src_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; } 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 ndims = 2; - int src_arr_sizes[2] = { XDIM, YDIM }; - int src_arr_subsizes[2] = { SUB_XDIM, SUB_YDIM }; - int src_arr_starts[2] = { 0, 0 }; - int dst_arr_sizes[2] = { XDIM, YDIM }; - int dst_arr_subsizes[2] = { SUB_XDIM, SUB_YDIM }; - int dst_arr_starts[2] = { 0, 0 }; - MPI_Datatype src_type, dst_type; + int ndims = 2; + int src_arr_sizes[2] = { XDIM, YDIM }; + int src_arr_subsizes[2] = { SUB_XDIM, SUB_YDIM }; + int src_arr_starts[2] = { 0, 0 }; + int dst_arr_sizes[2] = { XDIM, YDIM }; + int dst_arr_subsizes[2] = { SUB_XDIM, SUB_YDIM }; + int dst_arr_starts[2] = { 0, 0 }; + MPI_Datatype src_type, dst_type; - MPI_Type_create_subarray(ndims, src_arr_sizes, src_arr_subsizes, src_arr_starts, - MPI_ORDER_C, MPI_DOUBLE, &src_type); + MPI_Type_create_subarray(ndims, src_arr_sizes, src_arr_subsizes, src_arr_starts, + MPI_ORDER_C, MPI_DOUBLE, &src_type); - MPI_Type_create_subarray(ndims, dst_arr_sizes, dst_arr_subsizes, dst_arr_starts, - MPI_ORDER_C, MPI_DOUBLE, &dst_type); + MPI_Type_create_subarray(ndims, dst_arr_sizes, dst_arr_subsizes, dst_arr_starts, + MPI_ORDER_C, MPI_DOUBLE, &dst_type); - MPI_Type_commit(&src_type); - MPI_Type_commit(&dst_type); + MPI_Type_commit(&src_type); + MPI_Type_commit(&dst_type); - MPI_Win_lock(MPI_LOCK_EXCLUSIVE, peer, 0, buf_win); + MPI_Win_lock(MPI_LOCK_EXCLUSIVE, peer, 0, buf_win); - MPI_Accumulate(src_buf, 1, src_type, peer, 0, 1, dst_type, MPI_SUM, buf_win); + MPI_Accumulate(src_buf, 1, src_type, peer, 0, 1, dst_type, MPI_SUM, buf_win); - MPI_Win_unlock(peer, buf_win); + MPI_Win_unlock(peer, buf_win); - MPI_Type_free(&src_type); - MPI_Type_free(&dst_type); + MPI_Type_free(&src_type); + MPI_Type_free(&dst_type); } MPI_Barrier(MPI_COMM_WORLD); @@ -89,40 +90,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 + (1.0 + ((rank+nranks-1)%nranks)) * (ITERATIONS); - 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 + (1.0 + ((rank + nranks - 1) % nranks)) * (ITERATIONS); + 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); @@ -130,7 +131,7 @@ int main(int argc, char **argv) { MPI_Free_mem(win_buf); MPI_Free_mem(src_buf); - MTest_Finalize( errors ); + MTest_Finalize(errors); MPI_Finalize(); - return MTestReturnValue( errors ); + return MTestReturnValue(errors); }