Blahpack Translation Progress

Fortran BLAS/LAPACK → JavaScript — 1111 routines across 421 algorithms — generated 2026-04-13

769/1111
Routines Implemented
12230
Tests
74/79
BLAS
695/1043
LAPACK
421
Algorithms
Key: type prefixes & storage codes

Type prefixes

ddouble-precision real
zdouble-precision complex

Storage codes

bbblock bidiagonal
bdbidiagonal
didiagonal
gbgeneral band
gegeneral
gggeneral matrices (generalized problem)
gsgeneralized SVD
gtgeneral tridiagonal
hbHermitian (band)
heHermitian
hfHermitian (rectangular full packed)
hgupper Hessenberg (generalized)
hpHermitian (packed)
hsupper Hessenberg
laauxiliary
oporthogonal (packed)
ororthogonal
pbsymmetric/Hermitian positive definite (band)
pfsymmetric/Hermitian positive definite (rectangular full packed)
posymmetric/Hermitian positive definite
ppsymmetric/Hermitian positive definite (packed)
pspositive semi-definite
ptsymmetric/Hermitian positive definite (tridiagonal)
sbsymmetric (band)
sfsymmetric (rectangular full packed)
spsymmetric (packed)
stsymmetric tridiagonal
sysymmetric
tbtriangular (band)
tftriangular (rectangular full packed)
tgtriangular (generalized)
tptriangular (packed)
trtriangular
tztrapezoidal
ununitary
upunitary (packed)
Storage:
Search:

BLAS 25/28 algorithms (73/76 routines)

AlgorithmVariantsCoverageDescription
*amax
dz
100%Finds the index of the first element having maximum absolute value.
*asum
dz
100%Takes the sum of the absolute values.
*axpy
dz
100%Constant times a vector plus a vector.
dcabs1
z
100%Computes |Re(.)| + |Im(.)| of a double complex number
*copy
dz
100%Copies a vector, x, to a vector, y.
ddot
d
100%Forms the dot product of two vectors.
zdotc
z
100%Computes the dot product of two complex vectors, conjugating the first.
zdotu
z
100%Computes the dot product of two complex vectors, without conjugation.
xerbla
x
Is an error handler for the LAPACK routines.
xerbla_array
x
Assists other languages in calling XERBLA, the LAPACK and BLAS error handler.
*mm
gedz
hez
sydz
trdz
88%Performs one of the matrix-matrix operations \(C \coloneqq \alpha\operatorname{op}(A)\operatorname{op}(B) + \beta\,C\), where \(\operatorname{op}(X)\) is one of \(\operatorname{op}(X)\) = X or \(\operatorname{op}(X)\) = \(X^T\), \(\alpha\) and \(\beta\) are scalars, and A, B and C are matrices, with \(\operatorname{op}(A)\) an m by k matrix, \(\operatorname{op}(B)\) a k by n matrix and C an m by n matrix.
*mv
gbdz
gedz
hbz
hez
hpz
sbd
spd
syd
tbdz
tpdz
trdz
100%Performs one of the matrix-vector operations \(y \coloneqq \alpha\,Ax + \beta\,y\), or \(y \coloneqq \alpha\,A^Tx + \beta\,y\), where \(\alpha\) and \(\beta\) are scalars, x and y are vectors and A is an m by n band matrix, with kl sub-diagonals and ku super-diagonals.
*nrm2
dz
70%Returns the euclidean norm of a vector via the function name, so that DNRM2 := \(\sqrt{x'x }\)
*r
ged
hez
hpz
spd
syd
100%Performs the rank 1 operation \(A \coloneqq \alpha\,xy^T + A\), where \(\alpha\) is a scalar, x is an m element vector, y is an n element vector and A is an m by n matrix.
*r2
hez
hpz
spd
syd
100%Performs the symmetric rank 2 operation \(A \coloneqq \alpha\,xy^T + \alpha\,yx^T + A\), where \(\alpha\) is a scalar, x and y are n element vectors and A is an n by n symmetric matrix, supplied in packed form.
*r2k
hez
sydz
100%Performs one of the symmetric rank 2k operations \(C \coloneqq \alpha\,AB^T + \alpha\,BA^T + \beta\,C\), or \(C \coloneqq \alpha\,A^TB + \alpha\,B^TA + \beta\,C\), where \(\alpha\) and \(\beta\) are scalars, C is an n by n symmetric matrix and A and B are n by k matrices in the first case and k by n matrices in the second case.
zgerc
gez
100%Performs the rank 1 operation \(A \coloneqq \alpha\,xy^H + A\), where \(\alpha\) is a scalar, x is an m element vector, y is an n element vector and A is an m by n matrix.
*rk
hez
sydz
100%Performs one of the symmetric rank k operations \(C \coloneqq \alpha\,AA^T + \beta\,C\), or \(C \coloneqq \alpha\,A^TA + \beta\,C\), where \(\alpha\) and \(\beta\) are scalars, C is an n by n symmetric matrix and A is an n by k matrix in the first case and a k by n matrix in the second case.
*rot
dz
100%Applies a plane rotation.
*rotg
dz
79%Constructs a Givens plane rotation.
drotm
d
100%Applies a modified Givens rotation.
drotmg
d
92%Constructs a modified Givens rotation.
zgeru
gez
100%Performs the rank 1 operation \(A \coloneqq \alpha\,xy^T + A\), where \(\alpha\) is a scalar, x is an m element vector, y is an n element vector and A is an m by n matrix.
lsame
l
Returns .TRUE.
*scal
dz
100%Scales a vector by a constant.
*sm
trdz
97%Solves one of the matrix equations \(\operatorname{op}(A)\)*X = \(\alpha\)*B, or X*\(\operatorname{op}(A)\) = \(\alpha\)*B, where \(\alpha\) is a scalar, X and B are m by n matrices, A is a unit, or non-unit, upper or lower triangular matrix and \(\operatorname{op}(A)\) is one of \(\operatorname{op}(A)\) = A or \(\operatorname{op}(A)\) = \(A^T\).
*sv
tbdz
tpdz
trdz
100%Solves one of the systems of equations \(Ax = b\), or \(A^Tx = b\), where b and x are n element vectors and A is an n by n unit, or non-unit, upper or lower triangular band matrix, with ( k + 1 ) diagonals.
*swap
dz
100%Interchanges two vectors.

LAPACK 240/393 algorithms (693/1035 routines)

AlgorithmVariantsCoverageDescription
*2st_kernels
hbz
sbd
Is an internal routine used by the DSYTRD_SB2ST subroutine.
*_gbamv
ladz
100%Performs a matrix-vector operation to calculate error bounds.
dla_gbrcond
lad
100%Estimates the Skeel condition number for a general banded matrix.
zla_gbrcond_c
laz
100%Computes the infinity norm condition number of \(\operatorname{op}(A)\)*inv(diag(c)) for general banded matrices.
zla_gbrcond_x
laz
100%Computes the infinity norm condition number of \(\operatorname{op}(A)\)*diag(x) for general banded matrices.
*_gbrfsx_extended
ladz
95%Improves the computed solution to a general banded system using extra-precise iterative refinement.
*_gbrpvgrw
ladz
62%Computes the reciprocal pivot growth factor norm(A)/norm(U) for a general banded matrix.
*_geamv
ladz
99%Computes a matrix-vector product using a general matrix to calculate error bounds.
dla_gercond
lad
100%Estimates the Skeel condition number for a general matrix.
zla_gercond_c
laz
100%Computes the infinity norm condition number of \(\operatorname{op}(A)\)*inv(diag(c)) for general matrices.
zla_gercond_x
laz
100%Computes the infinity norm condition number of \(\operatorname{op}(A)\)*diag(x) for general matrices.
*_gerfsx_extended
ladz
97%Improves the computed solution to a general system using extra-precise iterative refinement.
*_gerpvgrw
ladz
57%Computes the reciprocal pivot growth factor norm(A)/norm(U).
zla_heamv
laz
100%Computes a matrix-vector product using a Hermitian indefinite matrix to calculate error bounds.
zla_hercond_c
laz
100%Computes the infinity norm condition number of \(\operatorname{op}(A)\)*inv(diag(c)) for Hermitian indefinite matrices.
zla_hercond_x
laz
100%Computes the infinity norm condition number of \(\operatorname{op}(A)\)*diag(x) for Hermitian indefinite matrices.
zla_herfsx_extended
laz
97%Improves the computed solution to a Hermitian indefinite system using extra-precise iterative refinement.
zla_herpvgrw
laz
20%Computes the reciprocal pivot growth factor norm(A)/norm(U).
*_lin_berr
ladz
100%Computes a component-wise relative backward error.
dla_porcond
lad
100%Estimates the Skeel condition number for a symmetric positive-definite matrix.
zla_porcond_c
laz
100%Computes the infinity norm condition number of \(\operatorname{op}(A)\)*inv(diag(c)) for Hermitian positive-definite matrices.
zla_porcond_x
laz
100%Computes the infinity norm condition number of \(\operatorname{op}(A)\)*diag(x) for Hermitian positive-definite matrices.
*_porfsx_extended
ladz
96%Improves the computed solution to a positive definite system using extra-precise iterative refinement.
*_porpvgrw
ladz
30%Computes the reciprocal pivot growth factor norm(A)/norm(U) for a symmetric or Hermitian positive-definite matrix.
*_syamv
ladz
99%Computes a matrix-vector product using a symmetric indefinite matrix to calculate error bounds.
dla_syrcond
lad
100%Estimates the Skeel condition number for a symmetric indefinite matrix.
zla_syrcond_c
laz
100%Computes the infinity norm condition number of \(\operatorname{op}(A)\)*inv(diag(c)) for symmetric indefinite matrices.
zla_syrcond_x
laz
100%Computes the infinity norm condition number of \(\operatorname{op}(A)\)*diag(x) for symmetric indefinite matrices.
*_syrfsx_extended
ladz
Improves the computed solution to a symmetric indefinite system using extra-precise iterative refinement.
*_syrpvgrw
ladz
18%Computes the reciprocal pivot growth factor norm(A)/norm(U) for a symmetric indefinite matrix.
*_wwaddw
ladz
100%Adds a vector into a doubled-single vector.
izmax1
z
100%Finds the index of the first vector element of maximum absolute value.
dlabad
lad
100%Is a no-op and kept for compatibility reasons.
*bak
gedz
ggdz
95%Forms the right or left eigenvectors of a real general matrix by backward transformation on the computed eigenvectors of the balanced matrix output by DGEBAL.
*bal
gedz
ggdz
97%Balances a general real matrix A.
*bd2
gedz
100%Reduces a general matrix to bidiagonal form using an unblocked algorithm.
*bdb
ord
unz
Simultaneously bidiagonalizes the blocks of a partitioned orthogonal/unitary matrix (CS decomposition driver).
*bdb1
ord
unz
Simultaneously bidiagonalizes the blocks of a tall and skinny matrix with orthonormal columns (case 1).
*bdb2
ord
unz
Simultaneously bidiagonalizes the blocks of a tall and skinny matrix with orthonormal columns (case 2).
*bdb3
ord
unz
Simultaneously bidiagonalizes the blocks of a tall and skinny matrix with orthonormal columns (case 3).
*bdb4
ord
unz
Simultaneously bidiagonalizes the blocks of a tall and skinny matrix with orthonormal columns (case 4).
*bdb5
ord
unz
Orthogonalizes the column vector X = [ X1 ] [ X2 ] with respect to the columns of Q = [ Q1 ] .
*bdb6
ord
unz
Orthogonalizes the column vector X = [ X1 ] [ X2 ] with respect to the columns of Q = [ Q1 ] .
*brd
gbdz
gedz
ladz
74%Reduces a real general m-by-n band matrix A to upper bidiagonal form B by an orthogonal transformation: \(Q^TAP = B\).
zcgesv
z
Computes the solution to a system of linear equations using mixed-precision iterative refinement.
zlacgv
laz
100%Conjugates a complex vector.
ilaclc
lai
Scans a matrix for its last non-zero column.
ilaclr
lai
Scans a matrix for its last non-zero row.
*cn2
ladz
60%Estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vector products.
*con
gbdz
gedz
gtdz
hez
hpz
ladz
pbdz
podz
ppdz
ptdz
spdz
sydz
tbdz
tpdz
trdz
63%Estimates the reciprocal of the condition number of a real general band matrix A, in either the 1-norm or the infinity-norm, using the LU factorization computed by DGBTRF.
*con_3
hez
sydz
Estimates the reciprocal of the condition number using the bounded Bunch-Kaufman factorization.
*con_rook
hez
sydz
Estimates the reciprocal of the condition number using rook pivoting.
*conv
sydz
95%Convert A given by TRF into L and D and vice-versa.
*convf
sydz
100%Converts between factorization storage formats for symmetric/Hermitian indefinite factorizations.
*convf_rook
sydz
98%Converts between factorization storage formats for rook-pivoted symmetric/Hermitian indefinite factorizations.
zlacp2
laz
100%Copies all or part of a real two-dimensional array to a complex array.
zcposv
z
Computes the solution to a positive definite system using mixed-precision iterative refinement.
*cpy
ladz
100%Copies all or part of one two-dimensional array to another.
zlacrm
laz
Multiplies a complex matrix by a square real matrix.
zlacrt
laz
100%Performs a linear transformation of a pair of complex vectors.
*csd
bbdz
ord
unz
Computes the CS decomposition of a partitioned orthogonal/unitary matrix.
*csd2by1
ord
unz
Computes the CS decomposition of an M-by-Q matrix with orthonormal columns partitioned into a 2-by-1 block structure.
iladiag
lai
This subroutine translated from a character string specifying if a matrix has unit diagonal or not to the relevant BLAST-specified integer constant.
*div
ladz
100%Performs complex division in real arithmetic, avoiding unnecessary overflow.
iladlc
lai
100%Scans a matrix for its last non-zero column.
iladlr
lai
100%Scans a matrix for its last non-zero row.
*dmd
gedz
Computes the Dynamic Mode Decomposition (DMD) for a pair of data snapshot matrices.
*dmdq
gedz
Computes the Dynamic Mode Decomposition (DMD) for a pair of data snapshot matrices.
dlae2
lad
100%Computes the eigenvalues of a 2-by-2 symmetric matrix.
*ebz
lad
std
56%Computes the number of eigenvalues of a real symmetric tridiagonal matrix which are less than or equal to a given value, and performs other tasks required by the routine sstebz.
*ed0
ladz
Used by DSTEDC. Computes all eigenvalues and corresponding eigenvectors of an unreduced symmetric tridiagonal matrix using the divide and conquer method.
dlaed1
lad
Used by DSTEDC. Computes the updated eigensystem of a diagonal matrix after modification by a rank-one symmetric matrix. Used when the original matrix is tridiagonal.
dlaed2
lad
95%Used by DSTEDC. Merges eigenvalues and deflates secular equation. Used when the original matrix is tridiagonal.
dlaed3
lad
Used by DSTEDC. Finds the roots of the secular equation and updates the eigenvectors. Used when the original matrix is tridiagonal.
dlaed4
lad
Used by DSTEDC. Finds a single root of the secular equation.
dlaed5
lad
100%Used by DSTEDC. Solves the 2-by-2 secular equation.
dlaed6
lad
94%Used by DSTEDC. Computes one Newton step in solution of the secular equation.
*ed7
ladz
Used by DSTEDC. Computes the updated eigensystem of a diagonal matrix after modification by a rank-one symmetric matrix. Used when the original matrix is dense.
*ed8
ladz
Used by DSTEDC. Merges eigenvalues and deflates secular equation. Used when the original matrix is dense.
dlaed9
lad
Used by DSTEDC. Finds the roots of the secular equation and updates the eigenvectors. Used when the original matrix is dense.
dlaeda
lad
Used by DSTEDC. Computes the Z vector determining the rank-one modification of the diagonal matrix. Used when the original matrix is dense.
*edc
stdz
Computes all eigenvalues and, optionally, eigenvectors of a symmetric tridiagonal matrix using the divide and conquer method.
ieeeck
i
Is called from the ILAENV to verify that Infinity and possibly NaN arithmetic is safe (i.e. will not trap).
*egr
stdz
Computes selected eigenvalues and, optionally, eigenvectors of a real symmetric tridiagonal matrix T.
*ein
hsdz
ladz
stdz
52%Uses inverse iteration to find specified right and/or left eigenvectors of a real upper Hessenberg matrix H.
*emr
stdz
Computes selected eigenvalues and, optionally, eigenvectors of a real symmetric tridiagonal matrix T.
ilaenv
lai
Is called from the LAPACK routines to choose problem-dependent parameters for the local environment.
ilaenv2stage
lai
Is called from the LAPACK routines to choose problem-dependent parameters for the local environment.
*eqr
hsdz
ptdz
stdz
93%Computes the eigenvalues of a Hessenberg matrix H and, optionally, the matrices T and Z from the Schur decomposition H = Z T \(Z^T\), where T is an upper quasi-triangular matrix (the Schur form), and Z is the orthogonal matrix of Schur vectors.
*equ
gbdz
gedz
pbdz
podz
ppdz
99%Computes row and column scalings intended to equilibrate an M-by-N band matrix A and reduce its condition number.
*equb
gbdz
gedz
hez
podz
sydz
Computes row and column scalings intended to equilibrate an M-by-N matrix A and reduce its condition number.
*eqz
hgdz
86%Computes the eigenvalues of a real matrix pair (H,T), where H is an upper Hessenberg matrix and T is upper triangular, using the double-shift QZ method.
xerbla
x
Is an error handler for the LAPACK routines.
xerbla_array
x
Assists other languages in calling XERBLA, the LAPACK and BLAS error handler.
dsterf
std
93%Computes all eigenvalues of a symmetric tridiagonal matrix using the Pal-Walker-Kahan variant of the QL or QR algorithm.
*es
gedz
ggdz
73%Computes the eigenvalues, the Schur form, and, optionally, the matrix of Schur vectors for GE matrices
*es3
ggdz
Computes the eigenvalues, the Schur form, and, optionally, the matrix of Schur vectors for GE matrices (blocked algorithm)
*esx
gedz
ggdz
Computes the eigenvalues, the Schur form, and, optionally, the matrix of Schur vectors for GE matrices
zlaesy
laz
100%Computes the eigenvalues and eigenvectors of a 2-by-2 complex symmetric matrix.
*ev
gedz
ggdz
hbz
hez
hpz
sbd
spd
std
syd
95%Computes the eigenvalues and, optionally, the left and/or right eigenvectors for GE matrices
*ev2
ladz
100%Computes the eigenvalues and eigenvectors of a 2-by-2 symmetric/Hermitian matrix.
*ev3
ggdz
Computes the eigenvalues and, optionally, the left and/or right eigenvectors for GE matrices (blocked algorithm)
*ev_2stage
hbz
hez
sbd
syd
Computes eigenvalues and optionally eigenvectors (2-stage algorithm).
*evc
tgdz
trdz
92%Computes some or all of the right and/or left eigenvectors of a pair of real matrices (S,P), where S is a quasi-triangular matrix and P is upper triangular.
*evc3
trdz
92%Computes some or all of the right and/or left eigenvectors of a real upper quasi-triangular matrix T.
*evd
hbz
hez
hpz
sbd
spd
std
syd
Computes eigenvalues and optionally eigenvectors using divide and conquer.
*evd_2stage
hbz
hez
sbd
syd
Computes eigenvalues and optionally eigenvectors using divide and conquer (2-stage algorithm).
*evr
hez
std
syd
94%Computes eigenvalues and optionally eigenvectors using the MRRR algorithm.
*evr_2stage
hez
syd
Computes eigenvalues and optionally eigenvectors using the MRRR algorithm (2-stage).
*evx
gedz
ggdz
hbz
hez
hpz
sbd
spd
std
syd
85%Computes the eigenvalues and, optionally, the left and/or right eigenvectors for GE matrices
*evx_2stage
hbz
hez
sbd
syd
Computes selected eigenvalues and optionally eigenvectors (2-stage algorithm).
*ex2
tgdz
97%Swaps adjacent diagonal blocks in an upper (quasi) triangular matrix pair by an orthogonal equivalence transformation.
*exc
lad
tgdz
trdz
85%Swaps adjacent diagonal blocks of a real upper quasi-triangular matrix in Schur canonical form, by an orthogonal similarity transformation.
dlag2
lad
99%Computes the eigenvalues of a 2-by-2 generalized eigenvalue problem, with scaling as necessary to avoid over-/underflow.
zlag2c
laz
100%Converts a complex double precision matrix to a complex single precision matrix.
*g2l
ord
unz
100%Generates all or part of the orthogonal matrix Q from a QL factorization determined by sgeqlf (unblocked algorithm).
*g2r
ord
unz
100%Generates all or part of the orthogonal matrix Q from a QR factorization determined by sgeqrf (unblocked algorithm).
dlag2s
lad
100%Converts a double precision matrix to a single precision matrix.
*gbr
ord
unz
100%Generates one of the real orthogonal matrices Q or \(P^T\) determined by DGEBRD when reducing a real matrix A to bidiagonal form: A = Q * B * \(P^T\).
*ghr
ord
unz
100%Generates a real orthogonal matrix Q which is defined as the product of IHI-ILO elementary reflectors of order N, as returned by DGEHRD: Q = H(ilo) H(ilo+1) .
*gl2
ord
unz
100%Generates an m by n real matrix Q with orthonormal rows, which is defined as the first m rows of a product of k elementary reflectors of order n Q = H(k) .
*glm
ggdz
Solves a general Gauss-Markov linear model (GLM) problem.
*glq
ord
unz
100%Generates an M-by-N real matrix Q with orthonormal rows, which is defined as the first M rows of a product of K elementary reflectors of order N Q = H(k) .
*gql
ord
unz
81%Generates an M-by-N real matrix Q with orthonormal columns, which is defined as the last N columns of a product of K elementary reflectors of order M Q = H(k) .
*gqr
ord
unz
100%Generates an M-by-N real matrix Q with orthonormal columns, which is defined as the first N columns of a product of K elementary reflectors of order M Q = H(1) H(2) .
*gr2
ord
unz
100%Generates all or part of the orthogonal matrix Q from an RQ factorization determined by sgerqf (unblocked algorithm).
*grq
ord
unz
Generates an M-by-N real matrix Q with orthonormal rows, which is defined as the last M rows of a product of K elementary reflectors of order N Q = H(1) H(2) .
*gs2
hez
ladz
syd
99%Computes 2-by-2 orthogonal matrices U, V, and Q, and applies them to matrices A and B such that the rows of the transformed A and B are parallel.
*gst
hbz
hez
hpz
sbd
spd
syd
92%Reduces a real symmetric-definite banded generalized eigenproblem A*x = \(\lambda\)*B*x to standard form C*y = \(\lambda\)*y, such that C has the same bandwidth as A.
dlagtf
lad
93%Computes an LU factorization of a matrix T-λI, where T is a general tridiagonal matrix, and λ a scalar, using partial pivoting with row interchanges.
*gtm
ladz
99%Performs a matrix-matrix product of the form C = αAB+βC, where A is a tridiagonal matrix, B and C are rectangular matrices, and α and β are scalars, which may be 0, 1, or -1.
*gtr
opd
ord
unz
upz
100%Generates a real orthogonal matrix Q which is defined as the product of n-1 elementary reflectors H(i) of order n, as returned by DSPTRD using packed storage: if UPLO = 'U', Q = H(n-1) .
dlagts
lad
77%Solves the system of equations (T-λI)x = y
*gtsqr
ord
unz
Generates an M-by-N real matrix Q_out with orthonormal columns, which are the first N columns of a product of real orthogonal matrices of order M which are returned by DLATSQR Q_out = first_N_columns_of( Q(1)_in * Q(2)_in * ... * Q(k)_in ).
*gtsqr_row
ord
unz
Generates an M-by-N real matrix Q_out with orthonormal columns from the output of DLATSQR.
*gv
hbz
hez
hpz
sbd
spd
syd
98%Computes all the eigenvalues, and optionally, the eigenvectors of a real generalized symmetric-definite banded eigenproblem, of the form \(Ax=\lambda\,Bx\).
dlagv2
lad
100%Computes the Generalized Schur factorization of a real 2-by-2 matrix pencil (A,B) where B is upper triangular.
*gv_2stage
hez
syd
Computes all the eigenvalues, and optionally, the eigenvectors of a real generalized symmetric-definite eigenproblem, of the form \(Ax=\lambda\,Bx\), A*Bx=(\(\lambda\))*x, or B*A*x=(\(\lambda\))*x.
*gvd
hbz
hez
hpz
sbd
spd
syd
Computes all the eigenvalues, and optionally, the eigenvectors of a real generalized symmetric-definite banded eigenproblem, of the form \(Ax=\lambda\,Bx\).
*gvx
hbz
hez
hpz
sbd
spd
syd
96%Computes selected eigenvalues, and optionally, eigenvectors of a real generalized symmetric-definite banded eigenproblem, of the form \(Ax=\lambda\,Bx\).
*hd2
gedz
100%Reduces a general square matrix to upper Hessenberg form using an unblocked algorithm.
*hd3
ggdz
Reduces a pair of real matrices (A,B) to generalized upper Hessenberg form using orthogonal transformations, where A is a general matrix and B is upper triangular.
zlahef
laz
94%Computes a partial factorization of a complex Hermitian indefinite matrix using the Bunch-Kaufman diagonal pivoting method (blocked algorithm, calling Level 3 BLAS).
zlahef_aa
laz
DLAHEF_AA factorizes a panel of a complex hermitian matrix A using the Aasen's algorithm.
zlahef_rk
laz
Computes a partial factorization of a complex Hermitian indefinite matrix using bounded Bunch-Kaufman (rook) diagonal pivoting method.
zlahef_rook
laz
Computes a partial factorization of a complex Hermitian indefinite matrix using the bounded Bunch-Kaufman ("rook") diagonal pivoting method (blocked algorithm, calling Level 3 BLAS).
*hqr
ladz
95%Computes the eigenvalues and Schur factorization of an upper Hessenberg matrix, using the double-shift/single-shift QR algorithm.
*hr2
ladz
100%Reduces the specified number of first columns of a general rectangular matrix A so that elements below the specified subdiagonal are zero, and returns auxiliary matrices which are needed to apply the transformation to the unreduced part of A.
*hr_col
ord
unz
Takes an M-by-N real matrix Q_in with orthonormal columns as input, stored in A, and performs Householder Reconstruction (HR), i.e.
*hrd
gedz
ggdz
100%Reduces a real general matrix A to upper Hessenberg form H by an orthogonal similarity transformation: \(Q^TAQ = H\) .
*ic1
ladz
100%Applies one step of incremental condition estimation.
*isnan
lad
d
100%Tests input for NaN.
*jsv
gedz
Computes the singular value decomposition (SVD) of a real M-by-N matrix [A], where M >= N.
dlaln2
lad
84%Solves a 1-by-1 or 2-by-2 linear system of equations of the specified form.
*lq
gedz
Computes an LQ factorization of a real M-by-N matrix A: A = ( L 0 ) * Q where: Q is a N-by-N orthogonal matrix; L is a lower-triangular M-by-M matrix; 0 is a M-by-(N-M) zero matrix, if M < N.
*lq2
gedz
100%Computes the LQ factorization of a general rectangular matrix using an unblocked algorithm.
*lqf
gedz
100%Computes an LQ factorization of a real M-by-N matrix A: A = ( L 0 ) * Q where: Q is a N-by-N orthogonal matrix; L is a lower-triangular M-by-M matrix; 0 is a M-by-(N-M) zero matrix, if M < N.
*lqt
gedz
tpdz
Computes a blocked LQ factorization of a real M-by-N matrix A using the compact WY representation of Q.
*lqt2
tpdz
Computes a LQ factorization of a real or complex "triangular-pentagonal" matrix, which is composed of a triangular block and a pentagonal block, using the compact WY representation for Q.
*lqt3
gedz
Recursively computes a LQ factorization of a general real or complex matrix using the compact WY representation of Q.
*ls
gedz
98%Solves overdetermined or underdetermined systems for GE matrices
*ls0
ladz
Applies back multiplying factors in solving the least squares problem using divide and conquer SVD approach. Used by sgelsd.
*lsa
ladz
Computes the SVD of the coefficient matrix in compact form. Used by sgelsd.
*lsd
gedz
ladz
Computes the minimum-norm solution to a linear least squares problem for GE matrices
*lse
ggdz
Solves overdetermined or underdetermined systems for OTHER matrices
*lss
gedz
58%Solves overdetermined or underdetermined systems for GE matrices
*lst
gedz
Solves overdetermined or underdetermined systems for GE matrices using QR or LQ factorization with compact WY representation of Q.
*lsy
gedz
Solves overdetermined or underdetermined systems for GE matrices
*m22
ord
unz
100%Multiplies a general matrix by a banded orthogonal matrix.
*m2l
ord
unz
100%Multiplies a general matrix by the orthogonal matrix from a QL factorization determined by sgeqlf (unblocked algorithm).
*m2r
ord
unz
100%Multiplies a general matrix by the orthogonal matrix from a QR factorization determined by sgeqrf (unblocked algorithm).
*mbr
ord
unz
100%Multiplies a general matrix by the orthogonal/unitary matrix Q or P from a bidiagonal reduction.
*mhr
ord
unz
100%Multiplies a general matrix by the orthogonal/unitary matrix Q from a Hessenberg reduction.
*ml2
ord
unz
100%Multiplies a general matrix by the orthogonal matrix from a LQ factorization determined by sgelqf (unblocked algorithm).
*mlq
gedz
ord
unz
99%Multiplies a general matrix by the orthogonal/unitary matrix Q from an LQ factorization.
*mlqt
gedz
tpdz
Multiplies a general matrix by the orthogonal/unitary matrix Q from a blocked LQ factorization.
*mql
ord
unz
99%Multiplies a general matrix by the orthogonal/unitary matrix Q from a QL factorization.
*mqr
gedz
ord
unz
99%Multiplies a general matrix by the orthogonal/unitary matrix Q from a QR factorization.
*mqrt
gedz
tpdz
Multiplies a general matrix by the orthogonal/unitary matrix Q from a blocked QR factorization.
*mr2
ord
unz
100%Multiplies a general matrix by the orthogonal matrix from a RQ factorization determined by sgerqf (unblocked algorithm).
*mr3
ord
unz
Multiplies a general matrix by the orthogonal matrix from a RZ factorization determined by stzrzf (unblocked algorithm).
dlamrg
lad
100%Creates a permutation list to merge the entries of two independently sorted sets into a single set sorted in ascending order.
*mrq
ord
unz
98%Multiplies a general matrix by the orthogonal/unitary matrix Q from an RQ factorization.
*mrz
ord
unz
Multiplies a general matrix by the orthogonal/unitary matrix Q from an RZ factorization.
*mswlq
ladz
Multiplies a general matrix by the orthogonal/unitary matrix Q from a short-wide LQ factorization.
*mtr
opd
ord
unz
upz
100%Multiplies a general matrix by the orthogonal/unitary matrix Q stored in a triangular factorization.
*mtsqr
ladz
Multiplies a general matrix by the orthogonal/unitary matrix Q from a tall-skinny QR factorization.
*mv
spz
syz
100%Computes a matrix-vector product for complex vectors using a complex symmetric packed matrix
*n
gbdz
gedz
gtdz
hbz
hez
hfz
hpz
hsdz
sbdz
sfd
spdz
std
sydz
tbdz
tpdz
trdz
32%Returns the value of the 1-norm, Frobenius norm, infinity-norm, or the largest absolute value of any element of general band matrix.
dlaneg
lad
100%Computes the Sturm count.
zlanht
laz
37%Returns the value of the 1-norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex Hermitian tridiagonal matrix.
dlanv2
lad
93%Computes the Schur factorization of a real 2-by-2 nonsymmetric matrix in standard form.
dlaorhr_col_getrfnp
lad
Computes the modified LU factorization without pivoting of a real general M-by-N matrix A.
dlaorhr_col_getrfnp2
lad
Computes the modified LU factorization without pivoting of a real general M-by-N matrix A.
dsposv
spd
Computes the solution to system of linear equations \(AX = B\) for PO matrices
iparmq
i
This program sets problem and machine dependent parameters useful for xHSEQR and related subroutines for eigenvalue problems.
*pll
ladz
100%Measures the linear dependence of two vectors.
*pmr
ladz
100%Rearranges rows of a matrix as specified by a permutation vector.
*pmt
ladz
100%Performs a forward or backward permutation of the columns of a matrix.
ilaprec
lai
This subroutine translated from a character string specifying an intermediate precision to the relevant BLAST-specified integer constant.
dlapy2
lad
100%Returns \(\sqrt{x^2+y^2}\).
dlapy3
lad
100%Returns \(\sqrt{x^2+y^2+z^2}\).
*qgb
ladz
100%Scales a general band matrix, using row and column scaling factors computed by sgbequ.
*qge
ladz
100%Scales a general rectangular matrix, using row and column scaling factors computed by sgeequ.
zlaqhb
laz
100%Scales a Hermitian band matrix, using scaling factors computed by cpbequ.
zlaqhe
laz
100%Scales a Hermitian matrix.
zlaqhp
laz
100%Scales a Hermitian matrix stored in packed form.
*ql2
gedz
100%Computes the QL factorization of a general rectangular matrix using an unblocked algorithm.
*qlf
gedz
100%Computes a QL factorization of a real M-by-N matrix A: A = Q * L.
*qp2
ladz
99%Computes a QR factorization with column pivoting of the matrix block.
*qp2rk
ladz
Computes truncated QR factorization with column pivoting of a real matrix block using Level 2 BLAS and overwrites a real m-by-nrhs matrix B with \(Q^T\) * B.
*qp3
gedz
100%Computes a QR factorization with column pivoting of a matrix A: \(AP = Q\)*R using Level 3 BLAS.
*qp3rk
gedz
ladz
Computes a truncated Householder QR factorization with column pivoting of a real m-by-n matrix A by using Level 3 BLAS and overwrites a real m-by-nrhs matrix B with \(Q^T\) * B.
*qps
ladz
98%Computes a step of QR factorization with column pivoting of a real m-by-n matrix A by using BLAS level 3.
*qr
gedz
Computes a QR factorization of a real M-by-N matrix A: A = Q * ( R ), ( 0 ) where: Q is a M-by-M orthogonal matrix; R is an upper-triangular N-by-N matrix; 0 is a (M-N)-by-N zero matrix, if M > N.
*qr0
ladz
81%Computes the eigenvalues of a Hessenberg matrix, and optionally the matrices from the Schur decomposition.
*qr1
ladz
100%Sets a scalar multiple of the first column of the product of 2-by-2 or 3-by-3 matrix H and specified shifts.
*qr2
gedz
ladz
97%Computes the QR factorization of a general rectangular matrix using an unblocked algorithm.
*qr2p
gedz
100%Computes the QR factorization of a general rectangular matrix with non-negative diagonal elements using an unblocked algorithm.
*qr3
ladz
98%Performs the orthogonal similarity transformation of a Hessenberg matrix to detect and deflate fully converged eigenvalues from a trailing principal submatrix (aggressive early deflation).
*qr4
ladz
80%Computes the eigenvalues of a Hessenberg matrix, and optionally the matrices from the Schur decomposition.
*qr5
ladz
81%Performs a single small-bulge multi-shift QR sweep.
*qrf
gedz
ggdz
100%Computes a QR factorization of a real M-by-N matrix A: A = Q * ( R ), ( 0 ) where: Q is a M-by-M orthogonal matrix; R is an upper-triangular N-by-N matrix; 0 is a (M-N)-by-N zero matrix, if M > N.
*qrfp
gedz
98%Computes a QR factorization with non-negative diagonal entries.
*qrt
gedz
tpdz
Computes a blocked QR factorization of a real M-by-N matrix A using the compact WY representation of Q.
*qrt2
gedz
tpdz
Computes a QR factorization of a general real or complex matrix using the compact WY representation of Q.
*qrt3
gedz
Recursively computes a QR factorization of a general real or complex matrix using the compact WY representation of Q.
*qsb
ladz
100%Scales a symmetric/Hermitian band matrix, using scaling factors computed by spbequ.
*qsp
ladz
100%Scales a symmetric/Hermitian matrix in packed storage, using scaling factors computed by sppequ.
*qsy
ladz
100%Scales a symmetric/Hermitian matrix, using scaling factors computed by spoequ.
dlaqtr
lad
17%Solves a real quasi-triangular system of equations, or a complex quasi-triangular system of special form, in real arithmetic.
*qz0
ladz
Computes the eigenvalues of a real matrix pair (H,T), where H is an upper Hessenberg matrix and T is upper triangular, using the double-shift QZ method.
*qz1
ladz
Given a 3-by-3 matrix pencil (A,B), DLAQZ1 sets v to a scalar multiple of the first column of the product (*) K = (A - (beta2*sr2 - i*si)*B)*B^(-1)*(beta1*A - (sr2 + i*si2)*B)*B^(-1).
*qz2
ladz
Chases a 2x2 shift bulge in a matrix pencil down a single position
*qz3
ladz
Performs aggressive early deflation in the QZ algorithm.
dlaqz4
lad
Executes a single multishift QZ sweep
*r
spz
syz
100%Performs the symmetrical rank-1 update of a complex symmetric packed matrix.
*r1v
ladz
Computes the (scaled) r-th column of the inverse of the submatrix in rows b1 through bn of the tridiagonal matrix LDLT - λI.
*r2v
ladz
100%Applies a vector of plane rotations with real cosines and real sines from both sides to a sequence of 2-by-2 symmetric/Hermitian matrices.
zlarcm
laz
Copies all or part of a real two-dimensional array to a complex array.
*rf
ladz
100%Applies an elementary reflector to a general rectangular matrix.
*rfb
ladz
tpdz
61%Applies a block reflector or its transpose to a general rectangular matrix.
*rfb_gett
ladz
100%Applies a block Householder reflector to a triangular-pentagonal matrix.
*rfg
ladz
100%Generates an elementary reflector (Householder matrix).
*rfgp
ladz
91%Generates an elementary reflector (Householder matrix) with non-negative \(\beta\).
*rfs
gbdz
gedz
gtdz
hez
hpz
pbdz
podz
ppdz
ptdz
spdz
sydz
tbdz
tpdz
trdz
95%Improves the computed solution to a system of linear equations when the coefficient matrix is banded, and provides error bounds and backward error estimates for the solution.
*rfsx
gbdz
gedz
hez
podz
sydz
Improves the computed solution to a system of linear equations and provides error bounds and backward error estimates for the solution.
*rft
ladz
100%Forms the triangular factor T of a block reflector.
*rfx
ladz
98%Applies an elementary reflector to a general rectangular matrix, with loop unrolling when the reflector has order ≤ 10.
*rfy
ladz
100%Applies an elementary reflector, or Householder matrix, H, to an n x n symmetric matrix C, from both the left and the right.
*rgv
ladz
97%Generates a vector of plane rotations with real cosines and real sines.
*rk
hfz
sfd
100%Performs a symmetric rank-k operation for matrix in RFP format.
dlarmm
lad
100%Computes a scaling factor to prevent overflow in triangular matrix multiplication.
*rnv
ladz
100%Returns a vector of random numbers from a uniform or normal distribution.
zrot
z
100%Applies a plane rotation with real cosine and complex sine to a pair of complex vectors.
*rq2
gedz
100%Computes the RQ factorization of a general rectangular matrix using an unblocked algorithm.
*rqf
gedz
ggdz
100%Computes an RQ factorization of a real M-by-N matrix A: A = R * Q.
dlarra
lad
54%Computes the splitting points with the specified threshold.
dlarrb
lad
Provides limited bisection to locate eigenvalues for more accuracy.
dlarrc
lad
99%Computes the number of eigenvalues of the symmetric tridiagonal matrix.
dlarrd
lad
80%Computes the eigenvalues of a symmetric tridiagonal matrix to suitable accuracy.
dlarre
lad
Given the tridiagonal matrix T, sets small off-diagonal elements to zero and for each unreduced block Ti, finds base representations and eigenvalues.
dlarrf
lad
64%Finds a new relatively robust representation such that at least one of the eigenvalues is relatively isolated.
dlarrj
lad
100%Performs refinement of the initial estimates of the eigenvalues of the matrix T.
dlarrk
lad
98%Computes one eigenvalue of a symmetric tridiagonal matrix T to suitable accuracy.
dlarrr
lad
100%Performs tests to decide whether the symmetric tridiagonal matrix T warrants expensive computations which guarantee high relative accuracy in the eigenvalues.
*rrv
ladz
Computes the eigenvectors of the tridiagonal matrix T = L D LT given L, D and the eigenvalues of L D LT.
*rscl
dz
99%Multiplies a vector by the reciprocal of a real scalar.
*rscl2
ladz
100%Performs reciprocal diagonal scaling on a matrix.
*rtg
ladz
98%Generates a plane rotation with real cosine and real sine.
dlartgp
lad
97%Generates a plane rotation so that the diagonal is nonnegative.
dlartgs
lad
76%Generates a plane rotation designed to introduce a bulge in implicit QR iteration for the bidiagonal SVD problem.
*rtv
ladz
100%Applies a vector of plane rotations with real cosines and real sines to the elements of a pair of vectors.
dlaruv
lad
99%Returns a vector of n random real numbers from a uniform distribution.
*rz
ladz
99%Applies an elementary reflector (as returned by stzrzf) to a general matrix.
*rzb
ladz
Applies a block reflector or its transpose to a general matrix.
*rzf
tzdz
Reduces the M-by-N ( M<=N ) real upper trapezoidal matrix A to upper triangular form by means of orthogonal transformations.
*rzt
ladz
100%Forms the triangular factor T of a block reflector H = I - vtvH.
dlas2
lad
100%Computes singular values of a 2-by-2 triangular matrix.
lsamen
l
Tests if the first N letters of CA are the same as the first N letters of CB, regardless of case.
*sc2
gedz
97%Solves a system of linear equations using the LU factorization with complete pivoting computed by sgetc2.
*scl
ladz
100%Multiplies a general rectangular matrix by a real scalar defined as cto/cfrom.
*scl2
ladz
100%Performs diagonal scaling on a matrix.
dlasd0
lad
Computes the singular values of a real upper bidiagonal n-by-m matrix B with diagonal d and off-diagonal e. Used by sbdsdc.
dlasd1
lad
Computes the SVD of an upper bidiagonal matrix B of the specified size. Used by sbdsdc.
dlasd2
lad
98%Merges the two sets of singular values together into a single sorted set. Used by sbdsdc.
dlasd3
lad
Finds all square roots of the roots of the secular equation, as defined by the values in D and Z, and then updates the singular vectors by matrix multiplication. Used by sbdsdc.
dlasd4
lad
Computes the square root of the i-th updated eigenvalue of a positive symmetric rank-one modification to a positive diagonal matrix. Used by dbdsdc.
dlasd5
lad
100%Computes the square root of the i-th eigenvalue of a positive symmetric rank-one modification of a 2-by-2 diagonal matrix. Used by sbdsdc.
dlasd6
lad
Computes the SVD of an updated upper bidiagonal matrix obtained by merging two smaller ones by appending a row. Used by sbdsdc.
dlasd7
lad
98%Merges the two sets of singular values together into a single sorted set. Then it tries to deflate the size of the problem. Used by sbdsdc.
dlasd8
lad
Finds the square roots of the roots of the secular equation, and stores, for each element in D, the distance to its two nearest poles. Used by sbdsdc.
dlasda
lad
Computes the singular value decomposition (SVD) of a real upper bidiagonal matrix with diagonal d and off-diagonal e. Used by sbdsdc.
dbdsdc
bdd
Computes the SVD of a bidiagonal matrix using a divide and conquer method.
*sdd
gedz
Computes the singular value decomposition (SVD) of a real M-by-N matrix A, optionally computing the left and right singular vectors.
dlasdq
lad
Computes the SVD of a real bidiagonal matrix with diagonal d and off-diagonal e. Used by sbdsdc.
dlasdt
lad
100%Creates a tree of subproblems for bidiagonal divide and conquer. Used by sbdsdc.
*sen
tgdz
trdz
62%Reorders the generalized Schur decomposition and computes reciprocal condition numbers.
*set
ladz
100%Initializes the off-diagonal elements and the diagonal elements of a matrix to given values.
dsgesv
d
Computes the solution to a system of linear equations using mixed-precision iterative refinement.
*sja
tgdz
97%Computes the generalized singular value decomposition (GSVD) of two real upper triangular (or trapezoidal) matrices A and B.
ilaslc
lai
Scans a matrix for its last non-zero column.
ilaslr
lai
Scans a matrix for its last non-zero row.
*sm
tfdz
100%Solves a matrix equation (one operand is a triangular matrix in RFP format).
*sna
did
tgdz
trdz
45%Computes the reciprocal condition numbers for the eigenvectors of a real symmetric or complex Hermitian matrix or for the left or right singular vectors of a general m-by-n matrix.
dlasq1
lad
91%Computes the singular values of a real square bidiagonal matrix. Used by sbdsqr.
dlasq2
lad
87%Computes all the eigenvalues of the symmetric positive definite tridiagonal matrix associated with the qd Array Z to high relative accuracy. Used by sbdsqr and sstegr.
dlasq3
lad
90%Checks for deflation, computes a shift and calls dqds. Used by sbdsqr.
dlasq4
lad
98%Computes an approximation to the smallest eigenvalue using values of d from the previous transform. Used by sbdsqr.
dlasq5
lad
100%Computes one dqds transform in ping-pong form. Used by sbdsqr and sstegr.
dlasq6
lad
100%Computes one dqd transform in ping-pong form. Used by sbdsqr and sstegr.
*sqr
bddz
95%Computes singular values and vectors from the SVD of a bidiagonal matrix.
*sr
ladz
100%Applies a sequence of plane rotations to a general rectangular matrix.
dlasrt
lad
98%Sorts numbers in increasing or decreasing order.
*ssq
ladz
100%Updates a sum of squares represented in scaled form.
*stf
pbdz
97%Computes a split Cholesky factorization of a real symmetric positive definite band matrix A.
dzsum1
z
100%Forms the 1-norm of the complex vector using the true absolute value.
*sv
gbdz
gedz
gtdz
hez
hpz
pbdz
podz
ppdz
ptdz
spdz
sydz
100%Computes the solution to system of linear equations \(AX = B\) for GB matrices (simple driver)
dlasv2
lad
100%Computes the singular value decomposition of a 2-by-2 triangular matrix.
*sv_aa
hez
sydz
Computes the solution to system of linear equations \(AX = B\) for SY matrices
*sv_aa_2stage
hez
sydz
Computes the solution to system of linear equations \(AX = B\) for SY matrices
*sv_rk
hez
sydz
Computes the solution to system of linear equations \(AX = B\) for SY matrices
*sv_rook
hez
sydz
Computes the solution to system of linear equations \(AX = B\) for SY matrices
*svd
gedz
99%Computes the singular value decomposition (SVD) for GE matrices
*svd3
ggdz
98%Computes the singular value decomposition (SVD) for OTHER matrices
*svdq
gedz
Computes the singular value decomposition (SVD) with a QR-Preconditioned QR SVD Method for GE matrices
*svdx
bdd
gedz
Computes selected singular values and vectors of a bidiagonal matrix.
*svj
gedz
Computes the singular value decomposition (SVD) of a real M-by-N matrix A, where M >= N.
*svp3
ggdz
67%Computes the generalized singular value decomposition (GSVD).
*svx
gbdz
gedz
gtdz
hez
hpz
pbdz
podz
ppdz
ptdz
spdz
sydz
97%Computes the solution to system of linear equations \(AX = B\) for GB matrices
*svxx
gbdz
gedz
hez
podz
sydz
Computes the solution to system of linear equations \(AX = B\) for GB matrices
*swapr
hez
sydz
100%Applies an elementary permutation on the rows and columns of a symmetric matrix.
*swlq
ladz
Computes a blocked short-wide LQ factorization.
*swp
ladz
100%Performs a series of row interchanges on a general rectangular matrix.
*sy2
lad
tgdz
92%Solves the Sylvester matrix equation where the matrices are of order 1 or 2.
*syf
ladz
92%Computes a partial factorization of a real symmetric matrix using the Bunch-Kaufman diagonal pivoting method.
*syf_aa
ladz
DLATRF_AA factorizes a panel of a real symmetric matrix A using the Aasen's algorithm.
*syf_rk
ladz
Computes a partial factorization of a real symmetric indefinite matrix using bounded Bunch-Kaufman (rook) diagonal pivoting method.
*syf_rook
ladz
Computes a partial factorization of a symmetric matrix using rook pivoting (blocked).
*syl
tgdz
trdz
70%Solves the generalized Sylvester equation.
*syl3
trdz
Solves the Sylvester matrix equation (blocked algorithm).
zlat2c
laz
100%Converts a double complex triangular matrix to a complex triangular matrix.
dlat2s
lad
100%Converts a double-precision triangular matrix to a single-precision triangular matrix.
*tbs
ladz
54%Solves a triangular banded system of equations.
*tc2
gedz
97%Computes the LU factorization with complete pivoting of the general n-by-n matrix.
*td2
hez
syd
100%Reduces a symmetric matrix to real symmetric tridiagonal form by an orthogonal similarity transformation (unblocked algorithm).
*tdf
ladz
99%Uses the LU factorization of the n-by-n matrix computed by sgetc2 and computes a contribution to the reciprocal Dif-estimate.
*tf2
gbdz
gedz
hez
pbdz
podz
psdz
sydz
98%Computes the LU factorization of a general band matrix using the unblocked version of the algorithm.
*tf2_rk
hez
sydz
Computes the factorization of a real symmetric indefinite matrix using the bounded Bunch-Kaufman (rook) diagonal pivoting method (BLAS2 unblocked algorithm).
*tf2_rook
hez
sydz
Computes the factorization of a real symmetric indefinite matrix using the bounded Bunch-Kaufman ("rook") diagonal pivoting method (unblocked algorithm).
*ti2
trdz
100%Computes the inverse of a triangular matrix (unblocked algorithm).
*tps
ladz
18%Solves a triangular system of equations with the matrix held in packed storage.
ilatrans
lai
This subroutine translates from a character string specifying a transposition operation to the relevant BLAST-specified integer constant.
*trd
hbz
hez
hpz
ladz
sbd
spd
syd
97%Reduces the first nb rows and columns of a symmetric/Hermitian matrix A to real tridiagonal form by an orthogonal similarity transformation.
*trd_2stage
hez
syd
Reduces a real symmetric matrix A to real symmetric tridiagonal form T by a orthogonal similarity transformation: Q1**T Q2**T* A * Q2 * Q1 = T.
zhetrd_he2hb
hez
Reduces a complex Hermitian matrix A to complex Hermitian band-diagonal form AB by a unitary similarity transformation: \(Q^H\) * A * Q = AB.
dsytrd_sy2sb
syd
Reduces a real symmetric matrix A to real symmetric band-diagonal form AB by a orthogonal similarity transformation: \(Q^T\) * A * Q = AB.
*trf
gbdz
gedz
gtdz
hez
hpz
pbdz
pfdz
podz
ppdz
psdz
ptdz
spdz
sydz
96%Computes an LU factorization of a real m-by-n band matrix A using partial pivoting with row interchanges.
*trf2
gedz
podz
100%Computes an LU factorization of a general M-by-N matrix A using partial pivoting with row interchanges.
*trf_aa
hez
sydz
Computes the factorization of a real symmetric matrix A using the Aasen's algorithm.
*trf_aa_2stage
hez
sydz
Computes the factorization of a real symmetric matrix A using the Aasen's algorithm.
*trf_rk
hez
sydz
Computes the factorization of a real symmetric indefinite matrix using the bounded Bunch-Kaufman (rook) diagonal pivoting method (BLAS3 blocked algorithm).
*trf_rook
hez
sydz
Computes the factorization of a real symmetric matrix A using the bounded Bunch-Kaufman ("rook") diagonal pivoting method.
*tri
gedz
hez
hpz
pfdz
podz
ppdz
spdz
sydz
tfdz
tpdz
trdz
90%Computes the inverse of a matrix using the LU factorization computed by DGETRF.
*tri2
hez
sydz
Computes the inverse of a DOUBLE PRECISION symmetric indefinite matrix A using the factorization \(A = UDU^T\) or \(A = LDL^T\) computed by DSYTRF.
*tri2x
hez
sydz
Computes the inverse of a real symmetric indefinite matrix A using the factorization \(A = UDU^T\) or \(A = LDL^T\) computed by DSYTRF.
*tri_3
hez
sydz
Computes the inverse of a symmetric/Hermitian indefinite matrix using the bounded Bunch-Kaufman factorization.
*tri_3x
hez
sydz
Computes the inverse of a symmetric/Hermitian indefinite matrix using the bounded Bunch-Kaufman factorization (expert).
*tri_rook
hez
sydz
Computes the inverse of a real symmetric matrix A using the factorization \(A = UDU^T\) or \(A = LDL^T\) computed by DSYTRF_ROOK.
*trs
gbdz
gedz
gtdz
hez
hpz
ladz
pbdz
pfdz
podz
ppdz
ptdz
spdz
sydz
tbdz
tpdz
trdz
98%Solves a system of linear equations \(AX = B\) or \(A^TX = B\) with a general band matrix A using the LU factorization computed by DGBTRF.
*trs2
hez
sydz
97%Solves a system of linear equations \(AX = B\) with a real symmetric matrix A using the factorization \(A = UDU^T\) or \(A = LDL^T\) computed by DSYTRF and converted by DSYCONV.
*trs3
ladz
Solves a triangular system of equations with the scale factors set to prevent overflow.
*trs_3
hez
sydz
Solves a system of linear equations using the bounded Bunch-Kaufman factorization.
*trs_aa
hez
sydz
Solves a system of linear equations \(AX = B\) with a real symmetric matrix A using the factorization \(A = U^TTU\) or \(A = LTL^T\) computed by DSYTRF_AA.
*trs_aa_2stage
hez
sydz
Solves a system of linear equations \(AX = B\) with a real symmetric matrix A using the factorization \(A = U^TTU\) or \(A = LTL^T\) computed by DSYTRF_AA_2STAGE.
*trs_rook
hez
sydz
Solves a system of linear equations \(AX = B\) with a real symmetric matrix A using the factorization \(A = UDU^T\) or \(A = LDL^T\) computed by DSYTRF_ROOK.
*trz
ladz
Factors an upper trapezoidal matrix by means of orthogonal transformations.
*ts2
gtdz
ptdz
99%Solves a system of linear equations with a tridiagonal matrix using the LU factorization computed by sgttrf.
*tsls
gedz
Solves overdetermined or underdetermined real linear systems involving an M-by-N matrix A, using a tall skinny QR or short wide LQ factorization of A.
*tsqr
ladz
Computes a blocked tall-skinny QR factorization.
*tsqrhrt
gedz
Computes an upper-triangular factor of an orthogonal/unitary matrix from a tall-skinny QR factorization.
*ttf
tpdz
trdz
100%Copies a triangular matrix from the standard packed format (TP) to the rectangular full packed format (TF).
*ttp
tfdz
trdz
100%Copies a triangular matrix from the rectangular full packed format (TF) to the standard packed format (TP).
*ttr
tfdz
tpdz
100%Copies a triangular matrix from the rectangular full packed format (TF) to the standard full format (TR).
zlaunhr_col_getrfnp
laz
Computes the modified LU factorization without pivoting of a complex general M-by-N matrix A.
zlaunhr_col_getrfnp2
laz
Computes the modified LU factorization without pivoting of a complex general M-by-N matrix A.
ilauplo
lai
This subroutine translated from a character string specifying a upper- or lower-triangular matrix to the relevant BLAST-specified integer constant.
*uu2
ladz
100%Computes the product UUH or LHL, where U and L are upper or lower triangular matrices (unblocked algorithm).
*uum
ladz
100%Computes the product UUH or LHL, where U and L are upper or lower triangular matrices (blocked algorithm).
*vj0
gsdz
59%Pre-processor for the routine dgesvj.
*vj1
gsdz
Pre-processor for the routine dgesvj, applies Jacobi rotations targeting only particular pivots.
ilazlc
lai
100%Scans a matrix for its last non-zero column.
ilazlr
lai
100%Scans a matrix for its last non-zero row.