#include "redist.h"
#include <stdio.h>
#include <stdlib.h>
#include <assert.h>

Data Structures
struct	MDESC
struct	IDESC

Macros
#define	static2 static
#define	fortran_mr2d pigemr2do
#define	fortran_mr2dnew pigemr2d
#define	icopy_ icopy
#define	ilacpy_ ilacpy
#define	Clacpy Cigelacpy
#define	BLOCK_CYCLIC_2D 1
#define	SHIFT(row, sprow, nbrow)
#define	max(A, B)
#define	min(A, B)
#define	DIVUP(a, b)
#define	ROUNDUP(a, b)
#define	scanD0 igescanD0
#define	dispmat igedispmat
#define	setmemory igesetmemory
#define	freememory igefreememory
#define	scan_intervals igescan_intervals
#define	SENDBUFF 0
#define	RECVBUFF 1
#define	SIZEBUFF 2
#define	NDEBUG

Functions
void	Cblacs_pcoord ()
Int	Cblacs_pnum ()
void	Csetpvmtids ()
void	Cblacs_get ()
void	Cblacs_pinfo ()
void	Cblacs_gridinfo ()
void	Cblacs_gridinit ()
void	Cblacs_exit ()
void	Cblacs_gridexit ()
void	Cblacs_setup ()
void	Cigebs2d ()
void	Cigebr2d ()
void	Cigesd2d ()
void	Cigerv2d ()
void	Cigsum2d ()
void	Cigamn2d ()
void	Cigamx2d ()
Int	localindice ()
void *	mr2d_malloc ()
Int	ppcm ()
Int	localsize ()
Int	memoryblocksize ()
Int	changeorigin ()
void	paramcheck ()
void	Cpigemr2do ()
void	Cpigemr2d ()
void *	mr2d_malloc (Int n)
Int	pgcd (Int a, Int b)
Int	ppcm (Int a, Int b)
Int	localsize (Int myprow, Int p, Int nbrow, Int m)
Int	memoryblocksize (MDESC *a)
void	checkequal (Int ctxt, Int a)
void	paramcheck (MDESC *a, Int i, Int j, Int m, Int n, Int p, Int q, Int gcontext)
Int	changeorigin (Int myp, Int sp, Int p, Int bs, Int i, Int decal, Int newsp)
Int	localindice (Int ig, Int jg, Int templateheight, Int templatewidth, MDESC *a)

Macro Definition Documentation

◆ BLOCK_CYCLIC_2D

#define BLOCK_CYCLIC_2D 1

Definition at line 36 of file pgemraux.c.

◆ Clacpy

#define Clacpy Cigelacpy

Definition at line 23 of file pgemraux.c.

◆ dispmat

#define dispmat igedispmat

Definition at line 83 of file pgemraux.c.

◆ DIVUP

#define DIVUP	(		a,
			b )

Value:

( ((a)-1) /(b)+1)

Definition at line 44 of file pgemraux.c.

◆ fortran_mr2d

#define fortran_mr2d pigemr2do

Definition at line 18 of file pgemraux.c.

◆ fortran_mr2dnew

#define fortran_mr2dnew pigemr2d

Definition at line 19 of file pgemraux.c.

◆ freememory

#define freememory igefreememory

Definition at line 85 of file pgemraux.c.

◆ icopy_

#define icopy_ icopy

Definition at line 20 of file pgemraux.c.

◆ ilacpy_

void ilacpy_ ilacpy

Definition at line 21 of file pgemraux.c.

◆ max

#define max	(		A,
			B )

Value:

((A)>(B)?(A):(B))

Definition at line 42 of file pgemraux.c.

◆ min

#define min	(		A,
			B )

Value:

((A)>(B)?(B):(A))

Definition at line 43 of file pgemraux.c.

◆ NDEBUG

#define NDEBUG

Definition at line 102 of file pgemraux.c.

◆ RECVBUFF

#define RECVBUFF 1

Definition at line 96 of file pgemraux.c.

◆ ROUNDUP

#define ROUNDUP	(		a,
			b )

Value:

(DIVUP(a,b)*(b))

DIVUP

#define DIVUP(a, b)

Definition pcgemr.c:182

Definition at line 45 of file pgemraux.c.

◆ scan_intervals

#define scan_intervals igescan_intervals

Definition at line 86 of file pgemraux.c.

◆ scanD0

#define scanD0 igescanD0

Definition at line 82 of file pgemraux.c.

◆ SENDBUFF

#define SENDBUFF 0

Definition at line 95 of file pgemraux.c.

◆ setmemory

#define setmemory igesetmemory

Definition at line 84 of file pgemraux.c.

◆ SHIFT

#define SHIFT	(	row,
		sprow,
		nbrow )

Value:

((row)-(sprow)+ ((row) >= (sprow) ? 0 : (nbrow)))

Definition at line 41 of file pgemraux.c.

◆ SIZEBUFF

#define SIZEBUFF 2

Definition at line 97 of file pgemraux.c.

◆ static2

#define static2 static

Definition at line 8 of file pgemraux.c.

Function Documentation

◆ Cblacs_exit()

void Cblacs_exit ( )

extern

◆ Cblacs_get()

void Cblacs_get ( )

extern

◆ Cblacs_gridexit()

void Cblacs_gridexit ( )

extern

◆ Cblacs_gridinfo()

void Cblacs_gridinfo ( )

extern

◆ Cblacs_gridinit()

void Cblacs_gridinit ( )

extern

◆ Cblacs_pcoord()

void Cblacs_pcoord ( )

extern

◆ Cblacs_pinfo()

void Cblacs_pinfo ( )

extern

◆ Cblacs_pnum()

Int Cblacs_pnum ( )

extern

◆ Cblacs_setup()

void Cblacs_setup ( )

extern

◆ changeorigin() [1/2]

Int changeorigin ( )

extern

◆ changeorigin() [2/2]

Int changeorigin	(	Int	myp,
		Int	sp,
		Int	p,
		Int	bs,
		Int	i,
		Int *	decal,
		Int *	newsp )

Definition at line 250 of file pgemraux.c.

{
  Int   tempheight, firstblock, firsttemp;
  /* we begin by changing the parameters so that ia < templatewidth,... */
  tempheight = bs * p;
  firsttemp = i / tempheight;
  firstblock = (i / bs) % p;
  *newsp = (sp + firstblock) % p;
  if (myp >= 0)
    *decal = firsttemp * bs + (SHIFT(myp, sp, p) < firstblock ? bs : 0);
  else
    *decal = 0;
  return i % bs;
}

◆ checkequal()

void checkequal	(	Int	ctxt,
		Int	a )

Definition at line 186 of file pgemraux.c.

{
  Int   np, dummy, nbrow, myp, b;
  Cblacs_gridinfo(ctxt, &nbrow, &np, &dummy, &myp);
  assert(nbrow == 1);
  if (np == 1)
    return;
  if (myp == 0) {
    Cigesd2d(ctxt, (Int)1, (Int)1, &a, (Int)1, (Int)0, (Int)1);
    Cigerv2d(ctxt, (Int)1, (Int)1, &b, (Int)1, (Int)0, np - 1);
    assert(a == b);
  } else {
    Cigerv2d(ctxt, (Int)1, (Int)1, &b, (Int)1, (Int)0, myp - 1);
    assert(a == b);
    Cigesd2d(ctxt, (Int)1, (Int)1, &a, (Int)1, (Int)0, (myp + 1) % np);
  }
}

◆ Cigamn2d()

void Cigamn2d ( )

extern

◆ Cigamx2d()

void Cigamx2d ( )

extern

◆ Cigebr2d()

void Cigebr2d ( )

extern

◆ Cigebs2d()

void Cigebs2d ( )

extern

◆ Cigerv2d()

void Cigerv2d ( )

extern

◆ Cigesd2d()

void Cigesd2d ( )

extern

◆ Cigsum2d()

void Cigsum2d ( )

extern

◆ Cpigemr2d()

void Cpigemr2d ( )

extern

◆ Cpigemr2do()

void Cpigemr2do ( )

extern

◆ Csetpvmtids()

void Csetpvmtids ( )

extern

◆ localindice() [1/2]

Int localindice ( )

extern

◆ localindice() [2/2]

Int localindice	(	Int	ig,
		Int	jg,
		Int	templateheight,
		Int	templatewidth,
		MDESC *	a )

Definition at line 267 of file pgemraux.c.

{
  Int   vtemp, htemp, vsubtemp, hsubtemp, il, jl;
  assert(ig >= 0 && ig < a->m && jg >= 0 && jg < a->n);
  /* coordinates in global matrix with the tests in intersect, ig MUST BE in
   * [0..m] and jg in [0..n] */
  /* coordinates of the template that "owns" the element */
  vtemp = ig / templateheight;
  htemp = jg / templatewidth;
  /* coordinates of the element in the subblock of the (vtemp, htemp)
   * template */
  vsubtemp = ig % a->nbrow;
  hsubtemp = jg % a->nbcol;
  /* coordinates of the element in the local block of the processor */
  il = a->nbrow * vtemp + vsubtemp;
  jl = a->nbcol * htemp + hsubtemp;
  assert(il < a->lda);
#ifndef NDEBUG
  {
    Int   pr, pc, p, q, lp, lq;
    Cblacs_gridinfo(a->ctxt, &p, &q, &pr, &pc);
    p = templateheight / a->nbrow;
    q = templatewidth / a->nbcol;
    lp = ig % templateheight / a->nbrow;
    lq = jg % templatewidth / a->nbcol;
    assert(lp == SHIFT(pr, a->sprow, p));
    assert(lq == SHIFT(pc, a->spcol, q));
  }
#endif
  return (jl * a->lda + il);
}

◆ localsize() [1/2]

Int localsize ( )

extern

◆ localsize() [2/2]

Int localsize	(	Int	myprow,
		Int	p,
		Int	nbrow,
		Int	m )

Definition at line 145 of file pgemraux.c.

{
  Int   templateheight, blockheight;
  templateheight = p * nbrow;
  if (m % templateheight != 0) {    /* not an exact boundary */
    if ((m % templateheight) > (nbrow * myprow)) {  /* processor
                             * (myprow,mypcol) has
                             * some elements in that
                             * incomplete template */
      if ((m % templateheight) >= (nbrow * (myprow + 1))) { /* processor
                                 * (myprow,mypcol)'s
                                 * part is complete */
    blockheight = (m / templateheight) * nbrow + nbrow;
      } else {  /* processor (myprow,mypcol)'s part is not complete */
    blockheight = (m / templateheight) * nbrow + (m % nbrow);
      };    /* if ((m%templateheight) > (nbrow*(myprow+1))) */
    } else {    /* processor (myprow,mypcol) has no element in that
         * incomplete template */
      blockheight = (m / templateheight) * nbrow;
    };  /* if ((m%templateheight) > (nbrow*myprow)) */
  } else {  /* exact boundary */
    blockheight = m / p;    /* (m/templateheight) * nbrow */
  };    /* if (m%templateheight !=0) */
  return blockheight;
}

◆ memoryblocksize() [1/2]

Int memoryblocksize ( )

extern

◆ memoryblocksize() [2/2]

Int memoryblocksize ( MDESC * a )

Definition at line 173 of file pgemraux.c.

{
  Int   myprow, mypcol, p, q;
  /* Compute the (myprow,mypcol) indices of processor mypnum in P0xQ0 We
   * assume the row-major ordering of the BLACS */
  Cblacs_gridinfo(a->ctxt, &p, &q, &myprow, &mypcol);
  myprow = SHIFT(myprow, a->sprow, p);
  mypcol = SHIFT(mypcol, a->spcol, q);
  assert(myprow >= 0 && mypcol >= 0);
  return localsize(myprow, p, a->nbrow, a->m) *
    localsize(mypcol, q, a->nbcol, a->n);
}

◆ mr2d_malloc() [1/2]

void * mr2d_malloc ( )

extern

◆ mr2d_malloc() [2/2]

void * mr2d_malloc ( Int n )

Definition at line 108 of file pgemraux.c.

{
  void *ptr;
  assert(n > 0);
  ptr = (void *) malloc(n);
  if (ptr == NULL) {
    fprintf(stderr, "xxmr2d:out of memory\n");
    exit(2);
  }
  return ptr;
}

◆ paramcheck() [1/2]

void paramcheck ( )

extern

◆ paramcheck() [2/2]

void paramcheck	(	MDESC *	a,
		Int	i,
		Int	j,
		Int	m,
		Int	n,
		Int	p,
		Int	q,
		Int	gcontext )

Definition at line 204 of file pgemraux.c.

{
  Int   p2, q2, myprow, mypcol;
#ifndef NDEBUG
  checkequal(gcontext, p);
  checkequal(gcontext, q);
  checkequal(gcontext, a->sprow);
  checkequal(gcontext, a->spcol);
  checkequal(gcontext, a->m);
  checkequal(gcontext, a->n);
  checkequal(gcontext, i);
  checkequal(gcontext, j);
  checkequal(gcontext, a->nbrow);
  checkequal(gcontext, a->nbcol);
#endif
  Cblacs_gridinfo(a->ctxt, &p2, &q2, &myprow, &mypcol);
  /* compatibility T3D, must check myprow  and mypcol are within bounds */
  if (myprow >= p2 || mypcol >= q2)
    myprow = mypcol = -1;
  if ((myprow >= 0 || mypcol >= 0) && (p2 != p && q2 != q)) {
    fprintf(stderr, "??MR2D:incoherent p,q parameters\n");
    exit(1);
  }
  assert(myprow < p && mypcol < q);
  if (a->sprow < 0 || a->sprow >= p || a->spcol < 0 || a->spcol >= q) {
    fprintf(stderr, "??MR2D:Bad first processor coordinates\n");
    exit(1);
  }
  if (i < 0 || j < 0 || i + m > a->m || j + n > a->n) {
    fprintf(stderr, "??MR2D:Bad submatrix:i=%d,j=%d,\
m=%d,n=%d,M=%d,N=%d\n",
        i, j, m, n, a->m, a->n);
    exit(1);
  }
  if ((myprow >= 0 || mypcol >= 0) &&
      localsize(SHIFT(myprow, a->sprow, p), p, a->nbrow, a->m) > a->lda) {
    fprintf(stderr, "??MR2D:bad lda arg:row=%d,m=%d,p=%d,\
nbrow=%d,lda=%d,sprow=%d\n",
        myprow, a->m, p, a->nbrow, a->lda, a->sprow);
    exit(1);
  }
}

◆ pgcd()

Int pgcd	(	Int	a,
		Int	b )

Definition at line 120 of file pgemraux.c.

{
  Int   aux;
  if (a < b)
    return pgcd(b, a);
  else {
    aux = a % b;
    if (aux == 0)
      return b;
    else
      return pgcd(b, aux);
  }
}

◆ ppcm() [1/2]

Int ppcm ( )

extern

◆ ppcm() [2/2]

Int ppcm	(	Int	a,
		Int	b )

Definition at line 134 of file pgemraux.c.

{
  Int   pg;
  pg = pgcd(a, b);
  return a * (b / pg);
}

OpenRadioss 2025.1.11 OpenRadioss project

Data Structures

Macros

Functions

Macro Definition Documentation

◆ BLOCK_CYCLIC_2D

◆ Clacpy

◆ dispmat

◆ DIVUP

◆ fortran_mr2d

◆ fortran_mr2dnew

◆ freememory

◆ icopy_

◆ ilacpy_

◆ max

◆ min

◆ NDEBUG

◆ RECVBUFF

◆ ROUNDUP

◆ scan_intervals

◆ scanD0

◆ SENDBUFF

◆ setmemory

◆ SHIFT

◆ SIZEBUFF

◆ static2

Function Documentation

◆ Cblacs_exit()

◆ Cblacs_get()

◆ Cblacs_gridexit()

◆ Cblacs_gridinfo()

◆ Cblacs_gridinit()

◆ Cblacs_pcoord()

◆ Cblacs_pinfo()

◆ Cblacs_pnum()

◆ Cblacs_setup()

◆ changeorigin() [1/2]

◆ changeorigin() [2/2]

◆ checkequal()

◆ Cigamn2d()

◆ Cigamx2d()

◆ Cigebr2d()

◆ Cigebs2d()

◆ Cigerv2d()

◆ Cigesd2d()

◆ Cigsum2d()

◆ Cpigemr2d()

◆ Cpigemr2do()

◆ Csetpvmtids()

◆ localindice() [1/2]

◆ localindice() [2/2]

◆ localsize() [1/2]

◆ localsize() [2/2]

◆ memoryblocksize() [1/2]

◆ memoryblocksize() [2/2]

◆ mr2d_malloc() [1/2]

◆ mr2d_malloc() [2/2]

◆ paramcheck() [1/2]

◆ paramcheck() [2/2]

◆ pgcd()

◆ ppcm() [1/2]

◆ ppcm() [2/2]