257#define nbparameter 24
263 Int ia, ja, ib, jb, m, n;
264 Int gcontext, context0, context1;
265 Int myprow1, myprow0, mypcol0, mypcol1;
268 dcomplex *ptrmyblock, *ptrsavemyblock, *ptrmyblockcopy, *ptrmyblockvide;
270 MPI_Init(&argc, &argv);
272 setvbuf(stdout, NULL, _IOLBF, 0);
273 setvbuf(stderr, NULL, _IOLBF, 0);
275 free(malloc(14000000));
285 if ((fp = fopen(
"GEMR2D.dat",
"r")) == NULL) {
286 fprintf(stderr,
"Can't open GEMR2D.dat\n");
289 printf(
"\n// ZGEMR2D TESTER for COMPLEX*16 //\n");
290 getparam(fp, &nbre, NULL);
291 printf(
"////////// %d tests \n\n", nbre);
292 parameters[0] = nbre;
296 nbre = parameters[0];
299 printf(
"\n m n m0 n0 sr0 sc0 i0 j0 p0 q0 nbr0 nbc0 \
300m1 n1 sr1 sc1 i1 j1 p1 q1 nbr1 nbc1\n\n");
306 fprintf(stderr,
"bonjour,je suis le noeud %d\n", mypnum);
308 while (nbre-- != 0) {
321 printf(
"\t\t************* TEST # %d **********\n",
323 printf(
" %3d %3d %3d %3d %3d %3d %3d %3d \
324%3d %3d %3d %3d %3d %3d %3d %3d %3d %3d %3d %3d %3d %3d",
331 if (p0 * q0 > nprocs || p1 * q1 > nprocs) {
332 fprintf(stderr,
"not enough nodes:%d processors required\n",
333 max(p0 * q0, p1 * q1));
338 parameters[2] = ma.
nbrow;
339 parameters[3] = ma.
nbcol;
342 parameters[6] = mb.
nbrow;
343 parameters[7] = mb.
nbcol;
344 parameters[8] = ma.
m;
345 parameters[9] = ma.
n;
346 parameters[10] = ma.
sprow;
347 parameters[11] = ma.
spcol;
348 parameters[12] = mb.
sprow;
349 parameters[13] = mb.
spcol;
356 parameters[20] = mb.
m;
357 parameters[21] = mb.
n;
363 ma.
nbrow = parameters[2];
364 ma.
nbcol = parameters[3];
367 mb.
nbrow = parameters[6];
368 mb.
nbcol = parameters[7];
369 ma.
m = parameters[8];
370 ma.
n = parameters[9];
371 ma.
sprow = parameters[10];
372 ma.
spcol = parameters[11];
373 mb.
sprow = parameters[12];
374 mb.
spcol = parameters[13];
381 mb.
m = parameters[20];
382 mb.
n = parameters[21];
391 if (myprow0 >= p0 || mypcol0 >= q0)
392 myprow0 = mypcol0 = -1;
394 if (myprow1 >= p1 || mypcol1 >= q1)
395 myprow1 = mypcol1 = -1;
396 assert((myprow0 < p0 && mypcol0 < q0) || (myprow0 == -1 && mypcol0 == -1));
397 assert((myprow1 < p1 && mypcol1 < q1) || (myprow1 == -1 && mypcol1 == -1));
405 if (myprow0 >= 0 && mypcol0 >= 0) {
411 memcpy((
char *) ptrmyblockcopy, (
char *) ptrmyblock,
414 for (i = 0; i < blocksize0; i++)
415 ptrmyblockvide[i].r = -1;
417 if (myprow1 >= 0 && mypcol1 >= 0) {
424 ptrmyblock, ia, ja, &ma,
425 ptrsavemyblock, ib, jb, &mb, gcontext);
429 ptrsavemyblock, ib, jb, &mb,
430 ptrmyblockvide, ia, ja, &ma, gcontext);
433 if (myprow0 >= 0 && mypcol0 >= 0) {
435 for (i = 0; i < blocksize0; i++) {
445 assert(gi < ma.
m && gj < ma.
n);
448 if (gi < 0 || gj < 0 || gi >= m || gj >= n)
451 ptrmyblockcopy[i].
r = -1;
453 if (ptrmyblockvide[i].r != ptrmyblockcopy[i].r) {
455 printf(
"Proc %d : Error element number %d, value = %f , initvalue =%f \n"
457 ptrmyblockvide[i].r, ptrmyblockcopy[i].r);
461 printf(
"Processor %d, has tested %d COMPLEX*16 elements,\
462Number of redistribution errors = %d \n",
463 mypnum, blocksize0, nberrors);
470 printf(
" => Total number of redistribution errors = %d \n",
473 printf(
"TEST PASSED OK\n");
475 if (myprow0 >= 0 && mypcol0 >= 0) {
481 if (myprow1 >= 0 && mypcol1 >= 0) {