pbdtr2a1()

subroutine pbdtr2a1	(	integer	icontxt,
		integer	n,
		integer	nb,
		integer	nz,
		double precision, dimension( * )	x,
		integer	incx,
		double precision	beta,
		double precision, dimension( * )	y,
		integer	incy,
		integer	intv
	)

Definition at line 646 of file pbdtrnv.f.

*
*  -- PB-BLAS routine (version 2.1) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory.
*     April 28, 1996
*
*     .. Scalar Arguments ..
      INTEGER              ICONTXT, N, NB, NZ, INCX, INCY, INTV
      DOUBLE PRECISION     BETA
*     ..
*     .. Array Arguments ..
      DOUBLE PRECISION     X( * ), Y( * )
*     ..
*
*  Purpose
*  =======
*
*     y <== x
*     y is a scattered vector, copied from a condensed vector x.
*
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC            min
*     ..
*     .. External Functions ..
      INTEGER              ICEIL
      EXTERNAL             iceil
*     ..
*     .. External Subroutines ..
      EXTERNAL             pbdvecadd
*     ..
*     .. Parameters ..
      DOUBLE PRECISION     ONE
      parameter( one = 1.0d+0 )
*     ..
*     .. Local Variables ..
      INTEGER              IX, IY, JZ, K, ITER
*
      ix = 0
      iy = 0
      jz = nz
      iter = iceil( n+nz, intv )
*
      IF( iter.GT.1 ) THEN
         CALL pbdvecadd( icontxt, 'G', nb-jz, one, x(ix*incx+1), incx,
     $                   beta, y(iy*incy+1), incy )
         ix = ix + nb   - jz
         iy = iy + intv - jz
         jz = 0
*
         DO 10 k = 2, iter-1
            CALL pbdvecadd( icontxt, 'G', nb, one, x(ix*incx+1), incx,
     $                      beta, y(iy*incy+1), incy )
            ix = ix + nb
            iy = iy + intv
   10    CONTINUE
      END IF
*
      CALL pbdvecadd( icontxt, 'G', min( n-iy, nb-jz ), one,
     $                x(ix*incx+1), incx, beta, y(iy*incy+1), incy )
*
      RETURN
*
*     End of PBDTR2A1
*

Here is the call graph for this function:

Here is the caller graph for this function: