fortran - scalapack での一貫性のない行割り当て

Question

次の単純な fortran プログラムを考えてみましょう

program test_vec_allocation
    use mpi
    implicit none
    integer(kind=8)             :: N
    ! =========================BLACS and MPI=======================
    integer                     :: ierr, size, rank,dims(2)
    ! -------------------------------------------------------------
    integer, parameter          :: block_size = 100
    integer                     :: context, nprow, npcol, local_nprow, local_npcol
    integer                     :: numroc, indxl2g, descmat(9),descvec(9)
    integer                     :: mloc_mat ,nloc_mat ,mloc_vec ,nloc_vec

    call blacs_pinfo(rank,size)
    dims=0
    call MPI_Dims_create(size, 2, dims, ierr)
    nprow = dims(1);npcol = dims(2)
    call blacs_get(0,0,context)
    call blacs_gridinit(context, 'R', nprow, npcol)
    call blacs_gridinfo(context, nprow, npcol, local_nprow,local_npcol)

    N = 700

    mloc_vec = numroc(N,block_size,local_nprow,0, nprow)
    nloc_vec = numroc(1,block_size,local_npcol,0, npcol)
    print *,"Rank", rank, mloc_vec, nloc_vec

    call blacs_gridexit(context)
    call blacs_exit(0)

end program test_vec_allocation

11 mpi ランクで実行すると、

 Rank           0         100           1
 Rank           4         100           1
 Rank           2         100           1
 Rank           1         100           1
 Rank           3         100           1
 Rank          10           0           1
 Rank           6         100           1
 Rank           5         100           1
 Rank           9           0           1
 Rank           8           0           1
 Rank           7           0           1

これは、scalapackがこの配列を分割することを期待する方法ですが、ランクの数が偶数の場合は次のようになります。

 Rank           0         200           1
 Rank           8         200           0
 Rank           9         100           1
 Rank          10         100           0
 Rank           1         200           0
 Rank           6         200           1
 Rank          11         100           0
 Rank           3         200           1
 Rank           4         200           0
 Rank           2         200           0
 Rank           7         200           0
 Rank           5         200           0

これは意味がありません。ブロック サイズ 100 でランク 0 が 200 要素を取得し、ランク * ブロック サイズ > N になるのはなぜですか。このため、私のプログラムは mpi ランク 1、2、3、5、7、11 では機能しますが、ランク 4 では失敗します。 、6、8、9、10、12など（ランク9で失敗する理由はわかりません！）。私のアプローチで何が間違っているのか誰かが説明できますか?

GFortran バージョン: 6.1.0

SCALPACK バージョン: 2.1.0

macOS バージョン: 10.11