OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
pctrmrdrv.c File Reference
#include "redist.h"
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <ctype.h>
#include <assert.h>
#include <varargs.h>

Go to the source code of this file.

Data Structures

struct  complex
struct  MDESC
struct  IDESC

Macros

#define static2   static
#define fortran_mr2d   pctrmr2do
#define fortran_mr2dnew   pctrmr2d
#define ccopy_   ccopy
#define clacpy_   clacpy
#define Clacpy   Cctrlacpy
#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   ctrscanD0
#define dispmat   ctrdispmat
#define setmemory   ctrsetmemory
#define freememory   ctrfreememory
#define scan_intervals   ctrscan_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 Ccgesd2d ()
void Ccgerv2d ()
Int localindice ()
void * mr2d_malloc ()
Int ppcm ()
Int localsize ()
Int memoryblocksize ()
Int changeorigin ()
void paramcheck ()
void Cpctrmr2do ()
void Cpctrmr2d ()
static2 void initblock (complex *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 92 of file pctrmrdrv.c.

◆ ccopy_

#define ccopy_   ccopy

Definition at line 73 of file pctrmrdrv.c.

◆ Clacpy

#define Clacpy   Cctrlacpy

Definition at line 76 of file pctrmrdrv.c.

◆ clacpy_

#define clacpy_   clacpy

Definition at line 74 of file pctrmrdrv.c.

◆ dispmat

#define dispmat   ctrdispmat

Definition at line 139 of file pctrmrdrv.c.

◆ DIVUP

#define DIVUP ( a,
b )
Value:
( ((a)-1) /(b)+1)

Definition at line 100 of file pctrmrdrv.c.

◆ fortran_mr2d

#define fortran_mr2d   pctrmr2do

Definition at line 71 of file pctrmrdrv.c.

◆ fortran_mr2dnew

#define fortran_mr2dnew   pctrmr2d

Definition at line 72 of file pctrmrdrv.c.

◆ freememory

#define freememory   ctrfreememory

Definition at line 141 of file pctrmrdrv.c.

◆ max

#define max ( A,
B )
Value:
((A)>(B)?(A):(B))

Definition at line 98 of file pctrmrdrv.c.

◆ min

#define min ( A,
B )
Value:
((A)>(B)?(B):(A))

Definition at line 99 of file pctrmrdrv.c.

◆ nbparameter

#define nbparameter   24

◆ NDEBUG

#define NDEBUG

Definition at line 158 of file pctrmrdrv.c.

◆ RECVBUFF

#define RECVBUFF   1

Definition at line 152 of file pctrmrdrv.c.

◆ ROUNDUP

#define ROUNDUP ( a,
b )
Value:
(DIVUP(a,b)*(b))
DIVUP
#define DIVUP(a, b)
Definition pcgemr.c:182

Definition at line 101 of file pctrmrdrv.c.

◆ scan_intervals

#define scan_intervals   ctrscan_intervals

Definition at line 142 of file pctrmrdrv.c.

◆ scanD0

#define scanD0   ctrscanD0

Definition at line 138 of file pctrmrdrv.c.

◆ SENDBUFF

#define SENDBUFF   0

Definition at line 151 of file pctrmrdrv.c.

◆ setmemory

#define setmemory   ctrsetmemory

Definition at line 140 of file pctrmrdrv.c.

◆ SHIFT

#define SHIFT ( row,
sprow,
nbrow )
Value:
((row)-(sprow)+ ((row) >= (sprow) ? 0 : (nbrow)))

Definition at line 97 of file pctrmrdrv.c.

◆ SIZEBUFF

#define SIZEBUFF   2

Definition at line 153 of file pctrmrdrv.c.

◆ static2

#define static2   static

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

◆ Ccgerv2d()

void Ccgerv2d ( )
extern

◆ Ccgesd2d()

void Ccgesd2d ( )
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

◆ Cpctrmr2d()

void Cpctrmr2d ( )
extern

◆ Cpctrmr2do()

void Cpctrmr2do ( )
extern

◆ Csetpvmtids()

void Csetpvmtids ( )
extern

◆ getparam()

void getparam ( va_alist )
static

Definition at line 187 of file pctrmrdrv.c.

189{
190 FILE *f;
191#endif
192 va_list ap;
193 Int i;
194 static Int nbline;
195 char *ptr, *next;
196 Int *var;
197 static char buffer[200];
198#ifdef __STDC__
199 va_start(ap, f);
200#else
201 va_start(ap);
202 f = va_arg(ap, FILE *);
203#endif
204 do {
205 next = fgets(buffer, 200, f);
206 if (next == NULL) {
207 fprintf(stderr, "bad configuration driver file:after line %d\n", nbline);
208 exit(1);
209 }
210 nbline += 1;
211 } while (buffer[0] == '#');
212 ptr = buffer;
213 var = va_arg(ap, Int *);
214 while (var != NULL) {
215 *var = strtol(ptr, &next, 10);
216 if (ptr == next) {
217 fprintf(stderr, "bad configuration driver file:error line %d\n", nbline);
218 exit(1);
219 }
220 ptr = next;
221 var = va_arg(ap, Int *);
222 }
223 va_end(ap);
224}
Int
#define Int
Definition Bconfig.h:22
cmumps_buf::next
integer next
Definition cmumps_comm_buffer.F:39
lapack_testing.f
f
Definition lapack_testing.py:164

◆ initblock()

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

Definition at line 168 of file pctrmrdrv.c.

169{
170 complex *pdata;
171 Int i;
172 pdata = block;
173 for (i = 0; i < m * n; i++, pdata++) {
174 (*pdata).r = i;
175 };
176}
n
n
Definition schur_example.m:9

◆ initforpvm()

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

Definition at line 226 of file pctrmrdrv.c.

227{
228 Int pnum, nproc;
229 Cblacs_pinfo(&pnum, &nproc);
230 if (nproc < 1) { /* we are with PVM */
231 if (pnum == 0) {
232 if (argc < 2) {
233 fprintf(stderr, "usage with PVM:xctrmr nbproc\n\
234\t where nbproc is the number of nodes to initialize\n");
235 exit(1);
236 }
237 nproc = atoi(argv[1]);
238 }
239 Cblacs_setup(&pnum, &nproc);
240 }
241}
Cblacs_pinfo
void Cblacs_pinfo()
Cblacs_setup
void Cblacs_setup()

◆ localindice()

Int localindice ( )
extern

◆ localsize()

Int localsize ( )
extern

◆ main()

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

Definition at line 243 of file pctrmrdrv.c.

244{
245 /* We initialize the data-block on the current processor, then redistribute
246 * it, and perform the inverse redistribution to compare the local memory
247 * with the initial one. */
248 /* Data file */
249 FILE *fp;
250 Int nbre, nbremax;
251 /* Data distribution 0 parameters */
252 Int p0, /* # of rows in the processor grid */
253 q0; /* # of columns in the processor grid */
254 /* Data distribution 1 parameters */
255 Int p1, q1;
256 /* # of parameter to be read on the keyboard */
257#define nbparameter 24
258 /* General variables */
259 Int blocksize0;
260 Int mypnum, nprocs;
261 Int parameters[nbparameter], nberrors;
262 Int i;
263 Int ia, ja, ib, jb, m, n;
264 Int gcontext, context0, context1;
265 Int myprow1, myprow0, mypcol0, mypcol1;
266 Int dummy;
267 MDESC ma, mb;
268 char *uplo, *diag;
269 complex *ptrmyblock, *ptrsavemyblock, *ptrmyblockcopy, *ptrmyblockvide;
270#ifdef UsingMpiBlacs
271 MPI_Init(&argc, &argv);
272#endif
273 setvbuf(stdout, NULL, _IOLBF, 0);
274 setvbuf(stderr, NULL, _IOLBF, 0);
275#ifdef T3D
276 free(malloc(14000000));
277#endif
278 initforpvm(argc, argv);
279 /* Read physical parameters */
280 Cblacs_pinfo(&mypnum, &nprocs);
281 /* initialize BLACS for the parameter communication */
282 Cblacs_get((Int)0, (Int)0, &gcontext);
283 Cblacs_gridinit(&gcontext, "R", nprocs, (Int)1);
284 Cblacs_gridinfo(gcontext, &dummy, &dummy, &mypnum, &dummy);
285 if (mypnum == 0) {
286 if ((fp = fopen("TRMR2D.dat", "r")) == NULL) {
287 fprintf(stderr, "Can't open TRMR2D.dat\n");
288 exit(1);
289 };
290 printf("\n// CTRMR2D TESTER for COMPLEX //\n");
291 getparam(fp, &nbre, NULL);
292 printf("////////// %d tests \n\n", nbre);
293 parameters[0] = nbre;
294 Cigebs2d(gcontext, "All", "H", (Int)1, (Int)1, parameters, (Int)1);
295 } else {
296 Cigebr2d(gcontext, "All", "H", (Int)1, (Int)1, parameters, (Int)1, (Int)0, (Int)0);
297 nbre = parameters[0];
298 };
299 if (mypnum == 0) {
300 printf("\n m n m0 n0 sr0 sc0 i0 j0 p0 q0 nbr0 nbc0 \
301m1 n1 sr1 sc1 i1 j1 p1 q1 nbr1 nbc1\n\n");
302 };
303 /****** TEST LOOP *****/
304 /* Here we are in grip 1xnprocs */
305 nbremax = nbre;
306#ifdef DEBUG
307 fprintf(stderr, "bonjour,je suis le noeud %d\n", mypnum);
308#endif
309 while (nbre-- != 0) { /* Loop on the serie of tests */
310 /* All the processors read the parameters so we have to be in a 1xnprocs
311 * grid at each iteration */
312 /* Read processors grid and matrices parameters */
313 if (mypnum == 0) {
314 Int u, d;
315 getparam(fp,
316 &m, &n,
317 &ma.m, &ma.n, &ma.sprow, &ma.spcol,
318 &ia, &ja, &p0, &q0, &ma.nbrow, &ma.nbcol,
319 &mb.m, &mb.n, &mb.sprow, &mb.spcol,
320 &ib, &jb, &p1, &q1, &mb.nbrow, &mb.nbcol,
321 &u, &d,
322 NULL);
323 uplo = u ? "UPPER" : "LOWER";
324 diag = d ? "UNIT" : "NONUNIT";
325 printf("\t\t************* TEST # %d **********\n",
326 nbremax - nbre);
327 printf(" %3d %3d %3d %3d %3d %3d %3d %3d \
328%3d %3d %3d %3d %3d %3d %3d %3d %3d %3d %3d %3d %3d %3d",
329 m, n,
330 ma.m, ma.n, ma.sprow, ma.spcol,
331 ia, ja, p0, q0, ma.nbrow, ma.nbcol,
332 mb.m, mb.n, mb.sprow, mb.spcol,
333 ib, jb, p1, q1, mb.nbrow, mb.nbcol);
334 printf(" %s %s", toupper(*uplo) == 'U' ? "up" : "low",
335 toupper(*diag) == 'U' ? "unit" : "nonunit");
336 printf("\n");
337 if (p0 * q0 > nprocs || p1 * q1 > nprocs) {
338 fprintf(stderr, "not enough nodes:%d processors required\n",
339 max(p0 * q0, p1 * q1));
340 exit(1);
341 }
342 parameters[0] = p0;
343 parameters[1] = q0;
344 parameters[2] = ma.nbrow;
345 parameters[3] = ma.nbcol;
346 parameters[4] = p1;
347 parameters[5] = q1;
348 parameters[6] = mb.nbrow;
349 parameters[7] = mb.nbcol;
350 parameters[8] = ma.m;
351 parameters[9] = ma.n;
352 parameters[10] = ma.sprow;
353 parameters[11] = ma.spcol;
354 parameters[12] = mb.sprow;
355 parameters[13] = mb.spcol;
356 parameters[14] = ia;
357 parameters[15] = ja;
358 parameters[16] = ib;
359 parameters[17] = jb;
360 parameters[18] = m;
361 parameters[19] = n;
362 parameters[20] = mb.m;
363 parameters[21] = mb.n;
364 parameters[22] = *uplo == 'U';
365 parameters[23] = *diag == 'U';
366 Cigebs2d(gcontext, "All", "H", (Int)1, nbparameter, parameters, (Int)1);
367 } else {
368 Cigebr2d(gcontext, "All", "H", (Int)1, nbparameter, parameters, (Int)1, (Int)0, (Int)0);
369 p0 = parameters[0];
370 q0 = parameters[1];
371 ma.nbrow = parameters[2];
372 ma.nbcol = parameters[3];
373 p1 = parameters[4];
374 q1 = parameters[5];
375 mb.nbrow = parameters[6];
376 mb.nbcol = parameters[7];
377 ma.m = parameters[8];
378 ma.n = parameters[9];
379 ma.sprow = parameters[10];
380 ma.spcol = parameters[11];
381 mb.sprow = parameters[12];
382 mb.spcol = parameters[13];
383 ia = parameters[14];
384 ja = parameters[15];
385 ib = parameters[16];
386 jb = parameters[17];
387 m = parameters[18];
388 n = parameters[19];
389 mb.m = parameters[20];
390 mb.n = parameters[21];
391 ma.desctype = BLOCK_CYCLIC_2D;
392 mb.desctype = BLOCK_CYCLIC_2D;
393 uplo = parameters[22] ? "UPPER" : "LOWER";
394 diag = parameters[23] ? "UNIT" : "NONUNIT";
395 };
396 Cblacs_get((Int)0, (Int)0, &context0);
397 Cblacs_gridinit(&context0, "R", p0, q0);
398 Cblacs_get((Int)0, (Int)0, &context1);
399 Cblacs_gridinit(&context1, "R", p1, q1);
400 Cblacs_gridinfo(context0, &dummy, &dummy, &myprow0, &mypcol0);
401 if (myprow0 >= p0 || mypcol0 >= q0)
402 myprow0 = mypcol0 = -1;
403 Cblacs_gridinfo(context1, &dummy, &dummy, &myprow1, &mypcol1);
404 if (myprow1 >= p1 || mypcol1 >= q1)
405 myprow1 = mypcol1 = -1;
406 assert((myprow0 < p0 && mypcol0 < q0) || (myprow0 == -1 && mypcol0 == -1));
407 assert((myprow1 < p1 && mypcol1 < q1) || (myprow1 == -1 && mypcol1 == -1));
408 ma.ctxt = context0;
409 mb.ctxt = context1;
410 /* From here, we are not assuming that only the processors working in the
411 * redistribution are calling xxMR2D, but the ones not concerned will do
412 * nothing. */
413 /* We compute the exact size of the local memory block for the memory
414 * allocations */
415 if (myprow0 >= 0 && mypcol0 >= 0) {
416 blocksize0 = memoryblocksize(&ma);
417 ma.lda = localsize(SHIFT(myprow0, ma.sprow, p0), p0, ma.nbrow, ma.m);
418 setmemory(&ptrmyblock, blocksize0);
419 initblock(ptrmyblock, 1, blocksize0);
420 setmemory(&ptrmyblockcopy, blocksize0);
421 memcpy((char *) ptrmyblockcopy, (char *) ptrmyblock,
422 blocksize0 * sizeof(complex));
423 setmemory(&ptrmyblockvide, blocksize0);
424 for (i = 0; i < blocksize0; i++)
425 ptrmyblockvide[i].r = -1;
426 }; /* if (mypnum < p0 * q0) */
427 if (myprow1 >= 0 && mypcol1 >= 0) {
428 setmemory(&ptrsavemyblock, memoryblocksize(&mb));
429 mb.lda = localsize(SHIFT(myprow1, mb.sprow, p1), p1, mb.nbrow, mb.m);
430 }; /* if (mypnum < p1 * q1) */
431 /* Redistribute the matrix from grid 0 to grid 1 (memory location
432 * ptrmyblock to ptrsavemyblock) */
433 Cpctrmr2d(uplo, diag, m, n,
434 ptrmyblock, ia, ja, &ma,
435 ptrsavemyblock, ib, jb, &mb, gcontext);
436 /* Perform the inverse redistribution of the matrix from grid 1 to grid 0
437 * (memory location ptrsavemyblock to ptrmyblockvide) */
438 Cpctrmr2d(uplo, diag, m, n,
439 ptrsavemyblock, ib, jb, &mb,
440 ptrmyblockvide, ia, ja, &ma, gcontext);
441 /* Check the differences */
442 nberrors = 0;
443 if (myprow0 >= 0 && mypcol0 >= 0) {
444 /* only for the processors that do have data at the begining */
445 for (i = 0; i < blocksize0; i++) {
446 Int li, lj, gi, gj;
447 Int in;
448 in = 1;
449 li = i % ma.lda;
450 lj = i / ma.lda;
451 gi = (li / ma.nbrow) * p0 * ma.nbrow +
452 SHIFT(myprow0, ma.sprow, p0) * ma.nbrow + li % ma.nbrow;
453 gj = (lj / ma.nbcol) * q0 * ma.nbcol +
454 SHIFT(mypcol0, ma.spcol, q0) * ma.nbcol + lj % ma.nbcol;
455 assert(gi < ma.m && gj < ma.n);
456 gi -= (ia - 1);
457 gj -= (ja - 1);
458 if (gi < 0 || gj < 0 || gi >= m || gj >= n)
459 in = 0;
460 else if (toupper(*uplo) == 'U')
461 in = (gi <= gj + max(0, m - n) - (toupper(*diag) == 'U'));
462 else
463 in = (gi >= gj - max(0, n - m) + (toupper(*diag) == 'U'));
464 if (!in) {
465 ptrmyblockcopy[i].r = -1;
466 }
467 if (ptrmyblockvide[i].r != ptrmyblockcopy[i].r) {
468 nberrors++;
469 printf("Proc %d : Error element number %d, value = %f , initvalue =%f \n"
470 ,mypnum, i,
471 ptrmyblockvide[i].r, ptrmyblockcopy[i].r);
472 };
473 };
474 if (nberrors > 0) {
475 printf("Processor %d, has tested %d COMPLEX elements,\
476Number of redistribution errors = %d \n",
477 mypnum, blocksize0, nberrors);
478 }
479 }
480 /* Look at the errors on all the processors at this point. */
481 Cigsum2d(gcontext, "All", "H", (Int)1, (Int)1, &nberrors, (Int)1, (Int)0, (Int)0);
482 if (mypnum == 0)
483 if (nberrors)
484 printf(" => Total number of redistribution errors = %d \n",
485 nberrors);
486 else
487 printf("TEST PASSED OK\n");
488 /* release memory for the next iteration */
489 if (myprow0 >= 0 && mypcol0 >= 0) {
490 freememory((char *) ptrmyblock);
491 freememory((char *) ptrmyblockvide);
492 freememory((char *) ptrmyblockcopy);
493 }; /* if (mypnum < p0 * q0) */
494 /* release memory for the next iteration */
495 if (myprow1 >= 0 && mypcol1 >= 0) {
496 freememory((char *) ptrsavemyblock);
497 };
498 if (myprow0 >= 0)
499 Cblacs_gridexit(context0);
500 if (myprow1 >= 0)
501 Cblacs_gridexit(context1);
502 }; /* while nbre != 0 */
503 if (mypnum == 0) {
504 fclose(fp);
505 };
506 Cblacs_exit((Int)0);
507 return 0;
508}/* main */
BLOCK_CYCLIC_2D
#define BLOCK_CYCLIC_2D
Definition PBtools.h:20
MPI_Init
LIBSEQ_INT LIBSEQ_CALL MPI_Init(LIBSEQ_INT *pargc, char ***pargv)
Definition mpic.c:18
cmumps_load::nprocs
integer, save, private nprocs
Definition cmumps_load.F:57
SHIFT
#define SHIFT(row, sprow, nbrow)
Definition pcgemr.c:179
nbparameter
#define nbparameter
memoryblocksize
Int memoryblocksize()
initblock
static2 void initblock(complex *block, Int m, Int n)
Definition pctrmrdrv.c:168
freememory
#define freememory
Definition pctrmrdrv.c:141
SHIFT
#define SHIFT(row, sprow, nbrow)
Definition pctrmrdrv.c:97
Cblacs_gridexit
void Cblacs_gridexit()
max
#define max(A, B)
Definition pctrmrdrv.c:98
Cigsum2d
void Cigsum2d()
Cpctrmr2d
void Cpctrmr2d()
localsize
Int localsize()
setmemory
#define setmemory
Definition pctrmrdrv.c:140
Cblacs_get
void Cblacs_get()
Cigebs2d
void Cigebs2d()
getparam
static void getparam(va_alist)
Definition pctrmrdrv.c:187
Cblacs_gridinit
void Cblacs_gridinit()
Cblacs_gridinfo
void Cblacs_gridinfo()
initforpvm
void initforpvm(Int argc, char *argv[])
Definition pctrmrdrv.c:226
Cigebr2d
void Cigebr2d()
Cblacs_exit
void Cblacs_exit()
r
r
Definition schur_example.m:55
MDESC::m
Int m
Definition pcgemr.c:166
MDESC::spcol
Int spcol
Definition pcgemr.c:171
MDESC::nbcol
Int nbcol
Definition pcgemr.c:169
MDESC::sprow
Int sprow
Definition pcgemr.c:170
MDESC::nbrow
Int nbrow
Definition pcgemr.c:168
MDESC::ctxt
Int ctxt
Definition pcgemr.c:165
MDESC::desctype
Int desctype
Definition pcgemr.c:164
MDESC::n
Int n
Definition pcgemr.c:167
MDESC::lda
Int lda
Definition pcgemr.c:172

◆ memoryblocksize()

Int memoryblocksize ( )
extern

◆ mr2d_malloc()

void * mr2d_malloc ( )
extern

◆ paramcheck()

void paramcheck ( )
extern

◆ ppcm()

Int ppcm ( )
extern