File: Cinfog2l.c

package info (click to toggle)
scalapack 1.6-13
  • links: PTS
  • area: main
  • in suites: potato
  • size: 30,476 kB
  • ctags: 25,789
  • sloc: fortran: 296,718; ansic: 51,265; makefile: 1,541; sh: 4
file content (42 lines) | stat: -rw-r--r-- 1,027 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
#include "tools.h"
void Cinfog2l(I, J, desc, nprow, npcol, myrow, mycol, ii, jj, prow, pcol)
int I;
int J;
int *desc;
int nprow;
int npcol;
int myrow;
int mycol;
int *ii;
int *jj;
int *prow;
int  *pcol;
{
   int nblocks, dist, extrablocks;
/*
 * Figure values for row
 */
   nblocks = I / desc[MB_];
   *prow = (desc[RSRC_] + nblocks) % nprow;
   dist = (nprow + myrow - desc[RSRC_]) % nprow; /* distance from RSRC */

   *ii = (nblocks/nprow) * desc[MB_];  /* all processes get >= this many */

   extrablocks = nblocks % nprow;
   if (dist < extrablocks) *ii += desc[MB_];
   else if (dist == extrablocks) *ii += I % desc[MB_];

/*
 * Figure values for column
 */
   nblocks = J / desc[NB_];
   *pcol = (desc[CSRC_] + nblocks) % npcol;
   dist = (npcol + mycol - desc[CSRC_]) % npcol; /* distance from CSRC */

   *jj = (nblocks/npcol) * desc[NB_];  /* all processes get >= this many */

   extrablocks = nblocks % npcol;
   if (dist < extrablocks) *jj += desc[NB_];
   else if (dist == extrablocks) *jj += J % desc[NB_];

}