pbctr2at()

subroutine pbctr2at	(	integer	icontxt,
		character*1	adist,
		character*1	trans,
		integer	m,
		integer	n,
		integer	nb,
		complex, dimension( lda, * )	a,
		integer	lda,
		complex	beta,
		complex, dimension( ldb, * )	b,
		integer	ldb,
		integer	lcmp,
		integer	lcmq
	)

Definition at line 613 of file pbctran.f.

*
*  -- PB-BLAS routine (version 2.1) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory.
*     April 28, 1996
*
*     .. Scalar Arguments ..
      CHARACTER*1        ADIST, TRANS
      INTEGER            ICONTXT, LCMP, LCMQ, LDA, LDB, M, N, NB
      COMPLEX            BETA
*     ..
*     .. Array Arguments ..
      COMPLEX            A( LDA, * ), B( LDB, * )
*     ..
*
*  Purpose
*  =======
*
*  PBCTR2AT forms   B <== A^T + beta*B, or A^C + beta*B
*  B is a ((conjugate) transposed) scattered block row (or column),
*  copied from a scattered block column (or row) of A
*
*  =====================================================================
*
*     .. Parameters ..
      COMPLEX            ONE
      parameter( one  = ( 1.0e+0, 0.0e+0 ) )
*     ..
*     .. Local Scalars ..
      INTEGER            IA, IB, K, INTV, JNTV
*     ..
*     .. External Subroutines ..
      EXTERNAL           pbcmatadd
*     ..
*     .. External Functions ..
      LOGICAL            LSAME
      INTEGER            ICEIL
      EXTERNAL           lsame, iceil
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC          min
*     ..
*     .. Excutable Statements ..
*
      IF( lcmp.EQ.lcmq ) THEN
         CALL pbcmatadd( icontxt, trans, n, m, one, a, lda, beta, b,
     $                   ldb )
*
      ELSE
*
*        If A is a column block ( ADIST = 'C' ),
*
         IF( lsame( adist, 'C' ) ) THEN
            intv = lcmp * nb
            jntv = lcmq * nb
            ia = 1
            ib = 1
            DO 10 k = 1, iceil( m, intv )
               CALL pbcmatadd( icontxt, trans, n, min( m-ia+1, nb ),
     $                         one, a(ia,1), lda, beta, b(1,ib), ldb )
               ia = ia + intv
               ib = ib + jntv
   10       CONTINUE
*
*        If A is a row block ( ADIST = 'R' ),
*
         ELSE
            intv = lcmp * nb
            jntv = lcmq * nb
            ia = 1
            ib = 1
            DO 20 k = 1, iceil( n, jntv )
               CALL pbcmatadd( icontxt, trans, min( n-ia+1, nb ), m,
     $                         one, a(1,ia), lda, beta, b(ib,1), ldb )
               ia = ia + jntv
               ib = ib + intv
   20       CONTINUE
         END IF
      END IF
*
      RETURN
*
*     End of PBCTR2AT
*

Here is the call graph for this function:

Here is the caller graph for this function: