在MPI中使用Cblacs初始化网格时出现问题

如何解决在MPI中使用Cblacs初始化网格时出现问题

我正在尝试使用以下代码建立一个非常简单的1 * 2网格：

        int nprow,npcol,myrow,mycol,myid;
        char rowcol[1] = "R";

        nprow = 1;
        npcol = size / nprow;
        if(npcol * nprow != size){
                printf("Error");
                MPI_Finalize();
                exit(1);
        }

        Cblacs_pinfo(&myid,&size);
        Cblacs_get(0,&ictxt);
        Cblacs_gridinit(&ictxt,rowcol,nprow,npcol);
        Cblacs_pcoord(ictxt,myid,&myrow,&mycol);

        printf("rank = %d,nprow = %d,npcol = %d,myrow = %d,mycol = %d\n",rank,mycol);                                                                    }

问题在于Cblacs_pcoord函数似乎将nprow更改为0，无论其初始设置是什么，这反过来又为每个myrow赋予了0，而npcol和mycol变量始终正确使用的任何数量的处理器。我很困惑，因为该函数不应碰到nprow，但是我在每一行代码之后都打印了nprow，它是正确的值，直到调用该函数之后。

如果我缺少任何可以帮助您回答问题的信息，请告诉我，我会相应地进行更新。