[mpich-discuss] mpi windows with groups mpich2

[Helen Kershaw]

I've just built 3.1 on my laptop and the program runs fine. I'll have to 
bug the sys admins to upgrade to MPICH3.1 on the machine I want to run on.

Thanks for you time,
Helen

On 4/28/14 11:27 AM, Helen Kershaw wrote:
> Will do.
>
> On 4/28/14 11:24 AM, Rajeev Thakur wrote:
>> Can you try the latest, 3.1.
>>
>> Rajeev
>>
>> On Apr 28, 2014, at 12:21 PM, Helen Kershaw <hkershaw at ucar.edu>
>>   wrote:
>>
>>> Version 3.0.4.
>>>
>>> Helen
>>>
>>> On 4/28/14 10:38 AM, Rajeev Thakur wrote:
>>>> The program runs fine on my Mac laptop. Which version of MPICH are
>>>> you using?
>>>>
>>>> Rajeev
>>>>
>>>> On Apr 28, 2014, at 10:41 AM, Helen Kershaw <hkershaw at ucar.edu> wrote:
>>>>
>>>>> Hi,
>>>>>
>>>>> I'm trying to use MPI one-sided communication with sub-communicators.
>>>>> I've put my test code at the end of the email.  It runs on 8
>>>>> processors, creates two communicators, makes a window on the
>>>>> communicators, then tries to do some one-sided communication.
>>>>>
>>>>> I can use groups OR windows without a problem, but when I use them
>>>>> together (in the test code) I get a coredump from the first 4
>>>>> tasks.  I
>>>>> think the one-sided communication may be using the local rank in
>>>>> the sub-communicator (mpi_comm_grid) as if it is the global rank in
>>>>> mpi_comm_world.
>>>>>
>>>>> I can run the test code ok using open-mpi, but not mpich2.  I might be
>>>>> missing something obvious.
>>>>>
>>>>> Thanks for your help,
>>>>> Helen
>>>>>
>>>>> -----
>>>>> program test_group_win
>>>>>
>>>>> use mpi
>>>>>
>>>>> implicit none
>>>>>
>>>>> ! mpi
>>>>> integer            :: rank, size, ierr, status(mpi_status_size)
>>>>> ! groups
>>>>> integer, parameter :: group_size = 4
>>>>> integer            :: group_members(group_size)
>>>>> integer            :: group_all ! mpi_comm_world group
>>>>> integer            :: subgroup
>>>>> integer            :: mpi_comm_grid
>>>>> integer            :: local_rank ! rank within subgroup
>>>>>
>>>>>
>>>>> ! window
>>>>> integer                        :: win
>>>>> integer                        ::bytesize !> size in bytes of each
>>>>> element in the window
>>>>> integer(KIND=MPI_ADDRESS_KIND) :: window_size
>>>>> pointer(aa, duplicate)
>>>>> real                           :: duplicate(*)
>>>>> integer(KIND=MPI_ADDRESS_KIND) :: target_disp !> displacement in
>>>>> window
>>>>> integer                        :: owner
>>>>> real                           :: result
>>>>>
>>>>> call mpi_init(ierr)
>>>>> call mpi_comm_size(mpi_comm_world, size, ierr)
>>>>> call mpi_comm_rank(mpi_comm_world, rank, ierr)
>>>>>
>>>>> ! create groups
>>>>> call mpi_comm_group(mpi_comm_world, group_all, ierr)
>>>>>
>>>>> if (rank < 4 ) then
>>>>>    group_members = (/0,1,2,3/)
>>>>> else
>>>>>    group_members = (/4,5,6,7/)
>>>>> endif
>>>>>
>>>>> call mpi_group_incl(group_all, group_size, group_members, subgroup,
>>>>> ierr)
>>>>> call mpi_comm_create(mpi_comm_world, subgroup, mpi_comm_grid, ierr)
>>>>> call mpi_comm_rank(mpi_comm_grid, local_rank, ierr) ! rank within
>>>>> group
>>>>>
>>>>> ! create window
>>>>> call mpi_type_size(mpi_real, bytesize, ierr)
>>>>> window_size = bytesize ! one element in the window
>>>>> aa = malloc(1)
>>>>> call MPI_ALLOC_MEM(window_size, mpi_info_null, aa, ierr)
>>>>>
>>>>> duplicate(1) = rank
>>>>>
>>>>> call mpi_win_create(duplicate,  window_size, bytesize, MPI_INFO_NULL,
>>>>> mpi_comm_grid, win, ierr)
>>>>>
>>>>> ! grabbing data from rank above
>>>>> if (rank < 3 ) then
>>>>>    owner = rank + 1
>>>>> elseif (rank < 7) then
>>>>>    owner = rank + 1 - 4
>>>>> else
>>>>>    owner = 0
>>>>> endif
>>>>>
>>>>> print*, 'rank owner', rank, owner
>>>>>
>>>>> target_disp = 0
>>>>>
>>>>> call mpi_win_lock(MPI_LOCK_SHARED, owner, 0,  win, ierr)
>>>>> call mpi_get(result, 1, mpi_real, owner, target_disp, 1, mpi_real,
>>>>> win,
>>>>> ierr)
>>>>> call mpi_win_unlock(owner, win, ierr)
>>>>>
>>>>> print*, 'rank, result ', rank, result
>>>>>
>>>>> ! free window
>>>>> call mpi_win_free(win, ierr)
>>>>> call mpi_free_mem(duplicate, ierr)
>>>>>
>>>>> call mpi_finalize(ierr)
>>>>>
>>>>> end program test_group_win
>>>>> ----
>>>>>
>>>>>
>>>>> _______________________________________________
>>>>> discuss mailing list     discuss at mpich.org
>>>>> To manage subscription options or unsubscribe:
>>>>> https://lists.mpich.org/mailman/listinfo/discuss
>>>>
>>>> _______________________________________________
>>>> discuss mailing list     discuss at mpich.org
>>>> To manage subscription options or unsubscribe:
>>>> https://lists.mpich.org/mailman/listinfo/discuss
>>>>
>>> _______________________________________________
>>> discuss mailing list     discuss at mpich.org
>>> To manage subscription options or unsubscribe:
>>> https://lists.mpich.org/mailman/listinfo/discuss
>>
>> _______________________________________________
>> discuss mailing list     discuss at mpich.org
>> To manage subscription options or unsubscribe:
>> https://lists.mpich.org/mailman/listinfo/discuss
>>
> _______________________________________________
> discuss mailing list     discuss at mpich.org
> To manage subscription options or unsubscribe:
> https://lists.mpich.org/mailman/listinfo/discuss