#include "redist.h"
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <ctype.h>
#include <assert.h>
#include <varargs.h>

Data Structures
struct	MDESC
struct	IDESC

Macros
#define	static2 static
#define	fortran_mr2d pdtrmr2do
#define	fortran_mr2dnew pdtrmr2d
#define	dcopy_ dcopy
#define	dlacpy_ dlacpy
#define	Clacpy Cdtrlacpy
#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 dtrscanD0
#define	dispmat dtrdispmat
#define	setmemory dtrsetmemory
#define	freememory dtrfreememory
#define	scan_intervals dtrscan_intervals
#define	SENDBUFF 0
#define	RECVBUFF 1
#define	SIZEBUFF 2
#define	NDEBUG
#define	nbparameter 24

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 ()
void	Cdgesd2d ()
void	Cdgerv2d ()
Int	localindice ()
void *	mr2d_malloc ()
Int	ppcm ()
Int	localsize ()
Int	memoryblocksize ()
Int	changeorigin ()
void	paramcheck ()
void	Cpdtrmr2do ()
void	Cpdtrmr2d ()
static2 void	initblock (double *block, Int m, Int n)
static void	getparam (va_alist)
void	initforpvm (Int argc, char *argv[])
int	main (int argc, char *argv[])

Macro Definition Documentation

◆ BLOCK_CYCLIC_2D

#define BLOCK_CYCLIC_2D 1

Definition at line 89 of file pdtrmrdrv.c.

◆ Clacpy

#define Clacpy Cdtrlacpy

Definition at line 76 of file pdtrmrdrv.c.

◆ dcopy_

#define dcopy_ dcopy

Definition at line 73 of file pdtrmrdrv.c.

◆ dispmat

#define dispmat dtrdispmat

Definition at line 136 of file pdtrmrdrv.c.

◆ DIVUP

#define DIVUP	(		a,
			b )

Value:

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

Definition at line 97 of file pdtrmrdrv.c.

◆ dlacpy_

#define dlacpy_ dlacpy

Definition at line 74 of file pdtrmrdrv.c.

◆ fortran_mr2d

#define fortran_mr2d pdtrmr2do

Definition at line 71 of file pdtrmrdrv.c.

◆ fortran_mr2dnew

#define fortran_mr2dnew pdtrmr2d

Definition at line 72 of file pdtrmrdrv.c.

◆ freememory

#define freememory dtrfreememory

Definition at line 138 of file pdtrmrdrv.c.

◆ max

#define max	(		A,
			B )

Value:

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

Definition at line 95 of file pdtrmrdrv.c.

◆ min

#define min	(		A,
			B )

Value:

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

Definition at line 96 of file pdtrmrdrv.c.

◆ nbparameter

#define nbparameter 24

◆ NDEBUG

#define NDEBUG

Definition at line 155 of file pdtrmrdrv.c.

◆ RECVBUFF

#define RECVBUFF 1

Definition at line 149 of file pdtrmrdrv.c.

◆ ROUNDUP

#define ROUNDUP	(		a,
			b )

Value:

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

DIVUP

#define DIVUP(a, b)

Definition pcgemr.c:182

Definition at line 98 of file pdtrmrdrv.c.

◆ scan_intervals

#define scan_intervals dtrscan_intervals

Definition at line 139 of file pdtrmrdrv.c.

◆ scanD0

#define scanD0 dtrscanD0

Definition at line 135 of file pdtrmrdrv.c.

◆ SENDBUFF

#define SENDBUFF 0

Definition at line 148 of file pdtrmrdrv.c.

◆ setmemory

#define setmemory dtrsetmemory

Definition at line 137 of file pdtrmrdrv.c.

◆ SHIFT

#define SHIFT	(	row,
		sprow,
		nbrow )

Value:

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

Definition at line 94 of file pdtrmrdrv.c.

◆ SIZEBUFF

#define SIZEBUFF 2

Definition at line 150 of file pdtrmrdrv.c.

◆ static2

#define static2 static

Definition at line 61 of file pdtrmrdrv.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

◆ Cdgerv2d()

void Cdgerv2d ( )

extern

◆ Cdgesd2d()

void Cdgesd2d ( )

extern

◆ changeorigin()

Int changeorigin ( )

extern

◆ 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

◆ Cpdtrmr2d()

void Cpdtrmr2d ( )

extern

◆ Cpdtrmr2do()

void Cpdtrmr2do ( )

extern

◆ Csetpvmtids()

void Csetpvmtids ( )

extern

◆ getparam()

void getparam ( va_alist )

static

Definition at line 184 of file pdtrmrdrv.c.

{
  FILE *f;
#endif
  va_list ap;
  Int   i;
  static Int nbline;
  char *ptr, *next;
  Int  *var;
  static char buffer[200];
#ifdef __STDC__
  va_start(ap, f);
#else
  va_start(ap);
  f = va_arg(ap, FILE *);
#endif
  do {
    next = fgets(buffer, 200, f);
    if (next == NULL) {
      fprintf(stderr, "bad configuration driver file:after line %d\n", nbline);
      exit(1);
    }
    nbline += 1;
  } while (buffer[0] == '#');
  ptr = buffer;
  var = va_arg(ap, Int *);
  while (var != NULL) {
    *var = strtol(ptr, &next, 10);
    if (ptr == next) {
      fprintf(stderr, "bad configuration driver file:error line %d\n", nbline);
      exit(1);
    }
    ptr = next;
    var = va_arg(ap, Int *);
  }
  va_end(ap);
}

◆ initblock()

static2 void initblock	(	double *	block,
		Int	m,
		Int	n )

Definition at line 165 of file pdtrmrdrv.c.

{
  double *pdata;
  Int   i;
  pdata = block;
  for (i = 0; i < m * n; i++, pdata++) {
    (*pdata) = i;
  };
}

◆ initforpvm()

void initforpvm	(	Int	argc,
		char *	argv[] )

Definition at line 223 of file pdtrmrdrv.c.

{
  Int   pnum, nproc;
  Cblacs_pinfo(&pnum, &nproc);
  if (nproc < 1) {  /* we are with PVM */
    if (pnum == 0) {
      if (argc < 2) {
    fprintf(stderr, "usage with PVM:xdtrmr nbproc\n\
\t where nbproc is the number of nodes to initialize\n");
    exit(1);
      }
      nproc = atoi(argv[1]);
    }
    Cblacs_setup(&pnum, &nproc);
  }
}

◆ localindice()

Int localindice ( )

extern

◆ localsize()

Int localsize ( )

extern

◆ main()

int main	(	int	argc,
		char *	argv[] )

Definition at line 240 of file pdtrmrdrv.c.

{
  /* We initialize the data-block on the current processor, then redistribute
   * it, and perform the inverse redistribution  to compare the local memory
   * with the initial one. */
  /* Data file */
  FILE *fp;
  Int   nbre, nbremax;
  /* Data distribution 0 parameters */
  Int   p0,  /* # of rows in the processor grid */
        q0; /* # of columns in the processor grid */
  /* Data distribution 1 parameters */
  Int   p1, q1;
  /* # of parameter to be read on the keyboard */
#define nbparameter 24
  /* General variables */
  Int   blocksize0;
  Int   mypnum, nprocs;
  Int   parameters[nbparameter], nberrors;
  Int   i;
  Int   ia, ja, ib, jb, m, n;
  Int   gcontext, context0, context1;
  Int   myprow1, myprow0, mypcol0, mypcol1;
  Int   dummy;
  MDESC ma, mb;
  char *uplo, *diag;
  double *ptrmyblock, *ptrsavemyblock, *ptrmyblockcopy, *ptrmyblockvide;
#ifdef UsingMpiBlacs
   MPI_Init(&argc, &argv);
#endif
  setvbuf(stdout, NULL, _IOLBF, 0);
  setvbuf(stderr, NULL, _IOLBF, 0);
#ifdef T3D
  free(malloc(14000000));
#endif
  initforpvm(argc, argv);
  /* Read physical parameters */
  Cblacs_pinfo(&mypnum, &nprocs);
  /* initialize BLACS for the parameter communication */
  Cblacs_get((Int)0, (Int)0, &gcontext);
  Cblacs_gridinit(&gcontext, "R", nprocs, (Int)1);
  Cblacs_gridinfo(gcontext, &dummy, &dummy, &mypnum, &dummy);
  if (mypnum == 0) {
    if ((fp = fopen("TRMR2D.dat", "r")) == NULL) {
      fprintf(stderr, "Can't open TRMR2D.dat\n");
      exit(1);
    };
    printf("\n// DTRMR2D TESTER for DOUBLE PRECISION //\n");
    getparam(fp, &nbre, NULL);
    printf("////////// %d tests \n\n", nbre);
    parameters[0] = nbre;
    Cigebs2d(gcontext, "All", "H", (Int)1, (Int)1, parameters, (Int)1);
  } else {
    Cigebr2d(gcontext, "All", "H", (Int)1, (Int)1, parameters, (Int)1, (Int)0, (Int)0);
    nbre = parameters[0];
  };
  if (mypnum == 0) {
    printf("\n  m   n  m0  n0  sr0 sc0 i0  j0  p0  q0 nbr0 nbc0 \
m1  n1  sr1 sc1 i1  j1  p1  q1 nbr1 nbc1\n\n");
  };
  /****** TEST LOOP *****/
  /* Here we are in grip 1xnprocs */
  nbremax = nbre;
#ifdef DEBUG
  fprintf(stderr, "bonjour,je suis le noeud %d\n", mypnum);
#endif
  while (nbre-- != 0) { /* Loop on the serie of tests */
    /* All the processors read the parameters so we have to be in a 1xnprocs
     * grid at each iteration */
    /* Read processors grid and matrices parameters */
    if (mypnum == 0) {
      Int   u, d;
      getparam(fp,
           &m, &n,
           &ma.m, &ma.n, &ma.sprow, &ma.spcol,
           &ia, &ja, &p0, &q0, &ma.nbrow, &ma.nbcol,
           &mb.m, &mb.n, &mb.sprow, &mb.spcol,
           &ib, &jb, &p1, &q1, &mb.nbrow, &mb.nbcol,
           &u, &d,
           NULL);
      uplo = u ? "UPPER" : "LOWER";
      diag = d ? "UNIT" : "NONUNIT";
      printf("\t\t************* TEST # %d **********\n",
         nbremax - nbre);
      printf(" %3d %3d %3d %3d %3d %3d %3d %3d \
%3d %3d %3d %3d %3d %3d %3d %3d %3d %3d %3d %3d %3d %3d",
         m, n,
         ma.m, ma.n, ma.sprow, ma.spcol,
         ia, ja, p0, q0, ma.nbrow, ma.nbcol,
         mb.m, mb.n, mb.sprow, mb.spcol,
         ib, jb, p1, q1, mb.nbrow, mb.nbcol);
      printf(" %s %s", toupper(*uplo) == 'U' ? "up" : "low",
         toupper(*diag) == 'U' ? "unit" : "nonunit");
      printf("\n");
      if (p0 * q0 > nprocs || p1 * q1 > nprocs) {
    fprintf(stderr, "not enough nodes:%d processors required\n",
        max(p0 * q0, p1 * q1));
    exit(1);
      }
      parameters[0] = p0;
      parameters[1] = q0;
      parameters[2] = ma.nbrow;
      parameters[3] = ma.nbcol;
      parameters[4] = p1;
      parameters[5] = q1;
      parameters[6] = mb.nbrow;
      parameters[7] = mb.nbcol;
      parameters[8] = ma.m;
      parameters[9] = ma.n;
      parameters[10] = ma.sprow;
      parameters[11] = ma.spcol;
      parameters[12] = mb.sprow;
      parameters[13] = mb.spcol;
      parameters[14] = ia;
      parameters[15] = ja;
      parameters[16] = ib;
      parameters[17] = jb;
      parameters[18] = m;
      parameters[19] = n;
      parameters[20] = mb.m;
      parameters[21] = mb.n;
      parameters[22] = *uplo == 'U';
      parameters[23] = *diag == 'U';
      Cigebs2d(gcontext, "All", "H", (Int)1, nbparameter, parameters, (Int)1);
    } else {
      Cigebr2d(gcontext, "All", "H", (Int)1, nbparameter, parameters, (Int)1, (Int)0, (Int)0);
      p0 = parameters[0];
      q0 = parameters[1];
      ma.nbrow = parameters[2];
      ma.nbcol = parameters[3];
      p1 = parameters[4];
      q1 = parameters[5];
      mb.nbrow = parameters[6];
      mb.nbcol = parameters[7];
      ma.m = parameters[8];
      ma.n = parameters[9];
      ma.sprow = parameters[10];
      ma.spcol = parameters[11];
      mb.sprow = parameters[12];
      mb.spcol = parameters[13];
      ia = parameters[14];
      ja = parameters[15];
      ib = parameters[16];
      jb = parameters[17];
      m = parameters[18];
      n = parameters[19];
      mb.m = parameters[20];
      mb.n = parameters[21];
      ma.desctype = BLOCK_CYCLIC_2D;
      mb.desctype = BLOCK_CYCLIC_2D;
      uplo = parameters[22] ? "UPPER" : "LOWER";
      diag = parameters[23] ? "UNIT" : "NONUNIT";
    };
    Cblacs_get((Int)0, (Int)0, &context0);
    Cblacs_gridinit(&context0, "R", p0, q0);
    Cblacs_get((Int)0, (Int)0, &context1);
    Cblacs_gridinit(&context1, "R", p1, q1);
    Cblacs_gridinfo(context0, &dummy, &dummy, &myprow0, &mypcol0);
    if (myprow0 >= p0 || mypcol0 >= q0)
      myprow0 = mypcol0 = -1;
    Cblacs_gridinfo(context1, &dummy, &dummy, &myprow1, &mypcol1);
    if (myprow1 >= p1 || mypcol1 >= q1)
      myprow1 = mypcol1 = -1;
    assert((myprow0 < p0 && mypcol0 < q0) || (myprow0 == -1 && mypcol0 == -1));
    assert((myprow1 < p1 && mypcol1 < q1) || (myprow1 == -1 && mypcol1 == -1));
    ma.ctxt = context0;
    mb.ctxt = context1;
    /* From here, we are not assuming that only the processors working in the
     * redistribution are calling  xxMR2D, but the ones not concerned will do
     * nothing. */
    /* We compute the exact size of the local memory block for the memory
     * allocations */
    if (myprow0 >= 0 && mypcol0 >= 0) {
      blocksize0 = memoryblocksize(&ma);
      ma.lda = localsize(SHIFT(myprow0, ma.sprow, p0), p0, ma.nbrow, ma.m);
      setmemory(&ptrmyblock, blocksize0);
      initblock(ptrmyblock, 1, blocksize0);
      setmemory(&ptrmyblockcopy, blocksize0);
      memcpy((char *) ptrmyblockcopy, (char *) ptrmyblock,
         blocksize0 * sizeof(double));
      setmemory(&ptrmyblockvide, blocksize0);
      for (i = 0; i < blocksize0; i++)
    ptrmyblockvide[i] = -1;
    };  /* if (mypnum < p0 * q0) */
    if (myprow1 >= 0 && mypcol1 >= 0) {
      setmemory(&ptrsavemyblock, memoryblocksize(&mb));
      mb.lda = localsize(SHIFT(myprow1, mb.sprow, p1), p1, mb.nbrow, mb.m);
    };  /* if (mypnum < p1 * q1)  */
    /* Redistribute the matrix from grid 0 to grid 1 (memory location
     * ptrmyblock to ptrsavemyblock) */
    Cpdtrmr2d(uplo, diag, m, n,
          ptrmyblock, ia, ja, &ma,
          ptrsavemyblock, ib, jb, &mb, gcontext);
    /* Perform the inverse redistribution of the matrix from grid 1 to grid 0
     * (memory location ptrsavemyblock to ptrmyblockvide) */
    Cpdtrmr2d(uplo, diag, m, n,
          ptrsavemyblock, ib, jb, &mb,
          ptrmyblockvide, ia, ja, &ma, gcontext);
    /* Check the differences */
    nberrors = 0;
    if (myprow0 >= 0 && mypcol0 >= 0) {
      /* only for the processors that do have data at the begining */
      for (i = 0; i < blocksize0; i++) {
    Int   li, lj, gi, gj;
    Int   in;
    in = 1;
    li = i % ma.lda;
    lj = i / ma.lda;
    gi = (li / ma.nbrow) * p0 * ma.nbrow +
          SHIFT(myprow0, ma.sprow, p0) * ma.nbrow + li % ma.nbrow;
    gj = (lj / ma.nbcol) * q0 * ma.nbcol +
          SHIFT(mypcol0, ma.spcol, q0) * ma.nbcol + lj % ma.nbcol;
    assert(gi < ma.m && gj < ma.n);
    gi -= (ia - 1);
    gj -= (ja - 1);
    if (gi < 0 || gj < 0 || gi >= m || gj >= n)
      in = 0;
    else if (toupper(*uplo) == 'U')
      in = (gi <= gj + max(0, m - n) - (toupper(*diag) == 'U'));
    else
      in = (gi >= gj - max(0, n - m) + (toupper(*diag) == 'U'));
    if (!in) {
      ptrmyblockcopy[i] = -1;
    }
    if (ptrmyblockvide[i] != ptrmyblockcopy[i]) {
      nberrors++;
      printf("Proc %d : Error element number %d, value = %f , initvalue =%f \n"
         ,mypnum, i,
         ptrmyblockvide[i], ptrmyblockcopy[i]);
    };
      };
      if (nberrors > 0) {
    printf("Processor %d, has tested  %d DOUBLE PRECISION elements,\
Number of redistribution errors = %d \n",
           mypnum, blocksize0, nberrors);
      }
    }
    /* Look at the errors on all the processors at this point. */
    Cigsum2d(gcontext, "All", "H", (Int)1, (Int)1, &nberrors, (Int)1, (Int)0, (Int)0);
    if (mypnum == 0)
      if (nberrors)
    printf("  => Total number of redistribution errors = %d \n",
           nberrors);
      else
    printf("TEST PASSED OK\n");
    /* release memory for the next iteration */
    if (myprow0 >= 0 && mypcol0 >= 0) {
      freememory((char *) ptrmyblock);
      freememory((char *) ptrmyblockvide);
      freememory((char *) ptrmyblockcopy);
    };  /* if (mypnum < p0 * q0) */
    /* release memory for the next iteration */
    if (myprow1 >= 0 && mypcol1 >= 0) {
      freememory((char *) ptrsavemyblock);
    };
    if (myprow0 >= 0)
      Cblacs_gridexit(context0);
    if (myprow1 >= 0)
      Cblacs_gridexit(context1);
  };    /* while nbre != 0 */
  if (mypnum == 0) {
    fclose(fp);
  };
  Cblacs_exit((Int)0);
  return 0;
}/* main */

◆ memoryblocksize()

Int memoryblocksize ( )

extern

◆ mr2d_malloc()

void * mr2d_malloc ( )

extern

◆ paramcheck()

void paramcheck ( )

extern

◆ ppcm()

Int ppcm ( )

extern

OpenRadioss 2025.1.11 OpenRadioss project

Data Structures

Macros

Functions

Macro Definition Documentation

◆ BLOCK_CYCLIC_2D

◆ Clacpy

◆ dcopy_

◆ dispmat

◆ DIVUP

◆ dlacpy_

◆ fortran_mr2d

◆ fortran_mr2dnew

◆ freememory

◆ max

◆ min

◆ nbparameter

◆ 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()

◆ Cdgerv2d()

◆ Cdgesd2d()

◆ changeorigin()

◆ Cigamn2d()

◆ Cigamx2d()

◆ Cigebr2d()

◆ Cigebs2d()

◆ Cigerv2d()

◆ Cigesd2d()

◆ Cigsum2d()

◆ Cpdtrmr2d()

◆ Cpdtrmr2do()

◆ Csetpvmtids()

◆ getparam()

◆ initblock()

◆ initforpvm()

◆ localindice()

◆ localsize()

◆ main()

◆ memoryblocksize()

◆ mr2d_malloc()

◆ paramcheck()

◆ ppcm()