pbstr2b1()

subroutine pbstr2b1	(	integer	icontxt,
		character*1	trans,
		integer	n,
		integer	nb,
		integer	nz,
		real, dimension( * )	x,
		integer	incx,
		real	beta,
		real, dimension( * )	y,
		integer	incy,
		integer	jinx,
		integer	jiny
	)

Definition at line 718 of file pbstrnv.f.

*
*  -- PB-BLAS routine (version 2.1) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory.
*     April 28, 1996
*
*     .. Scalar Arguments ..
      CHARACTER*1          TRANS
      INTEGER              ICONTXT, N, NB, NZ, INCX, INCY, JINX, JINY
      REAL                 BETA
*     ..
*     .. Array Arguments ..
      REAL                 X( * ), Y( * )
*     ..
*
*  Purpose
*  =======
*
*     y <== x + beta * y
*     y is a condensed vector, copied from a scattered vector x
*
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC            min
*     ..
*     .. External Functions ..
      INTEGER              ICEIL
      EXTERNAL             iceil
*     ..
*     .. External Subroutines ..
      EXTERNAL             pbsvecadd
*     ..
*     .. Parameters ..
      REAL                 ONE
      parameter( one = 1.0e+0 )
*     ..
*     .. Local Variables ..
      INTEGER              IX, IY, JZ, K, ITER, LENX, LENY
*
      IF( jinx.EQ.1 .AND. jiny.EQ.1 ) THEN
         CALL pbsvecadd( icontxt, trans, n, one, x, incx, beta,
     $                   y, incy )
*
      ELSE
         ix   = 0
         iy   = 0
         jz   = nz
         lenx = nb * jinx
         leny = nb * jiny
         iter = iceil( n+nz, lenx )
*
         IF( iter.GT.1 ) THEN
            CALL pbsvecadd( icontxt, trans, nb-jz, one, x(ix*incx+1),
     $                      incx, beta, y(iy*incy+1), incy )
            ix = ix + lenx - jz
            iy = iy + leny - jz
            jz = 0
*
            DO 10 k = 2, iter-1
               CALL pbsvecadd( icontxt, trans, nb, one, x(ix*incx+1),
     $                         incx, beta, y(iy*incy+1), incy )
               ix = ix + lenx
               iy = iy + leny
   10       CONTINUE
         END IF
*
         CALL pbsvecadd( icontxt, trans, min( n-ix, nb-jz ), one,
     $                   x(ix*incx+1), incx, beta, y(iy*incy+1), incy )
      END IF
*
      RETURN
*
*     End of PBSTR2B1
*

Here is the call graph for this function:

Here is the caller graph for this function: