SRC/DEPRECATED/clahrd.f(3)
π NAME π
Ok bestie, this is the page for SRC/DEPRECATED/clahrd.f(3) from the Library Functions Manual. The name of the game is SRC/DEPRECATED/clahrd.f β yes, that's literally the whole NAME line, iconic. π
β¨ SYNOPSIS β¨
πΈ Functions/Subroutines πΈ
subroutine clahrd (n, k, nb, a, lda, tau, t, ldt, y, ldy)βCLAHRDreduces the firstnbcolumns of a general rectangular matrixAso that elements below thek-th subdiagonal are zero, and returns auxiliary matrices which are needed to apply the transformation to the unreduced part ofA.
π«Ά Function/Subroutine Documentation π«Ά
β subroutine clahrd β
subroutine clahrd (integer n, integer k, integer nb, complex, dimension( lda, * ) a, integer lda, complex, dimension( nb ) tau, complex, dimension( ldt, nb ) t, integer ldt, complex, dimension( ldy, nb ) y, integer ldy)
CLAHRD reduces the first nb columns of a general rectangular matrix A so that elements below the k-th subdiagonal are zero, and returns auxiliary matrices which are needed to apply the transformation to the unreduced part of A.
π Purpose: π
Deprecation tea: This routine is deprecated and has been replaced by routine CLAHR2. π
CLAHRD reduces the first NB columns of a complex general n-by-(n-k+1) matrix A so that elements below the k-th subdiagonal are zero. The reduction is performed by a unitary similarity transformation Q**H * A * Q. The routine returns the matrices V and T which determine Q as a block reflector I - V*T*V**H, and also the matrix Y = A * V * T.
π Parameters π
NβN is INTEGER. The order of the matrixA.KβK is INTEGER. The offset for the reduction. Elements below thek-th subdiagonal in the firstNBcolumns are reduced to zero.NBβNB is INTEGER. The number of columns to be reduced.AβA is COMPLEX array, dimension (LDA,N-K+1). On entry, then-by-(n-k+1)general matrixA. On exit, the elements on and above thek-th subdiagonal in the firstNBcolumns are overwritten with the corresponding elements of the reduced matrix; the elements below thek-th subdiagonal, with the arrayTAU, represent the matrixQas a product of elementary reflectors. The other columns ofAare unchanged. See Further Details.LDAβLDA is INTEGER. The leading dimension of the arrayA.LDA >= max(1,N).TAUβTAU is COMPLEX array, dimension (NB). The scalar factors of the elementary reflectors. See Further Details.TβT is COMPLEX array, dimension (LDT,NB). The upper triangular matrixT.LDTβLDT is INTEGER. The leading dimension of the arrayT.LDT >= NB.YβY is COMPLEX array, dimension (LDY,NB). Then-by-nbmatrixY.LDYβLDY is INTEGER. The leading dimension of the arrayY.LDY >= max(1,N).
π Author π
Shoutout to the authors, periodt:
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π Further Details: π
The matrix Q is represented as a product of nb elementary reflectors:
Q = H(1) H(2) . . . H(nb).
Each H(i) has the form:
H(i) = I - tau * v * v**H
where tau is a complex scalar, and v is a complex vector with v(1:i+k-1) = 0, v(i+k) = 1; v(i+k+1:n) is stored on exit in A(i+k+1:n,i), and tau in TAU(i).
The elements of the vectors v together form the (n-k+1)-by-nb matrix V which is needed, with T and Y, to apply the transformation to the unreduced part of the matrix, using an update of the form:
A := (I - V*T*V**H) * (A - Y*V**H).
The contents of A on exit are illustrated by the following example with n = 7, k = 3 and nb = 2:
( a h a a a )
( a h a a a )
( a h a a a )
( h h a a a )
( v1 h a a a )
( v1 v2 a a a )
( v1 v2 a a a )
where a denotes an element of the original matrix A, h denotes a modified element of the upper Hessenberg matrix H, and vi denotes an element of the vector defining H(i).
Definition at line 166 of file clahrd.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code. π€
okay bestie, that's the whole vibe. This page is from LAPACK Version 3.12.0 and the file is SRC/DEPRECATED/clahrd.f(3) β go slay your LAPACK stuff, mwah! π
β¨π
π¦ where this page comes from
- arch rolling: lapack-doc 3.12.1-2 Β· custom
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms