Museum

Home

Lab Overview

Retrotechnology Articles

Online Manuals

⇒ Extended Math Library 3.4

Media Vault

Software Library

Restoration Projects

Artifacts Sought

Anchors

Section 3DXML

Section 3dxml

Section 3lapack

Section l

Manual — Extended Math Library 3.4

1521 entries

Section 3DXML

array-mathA library of linear algebra routines
blas1A library of linear algebra routines
blas1eA library of linear algebra routines
blas1sA library of linear algebra routines
blas2A library of linear algebra routines
blas3A library of linear algebra routines
caxpyi
sdoti, ddoti, cdotui, zdotui, zdotci, cdotciInner product of a vector and a sparse vector
sdoti, ddoti, zdotui, cdotci, zdotci, cdotuiInner product of a vector and a sparse vector
sgthr, dgthr, zgthr, cgthrGathers the specified elements of a vector
cgthrs
sgthrz, dgthrz, zgthrz, cgthrzGathers and zeros specified elements of a vector
ssctr, dsctr, zsctr, csctrScatters the elements of a sparse vector
ssctrs, dsctrs, zsctrs, csctrsScales and scatters the elements of a sparse vector
ssumi, dsumi, zsumi, csumiSum of a vector and a sparse vector
dapply_diag_allApply diagonal preconditioner for any storage scheme (Serial and Parallel Versions)
dapply_ilu_genr_lApply incomplete LU preconditioner for general storage by rows
dapply_ilu_genr_uApply incomplete LU preconditioner for general storage by rows
dapply_ilu_sdiaApply ILU preconditioner for symmetric diagonal storage
dapply_ilu_udia_lApply ILU preconditioner for unsymmetric diagonal storage
dapply_ilu_udia_uApply ILU preconditioner for unsymmetric diagonal storage
dapply_poly_genrApply polynomial preconditioner for general storage by rows (Serial and Parallel Versions)
dapply_poly_sdiaApply polynomial preconditioner for symmetric diagonal storage (Serial and Parallel Versions)
dapply_poly_udiaApply polynomial preconditioner for unsymmetric diagonal storage (Serial and Parallel Versions)
daxpyi
dcreate_diag_genrGenerate diagonal preconditioner for general storage by rows (Serial and Parallel Versions)
dcreate_diag_sdiaGenerate diagonal preconditioner for symmetric diagonal storage (Serial and Parallel Versions)
dcreate_diag_udiaGenerate diagonal preconditioner for unsymmetric diagonal storage (Serial and Parallel Versions)
dcreate_ilu_genrGenerate incomplete LU preconditioner for general storage by rows
dcreate_ilu_sdiaGenerate incomplete Cholesky preconditioner for symmetric diagonal storage
dcreate_ilu_udiaGenerate incomplete LU preconditioner for unsymmetric diagonal storage
dcreate_poly_genrGenerate polynomial preconditioner for general storage by rows (Serial and Parallel Versions)
dcreate_poly_sdiaGenerate polynomial preconditioner for symmetric diagonal storage (Serial and Parallel Versions)
dcreate_poly_udiaGenerate polynomial preconditioner for unsymmetric diagonal storage (Serial and Parallel Versions)
sdoti, cdotui, zdotui, cdotci, zdotci, ddotiInner product of a vector and a sparse vector
sgthr, cgthr, zgthr, dgthrGathers the specified elements of a vector
dgthrs
sgthrz, cgthrz, zgthrz, dgthrzGathers and zeros specified elements of a vector
ditsol_defaultsSet default values for iterative solver
ditsol_driverDriver for sparse iterative solvers (Serial and Parallel Versions)
ditsol_pbcgPreconditioned bi-conjugate gradient method (Serial and Parallel Versions)
ditsol_pcgPreconditioned conjugate gradient method (Serial and Parallel Versions)
ditsol_pcgsPreconditioned conjugate gradient squared method (Serial and Parallel Versions)
ditsol_pgmresPreconditioned generalized minimum residual method (Serial and Parallel Versions)
ditsol_plscgPreconditioned least square conjugate gradient method (Serial and Parallel Versions)
ditsol_ptfqmrPreconditioned transpose-free quasi-minimal method (Serial and Parallel Versions)
dmatvec_genrMatrix-vector product for general storage by rows (Serial and Parallel Versions)
dmatvec_sdiaMatrix-vector product for symmetric diagonal storage (Serial and Parallel Versions)
dmatvec_udiaMatrix-vector product for unsymmetric diagonal storage (Serial and Parallel Versions)
dpcgsPreconditioned conjugate gradient squared method
sroti, drotiReal givens plane rotation applied to sparse vector
ssctr, csctr, zsctr, dsctrScatters the elements of a sparse vector
ssctrs, csctrs, zsctrs, dsctrsScales and scatters the elements of a sparse vector
dsskycSymmetric sparse matrix condition number estimator using skyline storage scheme
dsskyd
dsskyfSymmetric sparse matrix factorization using skyline storage scheme (Serial and Parallel Versions)
dsskynSymmetric sparse matrix norm evaluation using skyline storage scheme
dsskyrSymmetric sparse iterative refinement using skyline storage scheme
dsskysSymmetric sparse matrix solve using skyline storage scheme
dsskyxSymmetric sparse expert driver using skyline storage scheme
ssumi, csumi, zsumi, dsumiSum of a vector and a sparse vector
duskycUnsymmetric sparse matrix condition number estimation using skyline storage scheme
duskydUnsymmetric sparse simple driver using skyline storage scheme
duskyf
duskynUnsymmetric sparse matrix norm evaluation using skyline storage scheme
duskyrUnsymmetric sparse iterative refinement using skyline storage scheme
duskysUnsymmetric sparse matrix solve using skyline storage scheme
duskyxUnsymmetric sparse expert driver using skyline storage scheme
iterative-solversA library of sparse linear solvers (iterative)
lapackA library of linear algebra routines
random-numbersRandom number generator subprograms
saxpyi
ddoti, cdotui, zdotui, cdotci, zdotci, sdotiInner product of a vector and a sparse vector
dgthr, cgthr, zgthr, sgthrGathers the specified elements of a vector
sgthrs
dgthrz, cgthrz, zgthrz, sgthrzGathers and zeros specified elements of a vector
signal-processingA library of signal processing routines
skyline-solversA library of sparse linear solvers (direct)
droti, srotiReal givens plane rotation applied to sparse vector
dsctr, csctr, zsctr, ssctrScatters the elements of a sparse vector
dsctrs, csctrs, zsctrs, ssctrsScales and scatters the elements of a sparse vector
dsumi, csumi, zsumi, ssumiSum of a vector and a sparse vector
vcosVector cosine
vcos_sinVector cosine and sine
vexpVector exponential
vlogVector logarithm
vrecipVector reciprocal
vsinVector sine
vsqrtVector square root
zaxpyi
sdoti, ddoti, cdotui, zdotui, cdotci, zdotciInner product of a vector and a sparse vector
sdoti, ddoti, cdotui, cdotci, zdotci, zdotuiInner product of a vector and a sparse vector
sgthr, dgthr, cgthr, zgthrGathers the specified elements of a vector
zgthrs
sgthrz, dgthrz, cgthrz, zgthrzGathers and zeros specified elements of a vector
ssctr, dsctr, csctr, zsctrScatters the elements of a sparse vector
ssctrs, dsctrs, csctrs, zsctrsScales and scatters the elements of a sparse vector
ssumi, dsumi, csumi, zsumiSum of a vector and a sparse vector

Section 3dxml

saxpy, daxpy, zaxpy, caxpyVector plus the product of a scalar and a vector
sconv_nonperiodic, dconv_nonperiodic, zconv_nonperiodic, cconv_nonperiodicNonperiodic convolution
sconv_nonperiodic_ext, dconv_nonperiodic_ext, zconv_nonperiodic_ext, cconv_nonperiodic_extExtended nonperiodic convolution
sconv_periodic, dconv_periodic, zconv_periodic, cconv_periodicPeriodic concolution
sconv_periodic_ext, dconv_periodic_ext, zconv_periodic_ext, cconv_periodic_extExtended periodic convolution
scopy, dcopy, zcopy, ccopyCopy of a vector
ccorr_nonperiodic
ccorr_nonperiodic_ext
ccorr_periodic
scorr_periodic_ext, dcorr_periodic_ext, zcorr_periodic_ext, ccorr_periodic_extExtended periodic correlation
sdot, ddot, dsdot, zdotc, cdotu, zdotu, cdotcINNER PRODUCT OF TWO VECTORS
sdot, ddot, dsdot, cdotc, zdotc, zdotu, cdotuINNER PRODUCT OF TWO VECTORS
sfft, dfft, zfft, cfftFast fourier transform in one dimension
sfft_2d, dfft_2d, zfft_2d, cfft_2dFast fourier transform in two dimensions
sfft_3d, dfft_3d, zfft_3d, cfft_3dFast fourier transform in three dimensions
sfft_apply, dfft_apply, zfft_apply, cfft_applyApplication step for fast fourier transform in one dimension
cfft_apply_2d
sfft_apply_3d, dfft_apply_3d, zfft_apply_3d, cfft_apply_3dApplication step for fast fourier transform in three dimensions
sfft_apply_grp, dfft_apply_grp, zfft_apply_grp, cfft_apply_grpApplication step for group fast fourier transform in one dimension
sfft_exit, dfft_exit, zfft_exit, cfft_exitFinal step for fast fourier transform in one dimension
sfft_exit_2d, dfft_exit_2d, zfft_exit_2d, cfft_exit_2dFinal step for fast fourier transform in two dimensions
sfft_exit_3d, dfft_exit_3d, zfft_exit_3d, cfft_exit_3dFinal step for fast fourier transfrom in three dimension
sfft_exit_grp, dfft_exit_grp, zfft_exit_grp, cfft_exit_grpExit step for group fast fourier transform in one dimension
sfft_grp, dfft_grp, zfft_grp, cfft_grpGroup fast fourier transform in one dimension
sfft_init, dfft_init, zfft_init, cfft_initInitialization step for fast fourier transform in one dimension
sfft_init_2d, dfft_init_2d, zfft_init_2d, cfft_init_2dInitialization step for fast fourier transform in two dimensions
sfft_init_3d, dfft_init_3d, zfft_init_3d, cfft_init_3dInitialization step for fast fourier transform in three dimension
sfft_init_grp, dfft_init_grp, zfft_init_grp, cfft_init_grpInitialization step for group fast fourier transform in one dimension
sgbmv, ddbmv, zgbmv, cgbmvMatrix-vector product for a general band matrix
sgema, dgema, zgema, cgemaMatrix-matrix addition
sgemm, dgemm, zgemm, cgemmMatrix-matrix product and addition
sgems, dgems, zgems, cgemsMatrix-matrix subtraction
sgemt, dgemt, zgemt, cgemtMatrix-matrix copy
sgemv, dgemv, zgemv, cgemvMatrix-vector product for a general matrix
sger, dger, zgerc, cgeru, zgeru, cgercRank-one update of a general matrix
sger, dger, cgerc, zgerc, zgeru, cgeruRank-one update of a general matrix
ssbmv, dsbmv, zhbmv, chbmvMatrix-vector product for a symmetric or hermitian band matrix
ssymm, dsymm, csymm, zsymm, zhemm, chemmMatrix-matrix product and addition for a symmetric or hermitian matrix
ssymv, dsymv, zhemv, chemvMatrix-vector product for a symmetric or hermitian matrix
ssyr, dsyr, zher, cherRank-one update of a symmetric or hermitian matrix
ssyr2, dsyr2, zher2, cher2Rank-two update of a symmetric or hermitian matrix
zher2k, cher2kRank-2k update of a complex hermitian matrix
zherk, cherkRank-k update of a complex hermitian matrix
sspmv, dspmv, zhpmv, chpmvMatrix-vector product for a symmetric or hermitian matrix stored in packed form
sspr, dspr, zhpr, chprRank-one update of a symmetric or hermitian matrix stored in packed form
sspr2, dspr2, zhpr2, chpr2Rank-two update of a symmetric or hermitian matrix stored in packed form
srot, drot, zrot, csrot, zdrot, crotApply givens plane rotation
crotg
sscal, dscal, zscal, csscal, zdscal, cscalProduct of a scalar and a vector
sset, dset, zset, csetSet all elements of a vector to a scalar
srot, drot, crot, zrot, zdrot, csrotApply givens plane rotation
sscal, dscal, cscal, zscal, zdscal, csscalProduct of a scalar and a vector
ssum, dsum, zsum, csumSum of the values of the elements of a vector
svcal, dvcal, cvcal, zvcal, zdvcal, csvcalProduct of a scalar and a vector
sswap, dswap, zswap, cswapExchange the elements of two vectors
ssymm, dsymm, zsymm, chemm, zhemm, csymmMatrix-matrix product and addition for a symmetric or hermitian matrix
ssyr2k, dsyr2k, zsyr2k, csyr2kRank-2k update of a symmetric matrix
ssyrk, dsyrk, zsyrk, csyrkRank-k update of a symmetric matrix
stbmv, dtbmv, ztbmv, ctbmvMatrix-vector product for a triangular band matrix
stbsv, dtbsv, ztbsv, ctbsvSolver of a system of linear equations with a triangular band matrix
stpmv, dtpmv, ztpmv, ctpmvMatrix-vector product for a triangular matrix in packed form
stpsv, dtpsv, ztpsv, ctpsvSolve a system of linear equations with a triangular matrix in packed form
strmm, dtrmm, ztrmm, ctrmmMatrix-matrix product for triangular matrix
strmv, dtrmv, ztrmv, ctrmvMarix-vector product for a triangular matrix
strsm, dtrsm, ztrsm, ctrsmSolve a triangular system of equations with a triangular coefficient matrix
strsv, dtrsv, ztrsv, ctrsvSolver of a system of linear equations with a triangular matrix
svcal, dvcal, zvcal, csvcal, zdvcal, cvcalProduct of a scalar and a vector
szaxpy, dzaxpy, zzaxpy, czaxpyVector plus the product of a scalar and a vector
samax, scamax, dzamax, damaxMaximum absolute value
damin
sasum, scasum, dzasum, dasumSum of the absolute value
saxpy, caxpy, zaxpy, daxpyVector plus the product of a scalar and a vector
sconv_nonperiodic, cconv_nonperiodic, zconv_nonperiodic, dconv_nonperiodicNonperiodic convolution
sconv_nonperiodic_ext, cconv_nonperiodic_ext, zconv_nonperiodic_ext, dconv_nonperiodic_extExtended nonperiodic convolution
sconv_periodic, cconv_periodic, zconv_periodic, dconv_periodicPeriodic concolution
sconv_periodic_ext, cconv_periodic_ext, zconv_periodic_ext, dconv_periodic_extExtended periodic convolution
scopy, ccopy, zcopy, dcopyCopy of a vector
dcorr_nonperiodic
dcorr_nonperiodic_ext
dcorr_periodic
scorr_periodic_ext, ccorr_periodic_ext, zcorr_periodic_ext, dcorr_periodic_extExtended periodic correlation
sdot, dsdot, cdotc, zdotc, cdotu, zdotu, ddotINNER PRODUCT OF TWO VECTORS
sfct, dfctFast cosine transform in one dimension
sfct_apply, dfct_applyApplication step for fast cosine transform in one dimension
sfct_exit, dfct_exitFinal step for fast cosine transform in one dimension
sfct_init, dfct_initInitialization step for fast cosine transform in one dimension
sfft, cfft, zfft, dfftFast fourier transform in one dimension
sfft_2d, cfft_2d, zfft_2d, dfft_2dFast fourier transform in two dimensions
sfft_3d, cfft_3d, zfft_3d, dfft_3dFast fourier transform in three dimensions
sfft_apply, cfft_apply, zfft_apply, dfft_applyApplication step for fast fourier transform in one dimension
dfft_apply_2d
sfft_apply_3d, cfft_apply_3d, zfft_apply_3d, dfft_apply_3dApplication step for fast fourier transform in three dimensions
sfft_apply_grp, cfft_apply_grp, zfft_apply_grp, dfft_apply_grpApplication step for group fast fourier transform in one dimension
sfft_exit, cfft_exit, zfft_exit, dfft_exitFinal step for fast fourier transform in one dimension
sfft_exit_2d, cfft_exit_2d, zfft_exit_2d, dfft_exit_2dFinal step for fast fourier transform in two dimensions
sfft_exit_3d, cfft_exit_3d, zfft_exit_3d, dfft_exit_3dFinal step for fast fourier transfrom in three dimension
sfft_exit_grp, cfft_exit_grp, zfft_exit_grp, dfft_exit_grpExit step for group fast fourier transform in one dimension
sfft_grp, cfft_grp, zfft_grp, dfft_grpGroup fast fourier transform in one dimension
sfft_init, cfft_init, zfft_init, dfft_initInitialization step for fast fourier transform in one dimension
sfft_init_2d, cfft_init_2d, zfft_init_2d, dfft_init_2dInitialization step for fast fourier transform in two dimensions
sfft_init_3d, cfft_init_3d, zfft_init_3d, dfft_init_3dInitialization step for fast fourier transform in three dimension
sfft_init_grp, cfft_init_grp, zfft_init_grp, dfft_init_grpInitialization step for group fast fourier transform in one dimension
sfst, dfstFast sine transform in one dimension
sfst_apply, dfst_applyApplication step for fast sine transform in one dimension
sfst_exit, dfst_exitFinal step for fast sine transform in one dimension
sfst_init, dfst_initInitialization step for fast sine transform in one dimension
sgbmv, ddbmv, cgbmv, zgbmv, dgbmvMatrix-vector product for a general band matrix
sgema, cgema, zgema, dgemaMatrix-matrix addition
sgemm, cgemm, zgemm, dgemmMatrix-matrix product and addition
sgems, cgems, zgems, dgemsMatrix-matrix subtraction
sgemt, cgemt, zgemt, dgemtMatrix-matrix copy
sgemv, cgemv, zgemv, dgemvMatrix-vector product for a general matrix
sger, cgerc, zgerc, cgeru, zgeru, dgerRank-one update of a general matrix
smax, dmaxLargest element in a real vector
smin, dminMinimum value of the elements of a real vector
snorm2, scnorm2, dznorm2, dnorm2Square root of sum of the squares of the elements of a vector
snrm2, scnrm2, dznrm2, dnrm2Square root of sum of the squares of the elements of a vector
snrsq, scnrsq, dznrsq, dnrsqSum of the squares of the elements of a vector
srot, crot, zrot, csrot, zdrot, drotApply givens plane rotation
drotg
srotm, drotmApply modified givens transformation
srotmg, drotmgGenerate elements for a modified Givens transform
ssbmv, chbmv, zhbmv, dsbmvMatrix-vector product for a symmetric or hermitian band matrix
sscal, cscal, zscal, csscal, zdscal, dscalProduct of a scalar and a vector
sdot, ddot, cdotc, zdotc, cdotu, zdotu, dsdotINNER PRODUCT OF TWO VECTORS
sset, cset, zset, dsetSet all elements of a vector to a scalar
isortq, ssortq, dsortqSort the elements of a vector
isortqx, ssortqx, dsortqxPerforms an indexed sort of a vector
sspmv, chpmv, zhpmv, dspmvMatrix-vector product for a symmetric or hermitian matrix stored in packed form
sspr, chpr, zhpr, dsprRank-one update of a symmetric or hermitian matrix stored in packed form
sspr2, chpr2, zhpr2, dspr2Rank-two update of a symmetric or hermitian matrix stored in packed form
ssum, csum, zsum, dsumSum of the values of the elements of a vector
sswap, cswap, zswap, dswapExchange the elements of two vectors
ssymm, csymm, zsymm, chemm, zhemm, dsymmMatrix-matrix product and addition for a symmetric or hermitian matrix
ssymv, chemv, zhemv, dsymvMatrix-vector product for a symmetric or hermitian matrix
ssyr, cher, zher, dsyrRank-one update of a symmetric or hermitian matrix
ssyr2, cher2, zher2, dsyr2Rank-two update of a symmetric or hermitian matrix
ssyr2k, csyr2k, zsyr2k, dsyr2kRank-2k update of a symmetric matrix
ssyrk, csyrk, zsyrk, dsyrkRank-k update of a symmetric matrix
stbmv, ctbmv, ztbmv, dtbmvMatrix-vector product for a triangular band matrix
stbsv, ctbsv, ztbsv, dtbsvSolver of a system of linear equations with a triangular band matrix
stpmv, ctpmv, ztpmv, dtpmvMatrix-vector product for a triangular matrix in packed form
stpsv, ctpsv, ztpsv, dtpsvSolve a system of linear equations with a triangular matrix in packed form
strmm, ctrmm, ztrmm, dtrmmMatrix-matrix product for triangular matrix
strmv, ctrmv, ztrmv, dtrmvMarix-vector product for a triangular matrix
strsm, ctrsm, ztrsm, dtrsmSolve a triangular system of equations with a triangular coefficient matrix
strsv, ctrsv, ztrsv, dtrsvSolver of a system of linear equations with a triangular matrix
svcal, cvcal, zvcal, csvcal, zdvcal, dvcalProduct of a scalar and a vector
dxmlA library of linear algebra and signal processing routines
samax, damax, scamax, dzamaxMaximum absolute value
dzamin
sasum, dasum, scasum, dzasumSum of the absolute value
szaxpy, czaxpy, zzaxpy, dzaxpyVector plus the product of a scalar and a vector
snorm2, dnorm2, scnorm2, dznorm2Square root of sum of the squares of the elements of a vector
snrm2, dnrm2, scnrm2, dznrm2Square root of sum of the squares of the elements of a vector
snrsq, dnrsq, scnrsq, dznrsqSum of the squares of the elements of a vector
gen_sortSort the elements of a vector
gen_sortxSort the elements of an indexed vector
isamax, idamax, izamax, icamaxIndex of the element of a vector with maximum absolute value
icamin
isamax, icamax, izamax, idamaxIndex of the element of a vector with maximum absolute value
idamin
ismax, idmaxIndex of the real vector element with maximum value
ismin, idminIndex of the real vector element with minimum value
idamax, icamax, izamax, isamaxIndex of the element of a vector with maximum absolute value
isamin
idmax, ismaxIndex of the real vector element with maximum value
idmin, isminIndex of the real vector element with minimum value
ssortq, dsortq, isortqSort the elements of a vector
ssortqx, dsortqx, isortqxPerforms an indexed sort of a vector
isamax, idamax, icamax, izamaxIndex of the element of a vector with maximum absolute value
izamin
ran16807Routine to generate single precision random numbers using a=16807 and m=2∗∗31-1
ran69069Routine to generate single precision random numbers using a=69069 and m=2∗∗32
ranlRandom number generator based on L’Ecuyer method
ranl_normalRoutine to generate normally distributed random numbers using summation of uniformly distributed random numbers
ranl_skip2Routine to skip forward 2∗∗d seeds for the RANL and RANL_NORMAL random number generators
ranl_skip64Routine to skip forward a given number, d, of seeds for the RANL and RANL_NORMAL random number generators
damax, scamax, dzamax, samaxMaximum absolute value
samin
dasum, scasum, dzasum, sasumSum of the absolute value
daxpy, caxpy, zaxpy, saxpyVector plus the product of a scalar and a vector
samax, damax, dzamax, scamaxMaximum absolute value
scamin
sasum, dasum, dzasum, scasumSum of the absolute value
snorm2, dnorm2, dznorm2, scnorm2Square root of sum of the squares of the elements of a vector
snrm2, dnrm2, dznrm2, scnrm2Square root of sum of the squares of the elements of a vector
snrsq, dnrsq, dznrsq, scnrsqSum of the squares of the elements of a vector
dconv_nonperiodic, cconv_nonperiodic, zconv_nonperiodic, sconv_nonperiodicNonperiodic convolution
dconv_nonperiodic_ext, cconv_nonperiodic_ext, zconv_nonperiodic_ext, sconv_nonperiodic_extExtended nonperiodic convolution
dconv_periodic, cconv_periodic, zconv_periodic, sconv_periodicPeriodic concolution
dconv_periodic_ext, cconv_periodic_ext, zconv_periodic_ext, sconv_periodic_extExtended periodic convolution
dcopy, ccopy, zcopy, scopyCopy of a vector
scorr_nonperiodic
scorr_nonperiodic_ext
scorr_periodic
dcorr_periodic_ext, ccorr_periodic_ext, zcorr_periodic_ext, scorr_periodic_extExtended periodic correlation
ddot, dsdot, cdotc, zdotc, cdotu, zdotu, sdotINNER PRODUCT OF TWO VECTORS
sdsdotProduct of scaled vector and vector
dfct, sfctFast cosine transform in one dimension
dfct_apply, sfct_applyApplication step for fast cosine transform in one dimension
dfct_exit, sfct_exitFinal step for fast cosine transform in one dimension
dfct_init, sfct_initInitialization step for fast cosine transform in one dimension
dfft, cfft, zfft, sfftFast fourier transform in one dimension
dfft_2d, cfft_2d, zfft_2d, sfft_2dFast fourier transform in two dimensions
dfft_3d, cfft_3d, zfft_3d, sfft_3dFast fourier transform in three dimensions
dfft_apply, cfft_apply, zfft_apply, sfft_applyApplication step for fast fourier transform in one dimension
sfft_apply_2d
dfft_apply_3d, cfft_apply_3d, zfft_apply_3d, sfft_apply_3dApplication step for fast fourier transform in three dimensions
dfft_apply_grp, cfft_apply_grp, zfft_apply_grp, sfft_apply_grpApplication step for group fast fourier transform in one dimension
dfft_exit, cfft_exit, zfft_exit, sfft_exitFinal step for fast fourier transform in one dimension
dfft_exit_2d, cfft_exit_2d, zfft_exit_2d, sfft_exit_2dFinal step for fast fourier transform in two dimensions
dfft_exit_3d, cfft_exit_3d, zfft_exit_3d, sfft_exit_3dFinal step for fast fourier transfrom in three dimension
dfft_exit_grp, cfft_exit_grp, zfft_exit_grp, sfft_exit_grpExit step for group fast fourier transform in one dimension
dfft_grp, cfft_grp, zfft_grp, sfft_grpGroup fast fourier transform in one dimension
dfft_init, cfft_init, zfft_init, sfft_initInitialization step for fast fourier transform in one dimension
dfft_init_2d, cfft_init_2d, zfft_init_2d, sfft_init_2dInitialization step for fast fourier transform in two dimensions
dfft_init_3d, cfft_init_3d, zfft_init_3d, sfft_init_3dInitialization step for fast fourier transform in three dimension
dfft_init_grp, cfft_init_grp, zfft_init_grp, sfft_init_grpInitialization step for group fast fourier transform in one dimension
sfilter_apply_nonrecPerforms filtering in lowpass, highpass, bandpass, or bandstop (notch) mode by using the working array that was computed by SFILTER_INIT_NONREC. 
sfilter_init_nonrecComputes a working array that is used by sfilter_apply_nonrec routine. 
sfilter_nonrecPerforms filtering in lowpass, highpass, bandpass, or bandstop (notch) mode. 
dfst, sfstFast sine transform in one dimension
dfst_apply, sfst_applyApplication step for fast sine transform in one dimension
dfst_exit, sfst_exitFinal step for fast sine transform in one dimension
dfst_init, sfst_initInitialization step for fast sine transform in one dimension
ddbmv, cgbmv, zgbmv, sgbmvMatrix-vector product for a general band matrix
dgema, cgema, zgema, sgemaMatrix-matrix addition
dgemm, cgemm, zgemm, sgemmMatrix-matrix product and addition
dgems, cgems, zgems, sgemsMatrix-matrix subtraction
dgemt, cgemt, zgemt, sgemtMatrix-matrix copy
dgemv, cgemv, zgemv, sgemvMatrix-vector product for a general matrix
dger, cgerc, zgerc, cgeru, zgeru, sgerRank-one update of a general matrix
dmax, smaxLargest element in a real vector
dmin, sminMinimum value of the elements of a real vector
dnorm2, scnorm2, dznorm2, snorm2Square root of sum of the squares of the elements of a vector
dnrm2, scnrm2, dznrm2, snrm2Square root of sum of the squares of the elements of a vector
dnrsq, scnrsq, dznrsq, snrsqSum of the squares of the elements of a vector
sortsA library of sort routines
drot, crot, zrot, csrot, zdrot, srotApply givens plane rotation
srotg
drotm, srotmApply modified givens transformation
drotmg, srotmgGenerate elements for a modified Givens transform
dsbmv, chbmv, zhbmv, ssbmvMatrix-vector product for a symmetric or hermitian band matrix
dscal, cscal, zscal, csscal, zdscal, sscalProduct of a scalar and a vector
dset, cset, zset, ssetSet all elements of a vector to a scalar
isortq, dsortq, ssortqSort the elements of a vector
isortqx, dsortqx, ssortqxPerforms an indexed sort of a vector
dspmv, chpmv, zhpmv, sspmvMatrix-vector product for a symmetric or hermitian matrix stored in packed form
dspr, chpr, zhpr, ssprRank-one update of a symmetric or hermitian matrix stored in packed form
dspr2, chpr2, zhpr2, sspr2Rank-two update of a symmetric or hermitian matrix stored in packed form
dsum, csum, zsum, ssumSum of the values of the elements of a vector
dswap, cswap, zswap, sswapExchange the elements of two vectors
dsymm, csymm, zsymm, chemm, zhemm, ssymmMatrix-matrix product and addition for a symmetric or hermitian matrix
dsymv, chemv, zhemv, ssymvMatrix-vector product for a symmetric or hermitian matrix
dsyr, cher, zher, ssyrRank-one update of a symmetric or hermitian matrix
dsyr2, cher2, zher2, ssyr2Rank-two update of a symmetric or hermitian matrix
dsyr2k, csyr2k, zsyr2k, ssyr2kRank-2k update of a symmetric matrix
dsyrk, csyrk, zsyrk, ssyrkRank-k update of a symmetric matrix
dtbmv, ctbmv, ztbmv, stbmvMatrix-vector product for a triangular band matrix
dtbsv, ctbsv, ztbsv, stbsvSolver of a system of linear equations with a triangular band matrix
dtpmv, ctpmv, ztpmv, stpmvMatrix-vector product for a triangular matrix in packed form
dtpsv, ctpsv, ztpsv, stpsvSolve a system of linear equations with a triangular matrix in packed form
dtrmm, ctrmm, ztrmm, strmmMatrix-matrix product for triangular matrix
dtrmv, ctrmv, ztrmv, strmvMarix-vector product for a triangular matrix
dtrsm, ctrsm, ztrsm, strsmSolve a triangular system of equations with a triangular coefficient matrix
dtrsv, ctrsv, ztrsv, strsvSolver of a system of linear equations with a triangular matrix
dvcal, cvcal, zvcal, csvcal, zdvcal, svcalProduct of a scalar and a vector
dzaxpy, czaxpy, zzaxpy, szaxpyVector plus the product of a scalar and a vector
vxworks_dxmlUsing DXML on VxWorks
saxpy, daxpy, caxpy, zaxpyVector plus the product of a scalar and a vector
sconv_nonperiodic, dconv_nonperiodic, cconv_nonperiodic, zconv_nonperiodicNonperiodic convolution
sconv_nonperiodic_ext, dconv_nonperiodic_ext, cconv_nonperiodic_ext, zconv_nonperiodic_extExtended nonperiodic convolution
sconv_periodic, dconv_periodic, cconv_periodic, zconv_periodicPeriodic concolution
sconv_periodic_ext, dconv_periodic_ext, cconv_periodic_ext, zconv_periodic_extExtended periodic convolution
scopy, dcopy, ccopy, zcopyCopy of a vector
zcorr_nonperiodic
zcorr_nonperiodic_ext
zcorr_periodic
scorr_periodic_ext, dcorr_periodic_ext, ccorr_periodic_ext, zcorr_periodic_extExtended periodic correlation
sdot, ddot, dsdot, cdotc, cdotu, zdotu, zdotcINNER PRODUCT OF TWO VECTORS
sdot, ddot, dsdot, cdotc, zdotc, cdotu, zdotuINNER PRODUCT OF TWO VECTORS
srot, drot, crot, zrot, csrot, zdrotApply givens plane rotation
sscal, dscal, cscal, zscal, csscal, zdscalProduct of a scalar and a vector
svcal, dvcal, cvcal, zvcal, csvcal, zdvcalProduct of a scalar and a vector
sfft, dfft, cfft, zfftFast fourier transform in one dimension
sfft_2d, dfft_2d, cfft_2d, zfft_2dFast fourier transform in two dimensions
sfft_3d, dfft_3d, cfft_3d, zfft_3dFast fourier transform in three dimensions
sfft_apply, dfft_apply, cfft_apply, zfft_applyApplication step for fast fourier transform in one dimension
zfft_apply_2d
sfft_apply_3d, dfft_apply_3d, cfft_apply_3d, zfft_apply_3dApplication step for fast fourier transform in three dimensions
sfft_apply_grp, dfft_apply_grp, cfft_apply_grp, zfft_apply_grpApplication step for group fast fourier transform in one dimension
sfft_exit, dfft_exit, cfft_exit, zfft_exitFinal step for fast fourier transform in one dimension
sfft_exit_2d, dfft_exit_2d, cfft_exit_2d, zfft_exit_2dFinal step for fast fourier transform in two dimensions
sfft_exit_3d, dfft_exit_3d, cfft_exit_3d, zfft_exit_3dFinal step for fast fourier transfrom in three dimension
sfft_exit_grp, dfft_exit_grp, cfft_exit_grp, zfft_exit_grpExit step for group fast fourier transform in one dimension
sfft_grp, dfft_grp, cfft_grp, zfft_grpGroup fast fourier transform in one dimension
sfft_init, dfft_init, cfft_init, zfft_initInitialization step for fast fourier transform in one dimension
sfft_init_2d, dfft_init_2d, cfft_init_2d, zfft_init_2dInitialization step for fast fourier transform in two dimensions
sfft_init_3d, dfft_init_3d, cfft_init_3d, zfft_init_3dInitialization step for fast fourier transform in three dimension
sfft_init_grp, dfft_init_grp, cfft_init_grp, zfft_init_grpInitialization step for group fast fourier transform in one dimension
sgbmv, ddbmv, cgbmv, zgbmvMatrix-vector product for a general band matrix
sgema, dgema, cgema, zgemaMatrix-matrix addition
sgemm, dgemm, cgemm, zgemmMatrix-matrix product and addition
sgems, dgems, cgems, zgemsMatrix-matrix subtraction
sgemt, dgemt, cgemt, zgemtMatrix-matrix copy
sgemv, dgemv, cgemv, zgemvMatrix-vector product for a general matrix
sger, dger, cgerc, cgeru, zgeru, zgercRank-one update of a general matrix
sger, dger, cgerc, zgerc, cgeru, zgeruRank-one update of a general matrix
ssbmv, dsbmv, chbmv, zhbmvMatrix-vector product for a symmetric or hermitian band matrix
ssymm, dsymm, csymm, zsymm, chemm, zhemmMatrix-matrix product and addition for a symmetric or hermitian matrix
ssymv, dsymv, chemv, zhemvMatrix-vector product for a symmetric or hermitian matrix
ssyr, dsyr, cher, zherRank-one update of a symmetric or hermitian matrix
ssyr2, dsyr2, cher2, zher2Rank-two update of a symmetric or hermitian matrix
cher2k, zher2kRank-2k update of a complex hermitian matrix
cherk, zherkRank-k update of a complex hermitian matrix
sspmv, dspmv, chpmv, zhpmvMatrix-vector product for a symmetric or hermitian matrix stored in packed form
sspr, dspr, chpr, zhprRank-one update of a symmetric or hermitian matrix stored in packed form
sspr2, dspr2, chpr2, zhpr2Rank-two update of a symmetric or hermitian matrix stored in packed form
srot, drot, crot, csrot, zdrot, zrotApply givens plane rotation
zrotg
sscal, dscal, cscal, csscal, zdscal, zscalProduct of a scalar and a vector
sset, dset, cset, zsetSet all elements of a vector to a scalar
ssum, dsum, csum, zsumSum of the values of the elements of a vector
sswap, dswap, cswap, zswapExchange the elements of two vectors
ssymm, dsymm, csymm, chemm, zhemm, zsymmMatrix-matrix product and addition for a symmetric or hermitian matrix
ssyr2k, dsyr2k, csyr2k, zsyr2kRank-2k update of a symmetric matrix
ssyrk, dsyrk, csyrk, zsyrkRank-k update of a symmetric matrix
stbmv, dtbmv, ctbmv, ztbmvMatrix-vector product for a triangular band matrix
stbsv, dtbsv, ctbsv, ztbsvSolver of a system of linear equations with a triangular band matrix
stpmv, dtpmv, ctpmv, ztpmvMatrix-vector product for a triangular matrix in packed form
stpsv, dtpsv, ctpsv, ztpsvSolve a system of linear equations with a triangular matrix in packed form
strmm, dtrmm, ctrmm, ztrmmMatrix-matrix product for triangular matrix
strmv, dtrmv, ctrmv, ztrmvMarix-vector product for a triangular matrix
strsm, dtrsm, ctrsm, ztrsmSolve a triangular system of equations with a triangular coefficient matrix
strsv, dtrsv, ctrsv, ztrsvSolver of a system of linear equations with a triangular matrix
svcal, dvcal, cvcal, csvcal, zdvcal, zvcalProduct of a scalar and a vector
szaxpy, dzaxpy, czaxpy, zzaxpyVector plus the product of a scalar and a vector

Section 3lapack

 , csrot
 , zdrot
 , zsrot

Section l

CBDSQR, cbdsqrcompute the singular value decomposition (SVD) of a real N-by-N (upper or lower) bidiagonal matrix B
CGBBRD, cgbbrdreduce a complex general m-by-n band matrix A to real upper bidiagonal form B by a unitary transformation
CGBCON, cgbconestimate the reciprocal of the condition number of a complex general band matrix A, in either the 1-norm or the infinity-norm,
CGBEQU, cgbequcompute row and column scalings intended to equilibrate an M-by-N band matrix A and reduce its condition number
CGBRFS, cgbrfsimprove 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
CGBSV, cgbsvcompute the solution to a complex system of linear equations A ∗ X = B, where A is a band matrix of order N with KL subdiagonals and KU superdiagonals, and X and B are N-by-NRHS matrices
CGBSVX, cgbsvxuse the LU factorization to compute the solution to a complex system of linear equations A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B,
CGBTF2, cgbtf2compute an LU factorization of a complex m-by-n band matrix A using partial pivoting with row interchanges
CGBTRF, cgbtrfcompute an LU factorization of a complex m-by-n band matrix A using partial pivoting with row interchanges
CGBTRS, cgbtrssolve a system of linear equations  A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B with a general band matrix A using the LU factorization computed by CGBTRF
CGEBAK, cgebakform the right or left eigenvectors of a complex general matrix by backward transformation on the computed eigenvectors of the balanced matrix output by CGEBAL
CGEBAL, cgebalbalance a general complex matrix A
CGEBD2, cgebd2reduce a complex general m by n matrix A to upper or lower real bidiagonal form B by a unitary transformation
CGEBRD, cgebrdreduce a general complex M-by-N matrix A to upper or lower bidiagonal form B by a unitary transformation
CGECON, cgeconestimate the reciprocal of the condition number of a general complex matrix A, in either the 1-norm or the infinity-norm, using the LU factorization computed by CGETRF
CGEEQU, cgeequcompute row and column scalings intended to equilibrate an M-by-N matrix A and reduce its condition number
CGEES, cgeescompute for an N-by-N complex nonsymmetric matrix A, the eigenvalues, the Schur form T, and, optionally, the matrix of Schur vectors Z
CGEESX, cgeesxcompute for an N-by-N complex nonsymmetric matrix A, the eigenvalues, the Schur form T, and, optionally, the matrix of Schur vectors Z
CGEEV, cgeevcompute for an N-by-N complex nonsymmetric matrix A, the eigenvalues and, optionally, the left and/or right eigenvectors
CGEEVX, cgeevxcompute for an N-by-N complex nonsymmetric matrix A, the eigenvalues and, optionally, the left and/or right eigenvectors
CGEGS, cgegscompute for a pair of N-by-N complex nonsymmetric matrices A,
CGEGV, cgegvcompute for a pair of N-by-N complex nonsymmetric matrices A and B, the generalized eigenvalues (alpha, beta), and optionally,
CGEHD2, cgehd2reduce a complex general matrix A to upper Hessenberg form H by a unitary similarity transformation
CGEHRD, cgehrdreduce a complex general matrix A to upper Hessenberg form H by a unitary similarity transformation
CGELQ2, cgelq2compute an LQ factorization of a complex m by n matrix A
CGELQF, cgelqfcompute an LQ factorization of a complex M-by-N matrix A
CGELS, cgelssolve overdetermined or underdetermined complex linear systems involving an M-by-N matrix A, or its conjugate-transpose, using a QR or LQ factorization of A
CGELSS, cgelsscompute the minimum norm solution to a complex linear least squares problem
CGELSX, cgelsxcompute the minimum-norm solution to a complex linear least squares problem
CGEQL2, cgeql2compute a QL factorization of a complex m by n matrix A
CGEQLF, cgeqlfcompute a QL factorization of a complex M-by-N matrix A
CGEQPF, cgeqpfcompute a QR factorization with column pivoting of a complex M-by-N matrix A
CGEQR2, cgeqr2compute a QR factorization of a complex m by n matrix A
CGEQRF, cgeqrfcompute a QR factorization of a complex M-by-N matrix A
CGERFS, cgerfsimprove the computed solution to a system of linear equations and provides error bounds and backward error estimates for the solution
CGERQ2, cgerq2compute an RQ factorization of a complex m by n matrix A
CGERQF, cgerqfcompute an RQ factorization of a complex M-by-N matrix A
CGESV, cgesvcompute the solution to a complex system of linear equations  A ∗ X = B,
CGESVD, cgesvdcompute the singular value decomposition (SVD) of a complex M-by-N matrix A, optionally computing the left and/or right singular vectors
CGESVX, cgesvxuse the LU factorization to compute the solution to a complex system of linear equations  A ∗ X = B,
CGETF2, cgetf2compute an LU factorization of a general m-by-n matrix A using partial pivoting with row interchanges
CGETRF, cgetrfcompute an LU factorization of a general M-by-N matrix A using partial pivoting with row interchanges
CGETRI, cgetricompute the inverse of a matrix using the LU factorization computed by CGETRF
CGETRS, cgetrssolve a system of linear equations  A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B with a general N-by-N matrix A using the LU factorization computed by CGETRF
CGGBAK, cggbakform the right or left eigenvectors of a complex generalized eigenvalue problem A∗x = lambda∗B∗x, by backward transformation on the computed eigenvectors of the balanced pair of matrices output by CGGBAL
CGGBAL, cggbalbalance a pair of general complex matrices (A,B)
CGGGLM, cggglmsolve a general Gauss-Markov linear model (GLM) problem
CGGHRD, cgghrdreduce a pair of complex matrices (A,B) to generalized upper Hessenberg form using unitary transformations, where A is a general matrix and B is upper triangular
CGGLSE, cgglsesolve the linear equality-constrained least squares (LSE) problem
CGGQRF, cggqrfcompute a generalized QR factorization of an N-by-M matrix A and an N-by-P matrix B
CGGRQF, cggrqfcompute a generalized RQ factorization of an M-by-N matrix A and a P-by-N matrix B
CGGSVD, cggsvdcompute the generalized singular value decomposition (GSVD) of an M-by-N complex matrix A and P-by-N complex matrix B
CGGSVP, cggsvpcompute unitary matrices U, V and Q such that   N-K-L K L  U’∗A∗Q = K ( 0 A12 A13 ) if M-K-L >= 0
CGTCON, cgtconestimate the reciprocal of the condition number of a complex tridiagonal matrix A using the LU factorization as computed by CGTTRF
CGTRFS, cgtrfsimprove the computed solution to a system of linear equations when the coefficient matrix is tridiagonal, and provides error bounds and backward error estimates for the solution
CGTSV, cgtsvsolve the equation   A∗X = B,
CGTSVX, cgtsvxuse the LU factorization to compute the solution to a complex system of linear equations A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B,
CGTTRF, cgttrfcompute an LU factorization of a complex tridiagonal matrix A using elimination with partial pivoting and row interchanges
CGTTRS, cgttrssolve one of the systems of equations  A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B,
CHBEV, chbevcompute all the eigenvalues and, optionally, eigenvectors of a complex Hermitian band matrix A
CHBEVD, chbevdcompute all the eigenvalues and, optionally, eigenvectors of a complex Hermitian band matrix A
CHBEVX, chbevxcompute selected eigenvalues and, optionally, eigenvectors of a complex Hermitian band matrix A
CHBGST, chbgstreduce a complex Hermitian-definite banded generalized eigenproblem A∗x = lambda∗B∗x to standard form C∗y = lambda∗y,
CHBGV, chbgvcompute all the eigenvalues, and optionally, the eigenvectors of a complex generalized Hermitian-definite banded eigenproblem, of the form A∗x=(lambda)∗B∗x
CHBTRD, chbtrdreduce a complex Hermitian band matrix A to real symmetric tridiagonal form T by a unitary similarity transformation
CHECON, checonestimate the reciprocal of the condition number of a complex Hermitian matrix A using the factorization A = U∗D∗U∗∗H or A = L∗D∗L∗∗H computed by CHETRF
CHEEV, cheevcompute all eigenvalues and, optionally, eigenvectors of a complex Hermitian matrix A
CHEEVD, cheevdcompute all eigenvalues and, optionally, eigenvectors of a complex Hermitian matrix A
CHEEVX, cheevxcompute selected eigenvalues and, optionally, eigenvectors of a complex Hermitian matrix A
CHEGS2, chegs2reduce a complex Hermitian-definite generalized eigenproblem to standard form
CHEGST, chegstreduce a complex Hermitian-definite generalized eigenproblem to standard form
CHEGV, chegvcompute all the eigenvalues, and optionally, the eigenvectors of a complex generalized Hermitian-definite eigenproblem, of the form A∗x=(lambda)∗B∗x, A∗Bx=(lambda)∗x, or B∗A∗x=(lambda)∗x
CHERFS, cherfsimprove the computed solution to a system of linear equations when the coefficient matrix is Hermitian indefinite, and provides error bounds and backward error estimates for the solution
CHESV, chesvcompute the solution to a complex system of linear equations  A ∗ X = B,
CHESVX, chesvxuse the diagonal pivoting factorization to compute the solution to a complex system of linear equations A ∗ X = B,
CHETD2, chetd2reduce a complex Hermitian matrix A to real symmetric tridiagonal form T by a unitary similarity transformation
CHETF2, chetf2compute the factorization of a complex Hermitian matrix A using the Bunch-Kaufman diagonal pivoting method
CHETRD, chetrdreduce a complex Hermitian matrix A to real symmetric tridiagonal form T by a unitary similarity transformation
CHETRF, chetrfcompute the factorization of a complex Hermitian matrix A using the Bunch-Kaufman diagonal pivoting method
CHETRI, chetricompute the inverse of a complex Hermitian indefinite matrix A using the factorization A = U∗D∗U∗∗H or A = L∗D∗L∗∗H computed by CHETRF
CHETRS, chetrssolve a system of linear equations A∗X = B with a complex Hermitian matrix A using the factorization A = U∗D∗U∗∗H or A = L∗D∗L∗∗H computed by CHETRF
CHGEQZ, chgeqzimplement a single-shift version of the QZ method for finding the generalized eigenvalues w(i)=ALPHA(i)/BETA(i) of the equation   det( A - w(i) B ) = 0  If JOB=’S’, then the pair (A,B) is simultaneously reduced to Schur form (i.e., A and B are both upper triangular) by applying one unitary tranformation (usually called Q) on the left and another (usually called Z) on the right
CHPCON, chpconestimate the reciprocal of the condition number of a complex Hermitian packed matrix A using the factorization A = U∗D∗U∗∗H or A = L∗D∗L∗∗H computed by CHPTRF
CHPEV, chpevcompute all the eigenvalues and, optionally, eigenvectors of a complex Hermitian matrix in packed storage
CHPEVD, chpevdcompute all the eigenvalues and, optionally, eigenvectors of a complex Hermitian matrix A in packed storage
CHPEVX, chpevxcompute selected eigenvalues and, optionally, eigenvectors of a complex Hermitian matrix A in packed storage
CHPGST, chpgstreduce a complex Hermitian-definite generalized eigenproblem to standard form, using packed storage
CHPGV, chpgvcompute all the eigenvalues and, optionally, the eigenvectors of a complex generalized Hermitian-definite eigenproblem, of the form A∗x=(lambda)∗B∗x, A∗Bx=(lambda)∗x, or B∗A∗x=(lambda)∗x
CHPRFS, chprfsimprove the computed solution to a system of linear equations when the coefficient matrix is Hermitian indefinite and packed, and provides error bounds and backward error estimates for the solution
CHPSV, chpsvcompute the solution to a complex system of linear equations  A ∗ X = B,
CHPSVX, chpsvxuse the diagonal pivoting factorization A = U∗D∗U∗∗H or A = L∗D∗L∗∗H to compute the solution to a complex system of linear equations A ∗ X = B, where A is an N-by-N Hermitian matrix stored in packed format and X and B are N-by-NRHS matrices
CHPTRD, chptrdreduce a complex Hermitian matrix A stored in packed form to real symmetric tridiagonal form T by a unitary similarity transformation
CHPTRF, chptrfcompute the factorization of a complex Hermitian packed matrix A using the Bunch-Kaufman diagonal pivoting method
CHPTRI, chptricompute the inverse of a complex Hermitian indefinite matrix A in packed storage using the factorization A = U∗D∗U∗∗H or A = L∗D∗L∗∗H computed by CHPTRF
CHPTRS, chptrssolve a system of linear equations A∗X = B with a complex Hermitian matrix A stored in packed format using the factorization A = U∗D∗U∗∗H or A = L∗D∗L∗∗H computed by CHPTRF
CHSEIN, chseinuse inverse iteration to find specified right and/or left eigenvectors of a complex upper Hessenberg matrix H
CHSEQR, chseqrcompute the eigenvalues of a complex upper Hessenberg matrix H, and, optionally, the matrices T and Z from the Schur decomposition H = Z T Z∗∗H, where T is an upper triangular matrix (the Schur form), and Z is the unitary matrix of Schur vectors
CLABRD, clabrdreduce the first NB rows and columns of a complex general m by n matrix A to upper or lower real bidiagonal form by a unitary transformation Q’ ∗ A ∗ P, and returns the matrices X and Y which are needed to apply the transformation to the unreduced part of A
CLACGV, clacgvconjugate a complex vector of length N
CLACON, claconestimate the 1-norm of a square, complex matrix A
CLACPY, clacpycopie all or part of a two-dimensional matrix A to another matrix B
CLACRM, clacrmperform a very simple matrix-matrix multiplication
CLACRT, clacrtapplie a plane rotation, where the cos and sin (C and S) are complex and the vectors CX and CY are complex
CLADIV, cladiv:= X / Y, where X and Y are complex
CLAED0, claed0the divide and conquer method, CLAED0 computes all eigenvalues of a symmetric tridiagonal matrix which is one diagonal block of those from reducing a dense or band Hermitian matrix and corresponding eigenvectors of the dense or band matrix
CLAED7, claed7compute the updated eigensystem of a diagonal matrix after modification by a rank-one symmetric matrix
CLAED8, claed8merge the two sets of eigenvalues together into a single sorted set
CLAEIN, claeinuse inverse iteration to find a right or left eigenvector corresponding to the eigenvalue W of a complex upper Hessenberg matrix H
CLAESY, claesycompute the eigendecomposition of a 2-by-2 symmetric matrix  ( ( A, B );( B, C ) ) provided the norm of the matrix of eigenvectors is larger than some threshold value
CLAEV2, claev2compute the eigendecomposition of a 2-by-2 Hermitian matrix  [ A B ]  [ CONJG(B) C ]
clags2
CLAGTM, clagtmperform a matrix-vector product of the form   B := alpha ∗ A ∗ X + beta ∗ B  where A is a tridiagonal matrix of order N, B and X are N by NRHS matrices, and alpha and beta are real scalars, each of which may be 0., 1., or -1
CLAHEF, clahefcompute a partial factorization of a complex Hermitian matrix A using the Bunch-Kaufman diagonal pivoting method
CLAHQR, clahqri an auxiliary routine called by CHSEQR to update the eigenvalues and Schur decomposition already computed by CHSEQR, by dealing with the Hessenberg submatrix in rows and columns ILO to IHI
CLAHRD, clahrdreduce the first NB columns of a complex general n-by-(n-k+1) matrix A so that elements below the k-th subdiagonal are zero
CLAIC1, claic1applie one step of incremental condition estimation in its simplest version
CLANGB, clangbreturn 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
CLANGE, clangereturn the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex matrix A
CLANGT, clangtreturn the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex tridiagonal matrix A
CLANHB, clanhbreturn 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 hermitian band matrix A, with k super-diagonals
CLANHE, clanhereturn the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex hermitian matrix A
CLANHP, clanhpreturn the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex hermitian matrix A, supplied in packed form
CLANHS, clanhsreturn 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
CLANHT, clanhtreturn the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex Hermitian tridiagonal matrix A
CLANSB, clansbreturn 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
CLANSP, clanspreturn the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex symmetric matrix A, supplied in packed form
CLANSY, clansyreturn the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex symmetric matrix A
CLANTB, clantbreturn 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
CLANTP, clantpreturn 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
CLANTR, clantrreturn 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
CLAPLL, claplltwo column vectors X and Y, let   A = ( X Y )
CLAPMT, clapmtrearrange 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
CLAQGB, claqgbequilibrate 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
CLAQGE, claqgeequilibrate a general M by N matrix A using the row and scaling factors in the vectors R and C
CLAQHB, claqhbequilibrate a symmetric band matrix A using the scaling factors in the vector S
CLAQHE, claqheequilibrate a Hermitian matrix A using the scaling factors in the vector S
CLAQHP, claqhpequilibrate a Hermitian matrix A using the scaling factors in the vector S
CLAQSB, claqsbequilibrate a symmetric band matrix A using the scaling factors in the vector S
CLAQSP, claqspequilibrate a symmetric matrix A using the scaling factors in the vector S
CLAQSY, claqsyequilibrate a symmetric matrix A using the scaling factors in the vector S
CLAR2V, clar2vapplie a vector of complex plane rotations with real cosines from both sides to a sequence of 2-by-2 complex Hermitian matrices,
CLARF, clarfapplie a complex elementary reflector H to a complex M-by-N matrix C, from either the left or the right
CLARFB, clarfbapplie a complex block reflector H or its transpose H’ to a complex M-by-N matrix C, from either the left or the right
CLARFG, clarfggenerate a complex elementary reflector H of order n, such that   H’ ∗ ( alpha ) = ( beta ), H’ ∗ H = I
CLARFT, clarftform the triangular factor T of a complex block reflector H of order n, which is defined as a product of k elementary reflectors
CLARFX, clarfxapplie a complex elementary reflector H to a complex m by n matrix C, from either the left or the right
CLARGV, clargvgenerate a vector of complex plane rotations with real cosines, determined by elements of the complex vectors x and y
CLARNV, clarnvreturn a vector of n random complex numbers from a uniform or normal distribution
CLARTG, clartggenerate a plane rotation so that   [ CS SN ] [ F ] [ R ]  [ __ ]
CLARTV, clartvapplie a vector of complex plane rotations with real cosines to elements of the complex vectors x and y
CLASCL, clasclmultiplie the M by N complex matrix A by the real scalar CTO/CFROM
CLASET, clasetinitialize a 2-D array A to BETA on the diagonal and ALPHA on the offdiagonals
CLASR, clasrperform the transformation   A := P∗A, when SIDE = ’L’ or ’l’ ( Left-hand side )   A := A∗P’, when SIDE = ’R’ or ’r’ ( Right-hand side )  where A is an m by n complex matrix and P is an orthogonal matrix,
CLASSQ, classqreturn the values scl and ssq such that   ( scl∗∗2 )∗ssq = x( 1 )∗∗2 +...+ x( n )∗∗2 + ( scale∗∗2 )∗sumsq,
CLASWP, claswpperform a series of row interchanges on the matrix A
CLASYF, clasyfcompute a partial factorization of a complex symmetric matrix A using the Bunch-Kaufman diagonal pivoting method
CLATBS, clatbssolve one of the triangular systems   A ∗ x = s∗b, A∗∗T ∗ x = s∗b, or A∗∗H ∗ x = s∗b,
CLATPS, clatpssolve one of the triangular systems   A ∗ x = s∗b, A∗∗T ∗ x = s∗b, or A∗∗H ∗ x = s∗b,
CLATRD, clatrdreduce NB rows and columns of a complex Hermitian matrix A to Hermitian tridiagonal form by a unitary similarity transformation Q’ ∗ A ∗ Q, and returns the matrices V and W which are needed to apply the transformation to the unreduced part of A
CLATRS, clatrssolve one of the triangular systems   A ∗ x = s∗b, A∗∗T ∗ x = s∗b, or A∗∗H ∗ x = s∗b,
CLATZM, clatzmapplie a Householder matrix generated by CTZRQF to a matrix
CLAUU2, clauu2compute the product U ∗ U’ or L’ ∗ L, where the triangular factor U or L is stored in the upper or lower triangular part of the array A
CLAUUM, clauumcompute the product U ∗ U’ or L’ ∗ L, where the triangular factor U or L is stored in the upper or lower triangular part of the array A
CLAZRO, clazroinitialize a 2-D array A to BETA on the diagonal and ALPHA on the offdiagonals
CPBCON, cpbconestimate the reciprocal of the condition number (in the 1-norm) of a complex Hermitian positive definite band matrix using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by CPBTRF
CPBEQU, cpbequcompute row and column scalings intended to equilibrate a Hermitian positive definite band matrix A and reduce its condition number (with respect to the two-norm)
CPBRFS, cpbrfsimprove the computed solution to a system of linear equations when the coefficient matrix is Hermitian positive definite and banded, and provides error bounds and backward error estimates for the solution
CPBSTF, cpbstfcompute a split Cholesky factorization of a complex Hermitian positive definite band matrix A
CPBSV, cpbsvcompute the solution to a complex system of linear equations  A ∗ X = B,
CPBSVX, cpbsvxuse the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H to compute the solution to a complex system of linear equations  A ∗ X = B,
CPBTF2, cpbtf2compute the Cholesky factorization of a complex Hermitian positive definite band matrix A
CPBTRF, cpbtrfcompute the Cholesky factorization of a complex Hermitian positive definite band matrix A
CPBTRS, cpbtrssolve a system of linear equations A∗X = B with a Hermitian positive definite band matrix A using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by CPBTRF
CPOCON, cpoconestimate the reciprocal of the condition number (in the 1-norm) of a complex Hermitian positive definite matrix using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by CPOTRF
CPOEQU, cpoequcompute row and column scalings intended to equilibrate a Hermitian positive definite matrix A and reduce its condition number (with respect to the two-norm)
CPORFS, cporfsimprove the computed solution to a system of linear equations when the coefficient matrix is Hermitian positive definite,
CPOSV, cposvcompute the solution to a complex system of linear equations  A ∗ X = B,
CPOSVX, cposvxuse the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H to compute the solution to a complex system of linear equations  A ∗ X = B,
CPOTF2, cpotf2compute the Cholesky factorization of a complex Hermitian positive definite matrix A
CPOTRF, cpotrfcompute the Cholesky factorization of a complex Hermitian positive definite matrix A
CPOTRI, cpotricompute the inverse of a complex Hermitian positive definite matrix A using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by CPOTRF
CPOTRS, cpotrssolve a system of linear equations A∗X = B with a Hermitian positive definite matrix A using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by CPOTRF
CPPCON, cppconestimate the reciprocal of the condition number (in the 1-norm) of a complex Hermitian positive definite packed matrix using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by CPPTRF
CPPEQU, cppequcompute row and column scalings intended to equilibrate a Hermitian positive definite matrix A in packed storage and reduce its condition number (with respect to the two-norm)
CPPRFS, cpprfsimprove the computed solution to a system of linear equations when the coefficient matrix is Hermitian positive definite and packed, and provides error bounds and backward error estimates for the solution
CPPSV, cppsvcompute the solution to a complex system of linear equations  A ∗ X = B,
CPPSVX, cppsvxuse the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H to compute the solution to a complex system of linear equations  A ∗ X = B,
CPPTRF, cpptrfcompute the Cholesky factorization of a complex Hermitian positive definite matrix A stored in packed format
CPPTRI, cpptricompute the inverse of a complex Hermitian positive definite matrix A using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by CPPTRF
CPPTRS, cpptrssolve a system of linear equations A∗X = B with a Hermitian positive definite matrix A in packed storage using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by CPPTRF
CPTCON, cptconcompute the reciprocal of the condition number (in the 1-norm) of a complex Hermitian positive definite tridiagonal matrix using the factorization A = L∗D∗L∗∗H or A = U∗∗H∗D∗U computed by CPTTRF
CPTEQR, cpteqrcompute all eigenvalues and, optionally, eigenvectors of a symmetric positive definite tridiagonal matrix by first factoring the matrix using SPTTRF and then calling CBDSQR to compute the singular values of the bidiagonal factor
CPTRFS, cptrfsimprove the computed solution to a system of linear equations when the coefficient matrix is Hermitian positive definite and tridiagonal, and provides error bounds and backward error estimates for the solution
CPTSV, cptsvcompute the solution to a complex system of linear equations A∗X = B, where A is an N-by-N Hermitian positive definite tridiagonal matrix, and X and B are N-by-NRHS matrices
CPTSVX, cptsvxuse the factorization A = L∗D∗L∗∗H to compute the solution to a complex system of linear equations A∗X = B, where A is an N-by-N Hermitian positive definite tridiagonal matrix and X and B are N-by-NRHS matrices
CPTTRF, cpttrfcompute the factorization of a complex Hermitian positive definite tridiagonal matrix A
CPTTRS, cpttrssolve a system of linear equations A ∗ X = B with a Hermitian positive definite tridiagonal matrix A using the factorization A = U∗∗H∗D∗U or A = L∗D∗L∗∗H computed by CPTTRF
CROT, crotapplie a plane rotation, where the cos (C) is real and the sin (S) is complex, and the vectors CX and CY are complex
CSPCON, cspconestimate the reciprocal of the condition number (in the 1-norm) of a complex symmetric packed matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by CSPTRF
CSPMV, cspmvperform the matrix-vector operation   y := alpha∗A∗x + beta∗y,
CSPR, csprperform the symmetric rank 1 operation   A := alpha∗x∗conjg( x’ ) + A,
CSPRFS, csprfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric indefinite and packed, and provides error bounds and backward error estimates for the solution
CSPSV, cspsvcompute the solution to a complex system of linear equations  A ∗ X = B,
CSPSVX, cspsvxuse the diagonal pivoting factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T to compute the solution to a complex system of linear equations A ∗ X = B, where A is an N-by-N symmetric matrix stored in packed format and X and B are N-by-NRHS matrices
CSPTRF, csptrfcompute the factorization of a complex symmetric matrix A stored in packed format using the Bunch-Kaufman diagonal pivoting method
CSPTRI, csptricompute the inverse of a complex symmetric indefinite matrix A in packed storage using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by CSPTRF
CSPTRS, csptrssolve a system of linear equations A∗X = B with a complex symmetric matrix A stored in packed format using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by CSPTRF
CSRSCL, csrsclmultiplie an n-element complex vector x by the real scalar 1/a
CSTEDC, cstedccompute all eigenvalues and, optionally, eigenvectors of a symmetric tridiagonal matrix using the divide and conquer method
CSTEIN, csteincompute the eigenvectors of a real symmetric tridiagonal matrix T corresponding to specified eigenvalues, using inverse iteration
CSTEQR, csteqrcompute all eigenvalues and, optionally, eigenvectors of a symmetric tridiagonal matrix using the implicit QL or QR method
CSYCON, csyconestimate the reciprocal of the condition number (in the 1-norm) of a complex symmetric matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by CSYTRF
CSYMV, csymvperform the matrix-vector operation   y := alpha∗A∗x + beta∗y,
CSYR, csyrperform the symmetric rank 1 operation   A := alpha∗x∗( x’ ) + A,
CSYRFS, csyrfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric indefinite, and provides error bounds and backward error estimates for the solution
CSYSV, csysvcompute the solution to a complex system of linear equations  A ∗ X = B,
CSYSVX, csysvxuse the diagonal pivoting factorization to compute the solution to a complex system of linear equations A ∗ X = B,
CSYTF2, csytf2compute the factorization of a complex symmetric matrix A using the Bunch-Kaufman diagonal pivoting method
CSYTRF, csytrfcompute the factorization of a complex symmetric matrix A using the Bunch-Kaufman diagonal pivoting method
CSYTRI, csytricompute the inverse of a complex symmetric indefinite matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by CSYTRF
CSYTRS, csytrssolve a system of linear equations A∗X = B with a complex symmetric matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by CSYTRF
CTBCON, ctbconestimate the reciprocal of the condition number of a triangular band matrix A, in either the 1-norm or the infinity-norm
CTBRFS, ctbrfsprovide error bounds and backward error estimates for the solution to a system of linear equations with a triangular band coefficient matrix
CTBTRS, ctbtrssolve a triangular system of the form   A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B,
CTGEVC, ctgevccompute some or all of the right and/or left generalized eigenvectors of a pair of complex upper triangular matrices (A,B)
CTGSJA, ctgsjacompute the generalized singular value decomposition (GSVD) of two complex upper triangular (or trapezoidal) matrices A and B
CTPCON, ctpconestimate the reciprocal of the condition number of a packed triangular matrix A, in either the 1-norm or the infinity-norm
CTPRFS, ctprfsprovide error bounds and backward error estimates for the solution to a system of linear equations with a triangular packed coefficient matrix
CTPTRI, ctptricompute the inverse of a complex upper or lower triangular matrix A stored in packed format
CTPTRS, ctptrssolve a triangular system of the form   A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B,
CTRCON, ctrconestimate the reciprocal of the condition number of a triangular matrix A, in either the 1-norm or the infinity-norm
CTREVC, ctrevccompute some or all of the right and/or left eigenvectors of a complex upper triangular matrix T
CTREXC, ctrexcreorder the Schur factorization of a complex matrix A = Q∗T∗Q∗∗H, so that the diagonal element of T with row index IFST is moved to row ILST
CTRRFS, ctrrfsprovide error bounds and backward error estimates for the solution to a system of linear equations with a triangular coefficient matrix
CTRSEN, ctrsenreorder the Schur factorization of a complex matrix A = Q∗T∗Q∗∗H, so that a selected cluster of eigenvalues appears in the leading positions on the diagonal of the upper triangular matrix T, and the leading columns of Q form an orthonormal basis of the corresponding right invariant subspace
CTRSNA, ctrsnaestimate reciprocal condition numbers for specified eigenvalues and/or right eigenvectors of a complex upper triangular matrix T (or of any matrix Q∗T∗Q∗∗H with Q unitary)
CTRSYL, ctrsylsolve the complex Sylvester matrix equation
CTRTI2, ctrti2compute the inverse of a complex upper or lower triangular matrix
CTRTRI, ctrtricompute the inverse of a complex upper or lower triangular matrix A
CTRTRS, ctrtrssolve a triangular system of the form   A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B,
CTZRQF, ctzrqfreduce the M-by-N ( M<=N ) complex upper trapezoidal matrix A to upper triangular form by means of unitary transformations
CUNG2L, cung2lgenerate an m by n complex matrix Q with orthonormal columns,
CUNG2R, cung2rgenerate an m by n complex matrix Q with orthonormal columns,
CUNGBR, cungbrgenerate one of the complex unitary matrices Q or P∗∗H determined by CGEBRD when reducing a complex matrix A to bidiagonal form
CUNGHR, cunghrgenerate a complex unitary matrix Q which is defined as the product of IHI-ILO elementary reflectors of order N, as returned by CGEHRD
CUNGL2, cungl2generate an m-by-n complex matrix Q with orthonormal rows,
CUNGLQ, cunglqgenerate an M-by-N complex matrix Q with orthonormal rows,
CUNGQL, cungqlgenerate an M-by-N complex matrix Q with orthonormal columns,
CUNGQR, cungqrgenerate an M-by-N complex matrix Q with orthonormal columns,
CUNGR2, cungr2generate an m by n complex matrix Q with orthonormal rows,
CUNGRQ, cungrqgenerate an M-by-N complex matrix Q with orthonormal rows,
CUNGTR, cungtrgenerate a complex unitary matrix Q which is defined as the product of n-1 elementary reflectors of order N, as returned by CHETRD
CUNM2L, cunm2loverwrite the general complex m-by-n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’C’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’C’,
CUNM2R, cunm2roverwrite the general complex m-by-n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’C’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’C’,
CUNMBR, cunmbrVECT = ’Q’, CUNMBR overwrites the general complex M-by-N matrix C with  SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
CUNMHR, cunmhroverwrite the general complex M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
CUNML2, cunml2overwrite the general complex m-by-n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’C’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’C’,
CUNMLQ, cunmlqoverwrite the general complex M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
CUNMQL, cunmqloverwrite the general complex M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
CUNMQR, cunmqroverwrite the general complex M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
CUNMR2, cunmr2overwrite the general complex m-by-n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’C’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’C’,
CUNMRQ, cunmrqoverwrite the general complex M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
CUNMTR, cunmtroverwrite the general complex M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
CUPGTR, cupgtrgenerate a complex unitary matrix Q which is defined as the product of n-1 elementary reflectors H(i) of order n, as returned by CHPTRD using packed storage
CUPMTR, cupmtroverwrite the general complex M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
DBDSQR, dbdsqrcompute the singular value decomposition (SVD) of a real N-by-N (upper or lower) bidiagonal matrix B
DDISNA, ddisnacompute 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
DGBBRD, dgbbrdreduce a real general m-by-n band matrix A to upper bidiagonal form B by an orthogonal transformation
DGBCON, dgbconestimate the reciprocal of the condition number of a real general band matrix A, in either the 1-norm or the infinity-norm,
DGBEQU, dgbequcompute row and column scalings intended to equilibrate an M-by-N band matrix A and reduce its condition number
DGBRFS, dgbrfsimprove 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
DGBSV, dgbsvcompute the solution to a real system of linear equations A ∗ X = B, where A is a band matrix of order N with KL subdiagonals and KU superdiagonals, and X and B are N-by-NRHS matrices
DGBSVX, dgbsvxuse the LU factorization to compute the solution to a real system of linear equations A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B,
DGBTF2, dgbtf2compute an LU factorization of a real m-by-n band matrix A using partial pivoting with row interchanges
DGBTRF, dgbtrfcompute an LU factorization of a real m-by-n band matrix A using partial pivoting with row interchanges
DGBTRS, dgbtrssolve a system of linear equations  A ∗ X = B or A’ ∗ X = B with a general band matrix A using the LU factorization computed by DGBTRF
DGEBAK, dgebakform the right or left eigenvectors of a real general matrix by backward transformation on the computed eigenvectors of the balanced matrix output by DGEBAL
DGEBAL, dgebalbalance a general real matrix A
DGEBD2, dgebd2reduce a real general m by n matrix A to upper or lower bidiagonal form B by an orthogonal transformation
DGEBRD, dgebrdreduce a general real M-by-N matrix A to upper or lower bidiagonal form B by an orthogonal transformation
DGECON, dgeconestimate the reciprocal of the condition number of a general real matrix A, in either the 1-norm or the infinity-norm, using the LU factorization computed by DGETRF
DGEEQU, dgeequcompute row and column scalings intended to equilibrate an M-by-N matrix A and reduce its condition number
DGEES, dgeescompute for an N-by-N real nonsymmetric matrix A, the eigenvalues, the real Schur form T, and, optionally, the matrix of Schur vectors Z
DGEESX, dgeesxcompute for an N-by-N real nonsymmetric matrix A, the eigenvalues, the real Schur form T, and, optionally, the matrix of Schur vectors Z
DGEEV, dgeevcompute for an N-by-N real nonsymmetric matrix A, the eigenvalues and, optionally, the left and/or right eigenvectors
DGEEVX, dgeevxcompute for an N-by-N real nonsymmetric matrix A, the eigenvalues and, optionally, the left and/or right eigenvectors
DGEGS, dgegscompute for a pair of N-by-N real nonsymmetric matrices A, B
DGEGV, dgegvcompute for a pair of n-by-n real nonsymmetric matrices A and B, the generalized eigenvalues (alphar +/- alphai∗i, beta), and optionally, the left and/or right generalized eigenvectors (VL and VR)
DGEHD2, dgehd2reduce a real general matrix A to upper Hessenberg form H by an orthogonal similarity transformation
DGEHRD, dgehrdreduce a real general matrix A to upper Hessenberg form H by an orthogonal similarity transformation
DGELQ2, dgelq2compute an LQ factorization of a real m by n matrix A
DGELQF, dgelqfcompute an LQ factorization of a real M-by-N matrix A
DGELS, dgelssolve overdetermined or underdetermined real linear systems involving an M-by-N matrix A, or its transpose, using a QR or LQ factorization of A
DGELSS, dgelsscompute the minimum norm solution to a real linear least squares problem
DGELSX, dgelsxcompute the minimum-norm solution to a real linear least squares problem
DGEQL2, dgeql2compute a QL factorization of a real m by n matrix A
DGEQLF, dgeqlfcompute a QL factorization of a real M-by-N matrix A
DGEQPF, dgeqpfcompute a QR factorization with column pivoting of a real M-by-N matrix A
DGEQR2, dgeqr2compute a QR factorization of a real m by n matrix A
DGEQRF, dgeqrfcompute a QR factorization of a real M-by-N matrix A
DGERFS, dgerfsimprove the computed solution to a system of linear equations and provides error bounds and backward error estimates for the solution
DGERQ2, dgerq2compute an RQ factorization of a real m by n matrix A
DGERQF, dgerqfcompute an RQ factorization of a real M-by-N matrix A
DGESV, dgesvcompute the solution to a real system of linear equations  A ∗ X = B,
DGESVD, dgesvdcompute the singular value decomposition (SVD) of a real M-by-N matrix A, optionally computing the left and/or right singular vectors
DGESVX, dgesvxuse the LU factorization to compute the solution to a real system of linear equations  A ∗ X = B,
DGETF2, dgetf2compute an LU factorization of a general m-by-n matrix A using partial pivoting with row interchanges
DGETRF, dgetrfcompute an LU factorization of a general M-by-N matrix A using partial pivoting with row interchanges
DGETRI, dgetricompute the inverse of a matrix using the LU factorization computed by DGETRF
DGETRS, dgetrssolve a system of linear equations  A ∗ X = B or A’ ∗ X = B with a general N-by-N matrix A using the LU factorization computed by DGETRF
DGGBAK, dggbakform the right or left eigenvectors of a real generalized eigenvalue problem A∗x = lambda∗B∗x, by backward transformation on the computed eigenvectors of the balanced pair of matrices output by DGGBAL
DGGBAL, dggbalbalance a pair of general real matrices (A,B)
DGGGLM, dggglmsolve a general Gauss-Markov linear model (GLM) problem
DGGHRD, dgghrdreduce 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
DGGLSE, dgglsesolve the linear equality-constrained least squares (LSE) problem
DGGQRF, dggqrfcompute a generalized QR factorization of an N-by-M matrix A and an N-by-P matrix B
DGGRQF, dggrqfcompute a generalized RQ factorization of an M-by-N matrix A and a P-by-N matrix B
DGGSVD, dggsvdcompute the generalized singular value decomposition (GSVD) of an M-by-N real matrix A and P-by-N real matrix B
DGGSVP, dggsvpcompute orthogonal matrices U, V and Q such that   N-K-L K L  U’∗A∗Q = K ( 0 A12 A13 ) if M-K-L >= 0
DGTCON, dgtconestimate the reciprocal of the condition number of a real tridiagonal matrix A using the LU factorization as computed by DGTTRF
DGTRFS, dgtrfsimprove the computed solution to a system of linear equations when the coefficient matrix is tridiagonal, and provides error bounds and backward error estimates for the solution
DGTSV, dgtsvsolve the equation   A∗X = B,
DGTSVX, dgtsvxuse the LU factorization to compute the solution to a real system of linear equations A ∗ X = B or A∗∗T ∗ X = B,
DGTTRF, dgttrfcompute an LU factorization of a real tridiagonal matrix A using elimination with partial pivoting and row interchanges
DGTTRS, dgttrssolve one of the systems of equations  A∗X = B or A’∗X = B,
DHGEQZ, dhgeqzimplement a single-/double-shift version of the QZ method for finding the generalized eigenvalues  w(j)=(ALPHAR(j) + i∗ALPHAI(j))/BETAR(j) of the equation   det( A - w(i) B ) = 0  In addition, the pair A,B may be reduced to generalized Schur form
DHSEIN, dhseinuse inverse iteration to find specified right and/or left eigenvectors of a real upper Hessenberg matrix H
DHSEQR, dhseqrcompute the eigenvalues of a real upper 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
DLABAD, dlabadtake as input the values computed by SLAMCH for underflow and overflow, and returns the square root of each of these values if the log of LARGE is sufficiently large
DLABRD, dlabrdreduce the first NB rows and columns of a real general m by n matrix A to upper or lower bidiagonal form by an orthogonal transformation Q’ ∗ A ∗ P, and returns the matrices X and Y which are needed to apply the transformation to the unreduced part of A
DLACON, dlaconestimate the 1-norm of a square, real matrix A
DLACPY, dlacpycopie all or part of a two-dimensional matrix A to another matrix B
DLADIV, dladivperform complex division in real arithmetic   a + i∗b  p + i∗q = ---------  c + i∗d  The algorithm is due to Robert L
DLAE2, dlae2compute the eigenvalues of a 2-by-2 symmetric matrix  [ A B ]  [ B C ]
DLAEBZ, dlaebzcontain the iteration loops which compute and use the function N(w), which is the count of eigenvalues of a symmetric tridiagonal matrix T less than or equal to its argument w
DLAED0, dlaed0compute all eigenvalues and corresponding eigenvectors of a symmetric tridiagonal matrix using the divide and conquer method
DLAED1, dlaed1compute the updated eigensystem of a diagonal matrix after modification by a rank-one symmetric matrix
DLAED2, dlaed2merge the two sets of eigenvalues together into a single sorted set
DLAED3, dlaed3find the roots of the secular equation, as defined by the values in D, W, and RHO, between KSTART and KSTOP
DLAED4, dlaed4subroutine 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
DLAED5, dlaed5subroutine 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
DLAED6, dlaed6compute the positive or negative root (closest to the origin) of  z(1) z(2) z(3) f(x) = rho + --------- + ---------- + ---------  d(1)-x d(2)-x d(3)-x  It is assumed that   if ORGATI = .true
DLAED7, dlaed7compute the updated eigensystem of a diagonal matrix after modification by a rank-one symmetric matrix
DLAED8, dlaed8merge the two sets of eigenvalues together into a single sorted set
DLAED9, dlaed9find the roots of the secular equation, as defined by the values in D, Z, and RHO, between KSTART and KSTOP
DLAEDA, dlaedacompute the Z vector corresponding to the merge step in the CURLVLth step of the merge process with TLVLS steps for the CURPBMth problem
DLAEIN, dlaeinuse inverse iteration to find a right or left eigenvector corresponding to the eigenvalue (WR,WI) of a real upper Hessenberg matrix H
DLAEV2, dlaev2compute the eigendecomposition of a 2-by-2 symmetric matrix  [ A B ]  [ B C ]
DLAEXC, dlaexcswap adjacent diagonal blocks T11 and T22 of order 1 or 2 in an upper quasi-triangular matrix T by an orthogonal similarity transformation
DLAG2, dlag2compute the eigenvalues of a 2 x 2 generalized eigenvalue problem A - w B, with scaling as necessary to avoid over-/underflow
DLAGTF, dlagtffactorize the matrix (T - lambda∗I), where T is an n by n tridiagonal matrix and lambda is a scalar, as   T - lambda∗I = PLU,
DLAGTM, dlagtmperform a matrix-vector product of the form   B := alpha ∗ A ∗ X + beta ∗ B  where A is a tridiagonal matrix of order N, B and X are N by NRHS matrices, and alpha and beta are real scalars, each of which may be 0., 1., or -1
DLAGTS, dlagtsmay be used to solve one of the systems of equations   (T - lambda∗I)∗x = y or (T - lambda∗I)’∗x = y,
DLAHQR, dlahqri an auxiliary routine called by DHSEQR to update the eigenvalues and Schur decomposition already computed by DHSEQR, by dealing with the Hessenberg submatrix in rows and columns ILO to IHI
DLAHRD, dlahrdreduce 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
DLAIC1, dlaic1applie one step of incremental condition estimation in its simplest version
DLALN2, dlaln2solve a system of the form (ca A - w D ) X = s B or (ca A’ - w D) X = s B with possible scaling ("s") and perturbation of A
DLAMCH, dlamchdetermine double precision machine parameters
DLAMRG, dlamrgwill 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
DLANGB, dlangbreturn 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
DLANGE, dlangereturn 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
DLANGT, dlangtreturn 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
DLANHS, dlanhsreturn 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
DLANSB, dlansbreturn 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
DLANSP, dlanspreturn 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
DLANST, dlanstreturn 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
DLANSY, dlansyreturn 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
DLANTB, dlantbreturn 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
DLANTP, dlantpreturn 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
DLANTR, dlantrreturn 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
DLANV2, dlanv2compute the Schur factorization of a real 2-by-2 nonsymmetric matrix in standard form
DLAPLL, dlaplltwo column vectors X and Y, let   A = ( X Y )
DLAPMT, dlapmtrearrange 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
DLAPY2, dlapy2return sqrt(x∗∗2+y∗∗2), taking care not to cause unnecessary overflow
DLAPY3, dlapy3return sqrt(x∗∗2+y∗∗2+z∗∗2), taking care not to cause unnecessary overflow
DLAQGB, dlaqgbequilibrate 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
DLAQGE, dlaqgeequilibrate a general M by N matrix A using the row and scaling factors in the vectors R and C
DLAQSB, dlaqsbequilibrate a symmetric band matrix A using the scaling factors in the vector S
DLAQSP, dlaqspequilibrate a symmetric matrix A using the scaling factors in the vector S
DLAQSY, dlaqsyequilibrate a symmetric matrix A using the scaling factors in the vector S
DLAQTR, dlaqtrsolve the real quasi-triangular system   op(T)∗p = scale∗c, if LREAL = .TRUE
DLAR2V, dlar2vapplie 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
DLARF, dlarfapplie a real elementary reflector H to a real m by n matrix C, from either the left or the right
DLARFB, dlarfbapplie a real block reflector H or its transpose H’ to a real m by n matrix C, from either the left or the right
DLARFG, dlarfggenerate a real elementary reflector H of order n, such that   H ∗ ( alpha ) = ( beta ), H’ ∗ H = I
DLARFT, dlarftform the triangular factor T of a real block reflector H of order n, which is defined as a product of k elementary reflectors
DLARFX, dlarfxapplie a real elementary reflector H to a real m by n matrix C, from either the left or the right
DLARGV, dlargvgenerate a vector of real plane rotations, determined by elements of the real vectors x and y
DLARNV, dlarnvreturn a vector of n random real numbers from a uniform or normal distribution
DLARTG, dlartggenerate a plane rotation so that   [ CS SN ]
DLARTV, dlartvapplie a vector of real plane rotations to elements of the real vectors x and y
DLARUV, dlaruvreturn a vector of n random real numbers from a uniform (0,1)
DLAS2, dlas2compute the singular values of the 2-by-2 matrix  [ F G ]  [ 0 H ]
DLASCL, dlasclmultiplie the M by N real matrix A by the real scalar CTO/CFROM
DLASET, dlasetinitialize an m-by-n matrix A to BETA on the diagonal and ALPHA on the offdiagonals
DLASQ1, dlasq1DLASQ1 computes the singular values of a real N-by-N bidiagonal  matrix with diagonal D and off-diagonal E
DLASQ2, dlasq2DLASQ2 computes the singular values of a real N-by-N unreduced  bidiagonal matrix with squared diagonal elements in Q and  squared off-diagonal elements in E
DLASQ3, dlasq3DLASQ3 is the workhorse of the whole bidiagonal SVD algorithm
DLASQ4, dlasq4DLASQ4 estimates TAU, the smallest eigenvalue of a matrix
DLASR, dlasrperform the transformation   A := P∗A, when SIDE = ’L’ or ’l’ ( Left-hand side )   A := A∗P’, when SIDE = ’R’ or ’r’ ( Right-hand side )  where A is an m by n real matrix and P is an orthogonal matrix,
DLASRT, dlasrtthe numbers in D in increasing order (if ID = ’I’) or in decreasing order (if ID = ’D’ )
DLASSQ, dlassqreturn the values scl and smsq such that   ( scl∗∗2 )∗smsq = x( 1 )∗∗2 +...+ x( n )∗∗2 + ( scale∗∗2 )∗sumsq,
DLASV2, dlasv2compute the singular value decomposition of a 2-by-2 triangular matrix  [ F G ]  [ 0 H ]
DLASWP, dlaswpperform a series of row interchanges on the matrix A
DLASY2, dlasy2solve for the N1 by N2 matrix X, 1 <= N1,N2 <= 2, in   op(TL)∗X + ISGN∗X∗op(TR) = SCALE∗B,
DLASYF, dlasyfcompute a partial factorization of a real symmetric matrix A using the Bunch-Kaufman diagonal pivoting method
DLATBS, dlatbssolve one of the triangular systems   A ∗x = s∗b or A’∗x = s∗b  with scaling to prevent overflow, where A is an upper or lower triangular band matrix
DLATPS, dlatpssolve one of the triangular systems   A ∗x = s∗b or A’∗x = s∗b  with scaling to prevent overflow, where A is an upper or lower triangular matrix stored in packed form
DLATRD, dlatrdreduce NB rows and columns of a real symmetric matrix A to symmetric tridiagonal form by an orthogonal similarity transformation Q’ ∗ A ∗ Q, and returns the matrices V and W which are needed to apply the transformation to the unreduced part of A
DLATRS, dlatrssolve one of the triangular systems   A ∗x = s∗b or A’∗x = s∗b  with scaling to prevent overflow
DLATZM, dlatzmapplie a Householder matrix generated by DTZRQF to a matrix
DLAUU2, dlauu2compute the product U ∗ U’ or L’ ∗ L, where the triangular factor U or L is stored in the upper or lower triangular part of the array A
DLAUUM, dlauumcompute the product U ∗ U’ or L’ ∗ L, where the triangular factor U or L is stored in the upper or lower triangular part of the array A
DLAZRO, dlazroinitialize a 2-D array A to BETA on the diagonal and ALPHA on the offdiagonals
DOPGTR, dopgtrgenerate 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
DOPMTR, dopmtroverwrite the general real M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
DORG2L, dorg2lgenerate an m by n real matrix Q with orthonormal columns,
DORG2R, dorg2rgenerate an m by n real matrix Q with orthonormal columns,
DORGBR, dorgbrgenerate one of the real orthogonal matrices Q or P∗∗T determined by DGEBRD when reducing a real matrix A to bidiagonal form
DORGHR, dorghrgenerate a real orthogonal matrix Q which is defined as the product of IHI-ILO elementary reflectors of order N, as returned by DGEHRD
DORGL2, dorgl2generate an m by n real matrix Q with orthonormal rows,
DORGLQ, dorglqgenerate an M-by-N real matrix Q with orthonormal rows,
DORGQL, dorgqlgenerate an M-by-N real matrix Q with orthonormal columns,
DORGQR, dorgqrgenerate an M-by-N real matrix Q with orthonormal columns,
DORGR2, dorgr2generate an m by n real matrix Q with orthonormal rows,
DORGRQ, dorgrqgenerate an M-by-N real matrix Q with orthonormal rows,
DORGTR, dorgtrgenerate a real orthogonal matrix Q which is defined as the product of n-1 elementary reflectors of order N, as returned by DSYTRD
DORM2L, dorm2loverwrite the general real m by n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’T’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’T’,
DORM2R, dorm2roverwrite the general real m by n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’T’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’T’,
DORMBR, dormbrVECT = ’Q’, DORMBR overwrites the general real M-by-N matrix C with  SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
DORMHR, dormhroverwrite the general real M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
DORML2, dorml2overwrite the general real m by n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’T’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’T’,
DORMLQ, dormlqoverwrite the general real M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
DORMQL, dormqloverwrite the general real M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
DORMQR, dormqroverwrite the general real M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
DORMR2, dormr2overwrite the general real m by n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’T’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’T’,
DORMRQ, dormrqoverwrite the general real M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
DORMTR, dormtroverwrite the general real M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
DPBCON, dpbconestimate the reciprocal of the condition number (in the 1-norm) of a real symmetric positive definite band matrix using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by DPBTRF
DPBEQU, dpbequcompute row and column scalings intended to equilibrate a symmetric positive definite band matrix A and reduce its condition number (with respect to the two-norm)
DPBRFS, dpbrfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric positive definite and banded, and provides error bounds and backward error estimates for the solution
DPBSTF, dpbstfcompute a split Cholesky factorization of a real symmetric positive definite band matrix A
DPBSV, dpbsvcompute the solution to a real system of linear equations  A ∗ X = B,
DPBSVX, dpbsvxuse the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T to compute the solution to a real system of linear equations  A ∗ X = B,
DPBTF2, dpbtf2compute the Cholesky factorization of a real symmetric positive definite band matrix A
DPBTRF, dpbtrfcompute the Cholesky factorization of a real symmetric positive definite band matrix A
DPBTRS, dpbtrssolve a system of linear equations A∗X = B with a symmetric positive definite band matrix A using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by DPBTRF
DPOCON, dpoconestimate the reciprocal of the condition number (in the 1-norm) of a real symmetric positive definite matrix using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by DPOTRF
DPOEQU, dpoequcompute row and column scalings intended to equilibrate a symmetric positive definite matrix A and reduce its condition number (with respect to the two-norm)
DPORFS, dporfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric positive definite,
DPOSV, dposvcompute the solution to a real system of linear equations  A ∗ X = B,
DPOSVX, dposvxuse the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T to compute the solution to a real system of linear equations  A ∗ X = B,
DPOTF2, dpotf2compute the Cholesky factorization of a real symmetric positive definite matrix A
DPOTRF, dpotrfcompute the Cholesky factorization of a real symmetric positive definite matrix A
DPOTRI, dpotricompute the inverse of a real symmetric positive definite matrix A using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by DPOTRF
DPOTRS, dpotrssolve a system of linear equations A∗X = B with a symmetric positive definite matrix A using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by DPOTRF
DPPCON, dppconestimate the reciprocal of the condition number (in the 1-norm) of a real symmetric positive definite packed matrix using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by DPPTRF
DPPEQU, dppequcompute row and column scalings intended to equilibrate a symmetric positive definite matrix A in packed storage and reduce its condition number (with respect to the two-norm)
DPPRFS, dpprfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric positive definite and packed, and provides error bounds and backward error estimates for the solution
DPPSV, dppsvcompute the solution to a real system of linear equations  A ∗ X = B,
DPPSVX, dppsvxuse the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T to compute the solution to a real system of linear equations  A ∗ X = B,
DPPTRF, dpptrfcompute the Cholesky factorization of a real symmetric positive definite matrix A stored in packed format
DPPTRI, dpptricompute the inverse of a real symmetric positive definite matrix A using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by DPPTRF
DPPTRS, dpptrssolve a system of linear equations A∗X = B with a symmetric positive definite matrix A in packed storage using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by DPPTRF
DPTCON, dptconcompute the reciprocal of the condition number (in the 1-norm) of a real symmetric positive definite tridiagonal matrix using the factorization A = L∗D∗L∗∗T or A = U∗∗T∗D∗U computed by DPTTRF
DPTEQR, dpteqrcompute all eigenvalues and, optionally, eigenvectors of a symmetric positive definite tridiagonal matrix by first factoring the matrix using DPTTRF, and then calling DBDSQR to compute the singular values of the bidiagonal factor
DPTRFS, dptrfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric positive definite and tridiagonal, and provides error bounds and backward error estimates for the solution
DPTSV, dptsvcompute the solution to a real system of linear equations A∗X = B, where A is an N-by-N symmetric positive definite tridiagonal matrix, and X and B are N-by-NRHS matrices
DPTSVX, dptsvxuse the factorization A = L∗D∗L∗∗T to compute the solution to a real system of linear equations A∗X = B, where A is an N-by-N symmetric positive definite tridiagonal matrix and X and B are N-by-NRHS matrices
DPTTRF, dpttrfcompute the factorization of a real symmetric positive definite tridiagonal matrix A
DPTTRS, dpttrssolve a system of linear equations A ∗ X = B with a symmetric positive definite tridiagonal matrix A using the factorization A = L∗D∗L∗∗T or A = U∗∗T∗D∗U computed by DPTTRF
DRSCL, drsclmultiplie an n-element real vector x by the real scalar 1/a
DSBEV, dsbevcompute all the eigenvalues and, optionally, eigenvectors of a real symmetric band matrix A
DSBEVD, dsbevdcompute all the eigenvalues and, optionally, eigenvectors of a real symmetric band matrix A
DSBEVX, dsbevxcompute selected eigenvalues and, optionally, eigenvectors of a real symmetric band matrix A
DSBGST, dsbgstreduce a real symmetric-definite banded generalized eigenproblem A∗x = lambda∗B∗x to standard form C∗y = lambda∗y,
DSBGV, dsbgvcompute all the eigenvalues, and optionally, the eigenvectors of a real generalized symmetric-definite banded eigenproblem, of the form A∗x=(lambda)∗B∗x
DSBTRD, dsbtrdreduce a real symmetric band matrix A to symmetric tridiagonal form T by an orthogonal similarity transformation
DSECND, dsecndreturn the user time for a process in seconds
DSPCON, dspconestimate the reciprocal of the condition number (in the 1-norm) of a real symmetric packed matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by DSPTRF
DSPEV, dspevcompute all the eigenvalues and, optionally, eigenvectors of a real symmetric matrix A in packed storage
DSPEVD, dspevdcompute all the eigenvalues and, optionally, eigenvectors of a real symmetric matrix A in packed storage
DSPEVX, dspevxcompute selected eigenvalues and, optionally, eigenvectors of a real symmetric matrix A in packed storage
DSPGST, dspgstreduce a real symmetric-definite generalized eigenproblem to standard form, using packed storage
DSPGV, dspgvcompute all the eigenvalues and, optionally, the eigenvectors of a real generalized symmetric-definite eigenproblem, of the form A∗x=(lambda)∗B∗x, A∗Bx=(lambda)∗x, or B∗A∗x=(lambda)∗x
DSPRFS, dsprfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric indefinite and packed, and provides error bounds and backward error estimates for the solution
DSPSV, dspsvcompute the solution to a real system of linear equations  A ∗ X = B,
DSPSVX, dspsvxuse the diagonal pivoting factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T to compute the solution to a real system of linear equations A ∗ X = B, where A is an N-by-N symmetric matrix stored in packed format and X and B are N-by-NRHS matrices
DSPTRD, dsptrdreduce a real symmetric matrix A stored in packed form to symmetric tridiagonal form T by an orthogonal similarity transformation
DSPTRF, dsptrfcompute the factorization of a real symmetric matrix A stored in packed format using the Bunch-Kaufman diagonal pivoting method
DSPTRI, dsptricompute the inverse of a real symmetric indefinite matrix A in packed storage using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by DSPTRF
DSPTRS, dsptrssolve a system of linear equations A∗X = B with a real symmetric matrix A stored in packed format using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by DSPTRF
DSTEBZ, dstebzcompute the eigenvalues of a symmetric tridiagonal matrix T
DSTEDC, dstedccompute all eigenvalues and, optionally, eigenvectors of a symmetric tridiagonal matrix using the divide and conquer method
DSTEIN, dsteincompute the eigenvectors of a real symmetric tridiagonal matrix T corresponding to specified eigenvalues, using inverse iteration
DSTEQR, dsteqrcompute all eigenvalues and, optionally, eigenvectors of a symmetric tridiagonal matrix using the implicit QL or QR method
DSTERF, dsterfcompute all eigenvalues of a symmetric tridiagonal matrix using the Pal-Walker-Kahan variant of the QL or QR algorithm
DSTEV, dstevcompute all eigenvalues and, optionally, eigenvectors of a real symmetric tridiagonal matrix A
DSTEVD, dstevdcompute all eigenvalues and, optionally, eigenvectors of a real symmetric tridiagonal matrix
DSTEVX, dstevxcompute selected eigenvalues and, optionally, eigenvectors of a real symmetric tridiagonal matrix A
DSYCON, dsyconestimate the reciprocal of the condition number (in the 1-norm) of a real symmetric matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by DSYTRF
DSYEV, dsyevcompute all eigenvalues and, optionally, eigenvectors of a real symmetric matrix A
DSYEVD, dsyevdcompute all eigenvalues and, optionally, eigenvectors of a real symmetric matrix A
DSYEVX, dsyevxcompute selected eigenvalues and, optionally, eigenvectors of a real symmetric matrix A
DSYGS2, dsygs2reduce a real symmetric-definite generalized eigenproblem to standard form
DSYGST, dsygstreduce a real symmetric-definite generalized eigenproblem to standard form
DSYGV, dsygvcompute all the eigenvalues, and optionally, the eigenvectors of a real generalized symmetric-definite eigenproblem, of the form A∗x=(lambda)∗B∗x, A∗Bx=(lambda)∗x, or B∗A∗x=(lambda)∗x
DSYRFS, dsyrfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric indefinite, and provides error bounds and backward error estimates for the solution
DSYSV, dsysvcompute the solution to a real system of linear equations  A ∗ X = B,
DSYSVX, dsysvxuse the diagonal pivoting factorization to compute the solution to a real system of linear equations A ∗ X = B,
DSYTD2, dsytd2reduce a real symmetric matrix A to symmetric tridiagonal form T by an orthogonal similarity transformation
DSYTF2, dsytf2compute the factorization of a real symmetric matrix A using the Bunch-Kaufman diagonal pivoting method
DSYTRD, dsytrdreduce a real symmetric matrix A to real symmetric tridiagonal form T by an orthogonal similarity transformation
DSYTRF, dsytrfcompute the factorization of a real symmetric matrix A using the Bunch-Kaufman diagonal pivoting method
DSYTRI, dsytricompute the inverse of a real symmetric indefinite matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by DSYTRF
DSYTRS, dsytrssolve a system of linear equations A∗X = B with a real symmetric matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by DSYTRF
DTBCON, dtbconestimate the reciprocal of the condition number of a triangular band matrix A, in either the 1-norm or the infinity-norm
DTBRFS, dtbrfsprovide error bounds and backward error estimates for the solution to a system of linear equations with a triangular band coefficient matrix
DTBTRS, dtbtrssolve a triangular system of the form   A ∗ X = B or A∗∗T ∗ X = B,
DTGEVC, dtgevccompute some or all of the right and/or left generalized eigenvectors of a pair of real upper triangular matrices (A,B)
DTGSJA, dtgsjacompute the generalized singular value decomposition (GSVD) of two real upper triangular (or trapezoidal) matrices A and B
DTPCON, dtpconestimate the reciprocal of the condition number of a packed triangular matrix A, in either the 1-norm or the infinity-norm
DTPRFS, dtprfsprovide error bounds and backward error estimates for the solution to a system of linear equations with a triangular packed coefficient matrix
DTPTRI, dtptricompute the inverse of a real upper or lower triangular matrix A stored in packed format
DTPTRS, dtptrssolve a triangular system of the form   A ∗ X = B or A∗∗T ∗ X = B,
DTRCON, dtrconestimate the reciprocal of the condition number of a triangular matrix A, in either the 1-norm or the infinity-norm
DTREVC, dtrevccompute some or all of the right and/or left eigenvectors of a real upper quasi-triangular matrix T
DTREXC, dtrexcreorder the real Schur factorization of a real matrix A = Q∗T∗Q∗∗T, so that the diagonal block of T with row index IFST is moved to row ILST
DTRRFS, dtrrfsprovide error bounds and backward error estimates for the solution to a system of linear equations with a triangular coefficient matrix
DTRSEN, dtrsenreorder the real Schur factorization of a real matrix A = Q∗T∗Q∗∗T, so that a selected cluster of eigenvalues appears in the leading diagonal blocks of the upper quasi-triangular matrix T,
DTRSNA, dtrsnaestimate reciprocal condition numbers for specified eigenvalues and/or right eigenvectors of a real upper quasi-triangular matrix T (or of any matrix Q∗T∗Q∗∗T with Q orthogonal)
DTRSYL, dtrsylsolve the real Sylvester matrix equation
DTRTI2, dtrti2compute the inverse of a real upper or lower triangular matrix
DTRTRI, dtrtricompute the inverse of a real upper or lower triangular matrix A
DTRTRS, dtrtrssolve a triangular system of the form   A ∗ X = B or A∗∗T ∗ X = B,
DTZRQF, dtzrqfreduce the M-by-N ( M<=N ) real upper trapezoidal matrix A to upper triangular form by means of orthogonal transformations
DZSUM1, dzsum1take the sum of the absolute values of a complex vector and returns a double precision result
ICMAX1, icmax1find the index of the element whose real part has maximum absolute value
ILAENV, ilaenvi called from the LAPACK routines to choose problem-dependent parameters for the local environment
IZMAX1, izmax1find the index of the element whose real part has maximum absolute value
lapack
LSAME, lsamereturn .TRUE
LSAMEN, lsamentest if the first N letters of CA are the same as the first N letters of CB, regardless of case
SBDSQR, sbdsqrcompute the singular value decomposition (SVD) of a real N-by-N (upper or lower) bidiagonal matrix B
SCSUM1, scsum1take the sum of the absolute values of a complex vector and returns a single precision result
SDISNA, sdisnacompute 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
SECOND, secondreturn the user time for a process in seconds
SGBBRD, sgbbrdreduce a real general m-by-n band matrix A to upper bidiagonal form B by an orthogonal transformation
SGBCON, sgbconestimate the reciprocal of the condition number of a real general band matrix A, in either the 1-norm or the infinity-norm,
SGBEQU, sgbequcompute row and column scalings intended to equilibrate an M-by-N band matrix A and reduce its condition number
SGBRFS, sgbrfsimprove 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
SGBSV, sgbsvcompute the solution to a real system of linear equations A ∗ X = B, where A is a band matrix of order N with KL subdiagonals and KU superdiagonals, and X and B are N-by-NRHS matrices
SGBSVX, sgbsvxuse the LU factorization to compute the solution to a real system of linear equations A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B,
SGBTF2, sgbtf2compute an LU factorization of a real m-by-n band matrix A using partial pivoting with row interchanges
SGBTRF, sgbtrfcompute an LU factorization of a real m-by-n band matrix A using partial pivoting with row interchanges
SGBTRS, sgbtrssolve a system of linear equations  A ∗ X = B or A’ ∗ X = B with a general band matrix A using the LU factorization computed by SGBTRF
SGEBAK, sgebakform the right or left eigenvectors of a real general matrix by backward transformation on the computed eigenvectors of the balanced matrix output by SGEBAL
SGEBAL, sgebalbalance a general real matrix A
SGEBD2, sgebd2reduce a real general m by n matrix A to upper or lower bidiagonal form B by an orthogonal transformation
SGEBRD, sgebrdreduce a general real M-by-N matrix A to upper or lower bidiagonal form B by an orthogonal transformation
SGECON, sgeconestimate the reciprocal of the condition number of a general real matrix A, in either the 1-norm or the infinity-norm, using the LU factorization computed by SGETRF
SGEEQU, sgeequcompute row and column scalings intended to equilibrate an M-by-N matrix A and reduce its condition number
SGEES, sgeescompute for an N-by-N real nonsymmetric matrix A, the eigenvalues, the real Schur form T, and, optionally, the matrix of Schur vectors Z
SGEESX, sgeesxcompute for an N-by-N real nonsymmetric matrix A, the eigenvalues, the real Schur form T, and, optionally, the matrix of Schur vectors Z
SGEEV, sgeevcompute for an N-by-N real nonsymmetric matrix A, the eigenvalues and, optionally, the left and/or right eigenvectors
SGEEVX, sgeevxcompute for an N-by-N real nonsymmetric matrix A, the eigenvalues and, optionally, the left and/or right eigenvectors
SGEGS, sgegscompute for a pair of N-by-N real nonsymmetric matrices A, B
SGEGV, sgegvcompute for a pair of n-by-n real nonsymmetric matrices A and B, the generalized eigenvalues (alphar +/- alphai∗i, beta), and optionally, the left and/or right generalized eigenvectors (VL and VR)
SGEHD2, sgehd2reduce a real general matrix A to upper Hessenberg form H by an orthogonal similarity transformation
SGEHRD, sgehrdreduce a real general matrix A to upper Hessenberg form H by an orthogonal similarity transformation
SGELQ2, sgelq2compute an LQ factorization of a real m by n matrix A
SGELQF, sgelqfcompute an LQ factorization of a real M-by-N matrix A
SGELS, sgelssolve overdetermined or underdetermined real linear systems involving an M-by-N matrix A, or its transpose, using a QR or LQ factorization of A
SGELSS, sgelsscompute the minimum norm solution to a real linear least squares problem
SGELSX, sgelsxcompute the minimum-norm solution to a real linear least squares problem
SGEQL2, sgeql2compute a QL factorization of a real m by n matrix A
SGEQLF, sgeqlfcompute a QL factorization of a real M-by-N matrix A
SGEQPF, sgeqpfcompute a QR factorization with column pivoting of a real M-by-N matrix A
SGEQR2, sgeqr2compute a QR factorization of a real m by n matrix A
SGEQRF, sgeqrfcompute a QR factorization of a real M-by-N matrix A
SGERFS, sgerfsimprove the computed solution to a system of linear equations and provides error bounds and backward error estimates for the solution
SGERQ2, sgerq2compute an RQ factorization of a real m by n matrix A
SGERQF, sgerqfcompute an RQ factorization of a real M-by-N matrix A
SGESV, sgesvcompute the solution to a real system of linear equations  A ∗ X = B,
SGESVD, sgesvdcompute the singular value decomposition (SVD) of a real M-by-N matrix A, optionally computing the left and/or right singular vectors
SGESVX, sgesvxuse the LU factorization to compute the solution to a real system of linear equations  A ∗ X = B,
SGETF2, sgetf2compute an LU factorization of a general m-by-n matrix A using partial pivoting with row interchanges
SGETRF, sgetrfcompute an LU factorization of a general M-by-N matrix A using partial pivoting with row interchanges
SGETRI, sgetricompute the inverse of a matrix using the LU factorization computed by SGETRF
SGETRS, sgetrssolve a system of linear equations  A ∗ X = B or A’ ∗ X = B with a general N-by-N matrix A using the LU factorization computed by SGETRF
SGGBAK, sggbakform the right or left eigenvectors of a real generalized eigenvalue problem A∗x = lambda∗B∗x, by backward transformation on the computed eigenvectors of the balanced pair of matrices output by SGGBAL
SGGBAL, sggbalbalance a pair of general real matrices (A,B)
SGGGLM, sggglmsolve a general Gauss-Markov linear model (GLM) problem
SGGHRD, sgghrdreduce 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
SGGLSE, sgglsesolve the linear equality-constrained least squares (LSE) problem
SGGQRF, sggqrfcompute a generalized QR factorization of an N-by-M matrix A and an N-by-P matrix B
SGGRQF, sggrqfcompute a generalized RQ factorization of an M-by-N matrix A and a P-by-N matrix B
SGGSVD, sggsvdcompute the generalized singular value decomposition (GSVD) of an M-by-N real matrix A and P-by-N real matrix B
SGGSVP, sggsvpcompute orthogonal matrices U, V and Q such that   N-K-L K L  U’∗A∗Q = K ( 0 A12 A13 ) if M-K-L >= 0
SGTCON, sgtconestimate the reciprocal of the condition number of a real tridiagonal matrix A using the LU factorization as computed by SGTTRF
SGTRFS, sgtrfsimprove the computed solution to a system of linear equations when the coefficient matrix is tridiagonal, and provides error bounds and backward error estimates for the solution
SGTSV, sgtsvsolve the equation   A∗X = B,
SGTSVX, sgtsvxuse the LU factorization to compute the solution to a real system of linear equations A ∗ X = B or A∗∗T ∗ X = B,
SGTTRF, sgttrfcompute an LU factorization of a real tridiagonal matrix A using elimination with partial pivoting and row interchanges
SGTTRS, sgttrssolve one of the systems of equations  A∗X = B or A’∗X = B,
SHGEQZ, shgeqzimplement a single-/double-shift version of the QZ method for finding the generalized eigenvalues  w(j)=(ALPHAR(j) + i∗ALPHAI(j))/BETAR(j) of the equation   det( A - w(i) B ) = 0  In addition, the pair A,B may be reduced to generalized Schur form
SHSEIN, shseinuse inverse iteration to find specified right and/or left eigenvectors of a real upper Hessenberg matrix H
SHSEQR, shseqrcompute the eigenvalues of a real upper 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
SLABAD, slabadtake as input the values computed by SLAMCH for underflow and overflow, and returns the square root of each of these values if the log of LARGE is sufficiently large
SLABRD, slabrdreduce the first NB rows and columns of a real general m by n matrix A to upper or lower bidiagonal form by an orthogonal transformation Q’ ∗ A ∗ P, and returns the matrices X and Y which are needed to apply the transformation to the unreduced part of A
SLACON, slaconestimate the 1-norm of a square, real matrix A
SLACPY, slacpycopie all or part of a two-dimensional matrix A to another matrix B
SLADIV, sladivperform complex division in real arithmetic   a + i∗b  p + i∗q = ---------  c + i∗d  The algorithm is due to Robert L
SLAE2, slae2compute the eigenvalues of a 2-by-2 symmetric matrix  [ A B ]  [ B C ]
SLAEBZ, slaebzcontain the iteration loops which compute and use the function N(w), which is the count of eigenvalues of a symmetric tridiagonal matrix T less than or equal to its argument w
SLAED0, slaed0compute all eigenvalues and corresponding eigenvectors of a symmetric tridiagonal matrix using the divide and conquer method
SLAED1, slaed1compute the updated eigensystem of a diagonal matrix after modification by a rank-one symmetric matrix
SLAED2, slaed2merge the two sets of eigenvalues together into a single sorted set
SLAED3, slaed3find the roots of the secular equation, as defined by the values in D, W, and RHO, between KSTART and KSTOP
SLAED4, slaed4subroutine 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, slaed5subroutine 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
SLAED6, slaed6compute the positive or negative root (closest to the origin) of  z(1) z(2) z(3) f(x) = rho + --------- + ---------- + ---------  d(1)-x d(2)-x d(3)-x  It is assumed that   if ORGATI = .true
SLAED7, slaed7compute the updated eigensystem of a diagonal matrix after modification by a rank-one symmetric matrix
SLAED8, slaed8merge the two sets of eigenvalues together into a single sorted set
SLAED9, slaed9find the roots of the secular equation, as defined by the values in D, Z, and RHO, between KSTART and KSTOP
SLAEDA, slaedacompute the Z vector corresponding to the merge step in the CURLVLth step of the merge process with TLVLS steps for the CURPBMth problem
SLAEIN, slaeinuse inverse iteration to find a right or left eigenvector corresponding to the eigenvalue (WR,WI) of a real upper Hessenberg matrix H
SLAEV2, slaev2compute the eigendecomposition of a 2-by-2 symmetric matrix  [ A B ]  [ B C ]
SLAEXC, slaexcswap adjacent diagonal blocks T11 and T22 of order 1 or 2 in an upper quasi-triangular matrix T by an orthogonal similarity transformation
SLAG2, slag2compute the eigenvalues of a 2 x 2 generalized eigenvalue problem A - w B, with scaling as necessary to avoid over-/underflow
SLAGTF, slagtffactorize the matrix (T - lambda∗I), where T is an n by n tridiagonal matrix and lambda is a scalar, as   T - lambda∗I = PLU,
SLAGTM, slagtmperform a matrix-vector product of the form   B := alpha ∗ A ∗ X + beta ∗ B  where A is a tridiagonal matrix of order N, B and X are N by NRHS matrices, and alpha and beta are real scalars, each of which may be 0., 1., or -1
SLAGTS, slagtsmay be used to solve one of the systems of equations   (T - lambda∗I)∗x = y or (T - lambda∗I)’∗x = y,
SLAHQR, slahqri 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, slahrdreduce 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, slaic1applie one step of incremental condition estimation in its simplest version
SLALN2, slaln2solve a system of the form (ca A - w D ) X = s B or (ca A’ - w D) X = s B with possible scaling ("s") and perturbation of A
SLAMCH, slamchdetermine single precision machine parameters
SLAMRG, slamrgwill 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, slangbreturn 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, slangereturn 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, slangtreturn 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, slanhsreturn 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, slansbreturn 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, slanspreturn 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, slanstreturn 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, slansyreturn 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, slantbreturn 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, slantpreturn 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, slantrreturn 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, slanv2compute the Schur factorization of a real 2-by-2 nonsymmetric matrix in standard form
SLAPLL, slaplltwo column vectors X and Y, let   A = ( X Y )
SLAPMT, slapmtrearrange 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, slapy2return sqrt(x∗∗2+y∗∗2), taking care not to cause unnecessary overflow
SLAPY3, slapy3return sqrt(x∗∗2+y∗∗2+z∗∗2), taking care not to cause unnecessary overflow
SLAQGB, slaqgbequilibrate 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, slaqgeequilibrate a general M by N matrix A using the row and scaling factors in the vectors R and C
SLAQSB, slaqsbequilibrate a symmetric band matrix A using the scaling factors in the vector S
SLAQSP, slaqspequilibrate a symmetric matrix A using the scaling factors in the vector S
SLAQSY, slaqsyequilibrate a symmetric matrix A using the scaling factors in the vector S
SLAQTR, slaqtrsolve the real quasi-triangular system   op(T)∗p = scale∗c, if LREAL = .TRUE
SLAR2V, slar2vapplie 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, slarfapplie a real elementary reflector H to a real m by n matrix C, from either the left or the right
SLARFB, slarfbapplie 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, slarfggenerate a real elementary reflector H of order n, such that   H ∗ ( alpha ) = ( beta ), H’ ∗ H = I
SLARFT, slarftform the triangular factor T of a real block reflector H of order n, which is defined as a product of k elementary reflectors
SLARFX, slarfxapplie a real elementary reflector H to a real m by n matrix C, from either the left or the right
SLARGV, slargvgenerate a vector of real plane rotations, determined by elements of the real vectors x and y
SLARNV, slarnvreturn a vector of n random real numbers from a uniform or normal distribution
SLARTG, slartggenerate a plane rotation so that   [ CS SN ]
SLARTV, slartvapplie a vector of real plane rotations to elements of the real vectors x and y
SLARUV, slaruvreturn a vector of n random real numbers from a uniform (0,1)
SLAS2, slas2compute the singular values of the 2-by-2 matrix  [ F G ]  [ 0 H ]
SLASCL, slasclmultiplie the M by N real matrix A by the real scalar CTO/CFROM
SLASET, slasetinitialize an m-by-n matrix A to BETA on the diagonal and ALPHA on the offdiagonals
SLASQ1, slasq1SLASQ1 computes the singular values of a real N-by-N bidiagonal  matrix with diagonal D and off-diagonal E
SLASQ2, slasq2SLASQ2 computes the singular values of a real N-by-N unreduced  bidiagonal matrix with squared diagonal elements in Q and  squared off-diagonal elements in E
SLASQ3, slasq3SLASQ3 is the workhorse of the whole bidiagonal SVD algorithm
SLASQ4, slasq4SLASQ4 estimates TAU, the smallest eigenvalue of a matrix
SLASR, slasrperform the transformation   A := P∗A, when SIDE = ’L’ or ’l’ ( Left-hand side )   A := A∗P’, when SIDE = ’R’ or ’r’ ( Right-hand side )  where A is an m by n real matrix and P is an orthogonal matrix,
SLASRT, slasrtthe numbers in D in increasing order (if ID = ’I’) or in decreasing order (if ID = ’D’ )
SLASSQ, slassqreturn the values scl and smsq such that   ( scl∗∗2 )∗smsq = x( 1 )∗∗2 +...+ x( n )∗∗2 + ( scale∗∗2 )∗sumsq,
SLASV2, slasv2compute the singular value decomposition of a 2-by-2 triangular matrix  [ F G ]  [ 0 H ]
SLASWP, slaswpperform a series of row interchanges on the matrix A
SLASY2, slasy2solve for the N1 by N2 matrix X, 1 <= N1,N2 <= 2, in   op(TL)∗X + ISGN∗X∗op(TR) = SCALE∗B,
SLASYF, slasyfcompute a partial factorization of a real symmetric matrix A using the Bunch-Kaufman diagonal pivoting method
SLATBS, slatbssolve one of the triangular systems   A ∗x = s∗b or A’∗x = s∗b  with scaling to prevent overflow, where A is an upper or lower triangular band matrix
SLATPS, slatpssolve one of the triangular systems   A ∗x = s∗b or A’∗x = s∗b  with scaling to prevent overflow, where A is an upper or lower triangular matrix stored in packed form
SLATRD, slatrdreduce NB rows and columns of a real symmetric matrix A to symmetric tridiagonal form by an orthogonal similarity transformation Q’ ∗ A ∗ Q, and returns the matrices V and W which are needed to apply the transformation to the unreduced part of A
SLATRS, slatrssolve one of the triangular systems   A ∗x = s∗b or A’∗x = s∗b  with scaling to prevent overflow
SLATZM, slatzmapplie a Householder matrix generated by STZRQF to a matrix
SLAUU2, slauu2compute the product U ∗ U’ or L’ ∗ L, where the triangular factor U or L is stored in the upper or lower triangular part of the array A
SLAUUM, slauumcompute the product U ∗ U’ or L’ ∗ L, where the triangular factor U or L is stored in the upper or lower triangular part of the array A
SLAZRO, slazroinitialize a 2-D array A to BETA on the diagonal and ALPHA on the offdiagonals
SOPGTR, sopgtrgenerate a real orthogonal matrix Q which is defined as the product of n-1 elementary reflectors H(i) of order n, as returned by SSPTRD using packed storage
SOPMTR, sopmtroverwrite the general real M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
SORG2L, sorg2lgenerate an m by n real matrix Q with orthonormal columns,
SORG2R, sorg2rgenerate an m by n real matrix Q with orthonormal columns,
SORGBR, sorgbrgenerate one of the real orthogonal matrices Q or P∗∗T determined by SGEBRD when reducing a real matrix A to bidiagonal form
SORGHR, sorghrgenerate a real orthogonal matrix Q which is defined as the product of IHI-ILO elementary reflectors of order N, as returned by SGEHRD
SORGL2, sorgl2generate an m by n real matrix Q with orthonormal rows,
SORGLQ, sorglqgenerate an M-by-N real matrix Q with orthonormal rows,
SORGQL, sorgqlgenerate an M-by-N real matrix Q with orthonormal columns,
SORGQR, sorgqrgenerate an M-by-N real matrix Q with orthonormal columns,
SORGR2, sorgr2generate an m by n real matrix Q with orthonormal rows,
SORGRQ, sorgrqgenerate an M-by-N real matrix Q with orthonormal rows,
SORGTR, sorgtrgenerate a real orthogonal matrix Q which is defined as the product of n-1 elementary reflectors of order N, as returned by SSYTRD
SORM2L, sorm2loverwrite the general real m by n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’T’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’T’,
SORM2R, sorm2roverwrite the general real m by n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’T’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’T’,
SORMBR, sormbrVECT = ’Q’, SORMBR overwrites the general real M-by-N matrix C with  SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
SORMHR, sormhroverwrite the general real M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
SORML2, sorml2overwrite the general real m by n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’T’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’T’,
SORMLQ, sormlqoverwrite the general real M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
SORMQL, sormqloverwrite the general real M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
SORMQR, sormqroverwrite the general real M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
SORMR2, sormr2overwrite the general real m by n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’T’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’T’,
SORMRQ, sormrqoverwrite the general real M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
SORMTR, sormtroverwrite the general real M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
SPBCON, spbconestimate the reciprocal of the condition number (in the 1-norm) of a real symmetric positive definite band matrix using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by SPBTRF
SPBEQU, spbequcompute row and column scalings intended to equilibrate a symmetric positive definite band matrix A and reduce its condition number (with respect to the two-norm)
SPBRFS, spbrfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric positive definite and banded, and provides error bounds and backward error estimates for the solution
SPBSTF, spbstfcompute a split Cholesky factorization of a real symmetric positive definite band matrix A
SPBSV, spbsvcompute the solution to a real system of linear equations  A ∗ X = B,
SPBSVX, spbsvxuse the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T to compute the solution to a real system of linear equations  A ∗ X = B,
SPBTF2, spbtf2compute the Cholesky factorization of a real symmetric positive definite band matrix A
SPBTRF, spbtrfcompute the Cholesky factorization of a real symmetric positive definite band matrix A
SPBTRS, spbtrssolve a system of linear equations A∗X = B with a symmetric positive definite band matrix A using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by SPBTRF
SPOCON, spoconestimate the reciprocal of the condition number (in the 1-norm) of a real symmetric positive definite matrix using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by SPOTRF
SPOEQU, spoequcompute row and column scalings intended to equilibrate a symmetric positive definite matrix A and reduce its condition number (with respect to the two-norm)
SPORFS, sporfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric positive definite,
SPOSV, sposvcompute the solution to a real system of linear equations  A ∗ X = B,
SPOSVX, sposvxuse the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T to compute the solution to a real system of linear equations  A ∗ X = B,
SPOTF2, spotf2compute the Cholesky factorization of a real symmetric positive definite matrix A
SPOTRF, spotrfcompute the Cholesky factorization of a real symmetric positive definite matrix A
SPOTRI, spotricompute the inverse of a real symmetric positive definite matrix A using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by SPOTRF
SPOTRS, spotrssolve a system of linear equations A∗X = B with a symmetric positive definite matrix A using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by SPOTRF
SPPCON, sppconestimate the reciprocal of the condition number (in the 1-norm) of a real symmetric positive definite packed matrix using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by SPPTRF
SPPEQU, sppequcompute row and column scalings intended to equilibrate a symmetric positive definite matrix A in packed storage and reduce its condition number (with respect to the two-norm)
SPPRFS, spprfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric positive definite and packed, and provides error bounds and backward error estimates for the solution
SPPSV, sppsvcompute the solution to a real system of linear equations  A ∗ X = B,
SPPSVX, sppsvxuse the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T to compute the solution to a real system of linear equations  A ∗ X = B,
SPPTRF, spptrfcompute the Cholesky factorization of a real symmetric positive definite matrix A stored in packed format
SPPTRI, spptricompute the inverse of a real symmetric positive definite matrix A using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by SPPTRF
SPPTRS, spptrssolve a system of linear equations A∗X = B with a symmetric positive definite matrix A in packed storage using the Cholesky factorization A = U∗∗T∗U or A = L∗L∗∗T computed by SPPTRF
SPTCON, sptconcompute the reciprocal of the condition number (in the 1-norm) of a real symmetric positive definite tridiagonal matrix using the factorization A = L∗D∗L∗∗T or A = U∗∗T∗D∗U computed by SPTTRF
SPTEQR, spteqrcompute all eigenvalues and, optionally, eigenvectors of a symmetric positive definite tridiagonal matrix by first factoring the matrix using SPTTRF, and then calling SBDSQR to compute the singular values of the bidiagonal factor
SPTRFS, sptrfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric positive definite and tridiagonal, and provides error bounds and backward error estimates for the solution
SPTSV, sptsvcompute the solution to a real system of linear equations A∗X = B, where A is an N-by-N symmetric positive definite tridiagonal matrix, and X and B are N-by-NRHS matrices
SPTSVX, sptsvxuse the factorization A = L∗D∗L∗∗T to compute the solution to a real system of linear equations A∗X = B, where A is an N-by-N symmetric positive definite tridiagonal matrix and X and B are N-by-NRHS matrices
SPTTRF, spttrfcompute the factorization of a real symmetric positive definite tridiagonal matrix A
SPTTRS, spttrssolve a system of linear equations A ∗ X = B with a symmetric positive definite tridiagonal matrix A using the factorization A = L∗D∗L∗∗T or A = U∗∗T∗D∗U computed by SPTTRF
SRSCL, srsclmultiplie an n-element real vector x by the real scalar 1/a
SSBEV, ssbevcompute all the eigenvalues and, optionally, eigenvectors of a real symmetric band matrix A
SSBEVD, ssbevdcompute all the eigenvalues and, optionally, eigenvectors of a real symmetric band matrix A
SSBEVX, ssbevxcompute selected eigenvalues and, optionally, eigenvectors of a real symmetric band matrix A
SSBGST, ssbgstreduce a real symmetric-definite banded generalized eigenproblem A∗x = lambda∗B∗x to standard form C∗y = lambda∗y,
SSBGV, ssbgvcompute all the eigenvalues, and optionally, the eigenvectors of a real generalized symmetric-definite banded eigenproblem, of the form A∗x=(lambda)∗B∗x
SSBTRD, ssbtrdreduce a real symmetric band matrix A to symmetric tridiagonal form T by an orthogonal similarity transformation
SSPCON, sspconestimate the reciprocal of the condition number (in the 1-norm) of a real symmetric packed matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by SSPTRF
SSPEV, sspevcompute all the eigenvalues and, optionally, eigenvectors of a real symmetric matrix A in packed storage
SSPEVD, sspevdcompute all the eigenvalues and, optionally, eigenvectors of a real symmetric matrix A in packed storage
SSPEVX, sspevxcompute selected eigenvalues and, optionally, eigenvectors of a real symmetric matrix A in packed storage
SSPGST, sspgstreduce a real symmetric-definite generalized eigenproblem to standard form, using packed storage
SSPGV, sspgvcompute all the eigenvalues and, optionally, the eigenvectors of a real generalized symmetric-definite eigenproblem, of the form A∗x=(lambda)∗B∗x, A∗Bx=(lambda)∗x, or B∗A∗x=(lambda)∗x
SSPRFS, ssprfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric indefinite and packed, and provides error bounds and backward error estimates for the solution
SSPSV, sspsvcompute the solution to a real system of linear equations  A ∗ X = B,
SSPSVX, sspsvxuse the diagonal pivoting factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T to compute the solution to a real system of linear equations A ∗ X = B, where A is an N-by-N symmetric matrix stored in packed format and X and B are N-by-NRHS matrices
SSPTRD, ssptrdreduce a real symmetric matrix A stored in packed form to symmetric tridiagonal form T by an orthogonal similarity transformation
SSPTRF, ssptrfcompute the factorization of a real symmetric matrix A stored in packed format using the Bunch-Kaufman diagonal pivoting method
SSPTRI, ssptricompute the inverse of a real symmetric indefinite matrix A in packed storage using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by SSPTRF
SSPTRS, ssptrssolve a system of linear equations A∗X = B with a real symmetric matrix A stored in packed format using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by SSPTRF
SSTEBZ, sstebzcompute the eigenvalues of a symmetric tridiagonal matrix T
SSTEDC, sstedccompute all eigenvalues and, optionally, eigenvectors of a symmetric tridiagonal matrix using the divide and conquer method
SSTEIN, ssteincompute the eigenvectors of a real symmetric tridiagonal matrix T corresponding to specified eigenvalues, using inverse iteration
SSTEQR, ssteqrcompute all eigenvalues and, optionally, eigenvectors of a symmetric tridiagonal matrix using the implicit QL or QR method
SSTERF, ssterfcompute all eigenvalues of a symmetric tridiagonal matrix using the Pal-Walker-Kahan variant of the QL or QR algorithm
SSTEV, sstevcompute all eigenvalues and, optionally, eigenvectors of a real symmetric tridiagonal matrix A
SSTEVD, sstevdcompute all eigenvalues and, optionally, eigenvectors of a real symmetric tridiagonal matrix
SSTEVX, sstevxcompute selected eigenvalues and, optionally, eigenvectors of a real symmetric tridiagonal matrix A
SSYCON, ssyconestimate the reciprocal of the condition number (in the 1-norm) of a real symmetric matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by SSYTRF
SSYEV, ssyevcompute all eigenvalues and, optionally, eigenvectors of a real symmetric matrix A
SSYEVD, ssyevdcompute all eigenvalues and, optionally, eigenvectors of a real symmetric matrix A
SSYEVX, ssyevxcompute selected eigenvalues and, optionally, eigenvectors of a real symmetric matrix A
SSYGS2, ssygs2reduce a real symmetric-definite generalized eigenproblem to standard form
SSYGST, ssygstreduce a real symmetric-definite generalized eigenproblem to standard form
SSYGV, ssygvcompute all the eigenvalues, and optionally, the eigenvectors of a real generalized symmetric-definite eigenproblem, of the form A∗x=(lambda)∗B∗x, A∗Bx=(lambda)∗x, or B∗A∗x=(lambda)∗x
SSYRFS, ssyrfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric indefinite, and provides error bounds and backward error estimates for the solution
SSYSV, ssysvcompute the solution to a real system of linear equations  A ∗ X = B,
SSYSVX, ssysvxuse the diagonal pivoting factorization to compute the solution to a real system of linear equations A ∗ X = B,
SSYTD2, ssytd2reduce a real symmetric matrix A to symmetric tridiagonal form T by an orthogonal similarity transformation
SSYTF2, ssytf2compute the factorization of a real symmetric matrix A using the Bunch-Kaufman diagonal pivoting method
SSYTRD, ssytrdreduce a real symmetric matrix A to real symmetric tridiagonal form T by an orthogonal similarity transformation
SSYTRF, ssytrfcompute the factorization of a real symmetric matrix A using the Bunch-Kaufman diagonal pivoting method
SSYTRI, ssytricompute the inverse of a real symmetric indefinite matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by SSYTRF
SSYTRS, ssytrssolve a system of linear equations A∗X = B with a real symmetric matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by SSYTRF
STBCON, stbconestimate the reciprocal of the condition number of a triangular band matrix A, in either the 1-norm or the infinity-norm
STBRFS, stbrfsprovide error bounds and backward error estimates for the solution to a system of linear equations with a triangular band coefficient matrix
STBTRS, stbtrssolve a triangular system of the form   A ∗ X = B or A∗∗T ∗ X = B,
STGEVC, stgevccompute some or all of the right and/or left generalized eigenvectors of a pair of real upper triangular matrices (A,B)
STGSJA, stgsjacompute the generalized singular value decomposition (GSVD) of two real upper triangular (or trapezoidal) matrices A and B
STPCON, stpconestimate the reciprocal of the condition number of a packed triangular matrix A, in either the 1-norm or the infinity-norm
STPRFS, stprfsprovide error bounds and backward error estimates for the solution to a system of linear equations with a triangular packed coefficient matrix
STPTRI, stptricompute the inverse of a real upper or lower triangular matrix A stored in packed format
STPTRS, stptrssolve a triangular system of the form   A ∗ X = B or A∗∗T ∗ X = B,
STRCON, strconestimate the reciprocal of the condition number of a triangular matrix A, in either the 1-norm or the infinity-norm
STREVC, strevccompute some or all of the right and/or left eigenvectors of a real upper quasi-triangular matrix T
STREXC, strexcreorder the real Schur factorization of a real matrix A = Q∗T∗Q∗∗T, so that the diagonal block of T with row index IFST is moved to row ILST
STRRFS, strrfsprovide error bounds and backward error estimates for the solution to a system of linear equations with a triangular coefficient matrix
STRSEN, strsenreorder the real Schur factorization of a real matrix A = Q∗T∗Q∗∗T, so that a selected cluster of eigenvalues appears in the leading diagonal blocks of the upper quasi-triangular matrix T,
STRSNA, strsnaestimate reciprocal condition numbers for specified eigenvalues and/or right eigenvectors of a real upper quasi-triangular matrix T (or of any matrix Q∗T∗Q∗∗T with Q orthogonal)
STRSYL, strsylsolve the real Sylvester matrix equation
STRTI2, strti2compute the inverse of a real upper or lower triangular matrix
STRTRI, strtricompute the inverse of a real upper or lower triangular matrix A
STRTRS, strtrssolve a triangular system of the form   A ∗ X = B or A∗∗T ∗ X = B,
STZRQF, stzrqfreduce the M-by-N ( M<=N ) real upper trapezoidal matrix A to upper triangular form by means of orthogonal transformations
XERBLA, xerblai an error handler for the LAPACK routines
ZBDSQR, zbdsqrcompute the singular value decomposition (SVD) of a real N-by-N (upper or lower) bidiagonal matrix B
ZDRSCL, zdrsclmultiplie an n-element complex vector x by the real scalar 1/a
ZGBBRD, zgbbrdreduce a complex general m-by-n band matrix A to real upper bidiagonal form B by a unitary transformation
ZGBCON, zgbconestimate the reciprocal of the condition number of a complex general band matrix A, in either the 1-norm or the infinity-norm,
ZGBEQU, zgbequcompute row and column scalings intended to equilibrate an M-by-N band matrix A and reduce its condition number
ZGBRFS, zgbrfsimprove 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
ZGBSV, zgbsvcompute the solution to a complex system of linear equations A ∗ X = B, where A is a band matrix of order N with KL subdiagonals and KU superdiagonals, and X and B are N-by-NRHS matrices
ZGBSVX, zgbsvxuse the LU factorization to compute the solution to a complex system of linear equations A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B,
ZGBTF2, zgbtf2compute an LU factorization of a complex m-by-n band matrix A using partial pivoting with row interchanges
ZGBTRF, zgbtrfcompute an LU factorization of a complex m-by-n band matrix A using partial pivoting with row interchanges
ZGBTRS, zgbtrssolve a system of linear equations  A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B with a general band matrix A using the LU factorization computed by ZGBTRF
ZGEBAK, zgebakform the right or left eigenvectors of a complex general matrix by backward transformation on the computed eigenvectors of the balanced matrix output by ZGEBAL
ZGEBAL, zgebalbalance a general complex matrix A
ZGEBD2, zgebd2reduce a complex general m by n matrix A to upper or lower real bidiagonal form B by a unitary transformation
ZGEBRD, zgebrdreduce a general complex M-by-N matrix A to upper or lower bidiagonal form B by a unitary transformation
ZGECON, zgeconestimate the reciprocal of the condition number of a general complex matrix A, in either the 1-norm or the infinity-norm, using the LU factorization computed by ZGETRF
ZGEEQU, zgeequcompute row and column scalings intended to equilibrate an M-by-N matrix A and reduce its condition number
ZGEES, zgeescompute for an N-by-N complex nonsymmetric matrix A, the eigenvalues, the Schur form T, and, optionally, the matrix of Schur vectors Z
ZGEESX, zgeesxcompute for an N-by-N complex nonsymmetric matrix A, the eigenvalues, the Schur form T, and, optionally, the matrix of Schur vectors Z
ZGEEV, zgeevcompute for an N-by-N complex nonsymmetric matrix A, the eigenvalues and, optionally, the left and/or right eigenvectors
ZGEEVX, zgeevxcompute for an N-by-N complex nonsymmetric matrix A, the eigenvalues and, optionally, the left and/or right eigenvectors
ZGEGS, zgegscompute for a pair of N-by-N complex nonsymmetric matrices A,
ZGEGV, zgegvcompute for a pair of N-by-N complex nonsymmetric matrices A and B, the generalized eigenvalues (alpha, beta), and optionally,
ZGEHD2, zgehd2reduce a complex general matrix A to upper Hessenberg form H by a unitary similarity transformation
ZGEHRD, zgehrdreduce a complex general matrix A to upper Hessenberg form H by a unitary similarity transformation
ZGELQ2, zgelq2compute an LQ factorization of a complex m by n matrix A
ZGELQF, zgelqfcompute an LQ factorization of a complex M-by-N matrix A
ZGELS, zgelssolve overdetermined or underdetermined complex linear systems involving an M-by-N matrix A, or its conjugate-transpose, using a QR or LQ factorization of A
ZGELSS, zgelsscompute the minimum norm solution to a complex linear least squares problem
ZGELSX, zgelsxcompute the minimum-norm solution to a complex linear least squares problem
ZGEQL2, zgeql2compute a QL factorization of a complex m by n matrix A
ZGEQLF, zgeqlfcompute a QL factorization of a complex M-by-N matrix A
ZGEQPF, zgeqpfcompute a QR factorization with column pivoting of a complex M-by-N matrix A
ZGEQR2, zgeqr2compute a QR factorization of a complex m by n matrix A
ZGEQRF, zgeqrfcompute a QR factorization of a complex M-by-N matrix A
ZGERFS, zgerfsimprove the computed solution to a system of linear equations and provides error bounds and backward error estimates for the solution
ZGERQ2, zgerq2compute an RQ factorization of a complex m by n matrix A
ZGERQF, zgerqfcompute an RQ factorization of a complex M-by-N matrix A
ZGESV, zgesvcompute the solution to a complex system of linear equations  A ∗ X = B,
ZGESVD, zgesvdcompute the singular value decomposition (SVD) of a complex M-by-N matrix A, optionally computing the left and/or right singular vectors
ZGESVX, zgesvxuse the LU factorization to compute the solution to a complex system of linear equations  A ∗ X = B,
ZGETF2, zgetf2compute an LU factorization of a general m-by-n matrix A using partial pivoting with row interchanges
ZGETRF, zgetrfcompute an LU factorization of a general M-by-N matrix A using partial pivoting with row interchanges
ZGETRI, zgetricompute the inverse of a matrix using the LU factorization computed by ZGETRF
ZGETRS, zgetrssolve a system of linear equations  A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B with a general N-by-N matrix A using the LU factorization computed by ZGETRF
ZGGBAK, zggbakform the right or left eigenvectors of a complex generalized eigenvalue problem A∗x = lambda∗B∗x, by backward transformation on the computed eigenvectors of the balanced pair of matrices output by ZGGBAL
ZGGBAL, zggbalbalance a pair of general complex matrices (A,B)
ZGGGLM, zggglmsolve a general Gauss-Markov linear model (GLM) problem
ZGGHRD, zgghrdreduce a pair of complex matrices (A,B) to generalized upper Hessenberg form using unitary transformations, where A is a general matrix and B is upper triangular
ZGGLSE, zgglsesolve the linear equality-constrained least squares (LSE) problem
ZGGQRF, zggqrfcompute a generalized QR factorization of an N-by-M matrix A and an N-by-P matrix B
ZGGRQF, zggrqfcompute a generalized RQ factorization of an M-by-N matrix A and a P-by-N matrix B
ZGGSVD, zggsvdcompute the generalized singular value decomposition (GSVD) of an M-by-N complex matrix A and P-by-N complex matrix B
ZGGSVP, zggsvpcompute unitary matrices U, V and Q such that   N-K-L K L  U’∗A∗Q = K ( 0 A12 A13 ) if M-K-L >= 0
ZGTCON, zgtconestimate the reciprocal of the condition number of a complex tridiagonal matrix A using the LU factorization as computed by ZGTTRF
ZGTRFS, zgtrfsimprove the computed solution to a system of linear equations when the coefficient matrix is tridiagonal, and provides error bounds and backward error estimates for the solution
ZGTSV, zgtsvsolve the equation   A∗X = B,
ZGTSVX, zgtsvxuse the LU factorization to compute the solution to a complex system of linear equations A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B,
ZGTTRF, zgttrfcompute an LU factorization of a complex tridiagonal matrix A using elimination with partial pivoting and row interchanges
ZGTTRS, zgttrssolve one of the systems of equations  A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B,
ZHBEV, zhbevcompute all the eigenvalues and, optionally, eigenvectors of a complex Hermitian band matrix A
ZHBEVD, zhbevdcompute all the eigenvalues and, optionally, eigenvectors of a complex Hermitian band matrix A
ZHBEVX, zhbevxcompute selected eigenvalues and, optionally, eigenvectors of a complex Hermitian band matrix A
ZHBGST, zhbgstreduce a complex Hermitian-definite banded generalized eigenproblem A∗x = lambda∗B∗x to standard form C∗y = lambda∗y,
ZHBGV, zhbgvcompute all the eigenvalues, and optionally, the eigenvectors of a complex generalized Hermitian-definite banded eigenproblem, of the form A∗x=(lambda)∗B∗x
ZHBTRD, zhbtrdreduce a complex Hermitian band matrix A to real symmetric tridiagonal form T by a unitary similarity transformation
ZHECON, zheconestimate the reciprocal of the condition number of a complex Hermitian matrix A using the factorization A = U∗D∗U∗∗H or A = L∗D∗L∗∗H computed by ZHETRF
ZHEEV, zheevcompute all eigenvalues and, optionally, eigenvectors of a complex Hermitian matrix A
ZHEEVD, zheevdcompute all eigenvalues and, optionally, eigenvectors of a complex Hermitian matrix A
ZHEEVX, zheevxcompute selected eigenvalues and, optionally, eigenvectors of a complex Hermitian matrix A
ZHEGS2, zhegs2reduce a complex Hermitian-definite generalized eigenproblem to standard form
ZHEGST, zhegstreduce a complex Hermitian-definite generalized eigenproblem to standard form
ZHEGV, zhegvcompute all the eigenvalues, and optionally, the eigenvectors of a complex generalized Hermitian-definite eigenproblem, of the form A∗x=(lambda)∗B∗x, A∗Bx=(lambda)∗x, or B∗A∗x=(lambda)∗x
ZHERFS, zherfsimprove the computed solution to a system of linear equations when the coefficient matrix is Hermitian indefinite, and provides error bounds and backward error estimates for the solution
ZHESV, zhesvcompute the solution to a complex system of linear equations  A ∗ X = B,
ZHESVX, zhesvxuse the diagonal pivoting factorization to compute the solution to a complex system of linear equations A ∗ X = B,
ZHETD2, zhetd2reduce a complex Hermitian matrix A to real symmetric tridiagonal form T by a unitary similarity transformation
ZHETF2, zhetf2compute the factorization of a complex Hermitian matrix A using the Bunch-Kaufman diagonal pivoting method
ZHETRD, zhetrdreduce a complex Hermitian matrix A to real symmetric tridiagonal form T by a unitary similarity transformation
ZHETRF, zhetrfcompute the factorization of a complex Hermitian matrix A using the Bunch-Kaufman diagonal pivoting method
ZHETRI, zhetricompute the inverse of a complex Hermitian indefinite matrix A using the factorization A = U∗D∗U∗∗H or A = L∗D∗L∗∗H computed by ZHETRF
ZHETRS, zhetrssolve a system of linear equations A∗X = B with a complex Hermitian matrix A using the factorization A = U∗D∗U∗∗H or A = L∗D∗L∗∗H computed by ZHETRF
ZHGEQZ, zhgeqzimplement a single-shift version of the QZ method for finding the generalized eigenvalues w(i)=ALPHA(i)/BETA(i) of the equation   det( A - w(i) B ) = 0  If JOB=’S’, then the pair (A,B) is simultaneously reduced to Schur form (i.e., A and B are both upper triangular) by applying one unitary tranformation (usually called Q) on the left and another (usually called Z) on the right
ZHPCON, zhpconestimate the reciprocal of the condition number of a complex Hermitian packed matrix A using the factorization A = U∗D∗U∗∗H or A = L∗D∗L∗∗H computed by ZHPTRF
ZHPEV, zhpevcompute all the eigenvalues and, optionally, eigenvectors of a complex Hermitian matrix in packed storage
ZHPEVD, zhpevdcompute all the eigenvalues and, optionally, eigenvectors of a complex Hermitian matrix A in packed storage
ZHPEVX, zhpevxcompute selected eigenvalues and, optionally, eigenvectors of a complex Hermitian matrix A in packed storage
ZHPGST, zhpgstreduce a complex Hermitian-definite generalized eigenproblem to standard form, using packed storage
ZHPGV, zhpgvcompute all the eigenvalues and, optionally, the eigenvectors of a complex generalized Hermitian-definite eigenproblem, of the form A∗x=(lambda)∗B∗x, A∗Bx=(lambda)∗x, or B∗A∗x=(lambda)∗x
ZHPRFS, zhprfsimprove the computed solution to a system of linear equations when the coefficient matrix is Hermitian indefinite and packed, and provides error bounds and backward error estimates for the solution
ZHPSV, zhpsvcompute the solution to a complex system of linear equations  A ∗ X = B,
ZHPSVX, zhpsvxuse the diagonal pivoting factorization A = U∗D∗U∗∗H or A = L∗D∗L∗∗H to compute the solution to a complex system of linear equations A ∗ X = B, where A is an N-by-N Hermitian matrix stored in packed format and X and B are N-by-NRHS matrices
ZHPTRD, zhptrdreduce a complex Hermitian matrix A stored in packed form to real symmetric tridiagonal form T by a unitary similarity transformation
ZHPTRF, zhptrfcompute the factorization of a complex Hermitian packed matrix A using the Bunch-Kaufman diagonal pivoting method
ZHPTRI, zhptricompute the inverse of a complex Hermitian indefinite matrix A in packed storage using the factorization A = U∗D∗U∗∗H or A = L∗D∗L∗∗H computed by ZHPTRF
ZHPTRS, zhptrssolve a system of linear equations A∗X = B with a complex Hermitian matrix A stored in packed format using the factorization A = U∗D∗U∗∗H or A = L∗D∗L∗∗H computed by ZHPTRF
ZHSEIN, zhseinuse inverse iteration to find specified right and/or left eigenvectors of a complex upper Hessenberg matrix H
ZHSEQR, zhseqrcompute the eigenvalues of a complex upper Hessenberg matrix H, and, optionally, the matrices T and Z from the Schur decomposition H = Z T Z∗∗H, where T is an upper triangular matrix (the Schur form), and Z is the unitary matrix of Schur vectors
ZLABRD, zlabrdreduce the first NB rows and columns of a complex general m by n matrix A to upper or lower real bidiagonal form by a unitary transformation Q’ ∗ A ∗ P, and returns the matrices X and Y which are needed to apply the transformation to the unreduced part of A
ZLACGV, zlacgvconjugate a complex vector of length N
ZLACON, zlaconestimate the 1-norm of a square, complex matrix A
ZLACPY, zlacpycopie all or part of a two-dimensional matrix A to another matrix B
ZLACRM, zlacrmperform a very simple matrix-matrix multiplication
ZLACRT, zlacrtapplie a plane rotation, where the cos and sin (C and S) are complex and the vectors CX and CY are complex
ZLADIV, zladiv:= X / Y, where X and Y are complex
ZLAED0, zlaed0the divide and conquer method, ZLAED0 computes all eigenvalues of a symmetric tridiagonal matrix which is one diagonal block of those from reducing a dense or band Hermitian matrix and corresponding eigenvectors of the dense or band matrix
ZLAED7, zlaed7compute the updated eigensystem of a diagonal matrix after modification by a rank-one symmetric matrix
ZLAED8, zlaed8merge the two sets of eigenvalues together into a single sorted set
ZLAEIN, zlaeinuse inverse iteration to find a right or left eigenvector corresponding to the eigenvalue W of a complex upper Hessenberg matrix H
ZLAESY, zlaesycompute the eigendecomposition of a 2-by-2 symmetric matrix  ( ( A, B );( B, C ) ) provided the norm of the matrix of eigenvectors is larger than some threshold value
ZLAEV2, zlaev2compute the eigendecomposition of a 2-by-2 Hermitian matrix  [ A B ]  [ CONJG(B) C ]
zlags2
ZLAGTM, zlagtmperform a matrix-vector product of the form   B := alpha ∗ A ∗ X + beta ∗ B  where A is a tridiagonal matrix of order N, B and X are N by NRHS matrices, and alpha and beta are real scalars, each of which may be 0., 1., or -1
ZLAHEF, zlahefcompute a partial factorization of a complex Hermitian matrix A using the Bunch-Kaufman diagonal pivoting method
ZLAHQR, zlahqri an auxiliary routine called by ZHSEQR to update the eigenvalues and Schur decomposition already computed by ZHSEQR, by dealing with the Hessenberg submatrix in rows and columns ILO to IHI
ZLAHRD, zlahrdreduce the first NB columns of a complex general n-by-(n-k+1) matrix A so that elements below the k-th subdiagonal are zero
ZLAIC1, zlaic1applie one step of incremental condition estimation in its simplest version
ZLANGB, zlangbreturn 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
ZLANGE, zlangereturn the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex matrix A
ZLANGT, zlangtreturn the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex tridiagonal matrix A
ZLANHB, zlanhbreturn 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 hermitian band matrix A, with k super-diagonals
ZLANHE, zlanhereturn the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex hermitian matrix A
ZLANHP, zlanhpreturn the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex hermitian matrix A, supplied in packed form
ZLANHS, zlanhsreturn 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
ZLANHT, zlanhtreturn the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex Hermitian tridiagonal matrix A
ZLANSB, zlansbreturn 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
ZLANSP, zlanspreturn the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex symmetric matrix A, supplied in packed form
ZLANSY, zlansyreturn the value of the one norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex symmetric matrix A
ZLANTB, zlantbreturn 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
ZLANTP, zlantpreturn 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
ZLANTR, zlantrreturn 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
ZLAPLL, zlaplltwo column vectors X and Y, let   A = ( X Y )
ZLAPMT, zlapmtrearrange 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
ZLAQGB, zlaqgbequilibrate 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
ZLAQGE, zlaqgeequilibrate a general M by N matrix A using the row and scaling factors in the vectors R and C
ZLAQHB, zlaqhbequilibrate a symmetric band matrix A using the scaling factors in the vector S
ZLAQHE, zlaqheequilibrate a Hermitian matrix A using the scaling factors in the vector S
ZLAQHP, zlaqhpequilibrate a Hermitian matrix A using the scaling factors in the vector S
ZLAQSB, zlaqsbequilibrate a symmetric band matrix A using the scaling factors in the vector S
ZLAQSP, zlaqspequilibrate a symmetric matrix A using the scaling factors in the vector S
ZLAQSY, zlaqsyequilibrate a symmetric matrix A using the scaling factors in the vector S
ZLAR2V, zlar2vapplie a vector of complex plane rotations with real cosines from both sides to a sequence of 2-by-2 complex Hermitian matrices,
ZLARF, zlarfapplie a complex elementary reflector H to a complex M-by-N matrix C, from either the left or the right
ZLARFB, zlarfbapplie a complex block reflector H or its transpose H’ to a complex M-by-N matrix C, from either the left or the right
ZLARFG, zlarfggenerate a complex elementary reflector H of order n, such that   H’ ∗ ( alpha ) = ( beta ), H’ ∗ H = I
ZLARFT, zlarftform the triangular factor T of a complex block reflector H of order n, which is defined as a product of k elementary reflectors
ZLARFX, zlarfxapplie a complex elementary reflector H to a complex m by n matrix C, from either the left or the right
ZLARGV, zlargvgenerate a vector of complex plane rotations with real cosines, determined by elements of the complex vectors x and y
ZLARNV, zlarnvreturn a vector of n random complex numbers from a uniform or normal distribution
ZLARTG, zlartggenerate a plane rotation so that   [ CS SN ] [ F ] [ R ]  [ __ ]
ZLARTV, zlartvapplie a vector of complex plane rotations with real cosines to elements of the complex vectors x and y
ZLASCL, zlasclmultiplie the M by N complex matrix A by the real scalar CTO/CFROM
ZLASET, zlasetinitialize a 2-D array A to BETA on the diagonal and ALPHA on the offdiagonals
ZLASR, zlasrperform the transformation   A := P∗A, when SIDE = ’L’ or ’l’ ( Left-hand side )   A := A∗P’, when SIDE = ’R’ or ’r’ ( Right-hand side )  where A is an m by n complex matrix and P is an orthogonal matrix,
ZLASSQ, zlassqreturn the values scl and ssq such that   ( scl∗∗2 )∗ssq = x( 1 )∗∗2 +...+ x( n )∗∗2 + ( scale∗∗2 )∗sumsq,
ZLASWP, zlaswpperform a series of row interchanges on the matrix A
ZLASYF, zlasyfcompute a partial factorization of a complex symmetric matrix A using the Bunch-Kaufman diagonal pivoting method
ZLATBS, zlatbssolve one of the triangular systems   A ∗ x = s∗b, A∗∗T ∗ x = s∗b, or A∗∗H ∗ x = s∗b,
ZLATPS, zlatpssolve one of the triangular systems   A ∗ x = s∗b, A∗∗T ∗ x = s∗b, or A∗∗H ∗ x = s∗b,
ZLATRD, zlatrdreduce NB rows and columns of a complex Hermitian matrix A to Hermitian tridiagonal form by a unitary similarity transformation Q’ ∗ A ∗ Q, and returns the matrices V and W which are needed to apply the transformation to the unreduced part of A
ZLATRS, zlatrssolve one of the triangular systems   A ∗ x = s∗b, A∗∗T ∗ x = s∗b, or A∗∗H ∗ x = s∗b,
ZLATZM, zlatzmapplie a Householder matrix generated by ZTZRQF to a matrix
ZLAUU2, zlauu2compute the product U ∗ U’ or L’ ∗ L, where the triangular factor U or L is stored in the upper or lower triangular part of the array A
ZLAUUM, zlauumcompute the product U ∗ U’ or L’ ∗ L, where the triangular factor U or L is stored in the upper or lower triangular part of the array A
ZLAZRO, zlazroinitialize a 2-D array A to BETA on the diagonal and ALPHA on the offdiagonals
ZPBCON, zpbconestimate the reciprocal of the condition number (in the 1-norm) of a complex Hermitian positive definite band matrix using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by ZPBTRF
ZPBEQU, zpbequcompute row and column scalings intended to equilibrate a Hermitian positive definite band matrix A and reduce its condition number (with respect to the two-norm)
ZPBRFS, zpbrfsimprove the computed solution to a system of linear equations when the coefficient matrix is Hermitian positive definite and banded, and provides error bounds and backward error estimates for the solution
ZPBSTF, zpbstfcompute a split Cholesky factorization of a complex Hermitian positive definite band matrix A
ZPBSV, zpbsvcompute the solution to a complex system of linear equations  A ∗ X = B,
ZPBSVX, zpbsvxuse the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H to compute the solution to a complex system of linear equations  A ∗ X = B,
ZPBTF2, zpbtf2compute the Cholesky factorization of a complex Hermitian positive definite band matrix A
ZPBTRF, zpbtrfcompute the Cholesky factorization of a complex Hermitian positive definite band matrix A
ZPBTRS, zpbtrssolve a system of linear equations A∗X = B with a Hermitian positive definite band matrix A using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by ZPBTRF
ZPOCON, zpoconestimate the reciprocal of the condition number (in the 1-norm) of a complex Hermitian positive definite matrix using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by ZPOTRF
ZPOEQU, zpoequcompute row and column scalings intended to equilibrate a Hermitian positive definite matrix A and reduce its condition number (with respect to the two-norm)
ZPORFS, zporfsimprove the computed solution to a system of linear equations when the coefficient matrix is Hermitian positive definite,
ZPOSV, zposvcompute the solution to a complex system of linear equations  A ∗ X = B,
ZPOSVX, zposvxuse the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H to compute the solution to a complex system of linear equations  A ∗ X = B,
ZPOTF2, zpotf2compute the Cholesky factorization of a complex Hermitian positive definite matrix A
ZPOTRF, zpotrfcompute the Cholesky factorization of a complex Hermitian positive definite matrix A
ZPOTRI, zpotricompute the inverse of a complex Hermitian positive definite matrix A using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by ZPOTRF
ZPOTRS, zpotrssolve a system of linear equations A∗X = B with a Hermitian positive definite matrix A using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by ZPOTRF
ZPPCON, zppconestimate the reciprocal of the condition number (in the 1-norm) of a complex Hermitian positive definite packed matrix using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by ZPPTRF
ZPPEQU, zppequcompute row and column scalings intended to equilibrate a Hermitian positive definite matrix A in packed storage and reduce its condition number (with respect to the two-norm)
ZPPRFS, zpprfsimprove the computed solution to a system of linear equations when the coefficient matrix is Hermitian positive definite and packed, and provides error bounds and backward error estimates for the solution
ZPPSV, zppsvcompute the solution to a complex system of linear equations  A ∗ X = B,
ZPPSVX, zppsvxuse the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H to compute the solution to a complex system of linear equations  A ∗ X = B,
ZPPTRF, zpptrfcompute the Cholesky factorization of a complex Hermitian positive definite matrix A stored in packed format
ZPPTRI, zpptricompute the inverse of a complex Hermitian positive definite matrix A using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by ZPPTRF
ZPPTRS, zpptrssolve a system of linear equations A∗X = B with a Hermitian positive definite matrix A in packed storage using the Cholesky factorization A = U∗∗H∗U or A = L∗L∗∗H computed by ZPPTRF
ZPTCON, zptconcompute the reciprocal of the condition number (in the 1-norm) of a complex Hermitian positive definite tridiagonal matrix using the factorization A = L∗D∗L∗∗H or A = U∗∗H∗D∗U computed by ZPTTRF
ZPTEQR, zpteqrcompute all eigenvalues and, optionally, eigenvectors of a symmetric positive definite tridiagonal matrix by first factoring the matrix using DPTTRF and then calling ZBDSQR to compute the singular values of the bidiagonal factor
ZPTRFS, zptrfsimprove the computed solution to a system of linear equations when the coefficient matrix is Hermitian positive definite and tridiagonal, and provides error bounds and backward error estimates for the solution
ZPTSV, zptsvcompute the solution to a complex system of linear equations A∗X = B, where A is an N-by-N Hermitian positive definite tridiagonal matrix, and X and B are N-by-NRHS matrices
ZPTSVX, zptsvxuse the factorization A = L∗D∗L∗∗H to compute the solution to a complex system of linear equations A∗X = B, where A is an N-by-N Hermitian positive definite tridiagonal matrix and X and B are N-by-NRHS matrices
ZPTTRF, zpttrfcompute the factorization of a complex Hermitian positive definite tridiagonal matrix A
ZPTTRS, zpttrssolve a system of linear equations A ∗ X = B with a Hermitian positive definite tridiagonal matrix A using the factorization A = U∗∗H∗D∗U or A = L∗D∗L∗∗H computed by ZPTTRF
ZROT, zrotapplie a plane rotation, where the cos (C) is real and the sin (S) is complex, and the vectors CX and CY are complex
ZSPCON, zspconestimate the reciprocal of the condition number (in the 1-norm) of a complex symmetric packed matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by ZSPTRF
ZSPMV, zspmvperform the matrix-vector operation   y := alpha∗A∗x + beta∗y,
ZSPR, zsprperform the symmetric rank 1 operation   A := alpha∗x∗conjg( x’ ) + A,
ZSPRFS, zsprfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric indefinite and packed, and provides error bounds and backward error estimates for the solution
ZSPSV, zspsvcompute the solution to a complex system of linear equations  A ∗ X = B,
ZSPSVX, zspsvxuse the diagonal pivoting factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T to compute the solution to a complex system of linear equations A ∗ X = B, where A is an N-by-N symmetric matrix stored in packed format and X and B are N-by-NRHS matrices
ZSPTRF, zsptrfcompute the factorization of a complex symmetric matrix A stored in packed format using the Bunch-Kaufman diagonal pivoting method
ZSPTRI, zsptricompute the inverse of a complex symmetric indefinite matrix A in packed storage using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by ZSPTRF
ZSPTRS, zsptrssolve a system of linear equations A∗X = B with a complex symmetric matrix A stored in packed format using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by ZSPTRF
ZSRSCL, zsrsclmultiplie an n-element complex vector x by the real scalar 1/a
ZSTEDC, zstedccompute all eigenvalues and, optionally, eigenvectors of a symmetric tridiagonal matrix using the divide and conquer method
ZSTEIN, zsteincompute the eigenvectors of a real symmetric tridiagonal matrix T corresponding to specified eigenvalues, using inverse iteration
ZSTEQR, zsteqrcompute all eigenvalues and, optionally, eigenvectors of a symmetric tridiagonal matrix using the implicit QL or QR method
ZSYCON, zsyconestimate the reciprocal of the condition number (in the 1-norm) of a complex symmetric matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by ZSYTRF
ZSYMV, zsymvperform the matrix-vector operation   y := alpha∗A∗x + beta∗y,
ZSYR, zsyrperform the symmetric rank 1 operation   A := alpha∗x∗( x’ ) + A,
ZSYRFS, zsyrfsimprove the computed solution to a system of linear equations when the coefficient matrix is symmetric indefinite, and provides error bounds and backward error estimates for the solution
ZSYSV, zsysvcompute the solution to a complex system of linear equations  A ∗ X = B,
ZSYSVX, zsysvxuse the diagonal pivoting factorization to compute the solution to a complex system of linear equations A ∗ X = B,
ZSYTF2, zsytf2compute the factorization of a complex symmetric matrix A using the Bunch-Kaufman diagonal pivoting method
ZSYTRF, zsytrfcompute the factorization of a complex symmetric matrix A using the Bunch-Kaufman diagonal pivoting method
ZSYTRI, zsytricompute the inverse of a complex symmetric indefinite matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by ZSYTRF
ZSYTRS, zsytrssolve a system of linear equations A∗X = B with a complex symmetric matrix A using the factorization A = U∗D∗U∗∗T or A = L∗D∗L∗∗T computed by ZSYTRF
ZTBCON, ztbconestimate the reciprocal of the condition number of a triangular band matrix A, in either the 1-norm or the infinity-norm
ZTBRFS, ztbrfsprovide error bounds and backward error estimates for the solution to a system of linear equations with a triangular band coefficient matrix
ZTBTRS, ztbtrssolve a triangular system of the form   A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B,
ZTGEVC, ztgevccompute some or all of the right and/or left generalized eigenvectors of a pair of complex upper triangular matrices (A,B)
ZTGSJA, ztgsjacompute the generalized singular value decomposition (GSVD) of two complex upper triangular (or trapezoidal) matrices A and B
ZTPCON, ztpconestimate the reciprocal of the condition number of a packed triangular matrix A, in either the 1-norm or the infinity-norm
ZTPRFS, ztprfsprovide error bounds and backward error estimates for the solution to a system of linear equations with a triangular packed coefficient matrix
ZTPTRI, ztptricompute the inverse of a complex upper or lower triangular matrix A stored in packed format
ZTPTRS, ztptrssolve a triangular system of the form   A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B,
ZTRCON, ztrconestimate the reciprocal of the condition number of a triangular matrix A, in either the 1-norm or the infinity-norm
ZTREVC, ztrevccompute some or all of the right and/or left eigenvectors of a complex upper triangular matrix T
ZTREXC, ztrexcreorder the Schur factorization of a complex matrix A = Q∗T∗Q∗∗H, so that the diagonal element of T with row index IFST is moved to row ILST
ZTRRFS, ztrrfsprovide error bounds and backward error estimates for the solution to a system of linear equations with a triangular coefficient matrix
ZTRSEN, ztrsenreorder the Schur factorization of a complex matrix A = Q∗T∗Q∗∗H, so that a selected cluster of eigenvalues appears in the leading positions on the diagonal of the upper triangular matrix T, and the leading columns of Q form an orthonormal basis of the corresponding right invariant subspace
ZTRSNA, ztrsnaestimate reciprocal condition numbers for specified eigenvalues and/or right eigenvectors of a complex upper triangular matrix T (or of any matrix Q∗T∗Q∗∗H with Q unitary)
ZTRSYL, ztrsylsolve the complex Sylvester matrix equation
ZTRTI2, ztrti2compute the inverse of a complex upper or lower triangular matrix
ZTRTRI, ztrtricompute the inverse of a complex upper or lower triangular matrix A
ZTRTRS, ztrtrssolve a triangular system of the form   A ∗ X = B, A∗∗T ∗ X = B, or A∗∗H ∗ X = B,
ZTZRQF, ztzrqfreduce the M-by-N ( M<=N ) complex upper trapezoidal matrix A to upper triangular form by means of unitary transformations
ZUNG2L, zung2lgenerate an m by n complex matrix Q with orthonormal columns,
ZUNG2R, zung2rgenerate an m by n complex matrix Q with orthonormal columns,
ZUNGBR, zungbrgenerate one of the complex unitary matrices Q or P∗∗H determined by ZGEBRD when reducing a complex matrix A to bidiagonal form
ZUNGHR, zunghrgenerate a complex unitary matrix Q which is defined as the product of IHI-ILO elementary reflectors of order N, as returned by ZGEHRD
ZUNGL2, zungl2generate an m-by-n complex matrix Q with orthonormal rows,
ZUNGLQ, zunglqgenerate an M-by-N complex matrix Q with orthonormal rows,
ZUNGQL, zungqlgenerate an M-by-N complex matrix Q with orthonormal columns,
ZUNGQR, zungqrgenerate an M-by-N complex matrix Q with orthonormal columns,
ZUNGR2, zungr2generate an m by n complex matrix Q with orthonormal rows,
ZUNGRQ, zungrqgenerate an M-by-N complex matrix Q with orthonormal rows,
ZUNGTR, zungtrgenerate a complex unitary matrix Q which is defined as the product of n-1 elementary reflectors of order N, as returned by ZHETRD
ZUNM2L, zunm2loverwrite the general complex m-by-n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’C’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’C’,
ZUNM2R, zunm2roverwrite the general complex m-by-n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’C’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’C’,
ZUNMBR, zunmbrVECT = ’Q’, ZUNMBR overwrites the general complex M-by-N matrix C with  SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
ZUNMHR, zunmhroverwrite the general complex M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
ZUNML2, zunml2overwrite the general complex m-by-n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’C’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’C’,
ZUNMLQ, zunmlqoverwrite the general complex M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
ZUNMQL, zunmqloverwrite the general complex M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
ZUNMQR, zunmqroverwrite the general complex M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
ZUNMR2, zunmr2overwrite the general complex m-by-n matrix C with   Q ∗ C if SIDE = ’L’ and TRANS = ’N’, or   Q’∗ C if SIDE = ’L’ and TRANS = ’C’, or   C ∗ Q if SIDE = ’R’ and TRANS = ’N’, or   C ∗ Q’ if SIDE = ’R’ and TRANS = ’C’,
ZUNMRQ, zunmrqoverwrite the general complex M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
ZUNMTR, zunmtroverwrite the general complex M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’
ZUPGTR, zupgtrgenerate a complex unitary matrix Q which is defined as the product of n-1 elementary reflectors H(i) of order n, as returned by ZHPTRD using packed storage
ZUPMTR, zupmtroverwrite the general complex M-by-N matrix C with   SIDE = ’L’ SIDE = ’R’ TRANS = ’N’

Typewritten Software • bear@typewritten.org • Edmonds, WA 98026