STGEXC(3)
π STGEXC(3) β Library Functions Manual π
okay so, bestie, this is the Library Functions Manual page for STGEXC(3). The header/footer path is /home/abui...RC/stgexc.f(3), and the full source path is chilling in NAME below. π
π NAME π
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/SRC/stgexc.f
β¨ SYNOPSIS β¨
π« Functions/Subroutines π«
subroutine STGEXC (wantq, wantz, n, a, lda, b, ldb, q, ldq, z, ldz,
ifst, ilst, work, lwork, info)
STGEXC
π¦ Function/Subroutine Documentation π¦
subroutine STGEXC (logical wantq, logical wantz, integer n, real, dimension(
lda, * ) a, integer lda, real, dimension( ldb, * ) b, integer ldb, real,
dimension( ldq, * ) q, integer ldq, real, dimension( ldz, * ) z, integer
ldz, integer ifst, integer ilst, real, dimension( * ) work, integer
lwork, integer info)
STGEXC
β¨ Purpose β¨
okay so hereβs the main-character job of this routine: STGEXC reorders the generalized real Schur decomposition of a real matrix pair (A,B) using an orthogonal equivalence transformation
(A, B) = Q * (A, B) * Z**T,
so that the diagonal block of (A, B) with row index IFST is moved to row ILST.
(A, B) must be in generalized real Schur canonical form (as returned by SGGES), i.e. A is block upper triangular with 1-by-1 and 2-by-2 diagonal blocks. B is upper triangular.
Optionally, the matrices Q and Z of generalized Schur vectors are updated.
Q(in) * A(in) * Z(in)**T = Q(out) * A(out) * Z(out)**T
Q(in) * B(in) * Z(in)**T = Q(out) * B(out) * Z(out)**T
π Parameters π
WANTQβWANTQisLOGICAL..TRUE.: update the left transformation matrixQ;.FALSE.: do not updateQ.WANTZβWANTZisLOGICAL..TRUE.: update the right transformation matrixZ;.FALSE.: do not updateZ.NβNisINTEGER. The order of the matricesAandB.N >= 0.AβAisREALarray, dimension(LDA,N). On entry, the matrixAin generalized real Schur canonical form. On exit, the updated matrixA, again in generalized real Schur canonical form.LDAβLDAisINTEGER. The leading dimension of the arrayA.LDA >= max(1,N).BβBisREALarray, dimension(LDB,N). On entry, the matrixBin generalized real Schur canonical form(A,B). On exit, the updated matrixB, again in generalized real Schur canonical form(A,B).LDBβLDBisINTEGER. The leading dimension of the arrayB.LDB >= max(1,N).QβQisREALarray, dimension(LDQ,N). On entry, ifWANTQ = .TRUE., the orthogonal matrixQ. On exit, the updated matrixQ. IfWANTQ = .FALSE.,Qis not referenced.LDQβLDQisINTEGER. The leading dimension of the arrayQ.LDQ >= 1. IfWANTQ = .TRUE.,LDQ >= N.ZβZisREALarray, dimension(LDZ,N). On entry, ifWANTZ = .TRUE., the orthogonal matrixZ. On exit, the updated matrixZ. IfWANTZ = .FALSE.,Zis not referenced.LDZβLDZisINTEGER. The leading dimension of the arrayZ.LDZ >= 1. IfWANTZ = .TRUE.,LDZ >= N.IFSTβIFSTisINTEGER.ILSTβILSTisINTEGER. Specify the reordering of the diagonal blocks of(A, B). The block with row indexIFSTis moved to rowILST, by a sequence of swapping between adjacent blocks. On exit, ifIFSTpointed on entry to the second row of a2-by-2block, it is changed to point to the first row;ILSTalways points to the first row of the block in its final position (which may differ from its input value by+1or-1).1 <= IFST, ILST <= N.WORKβWORKisREALarray, dimension(MAX(1,LWORK)). On exit, ifINFO = 0,WORK(1)returns the optimalLWORK.LWORKβLWORKisINTEGER. The dimension of the arrayWORK.LWORK >= 1whenN <= 1, otherwiseLWORK >= 4*N + 16.
If LWORK = -1, then a workspace query is assumed; the routine only calculates the optimal size of the WORK array, returns this value as the first entry of the WORK array, and no error message related to LWORK is issued by XERBLA.
INFOβINFOisINTEGER.=0: successful exit.<0: ifINFO = -i, thei-th argument had an illegal value.=1: The transformed matrix pair(A, B)would be too far from generalized Schur form; the problem is ill-conditioned.(A, B)may have been partially reordered, andILSTpoints to the first row of the current position of the block being moved.
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
πΈ Contributors: πΈ
Bo Kagstrom and Peter Poromaa, Department of Computing Science, Umea University, S-901 87 Umea, Sweden.
π References: π
[1] B. Kagstrom; A Direct Method for Reordering Eigenvalues in the Generalized Real Schur Form of a Regular Matrix Pair (A, B), in M.S. Moonen et al (eds), Linear Algebra for Large Scale and Real-Time Applications, Kluwer Academic Publ. 1993, pp 195-218.
Definition at line 218 of file stgexc.f.
πͺ Author πͺ
Generated automatically by Doxygen for LAPACK from the source code.
okay thatβs the whole manual page, bestie π«Ά β LAPACK Version 3.12.0, footer path /home/abui...RC/stgexc.f(3) included, Doxygen did its auto-generation thing, and yes, thatβs the cute little source/version sign-off. periodt. π
β¨
π¦ 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