DLAVSY_ROOK(3)
NAME
hey bestie! welcome to the TESTING/LI...avsy_rook.f(3) page from the Library Functions Manual π
this one is all about TESTING/LIN/dlavsy_rook.f, and the main character is the dlavsy_rook subroutine. no cap.
SYNOPSIS
Functions/Subroutines
subroutine dlavsy_rook (uplo, trans, diag, n, nrhs, a, lda, ipiv, b,
ldb, info)
DLAVSY_ROOK
Function/Subroutine Documentation
here's the full signature, literally exactly as the source has it π«Ά
subroutine dlavsy_rook (character uplo, character trans, character diag,
integer n, integer nrhs, double precision, dimension( lda, * ) a,
integer lda, integer, dimension( * ) ipiv, double precision, dimension(
ldb, * ) b, integer ldb, integer info)
DLAVSY_ROOK
Purpose
okay so DLAVSY_ROOK performs one of the matrix-vector operations
x := A*x or x := A'*x,
where x is an N element vector and A is one of the factors from the block U*D*U' or L*D*L' factorization computed by DSYTRF_ROOK.
- If
TRANS = 'N', multiplies byUorU * D(orLorL * D) - If
TRANS = 'T', multiplies byU'orD * U'(orL'orD * L') - If
TRANS = 'C', multiplies byU'orD * U'(orL'orD * L')
Parameters
UPLO
UPLO is CHARACTER*1. it specifies whether the factor stored in A is upper or lower triangular.
= 'U':Upper triangular= 'L':Lower triangular
TRANS
TRANS is CHARACTER*1. it specifies the operation to be performed:
= 'N':x := A*x= 'T':x := A'*x= 'C':x := A'*x
DIAG
DIAG is CHARACTER*1. it specifies whether or not the diagonal blocks are unit matrices. if the diagonal blocks are assumed to be unit, then A = U or A = L, otherwise A = U*D or A = L*D.
= 'U':Diagonal blocks are assumed to be unit matrices.= 'N':Diagonal blocks are assumed to be non-unit matrices.
N
N is INTEGER. the number of rows and columns of the matrix A. N >= 0.
NRHS
NRHS is INTEGER. the number of right hand sides, i.e., the number of vectors x to be multiplied by A. NRHS >= 0.
A
A is DOUBLE PRECISION array, dimension (LDA,N). the block diagonal matrix D and the multipliers used to obtain the factor U or L as computed by DSYTRF_ROOK. stored as a 2-D triangular matrix.
LDA
LDA is INTEGER. the leading dimension of the array A. LDA >= max(1,N).
IPIV
IPIV is INTEGER array, dimension (N). details of the interchanges and the block structure of D, as determined by DSYTRF_ROOK.
if UPLO = 'U':
- If
IPIV(k) > 0, then rows and columnskandIPIV(k)were interchanged andD(k,k)is a 1-by-1 diagonal block. (IfIPIV( k ) = k, no interchange was done). - If
IPIV(k) < 0andIPIV(k-1) < 0, then rows and columnskand-IPIV(k)were interchanged and rows and columnsk-1and-IPIV(k-1)were inerchaged,D(k-1:k,k-1:k)is a 2-by-2 diagonal block.
if UPLO = 'L':
- If
IPIV(k) > 0, then rows and columnskandIPIV(k)were interchanged andD(k,k)is a 1-by-1 diagonal block. (IfIPIV( k ) = k, no interchange was done). - If
IPIV(k) < 0andIPIV(k+1) < 0, then rows and columnskand-IPIV(k)were interchanged and rows and columnsk+1and-IPIV(k+1)were inerchaged,D(k:k+1,k:k+1)is a 2-by-2 diagonal block.
B
B is DOUBLE PRECISION array, dimension (LDB,NRHS). on entry, B contains NRHS vectors of length N. on exit, B is overwritten with the product A * B.
LDB
LDB is INTEGER. the leading dimension of the array B. LDB >= max(1,N).
INFO
INFO is INTEGER.
= 0:successful exit< 0:ifINFO = -k, thek-th argument had an illegal value
Author
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 155 of file dlavsy_rook.f.
Author
Generated automatically by Doxygen for LAPACK from the source code.
that's the whole page, bestie π the TESTING/LI...avsy_rook.f(3) page in the Library Functions Manual is brought to you by LAPACK Version 3.12.0 β¨ the footer has no date, so we're not inventing one β that would be a no-cap crime π«Άπ
π¦ 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