zlarfy()

subroutine zlarfy	(	character	uplo,
		integer	n,
		complex*16, dimension( * )	v,
		integer	incv,
		complex*16	tau,
		complex*16, dimension( ldc, * )	c,
		integer	ldc,
		complex*16, dimension( * )	work )

ZLARFY

Purpose:

!>
!> ZLARFY applies an elementary reflector, or Householder matrix, H,
!> to an n x n Hermitian matrix C, from both the left and the right.
!>
!> H is represented in the form
!>
!>    H = I - tau * v * v'
!>
!> where  tau  is a scalar and  v  is a vector.
!>
!> If  tau  is  zero, then  H  is taken to be the unit matrix.
!> 

Parameters

[in]	UPLO	!> UPLO is CHARACTER*1 !> Specifies whether the upper or lower triangular part of the !> Hermitian matrix C is stored. !> = 'U': Upper triangle !> = 'L': Lower triangle !>
[in]	N	!> N is INTEGER !> The number of rows and columns of the matrix C. N >= 0. !>
[in]	V	!> V is COMPLEX*16 array, dimension !> (1 + (N-1)*abs(INCV)) !> The vector v as described above. !>
[in]	INCV	!> INCV is INTEGER !> The increment between successive elements of v. INCV must !> not be zero. !>
[in]	TAU	!> TAU is COMPLEX*16 !> The value tau as described above. !>
[in,out]	C	!> C is COMPLEX*16 array, dimension (LDC, N) !> On entry, the matrix C. !> On exit, C is overwritten by H * C * H'. !>
[in]	LDC	!> LDC is INTEGER !> The leading dimension of the array C. LDC >= max( 1, N ). !>
[out]	WORK	!> WORK is COMPLEX*16 array, dimension (N) !>

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.

Definition at line 107 of file zlarfy.f.

*
*  -- LAPACK test routine --
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
*
*     .. Scalar Arguments ..
      CHARACTER          UPLO
      INTEGER            INCV, LDC, N
      COMPLEX*16         TAU
*     ..
*     .. Array Arguments ..
      COMPLEX*16         C( LDC, * ), V( * ), WORK( * )
*     ..
*
*  =====================================================================
*
*     .. Parameters ..
      COMPLEX*16         ONE, ZERO, HALF
      parameter( one = ( 1.0d+0, 0.0d+0 ),
     $                   zero = ( 0.0d+0, 0.0d+0 ),
     $                   half = ( 0.5d+0, 0.0d+0 ) )
*     ..
*     .. Local Scalars ..
      COMPLEX*16         ALPHA
*     ..
*     .. External Subroutines ..
      EXTERNAL           zaxpy, zhemv, zher2
*     ..
*     .. External Functions ..
      COMPLEX*16         ZDOTC
      EXTERNAL           zdotc
*     ..
*     .. Executable Statements ..
*
      IF( tau.EQ.zero )
     $   RETURN
*
*     Form  w:= C * v
*
      CALL zhemv( uplo, n, one, c, ldc, v, incv, zero, work, 1 )
*
      alpha = -half*tau*zdotc( n, work, 1, v, incv )
      CALL zaxpy( n, alpha, v, incv, work, 1 )
*
*     C := C - v * w' - w * v'
*
      CALL zher2( uplo, n, -tau, v, incv, work, 1, c, ldc )
*
      RETURN
*
*     End of ZLARFY
*

Here is the call graph for this function:

Here is the caller graph for this function: