/home/abui...N/zlarot.f(3)
π NAME π
ok so this is /home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/TESTING/MATGEN/zlarot.f β like, the whole entire file path is the name, bestie. it's from the MATGEN testing folder of LAPACK, and it's a Library Functions Manual page, no cap. π«Ά
β¨ SYNOPSIS β¨
Functions/Subroutines
subroutine ZLAROT (lrows, lleft, lright, nl, c, s, a, lda, xleft, xright)β that's the one and onlyZLAROTπ₯
π Function/Subroutine Documentation π
subroutine ZLAROT (logical lrows, logical lleft, logical lright, integer nl, complex*16 c, complex*16 s, complex*16, dimension( * ) a, integer lda, complex*16 xleft, complex*16 xright)
ok literally here's the whole signature, look at her, she's iconic π
subroutine ZLAROT (logical lrows, logical lleft, logical lright, integer nl,
complex*16 c, complex*16 s, complex*16, dimension( * ) a, integer lda,
complex*16 xleft, complex*16 xright)
π ZLAROT π
π¦ Purpose:
!>
!> ZLAROT applies a (Givens) rotation to two adjacent rows or
!> columns, where one element of the first and/or last column/row
!> for use on matrices stored in some format other than GE, so
!> that elements of the matrix may be used or modified for which
!> no array element is provided.
!>
!> One example is a symmetric matrix in SB format (bandwidth=4), for
!> which UPLO='L': Two adjacent rows will have the format:
!>
!> row j: C> C> C> C> C> . . . .
!> row j+1: C> C> C> C> C> . . . .
!>
!> '*' indicates elements for which storage is provided,
!> '.' indicates elements for which no storage is provided, but
!> are not necessarily zero; their values are determined by
!> symmetry. ' ' indicates elements which are necessarily zero,
!> and have no storage provided.
!>
!> Those columns which have two '*'s can be handled by DROT.
!> Those columns which have no '*'s can be ignored, since as long
!> as the Givens rotations are carefully applied to preserve
!> symmetry, their values are determined.
!> Those columns which have one '*' have to be handled separately,
!> by using separate variables and :
!>
!> row j: C> C> C> C> C> p . . .
!> row j+1: q C> C> C> C> C> . . . .
!>
!> The element p would have to be set correctly, then that column
!> is rotated, setting p to its new value. The next call to
!> ZLAROT would rotate columns j and j+1, using p, and restore
!> symmetry. The element q would start out being zero, and be
!> made non-zero by the rotation. Later, rotations would presumably
!> be chosen to zero q out.
!>
!> Typical Calling Sequences: rotating the i-th and (i+1)-st rows.
!> ------- ------- ---------
!>
!> General dense matrix:
!>
!> CALL ZLAROT(.TRUE.,.FALSE.,.FALSE., N, C,S,
!> A(i,1),LDA, DUMMY, DUMMY)
!>
!> General banded matrix in GB format:
!>
!> j = MAX(1, i-KL )
!> NL = MIN( N, i+KU+1 ) + 1-j
!> CALL ZLAROT( .TRUE., i-KL.GE.1, i+KU.LT.N, NL, C,S,
!> A(KU+i+1-j,j),LDA-1, XLEFT, XRIGHT )
!>
!> [ note that i+1-j is just MIN(i,KL+1) ]
!>
!> Symmetric banded matrix in SY format, bandwidth K,
!> lower triangle only:
!>
!> j = MAX(1, i-K )
!> NL = MIN( K+1, i ) + 1
!> CALL ZLAROT( .TRUE., i-K.GE.1, .TRUE., NL, C,S,
!> A(i,j), LDA, XLEFT, XRIGHT )
!>
!> Same, but upper triangle only:
!>
!> NL = MIN( K+1, N-i ) + 1
!> CALL ZLAROT( .TRUE., .TRUE., i+K.LT.N, NL, C,S,
!> A(i,i), LDA, XLEFT, XRIGHT )
!>
!> Symmetric banded matrix in SB format, bandwidth K,
!> lower triangle only:
!>
!> [ same as for SY, except:]
!> . . . .
!> A(i+1-j,j), LDA-1, XLEFT, XRIGHT )
!>
!> [ note that i+1-j is just MIN(i,K+1) ]
!>
!> Same, but upper triangle only:
!> . . .
!> A(K+1,i), LDA-1, XLEFT, XRIGHT )
!>
!> Rotating columns is just the transpose of rotating rows, except
!> for GB and SB: (rotating columns i and i+1)
!>
!> GB:
!> j = MAX(1, i-KU )
!> NL = MIN( N, i+KL+1 ) + 1-j
!> CALL ZLAROT( .TRUE., i-KU.GE.1, i+KL.LT.N, NL, C,S,
!> A(KU+j+1-i,i),LDA-1, XTOP, XBOTTM )
!>
!> [note that KU+j+1-i is just MAX(1,KU+2-i)]
!>
!> SB: (upper triangle)
!>
!> . . . . . .
!> A(K+j+1-i,i),LDA-1, XTOP, XBOTTM )
!>
!> SB: (lower triangle)
!>
!> . . . . . .
!> A(1,i),LDA-1, XTOP, XBOTTM )
!>
πΈ The Arguments (yes, all of them, we're doing this) πΈ
LROWSβ LOGICAL. if.TRUE., thenZLAROTrotates two rows; if.FALSE., it rotates two columns. not modified. π«LLEFTβ LOGICAL. if.TRUE., thenXLEFTgets used instead of the corresponding element ofAfor the first element in the second row (ifLROWS=.FALSE.) or column (ifLROWS=.TRUE.). if.FALSE., then the corresponding element ofAis what gets used. not modified. π«§LRIGHTβ LOGICAL. if.TRUE., thenXRIGHTgets used instead of the corresponding element ofAfor the last element in the first row (ifLROWS=.FALSE.) or column (ifLROWS=.TRUE.). if.FALSE., then the corresponding element ofAgets used. not modified. βNLβ INTEGER. the length of the rows (ifLROWS=.TRUE.) or columns (ifLROWS=.FALSE.) to be rotated. ifXLEFTand/orXRIGHTare used, the columns/rows they're in should be included inNLβ like, ifLLEFT = LRIGHT = .TRUE., thenNLmust be at least 2. the number of rows/columns to be rotated exclusive of those involvingXLEFTand/orXRIGHTmay not be negative, i.e.NLminus how many ofLLEFTandLRIGHTare.TRUE.must be at least zero; if it isn't,XERBLAwill be called π€. not modified.C,Sβ COMPLEX*16. these specify the Givens rotation to be applied. ifLROWSis true, then the matrix
( c s )
( _ _ )
(-s c )
gets applied from the left; if false, then the transpose (not conjugated) thereof gets applied from the right. bestie, note that in contrast to the output of ZROTG or to most versions of ZROT, both C and S are complex. for a Givens rotation, |C|**2 + |S|**2 should be 1, but this is not checked π. not modified.
Aβ COMPLEX*16 array. the array containing the rows/columns to be rotated. the first element ofAshould be the upper left element to be rotated. read and modified. πLDAβ INTEGER. the leading dimension ofA. ifAcontains a matrix stored in GE, HE, or SY format, then this is just the leading dimension ofAas dimensioned in the calling routine. ifAcontains a matrix stored in band (GB, HB, or SB) format, then this should be one less than the leading dimension used in the calling routine. so like, ifAwere dimensionedA(LDA,*)inZLAROT, thenA(1,j)would be the j-th element in the first of the two rows to be rotated, andA(2,j)would be the j-th in the second, regardless of how the array may be stored in the calling routine. (Acannot, however, actually be dimensioned thus, since for band format the row number may exceedLDA, which is not legal FORTRAN. π) ifLROWS=.TRUE., thenLDAmust be at least 1, otherwise it must be at leastNLminus the number of.TRUE.values inXLEFTandXRIGHT. not modified.XLEFTβ COMPLEX*16. ifLLEFTis.TRUE., thenXLEFTgets used and modified instead ofA(2,1)(ifLROWS=.TRUE.) orA(1,2)(ifLROWS=.FALSE.). read and modified. πXRIGHTβ COMPLEX*16. ifLRIGHTis.TRUE., thenXRIGHTgets used and modified instead ofA(1,NL)(ifLROWS=.TRUE.) orA(NL,1)(ifLROWS=.FALSE.). read and modified. π
π Author
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π₯ Definition
definition at line 227 of file zlarot.f. periodt.
π Author π
this page was generated automatically by Doxygen for LAPACK, straight from the source code. yeah, a robot wrote the original, but I wrote this version, so. π§
fr fr bestie, that's the whole page. every flag, every.TRUE./.FALSE., everyLDA-1β all of it. nothing skipped, no cap. π
that's a wrap on LAPACK, Version 3.12.0, <pink>/home/abui...N/zlarot.f(3)</pink> β mwah π
π¦ where this page comes from
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms