single sorted set
slaed3		slaed3 (3p)	- find the roots of the secular equation, as defined by the values in D, W, and RHO, between KSTART and KSTOP
slaed4		slaed4 (3p)	- subroutine computes the I-th updated eigenvalue of a symmetric rank-one modification to a diagonal matrix whose elements are given in the array d, and that	D(i) < D(j) for i < j  and that RHO > 0
slaed5		slaed5 (3p)	- subroutine computes the I-th eigenvalue of a symmetric rank-one modification of a 2-by-2 diagonal matrix   diag( D ) + RHO  The diagonal elements in the array D are assumed to satisfy   D(i) < D(j) for i < j 
slaed7		slaed7 (3p)	- compute the updated eigensystem of a diagonal matrix after modification by a rank-one symmetric matrix
slaed8		slaed8 (3p)	- merge the two sets of eigenvalues together into a single sorted set
slaed9		slaed9 (3p)	- find the roots of the secular equation, as defined by the values in D, Z, and RHO, between KSTART and KSTOP
slaeda		slaeda (3p)	- compute the Z vector corresponding to the merge step in the CURLVLth step of the merge process with TLVLS steps for the CURPBMth problem
slaein		slaein (3p)	- use inverse iteration to find a right or left eigenvector corresponding to the eigenvalue (WR,WI) of a real upper Hessenberg matrix H
slaev2		slaev2 (3p)	- compute the eigendecomposition of a 2-by-2 symmetric matrix  [ A B ]	[ B C ]
slaexc		slaexc (3p)	- swap adjacent diagonal blocks T11 and T22 of order 1 or 2 in an upper quasi-triangular matrix T by an orthogonal similarity transformation
slahqr		slahqr (3p)	- i an auxiliary routine called by SHSEQR to update the eigenvalues and Schur decomposition already computed by SHSEQR, by dealing with the Hessenberg submatrix in rows and columns ILO to IHI
slahrd		slahrd (3p)	- reduce the first NB columns of a real general n-by-(n-k+1) matrix A so that elements below the k-th subdiagonal are zero
slaic1		slaic1 (3p)	- applie one step of incremental condition estimation in its simplest version
slamch		slamch (3p)	- determine single precision machine parameters
slamrg		slamrg (3p)	- will create a permutation list which will merge the elements of A (which is composed of two independently sorted sets) into a single set which is sorted in ascending order
slangb		slangb (3p)	- return the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of an n by n band matrix A, with kl sub-diagonals and ku super-diagonals
slange		slange (3p)	- return the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a real matrix A
slangt		slangt (3p)	- return the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a real tridiagonal matrix A
slanhs		slanhs (3p)	- return the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a Hessenberg matrix A
slansb		slansb (3p)	- return the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of an n by n symmetric band matrix A, with k super-diagonals
slansp		slansp (3p)	- return the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a real symmetric matrix A, supplied in packed form
slanst		slanst (3p)	- return the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a real symmetric tridiagonal matrix A
slansy		slansy (3p)	- return the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a real symmetric matrix A
slantb		slantb (3p)	- return the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of an n by n triangular band matrix A, with ( k + 1 ) diagonals
slantp		slantp (3p)	- return the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a triangular matrix A, supplied in packed form
slantr		slantr (3p)	- return the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a trapezoidal or triangular matrix A
slanv2		slanv2 (3p)	- compute the Schur factorization of a real 2-by-2 nonsymmetric matrix in standard form
slapll		slapll (3p)	- two column vectors X and Y, let   A = ( X Y )
slapmt		slapmt (3p)	- rearrange the columns of the M by N matrix X as specified by the permutation K(1),K(2),...,K(N) of the integers 1,...,N
slapy2		slapy2 (3p)	- return sqrt(x**2+y**2), taking care not to cause unnecessary overflow
slapy3		slapy3 (3p)	- return sqrt(x**2+y**2+z**2), taking care not to cause unnecessary overflow
slaqgb		slaqgb (3p)	- equilibrate a general M by N band matrix A with KL subdiagonals and KU superdiagonals using the row and scaling factors in the vectors R and C
slaqge		slaqge (3p)	- equilibrate a general M by N matrix A using the row and scaling factors in the vectors R and C
slaqsb		slaqsb (3p)	- equilibrate a symmetric band matrix A using the scaling factors in the vector S
slaqsp		slaqsp (3p)	- equilibrate a symmetric matrix A using the scaling factors in the vector S
slaqsy		slaqsy (3p)	- equilibrate a symmetric matrix A using the scaling factors in the vector S
slaqtr		slaqtr (3p)	- solve the real quasi-triangular system   op(T)*p = scale*c, if LREAL = .TRUE
slar2v		slar2v (3p)	- applie a vector of real plane rotations from both sides to a sequence of 2-by-2 real symmetric matrices, defined by the elements of the vectors x, y and z
slarf		slarf (3p)	- applie a real elementary reflector H to a real m by n matrix C, from either the left or the right
slarfb		slarfb (3p)	- applie a real block reflector H or its transpose H' to a real m by n matrix C, from either the left or the right
slarfg		slarfg (3p)	- generate a real elementary reflector H of order n, such that	 H * ( alpha ) = ( beta ), H' * H = I
slarft		slarft (3p)	- form the triangular factor T of a real block reflector H of order n, which is defined as a product of k el