Blacs:不使用所有MPI核心的处理器网格

Question

当处理器网格没有使用所有可用的MPI核心时，我对Scalapack有问题。假设我知道设置Pc*Pr<=NP的网格

当我设置一个进程网格时

Cblacs_get(0, 0, &ctxt);
const char *order = ( column_major ? "Col" : "Row" );
Cblacs_gridinit(&ctxt, order, procrows, proccols);
Cblacs_gridinfo( ctxt, &procrows, &proccols, &myrow, &mycol );

Cblacs_gridinfo将输入网格大小从2x2更改为-1x-1。在本例中，我使用11 MPI内核进行测试。

这本身并不会引发任何错误，但是当我尝试设置描述符向量时。

int irsrc = 0, icsrc = 0;
descinit_(descA, &M, &N, &Mb, &Nb,&irsrc,&icsrc,&ctxt, &lda, &info);

对于i为4-10的进程，我得到了错误

{   -1,   -1}:  On entry to DESCINIT parameter number    6 had an illegal value

问题：

处理不在网格上的MPI核心的正确方法是什么？我要跳过所有的核心吗？

Cblacs_gridinfo( ctxt, &procrows, &proccols, &myrow, &mycol );

myrow=mycol=-1回来了？这是API的一部分吗？

Answer 1

除非有人有更好的源代码，否则http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=139中的第二篇文章提供了一个C示例，它只对以下进程调用ScaLAPACK函数

 if ((myrow>-1)&(mycol>-1)&(myrow<nprow)&(mycol<npcol)) {
    // do the job
    Cblacs_gridexit( 0 );
 }