SRC/zhetd2.f(3)
π
SRC/zhetd2.f(3) β Library Functions Manual π
okay so bestie, you found the page for SRC/zhetd2.f, and the manual name is Library Functions Manual. It's giving serious LAPACK energy, no cap, and we're gonna walk through the whole thing together. β¨
β¨ NAME β¨
SRC/zhetd2.f
πΈ SYNOPSIS πΈ
β¨ Functions/Subroutines β¨
subroutine zhetd2 (uplo, n, a, lda, d, e, tau, info)
ZHETD2 reduces a Hermitian matrix to real symmetric tridiagonal form
by an unitary similarity transformation (unblocked algorithm).
π Function/Subroutine Documentation π
subroutine zhetd2 (character uplo, integer n, complex*16, dimension( lda, *
) a, integer lda, double precision, dimension( * ) d, double precision,
dimension( * ) e, complex*16, dimension( * ) tau, integer info)
ZHETD2 reduces a Hermitian matrix to real symmetric tridiagonal form by an unitary similarity transformation (unblocked algorithm). π«
π Purpose π
Vibe check: ZHETD2 reduces a complex Hermitian matrix A to real symmetric tridiagonal form T by a unitary similarity transformation:
Q**H * A * Q = T.
π§ Parameters π§
UPLOβ isCHARACTER*1and specifies whether the upper or lower triangular part of the Hermitian matrixAis stored:= 'U': Upper triangular= 'L': Lower triangular
Nβ isINTEGER. The order of the matrixA.N >= 0.Aβ isCOMPLEX*16array, dimension(LDA,N). On entry, it's the Hermitian matrixA. IfUPLO = 'U', the leading n-by-n upper triangular part ofAcontains the upper triangular part of the matrixA, and the strictly lower triangular part ofAis not referenced. IfUPLO = 'L', the leading n-by-n lower triangular part ofAcontains the lower triangular part of the matrixA, and the strictly upper triangular part ofAis not referenced. On exit, ifUPLO = 'U', the diagonal and first superdiagonal ofAare overwritten by the corresponding elements of the tridiagonal matrixT, and the elements above the first superdiagonal, with the arrayTAU, represent the unitary matrixQas a product of elementary reflectors; ifUPLO = 'L', the diagonal and first subdiagonal ofAare overwritten by the corresponding elements of the tridiagonal matrixT, and the elements below the first subdiagonal, with the arrayTAU, represent the unitary matrixQas a product of elementary reflectors. See Further Details. πLDAβ isINTEGER. The leading dimension of the arrayA.LDA >= max(1,N).Dβ isDOUBLE PRECISIONarray, dimension(N). The diagonal elements of the tridiagonal matrixT:D(i) = A(i,i).Eβ isDOUBLE PRECISIONarray, dimension(N-1). The off-diagonal elements of the tridiagonal matrixT:E(i) = A(i,i+1)ifUPLO = 'U',E(i) = A(i+1,i)ifUPLO = 'L'.TAUβ isCOMPLEX*16array, dimension(N-1). The scalar factors of the elementary reflectors (see Further Details).INFOβ isINTEGER.= 0: successful exit< 0: ifINFO = -i, the i-th argument had an illegal value.
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π¦ Further Details π¦
If UPLO = 'U', the matrix Q is represented as a product of elementary reflectors:
Q = H(n-1) . . . H(2) H(1).
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(i+1:n) = 0 and v(i) = 1; v(1:i-1) is stored on exit in A(1:i-1,i+1), and tau in TAU(i).
If UPLO = 'L', the matrix Q is represented as a product of elementary reflectors:
Q = H(1) H(2) . . . H(n-1).
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) = 0 and v(i+1) = 1; v(i+2:n) is stored on exit in A(i+2:n,i), and tau in TAU(i).
The contents of A on exit are illustrated by the following examples with n = 5:
if UPLO = 'U': if UPLO = 'L':
( d e v2 v3 v4 ) ( d )
( d e v3 v4 ) ( e d )
( d e v4 ) ( v1 e d )
( d e ) ( v1 v2 e d )
( d ) ( v1 v2 v3 e d )
where d and e denote diagonal and off-diagonal elements of T, and vi denotes an element of the vector defining H(i).
Definition at line 174 of file zhetd2.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
That's the whole page, bestie. This footer is giving LAPACK, Version 3.12.0, and SRC/zhetd2.f(3) β there's no date in the footer, so that's the full cute sign-off. You ate, fr fr! ππβ¨
π¦ 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