- discuss - mailman3.cels.anl.gov

Re: [mpich-discuss] libmpich.so: undefined reference to ....
by Dave Goodell 04 Apr '13

04 Apr '13

It looks like Delft3D is building its libraries without using "mpicc" or any of the other compiler wrappers. These wrappers add numerous internal MPI libraries to the link line, including "-lmpl", only some of which are captured in the existing Delft3D build system. Please either use the compiler wrapper scripts to build your application or at least ensure that you are linking with all of the flags shown by "mpicc -show". I also strongly recommend moving to MPICH 3.0.3 instead of MPICH2 1.4.1. MPICH2 >=1.5 (and therefore MPICH >=3.0) have a newer build system that handles shared libraries in a substantially more robust manner. -Dave On Apr 4, 2013, at 9:44 AM CDT, "Boyington, Thomas Matthew" <tmboyington(a)tva.gov> wrote: > Dear mpich community, > > I’m getting numerous libmpich.so: undefined reference to XXX errors when building Delft3D (a CFD costal/river model). I’ve read the log files, Googled, and rtm’ed to no avail. I hope that the mpich mailing list can help me out. > > I am running Red Hat Enterprise 6.3 64-bit, and I am building mpich 1.4.1. I have been able to build mpich successfully (I think…the examples described in the README file work) using the commands shown below. I’m using the Intel 2013 Fortran compiler, and I’m using the –enabled-shared, CFLAGS=”-fPIC” and CXXFLAGS=”-fPIC” switches because of a previous Delft3D build error saying to build mpich with the –fPIC switch. > > > configure --prefix=/home/labuser/mpich2-install --enable-shared CFLAGS="-fPIC" CXXFLAGS="-fPIC" F77=ifort FC=ifort |& tee c.txt > > make 2>&1 | tee m.txt > > make install 2>&1 | tee mi.txt > > When building Delft3D (again with the Intel compiler) I get the following errors in the make log. I’ve included my log files as attachments. Any suggestions on how to address this will be greatly appreciated? > > libNefis.a ../../../../../engines_gpl/waq/packages/waq_utils_f/src/interact/.libs/libinteract.a ../../../../../engines_gpl/waq/packages/waq_utils_c/src/.libs/libwaq_utils_c.a ../../../../../engines_gpl/waq/packages/delwaq2_version_number/src/.libs/libdelwaq2_version_number.a -lfmpich -lmpich -lpthread -ldl -lrt -lexpat -lstdc++ -fopenmp -threads > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trid' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trvalid' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_env2int' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trrealloc' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trspace' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trDebugLevel' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_TrSetMaxMem' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trlevel' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trmalloc' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_env2bool' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_env2range' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trcalloc' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trfree' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_env2str' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trstrdup' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trdump' > > /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trinit' > > make[5]: *** [delwaq2] Error 1 > > make[5]: Leaving directory `/home/labuser/Delft3D/2367/src/engines_gpl/waq/packages/delwaq2/src' > > make[4]: *** [install-recursive] Error 1 > > make[4]: Leaving directory `/home/labuser/Delft3D/2367/src/engines_gpl/waq/packages/delwaq2' > > make[3]: *** [install-recursive] Error 1 > > make[3]: Leaving directory `/home/labuser/Delft3D/2367/src/engines_gpl/waq/packages' > > make[2]: *** [install-recursive] Error 1 > > make[2]: Leaving directory `/home/labuser/Delft3D/2367/src/engines_gpl/waq' > > make[1]: *** [install-recursive] Error 1 > > make[1]: Leaving directory `/home/labuser/Delft3D/2367/src/engines_gpl' > > make: *** [install-recursive] Error 1 > > > Thanks, > > Matt > > > -- > > T. Matthew Boyington > > Civil Engineer, Water Resources > > Hydrothermal Team > > River Operations & Renewables > > Tennessee Valley Authority > > 400 West Summit Hill Drive, Knoxville, TN 37902 > > Office: (865) 632-8365 > > Cell: (865) 599-8265 > > > <<mi.zip>> <<config_openpa.zip>> <<config.zip>> <<c.zip>> <<m.zip>> > > <mi.zip><config_openpa.zip><config.zip><c.zip><m.zip>_______________________________________________ > discuss mailing list discuss(a)mpich.org > To manage subscription options or unsubscribe: > https://lists.mpich.org/mailman/listinfo/discuss

1 0

libmpich.so: undefined reference to ....
by Boyington, Thomas Matthew 04 Apr '13

04 Apr '13

Dear mpich community, I'm getting numerous libmpich.so: undefined reference to XXX errors when building Delft3D (a CFD costal/river model). I've read the log files, Googled, and rtm'ed to no avail. I hope that the mpich mailing list can help me out. I am running Red Hat Enterprise 6.3 64-bit, and I am building mpich 1.4.1. I have been able to build mpich successfully (I think...the examples described in the README file work) using the commands shown below. I'm using the Intel 2013 Fortran compiler, and I'm using the -enabled-shared, CFLAGS="-fPIC" and CXXFLAGS="-fPIC" switches because of a previous Delft3D build error saying to build mpich with the -fPIC switch. configure --prefix=/home/labuser/mpich2-install --enable-shared CFLAGS="-fPIC" CXXFLAGS="-fPIC" F77=ifort FC=ifort |& tee c.txt make 2>&1 | tee m.txt make install 2>&1 | tee mi.txt When building Delft3D (again with the Intel compiler) I get the following errors in the make log. I've included my log files as attachments. Any suggestions on how to address this will be greatly appreciated? libNefis.a ../../../../../engines_gpl/waq/packages/waq_utils_f/src/interact/.libs/l ibinteract.a ../../../../../engines_gpl/waq/packages/waq_utils_c/src/.libs/libwaq_uti ls_c.a ../../../../../engines_gpl/waq/packages/delwaq2_version_number/src/.libs /libdelwaq2_version_number.a -lfmpich -lmpich -lpthread -ldl -lrt -lexpat -lstdc++ -fopenmp -threads /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trid' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trvalid' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_env2int' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trrealloc' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trspace' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trDebugLevel' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_TrSetMaxMem' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trlevel' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trmalloc' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_env2bool' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_env2range' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trcalloc' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trfree' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_env2str' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trstrdup' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trdump' /home/labuser/mpich2-install/lib/libmpich.so: undefined reference to `MPL_trinit' make[5]: *** [delwaq2] Error 1 make[5]: Leaving directory `/home/labuser/Delft3D/2367/src/engines_gpl/waq/packages/delwaq2/src' make[4]: *** [install-recursive] Error 1 make[4]: Leaving directory `/home/labuser/Delft3D/2367/src/engines_gpl/waq/packages/delwaq2' make[3]: *** [install-recursive] Error 1 make[3]: Leaving directory `/home/labuser/Delft3D/2367/src/engines_gpl/waq/packages' make[2]: *** [install-recursive] Error 1 make[2]: Leaving directory `/home/labuser/Delft3D/2367/src/engines_gpl/waq' make[1]: *** [install-recursive] Error 1 make[1]: Leaving directory `/home/labuser/Delft3D/2367/src/engines_gpl' make: *** [install-recursive] Error 1 Thanks, Matt -- T. Matthew Boyington Civil Engineer, Water Resources Hydrothermal Team River Operations & Renewables Tennessee Valley Authority 400 West Summit Hill Drive, Knoxville, TN 37902 Office: (865) 632-8365 Cell: (865) 599-8265 <<mi.zip>> <<config_openpa.zip>> <<config.zip>> <<c.zip>> <<m.zip>>

1 0

Re: [mpich-discuss] MPI IO, reading MPI Writes without MPI read
by Rob Latham 01 Apr '13

01 Apr '13

On Thu, Mar 28, 2013 at 12:53:42PM -0700, Ryan Crocker wrote: > I just went though the standard, i'm not doing anything to untoward. that i could see. > > I also would really not have to re-run these simulations. There has to be a way to get access to those MPI writes with an MPI read, and like i said, other programs have no trouble reading these in and using them (paravew, and ensight). I would just like to have my own stand alone in fortran. The writes come out of a fortran code compiled with mpich so there should be a way to read them back in with a fortran, mpich compiled code. Shouldn't there? If you have confirmed that the data really does live in that file, with a 244 byte header and three 'ncells_hexa*4'-byte regions of data, then heck yeah we can read it all back into fortran. The quoting's making it hard to read so here's the code again: ncells = parallel_sum(ncells_hexa_,ncells) allocate(buffer3(ncells,3)) openfile=trim(workdir)//'/'//'V/V.000002' call MPI_FILE_OPEN(comm,openfile,MPI_MODE_RDONLY,mpi_info,iunit,ierr) ! Read header bsize = 80 call MPI_FILE_READ(iunit,cbuffer,bsize,MPI_CHARACTER,status,ierr) print*,trim(cbuffer) bsize = 80 call MPI_FILE_READ(iunit,cbuffer,bsize,MPI_CHARACTER,status,ierr) print*,trim(cbuffer) bsize = 1 call MPI_FILE_READ(iunit,ibuffer,bsize,MPI_INTEGER,status,ierr) print*,ibuffer bsize = 80 call MPI_FILE_READ(iunit,cbuffer,bsize,MPI_CHARACTER,status,ierr) print*,trim(cbuffer),ncells ! Read the data disp = 3*80+4+0*ncells*4 call MPI_FILE_SET_VIEW(iunit,disp,MPI_REAL_SP,fileview_hexa,"native",mpi_info,ierr) call MPI_FILE_READ_ALL(iunit,buffer3(:,1),ncells,MPI_REAL_SP,status,ierr) disp = 3*80+4+1*ncells*4 call MPI_FILE_SET_VIEW(iunit,disp,MPI_REAL_SP,fileview_hexa,"native",mpi_info,ierr) call MPI_FILE_READ_ALL(iunit,buffer3(:,2),ncells,MPI_REAL_SP,status,ierr) disp = 3*80+4+2*ncells*4 call MPI_FILE_SET_VIEW(iunit,disp,MPI_REAL_SP,fileview_hexa,"native",mpi_info,ierr) call MPI_FILE_READ_ALL(iunit,buffer3(:,3),ncells,MPI_REAL_SP,status,ierr) ! Close the file call MPI_FILE_CLOSE(iunit,ierr) The symptom is htat 'buffer3' ends up with all the same values, right? I think what Rajeev was suggesting is that you pass not 'buffer3(:,1)' to MPI_File_read_all -- the entire array -- but rather pass in the first element of that (row? column? I have very little fortran experience, sorry). Your use of MPI_FILE_SET_VIEW to move the file pointer to the right place is uncommon but should be correct. One way to rationalize this is that many aspects of the fortran interface are designed to accommodate "speaking to C". ==rob -- Rob Latham Mathematics and Computer Science Division Argonne National Lab, IL USA

1 0

Re: [mpich-discuss] MPICH 3.0.3 Ubuntu, Debian and Slackware packages (x86 and amd64)
by Dave Goodell 01 Apr '13

01 Apr '13

Downloads page updated (http://www.mpich.org/downloads/) Thanks! -Dave On Apr 1, 2013, at 2:18 PM CDT, Marcelo Soares Souza <marcelo(a)cebacad.net> wrote: > MPICH 3.0.3 packages > > * Ubuntu 12.04.1 LTS (Precise Pangolin) > > - amd64 > http://www.cebacad.net/files/mpich/ubuntu/mpich-3.0.3/mpich_3.0.3-1ubuntu_a… > > - x86 > http://www.cebacad.net/files/mpich/ubuntu/mpich-3.0.3/mpich_3.0.3-1ubuntu_i… > > * Debian 7.0 (Wheezy) > > - amd64 > http://www.cebacad.net/files/mpich/debian/mpich-3.0.3/mpich_3.0.3-1debian_a… > > - x86 > http://www.cebacad.net/files/mpich/debian/mpich-3.0.3/mpich_3.0.3-1debian_i… > > * Slackware 14.0 > > - amd64 > http://www.cebacad.net/files/mpich/slackware/mpich-3.0.3/mpich-3.0.3-x86_64… > > - x86 > http://www.cebacad.net/files/mpich/slackware/mpich-3.0.3/mpich-3.0.3-i486-g… > > Enjoy! ;) > > -- > Abraços > Marcelo Soares Souza > http://marcelo.cebacad.net > _______________________________________________ > discuss mailing list discuss(a)mpich.org > To manage subscription options or unsubscribe: > https://lists.mpich.org/mailman/listinfo/discuss

1 0

MPICH 3.0.3 Ubuntu, Debian and Slackware packages (x86 and amd64)
by Marcelo Soares Souza 01 Apr '13

01 Apr '13

MPICH 3.0.3 packages * Ubuntu 12.04.1 LTS (Precise Pangolin) - amd64 http://www.cebacad.net/files/mpich/ubuntu/mpich-3.0.3/mpich_3.0.3-1ubuntu_a… - x86 http://www.cebacad.net/files/mpich/ubuntu/mpich-3.0.3/mpich_3.0.3-1ubuntu_i… * Debian 7.0 (Wheezy) - amd64 http://www.cebacad.net/files/mpich/debian/mpich-3.0.3/mpich_3.0.3-1debian_a… - x86 http://www.cebacad.net/files/mpich/debian/mpich-3.0.3/mpich_3.0.3-1debian_i… * Slackware 14.0 - amd64 http://www.cebacad.net/files/mpich/slackware/mpich-3.0.3/mpich-3.0.3-x86_64… - x86 http://www.cebacad.net/files/mpich/slackware/mpich-3.0.3/mpich-3.0.3-i486-g… Enjoy! ;) -- Abraços Marcelo Soares Souza http://marcelo.cebacad.net

1 0

Re: [mpich-discuss] jeangosselin
by Jean Gosselin 01 Apr '13

01 Apr '13

http://www.pf3paintsupply.com/nxxpwaa/xiue.dqzi?vck Jean Gosselin uxuzghet

1 0

building mpich
by Иващук Роман 01 Apr '13

01 Apr '13

There is no file libmpi_f77_so.0.0.1 in the directories after building. Please help to resolve the problem.

1 0

Re: [mpich-discuss] MPI_PUT Implementation Bug or Pilot Error?
by Timothy Stitt 29 Mar '13

29 Mar '13

Thanks Rajeev. That makes sense now...and a modified version using separate storage per iteration works perfectly. Tim. On Mar 29, 2013, at 1:55 PM, Rajeev Thakur <thakur(a)mcs.anl.gov<mailto:[email protected]>> wrote: MPI_Put is nonblocking. You are changing the contents of the source buffer before the put is complete. Rajeev On Mar 29, 2013, at 12:01 PM, Timothy Stitt wrote: Hi all, I've come across some strange behavior with a code that uses MPI_PUT(). I can confirm that this problem appears in older MPICH2 v1.4.x releases as well as the latest MPICH v3.0.3 release. I can also confirm that I receive the correct results when using Open MPI. I am attaching a reproducer code below that simulates the problem I am having (following the code will hopefully explain its purpose better than what I can describe here in a few lines). This case is hardcoded for a 4-process run. For some reason MPICH implementations are losing some data elements and duplicating others after the communication is complete. I would be grateful if someone can confirm whether this is indeed an implementation bug or a subtle problem with how I've implemented the one-sided communication. Thanks in advance, Tim. program reproducer use mpi implicit none integer :: error,dum_com(2),my_rank,i integer :: win_values_points, type_integer,NProc=4 INTEGER(KIND=MPI_ADDRESS_KIND) :: win_values_size, target_disp integer,allocatable :: neighbor_values(:), ncount_in(:) call MPI_INIT(error) call MPI_COMM_RANK(MPI_COMM_WORLD,my_rank,error) CALL MPI_Type_size(MPI_INTEGER,type_integer,error) ! Calculate Size Of Window (in bytes) win_values_size=2*type_integer*NProc ! Allocate Window Storage allocate(neighbor_values(0:2*Nproc-1)) neighbor_values=0 ! Define Input Values allocate(ncount_in(NProc)) if (my_rank .eq. 0) ncount_in = (/1,2,3,4/) if (my_rank .eq. 1) ncount_in = (/5,6,7,8/) if (my_rank .eq. 2) ncount_in = (/9,10,11,12/) if (my_rank .eq. 3) ncount_in = (/13,14,15,16/) ! Create Window CALL MPI_WIN_CREATE(neighbor_values, win_values_size, type_integer, MPI_INFO_NULL, MPI_COMM_WORLD, win_values_points, error) ! Create Neighbor Storage Win dow call MPI_Win_Fence(0,win_values_points,error) ! Print pre-communication data print*,my_rank,ncount_in call flush(5) call MPI_BARRIER(MPI_COMM_WORLD,error) ! Communicate Data Values to all proceses do i = 1,4 target_disp=2*my_rank dum_com(1)=ncount_in(i) dum_com(2)=99 CALL MPI_PUT(dum_com,2,MPI_INTEGER,Nproc-i,target_disp,2,MPI_INTEGER,win_values_points,error) end do call MPI_Win_Fence(0,win_values_points,error) ! Print post-communication values print * print *,my_rank,neighbor_values call flush(5) call MPI_BARRIER(MPI_COMM_WORLD,error) ! Free Window Storage call mpi_win_free(win_values_points,error) call MPI_FINALIZE(error) end program reproducer _______________________________________________ discuss mailing list discuss(a)mpich.org<mailto:[email protected]> To manage subscription options or unsubscribe: https://lists.mpich.org/mailman/listinfo/discuss _______________________________________________ discuss mailing list discuss(a)mpich.org<mailto:[email protected]> To manage subscription options or unsubscribe: https://lists.mpich.org/mailman/listinfo/discuss Tim Stitt PhD User Support Manager (CRC) Research Assistant Professor (Computer Science & Engineering) Room 108, Center for Research Computing | University of Notre Dame | P.O. Box 539, Notre Dame, IN 46556 | Phone: 574-631-5287 | Email: tstitt(a)nd.edu<mailto:[email protected]>

1 0

Re: [mpich-discuss] MPI_PUT Implementation Bug or Pilot Error?
by Rajeev Thakur 29 Mar '13

29 Mar '13

MPI_Put is nonblocking. You are changing the contents of the source buffer before the put is complete. Rajeev On Mar 29, 2013, at 12:01 PM, Timothy Stitt wrote: > Hi all, > > I've come across some strange behavior with a code that uses MPI_PUT(). I can confirm that this problem appears in older MPICH2 v1.4.x releases as well as the latest MPICH v3.0.3 release. I can also confirm that I receive the correct results when using Open MPI. > > I am attaching a reproducer code below that simulates the problem I am having (following the code will hopefully explain its purpose better than what I can describe here in a few lines). This case is hardcoded for a 4-process run. For some reason MPICH implementations are losing some data elements and duplicating others after the communication is complete. I would be grateful if someone can confirm whether this is indeed an implementation bug or a subtle problem with how I've implemented the one-sided communication. > > Thanks in advance, > > Tim. > > program reproducer > > use mpi > > implicit none > > integer :: error,dum_com(2),my_rank,i > integer :: win_values_points, type_integer,NProc=4 > INTEGER(KIND=MPI_ADDRESS_KIND) :: win_values_size, target_disp > integer,allocatable :: neighbor_values(:), ncount_in(:) > > call MPI_INIT(error) > call MPI_COMM_RANK(MPI_COMM_WORLD,my_rank,error) > > CALL MPI_Type_size(MPI_INTEGER,type_integer,error) > > ! Calculate Size Of Window (in bytes) > win_values_size=2*type_integer*NProc > > ! Allocate Window Storage > allocate(neighbor_values(0:2*Nproc-1)) > neighbor_values=0 > > ! Define Input Values > allocate(ncount_in(NProc)) > if (my_rank .eq. 0) ncount_in = (/1,2,3,4/) > if (my_rank .eq. 1) ncount_in = (/5,6,7,8/) > if (my_rank .eq. 2) ncount_in = (/9,10,11,12/) > if (my_rank .eq. 3) ncount_in = (/13,14,15,16/) > > ! Create Window > CALL MPI_WIN_CREATE(neighbor_values, win_values_size, type_integer, MPI_INFO_NULL, MPI_COMM_WORLD, win_values_points, error) ! Create Neighbor Storage Win > dow > call MPI_Win_Fence(0,win_values_points,error) > > ! Print pre-communication data > print*,my_rank,ncount_in > call flush(5) > call MPI_BARRIER(MPI_COMM_WORLD,error) > > ! Communicate Data Values to all proceses > do i = 1,4 > target_disp=2*my_rank > dum_com(1)=ncount_in(i) > dum_com(2)=99 > CALL MPI_PUT(dum_com,2,MPI_INTEGER,Nproc-i,target_disp,2,MPI_INTEGER,win_values_points,error) > end do > > call MPI_Win_Fence(0,win_values_points,error) > > ! Print post-communication values > print * > print *,my_rank,neighbor_values > call flush(5) > call MPI_BARRIER(MPI_COMM_WORLD,error) > > ! Free Window Storage > call mpi_win_free(win_values_points,error) > > call MPI_FINALIZE(error) > > end program reproducer > > > _______________________________________________ > discuss mailing list discuss(a)mpich.org > To manage subscription options or unsubscribe: > https://lists.mpich.org/mailman/listinfo/discuss

1 0

MPI_PUT Implementation Bug or Pilot Error?
by Timothy Stitt 29 Mar '13

29 Mar '13

Hi all, I've come across some strange behavior with a code that uses MPI_PUT(). I can confirm that this problem appears in older MPICH2 v1.4.x releases as well as the latest MPICH v3.0.3 release. I can also confirm that I receive the correct results when using Open MPI. I am attaching a reproducer code below that simulates the problem I am having (following the code will hopefully explain its purpose better than what I can describe here in a few lines). This case is hardcoded for a 4-process run. For some reason MPICH implementations are losing some data elements and duplicating others after the communication is complete. I would be grateful if someone can confirm whether this is indeed an implementation bug or a subtle problem with how I've implemented the one-sided communication. Thanks in advance, Tim. program reproducer use mpi implicit none integer :: error,dum_com(2),my_rank,i integer :: win_values_points, type_integer,NProc=4 INTEGER(KIND=MPI_ADDRESS_KIND) :: win_values_size, target_disp integer,allocatable :: neighbor_values(:), ncount_in(:) call MPI_INIT(error) call MPI_COMM_RANK(MPI_COMM_WORLD,my_rank,error) CALL MPI_Type_size(MPI_INTEGER,type_integer,error) ! Calculate Size Of Window (in bytes) win_values_size=2*type_integer*NProc ! Allocate Window Storage allocate(neighbor_values(0:2*Nproc-1)) neighbor_values=0 ! Define Input Values allocate(ncount_in(NProc)) if (my_rank .eq. 0) ncount_in = (/1,2,3,4/) if (my_rank .eq. 1) ncount_in = (/5,6,7,8/) if (my_rank .eq. 2) ncount_in = (/9,10,11,12/) if (my_rank .eq. 3) ncount_in = (/13,14,15,16/) ! Create Window CALL MPI_WIN_CREATE(neighbor_values, win_values_size, type_integer, MPI_INFO_NULL, MPI_COMM_WORLD, win_values_points, error) ! Create Neighbor Storage Win dow call MPI_Win_Fence(0,win_values_points,error) ! Print pre-communication data print*,my_rank,ncount_in call flush(5) call MPI_BARRIER(MPI_COMM_WORLD,error) ! Communicate Data Values to all proceses do i = 1,4 target_disp=2*my_rank dum_com(1)=ncount_in(i) dum_com(2)=99 CALL MPI_PUT(dum_com,2,MPI_INTEGER,Nproc-i,target_disp,2,MPI_INTEGER,win_values_points,error) end do call MPI_Win_Fence(0,win_values_points,error) ! Print post-communication values print * print *,my_rank,neighbor_values call flush(5) call MPI_BARRIER(MPI_COMM_WORLD,error) ! Free Window Storage call mpi_win_free(win_values_points,error) call MPI_FINALIZE(error) end program reproducer

1 0