101 SUBROUTINE clapll( N, X, INCX, Y, INCY, SSMIN )
109 INTEGER incx, incy, n
113 COMPLEX x( * ), y( * )
120 parameter( zero = 0.0e+0 )
122 parameter( cone = ( 1.0e+0, 0.0e+0 ) )
126 COMPLEX a11, a12, a22, c, tau
149 CALL
clarfg( n, x( 1 ), x( 1+incx ), incx, tau )
153 c = -conjg( tau )*
cdotc( n, x, incx, y, incy )
154 CALL
caxpy( n, c, x, incx, y, incy )
156 CALL
clarfg( n-1, y( 1+incy ), y( 1+2*incy ), incy, tau )
163 CALL
slas2( abs( a11 ), abs( a12 ), abs( a22 ), ssmin, ssmax )
subroutine caxpy(N, CA, CX, INCX, CY, INCY)
CAXPY
subroutine clapll(N, X, INCX, Y, INCY, SSMIN)
CLAPLL measures the linear dependence of two vectors.
complex function cdotc(N, CX, INCX, CY, INCY)
CDOTC
subroutine clarfg(N, ALPHA, X, INCX, TAU)
CLARFG generates an elementary reflector (Householder matrix).
subroutine slas2(F, G, H, SSMIN, SSMAX)
SLAS2 computes singular values of a 2-by-2 triangular matrix.