- /* To improve reporting of problems about operations, we
- change the error handler to errors return */
- MPI_Win_set_errhandler( win, MPI_ERRORS_RETURN );
-
- /* MPI_REPLACE on accumulate is almost the same
- as MPI_Put; the only difference is in the
- handling of overlapping accumulate operations,
- which are not tested here */
- err = MPI_Accumulate( sendtype.buf, sendtype.count,
- sendtype.datatype, dest, 0,
- recvtype.count, recvtype.datatype,
- MPI_REPLACE, win );
- if (err) {
- errs++;
- if (errs < 10) {
- printf( "Accumulate types: send %s, recv %s\n",
- MTestGetDatatypeName( &sendtype ),
- MTestGetDatatypeName( &recvtype ) );
- MTestPrintError( err );
- }
- }
- err = MPI_Win_fence( 0, win );
- if (err) {
- errs++;
- if (errs < 10) {
- MTestPrintError( err );
- }
- }
- }
- else if (rank == dest) {
- MPI_Win_fence( 0, win );
- /* This should have the same effect, in terms of
- transfering data, as a send/recv pair */
- err = MTestCheckRecv( 0, &recvtype );
- if (err) {
- errs += err;
- }
- }
- else {
- MPI_Win_fence( 0, win );
- }
- MPI_Win_free( &win );
- MTestFreeDatatype( &sendtype );
- MTestFreeDatatype( &recvtype );
- }
- }
+ MPI_Type_extent(recvtype.datatype, &extent);
+ MPI_Type_lb(recvtype.datatype, &lb);
+ MPI_Win_create(recvtype.buf, recvtype.count * extent + lb,
+ (int) extent, MPI_INFO_NULL, comm, &win);
+ MPI_Win_fence(0, win);
+ if (rank == source) {
+ sendtype.InitBuf(&sendtype);
+
+ /* To improve reporting of problems about operations, we
+ * change the error handler to errors return */
+ MPI_Win_set_errhandler(win, MPI_ERRORS_RETURN);
+
+ /* MPI_REPLACE on accumulate is almost the same
+ * as MPI_Put; the only difference is in the
+ * handling of overlapping accumulate operations,
+ * which are not tested here */
+ err = MPI_Accumulate(sendtype.buf, sendtype.count,
+ sendtype.datatype, dest, 0,
+ recvtype.count, recvtype.datatype, MPI_REPLACE, win);
+ if (err) {
+ errs++;
+ if (errs < 10) {
+ printf("Accumulate types: send %s, recv %s\n",
+ MTestGetDatatypeName(&sendtype),
+ MTestGetDatatypeName(&recvtype));
+ MTestPrintError(err);
+ }
+ }
+ err = MPI_Win_fence(0, win);
+ if (err) {
+ errs++;
+ if (errs < 10) {
+ MTestPrintError(err);
+ }
+ }
+ }
+ else if (rank == dest) {
+ MPI_Win_fence(0, win);
+ /* This should have the same effect, in terms of
+ * transfering data, as a send/recv pair */
+ err = MTestCheckRecv(0, &recvtype);
+ if (err) {
+ errs += err;
+ }
+ }
+ else {
+ MPI_Win_fence(0, win);
+ }
+ MPI_Win_free(&win);
+ MTestFreeDatatype(&sendtype);
+ MTestFreeDatatype(&recvtype);
+ }
+ }