1

当处理器网格不使用所有可用的 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将输入网格大小从2x更改2为-1x -1。在这种情况下，我正在使用11MPI 内核进行测试。

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

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

在带有 id 的进程上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 的一部分吗？

4

1 回答 1

0

除非有人有更好的来源，否则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 );
 }
于 2017-09-04T18:49:07.593 回答