PB_Cplaprn2()

void PB_Cplaprn2	(	PBTYP_T *	TYPE,
		Int	M,
		Int	N,
		char *	A,
		Int	IA,
		Int	JA,
		Int *	DESCA,
		Int	IRPRNT,
		Int	ICPRNT,
		char *	CMATNM,
		Int	PROW,
		Int	PCOL
	)

Definition at line 246 of file PB_Cplaprnt.c.

{
/*
*  .. Local Scalars ..
*/
   char           type;
   Int            Acol, Aii, AisColRep, AisRowRep, Ajj, Ald, Arow, ctxt, h, i,
                  ib, icurcol, icurrow, ii, in, j, jb, jj, jn, ldw, mycol,
                  myrow, npcol, nprow, size, usiz;
/*
*  .. Local Arrays ..
*/
   char           * buf = NULL;
/* ..
*  .. Executable Statements ..
*
*/
/*
*  Retrieve process grid information
*/
   Cblacs_gridinfo( ( ctxt = DESCA[CTXT_] ), &nprow, &npcol, &myrow, &mycol );
/*
*  Retrieve sub( A )'s local information: Aii, Ajj, Arow, Acol ...
*/
   Ald  = DESCA[LLD_];
   PB_Cinfog2l( IA, JA, DESCA, nprow, npcol, myrow, mycol, &Aii, &Ajj, &Arow,
                &Acol );
/*
*  Save the local first index of each row and column sub( A )
*/
   ii   = Aii;
   jj   = Ajj;
/*
*  When sub( A ) is row-replicated, print the copy in process row PROW.
*  Otherwise, print the distributed matrix rows starting in process row Arow.
*/
   if( Arow < 0 ) { AisRowRep = 1; icurrow = Arow = PROW; }
   else           { AisRowRep = 0; icurrow = Arow;        }
/*
*  When sub( A ) is column-replicated, print the copy in process column PCOL.
*  Otherwise, print the distributed matrix columns starting in process column
*  Acol.
*/
   if( Acol < 0 ) { AisColRep = 1; icurcol = Acol = PCOL; }
   else           { AisColRep = 0; icurcol = Acol;        }
 
   type = TYPE->type; usiz = TYPE->usiz;  size = TYPE->size;
/*
*  Allocate buffer in printing process
*/
   ldw  = MAX( DESCA[ IMB_ ], DESCA[ MB_ ] );
   if( ( myrow == IRPRNT ) && ( mycol == ICPRNT ) )
      buf = PB_Cmalloc( ldw * size );
/*
*  Handle the first block of column separately
*/
   jb = PB_Cfirstnb( N, JA, DESCA[INB_], DESCA[NB_] );
   jn = JA + jb - 1;
 
   for( h = 0; h < jb; h++ )
   {
      ib = PB_Cfirstnb( M, IA, DESCA[IMB_], DESCA[MB_] );
      in = IA + ib - 1;
 
      if( ( icurrow == IRPRNT ) && ( icurcol == ICPRNT ) )
      {
         if( ( myrow == IRPRNT ) && ( mycol == ICPRNT ) )
         {
            PB_Cprnt( type, size, usiz, ib, Mptr( A, ii, jj+h, Ald, size ),
                      IA+1, JA+h+1, CMATNM );
         }
      }
      else
      {
         if( ( myrow == icurrow ) && ( mycol == icurcol ) )
         {
            TYPE->Cgesd2d( ctxt, ib, 1, Mptr( A, ii, jj+h, Ald, size ), Ald,
                           IRPRNT, ICPRNT );
         }
         else if( ( myrow == IRPRNT ) && ( mycol == ICPRNT ) )
         {
            TYPE->Cgerv2d( ctxt, ib, 1, buf, ldw, icurrow, icurcol );
            PB_Cprnt( type, size, usiz, ib, buf, IA+1, JA+h+1, CMATNM );
         }
      }
/*
*  Go to next block of rows
*/
      if( myrow == icurrow ) ii += ib;
      if( !( AisRowRep ) ) icurrow = MModAdd1( icurrow, nprow );
 
      Cblacs_barrier( ctxt, ALL );
/*
*  Loop over remaining block of rows
*/
      for( i = in+1; i <= IA+M-1; i += DESCA[MB_] )
      {
         ib = MIN( DESCA[MB_], IA+M-i );
         if( ( icurrow == IRPRNT ) && ( icurcol == ICPRNT ) )
         {
            if( ( myrow == IRPRNT ) && ( mycol == ICPRNT ) )
            {
               PB_Cprnt( type, size, usiz, ib, Mptr( A, ii, jj+h, Ald, size ),
                         i+1, JA+h+1, CMATNM );
            }
         }
         else
         {
            if( ( myrow == icurrow ) && ( mycol == icurcol ) )
            {
               TYPE->Cgesd2d( ctxt, ib, 1, Mptr( A, ii, jj+h, Ald, size ), Ald,
                              IRPRNT, ICPRNT );
            }
            else if( ( myrow == IRPRNT ) && ( mycol == ICPRNT ) )
            {
               TYPE->Cgerv2d( ctxt, ib, 1, buf, ldw, icurrow, icurcol );
               PB_Cprnt( type, size, usiz, ib, buf, i+1, JA+h+1, CMATNM);
            }
         }
/*
*  Go to next block of rows
*/
         if( myrow == icurrow ) ii += ib;
         if( !( AisRowRep ) ) icurrow = MModAdd1( icurrow, nprow );
 
         Cblacs_barrier( ctxt, ALL );
      }
/*
*  Restart at the first row to be printed
*/
      ii = Aii;
      icurrow = Arow;
   }
/*
*  Go to next block of columns
*/
   if( mycol == icurcol ) jj += jb;
   if( !( AisColRep ) ) icurcol = MModAdd1( icurcol, npcol );
 
   Cblacs_barrier( ctxt, ALL );
/*
*  Loop over remaining column blocks
*/
   for( j = jn+1; j <= JA+N-1; j += DESCA[NB_] )
   {
      jb = MIN( DESCA[NB_], JA+N-j );
      for( h = 0; h < jb; h++ )
      {
         ib = PB_Cfirstnb( M, IA, DESCA[IMB_], DESCA[MB_] );
         in = IA + ib - 1;
 
         if( ( icurrow == IRPRNT ) && ( icurcol == ICPRNT ) )
         {
            if( ( myrow == IRPRNT ) && ( mycol == ICPRNT ) )
            {
               PB_Cprnt( type, size, usiz, ib, Mptr( A, ii, jj+h, Ald, size ),
                         IA+1, j+h+1, CMATNM );
            }
         }
         else
         {
            if( ( myrow == icurrow ) && ( mycol == icurcol ) )
            {
               TYPE->Cgesd2d( ctxt, ib, 1, Mptr( A, ii, jj+h, Ald, size ), Ald,
                              IRPRNT, ICPRNT );
            }
            else if( ( myrow == IRPRNT ) && ( mycol == ICPRNT ) )
            {
               TYPE->Cgerv2d( ctxt, ib, 1, buf, ldw, icurrow, icurcol );
               PB_Cprnt( type, size, usiz, ib, buf, IA+1, j+h+1, CMATNM );
            }
         }
/*
*  Go to next block of rows
*/
         if( myrow == icurrow ) ii += ib;
         if( !( AisRowRep ) ) icurrow = MModAdd1( icurrow, nprow );
 
         Cblacs_barrier( ctxt, ALL );
/*
*  Loop over remaining block of rows
*/
         for( i = in+1; i <= IA+M-1; i += DESCA[MB_] )
         {
            ib = MIN( DESCA[MB_], IA+M-i );
            if( ( icurrow == IRPRNT ) && ( icurcol == ICPRNT ) )
            {
               if( ( myrow == IRPRNT ) && ( mycol == ICPRNT ) )
               {
                  PB_Cprnt( type, size, usiz, ib, Mptr( A, ii, jj+h, Ald,
                            size ), i+1, j+h+1, CMATNM );
               }
            }
            else
            {
               if( ( myrow == icurrow ) && ( mycol == icurcol ) )
               {
                  TYPE->Cgesd2d( ctxt, ib, 1, Mptr( A, ii, jj+h, Ald, size ),
                                 Ald, IRPRNT, ICPRNT );
               }
               else if( ( myrow == IRPRNT ) && ( mycol == ICPRNT ) )
               {
                  TYPE->Cgerv2d( ctxt, ib, 1, buf, ldw, icurrow, icurcol );
                  PB_Cprnt( type, size, usiz, ib, buf, i+1, j+h+1, CMATNM );
               }
            }
/*
*  Go to next block of rows
*/
            if( myrow == icurrow ) ii += ib;
            if( !( AisRowRep ) ) icurrow = MModAdd1( icurrow, nprow );
 
            Cblacs_barrier( ctxt, ALL );
         }
/*
*  Restart at the first row to be printed
*/
         ii      = Aii;
         icurrow = Arow;
      }
/*
*  Go to next block of columns
*/
      if( mycol == icurcol ) jj += jb;
      if( !( AisColRep ) ) icurcol = MModAdd1( icurcol, npcol );
 
      Cblacs_barrier( ctxt, ALL );
   }
 
   if( ( myrow == IRPRNT ) && ( mycol == ICPRNT ) && ( buf ) ) free( buf );
/*
*  End of PB_Cplaprn2
*/
}

Here is the call graph for this function: