XLPack function list B: Basic functions M1: Linear computation(real), M2: Linear computation(complex), M3: Special functions, nonlinear computation,
(July 27, 2025) M4: Interpolation, differential/integral equations, random numbers, M5: Sparse matrix calculation.
Subject to change without notice B can be used without license. In-app purchase of addon licenses are required to use M1 to M5.  @ does not work on 32 bit Excel.
  Deprecated (to be deleted in the next version).
(*) shows the experimental version (may be changed in the future version). (#) shows that equivalent subroutines are also provided in version  7.0.10 or later. 
Category   VBA routine name (WS function name) Functions V6.1 V7.0
VBA WS Solver VBA WS Solver
A3. Real arithmetic A3. Real arithmetic D1num (#) IEEE 754 special numbers (double precision) B     B    
IsFinite (#) Determines if finite value (double precision) B     B    
IsInf (#) Determines if infinite value (double precision) B     B    
IsNan (#) Determines if NaN (not a number) (double precision) B     B    
IsNormal (#) Determines if nomal value (double precision) B     B    
Signbit (#) Determines if negative value (double precision) B     B    
A4. Complex arithmetic A4. Complex arithmetic Creal Real part of complex number B     B    
Cimag Imaginary part of complex number B     B    
Cabs Absolute value of complex number B     B    
Conj Conjugate number B     B    
Carg Argument of complex number B     B    
Cproj Projection of complex number on Riemann sphere B     B    
Cmplx Building complex number B     B    
Cpolar Building complex number (polar coordinate) B     B    
Cminus Sign inversion of complex number B     B    
Cadd Addition of complex numbers B     B    
Cadd3 Addition of three complex numbers       B    
Cadd4 Addition of four complex numbers       B    
Cadd5 Addition of five complex numbers       B    
Cdadd Addition of complex number and real number B     B    
Dcadd Addition of real number and complex number       B    
Csub Subtraction of complex number from complex number B     B    
Cdsub Subtraction of real number from complex number B     B    
Dcsub Subtraction of complex number from real number B     B    
Cmul Multiplication of complex numbers B     B    
Cmul3 Multiplication of three complex numbers       B    
Cmul4 Multiplication of four complex numbers       B    
Cmul5 Multiplication of five complex numbers       B    
Cdmul Multiplication of complex number and real number B     B    
Dcmul Multiplication of real number and complex number       B    
Cdiv Division of complex number by complex number B     B    
Cddiv Division of complex number by real number B     B    
Dcdiv Division of real number by complex number B     B    
Cpow Power of complex number B     B    
Cdpow Power of a complex number (real order) B     B    
Cipow Power of a complex number (integer order) B     B    
Ceq Comparison of complex numbers       B    
Cneq Negative comparison of complex numbers       B    
C. Elementary and special functions C1. Integer-valued functions Factorial (#) Factorial B   B  
C2. Powers, roots, reciprocals Fma (#) (WFma) (x*y)+z B B   B B  
Hypot (#) (WHypot) sqrt(x^2+y^2) B B   B B  
Cbrt (#) (WCbrt) Cube root B B   B B  
Csqrt Complex square root B     B    
Ccbrt Complex cube root B     B    
C3. Polynomials Laguerre (#) (WLaguerre) Laguerre polynomial Ln(x) M3     M3 M3  
Alaguerre (#) (WAlaguerre) Associated Laguerre polynomial Lnm(x) M3     M3 M3  
Legendre (#) (WLegendre) Legendre polynomial Pn(x) M3     M3 M3  
Legendred (#) (WLegendred) Derivative of Legendre polynomial Pn(x) M3     M3 M3  
Alegendre (#) (WAlegendre) Associated Legendre polynomial Pnm(x) M3     M3 M3  
Sharmonic (WSharmonic) Spherical harmonic Ylm(θ, φ) M3     M3 M3  
Sharmonicr (#) (WSharmonicr) Real part of spherical harmonic Ylm(θ, φ) M3     M3 M3  
Sharmonici (#) (WSharmonici) Imaginary part of spherical harmonic Ylm(θ, φ) M3     M3 M3  
Hermite (#) (WHermite) Hermite polynomial Hn(x) M3     M3 M3  
Chebt (#) (WChebt) Chebyshev polynomial of first kind Tn(x) M3     M3 M3  
Chebtd (#) (WChebtd) Derivative of Chebyshev polynomial of first kind Tn'(x) M3     M3 M3  
Chebu (#) (WChebu) Chebyshev polynomial of second kind Un(x) M3     M3 M3  
Chebs Evaluation of Chebyshev series B     B    
Gegenbauer (#) (WGegenbauer) Gegenbauer polynomial Cn(λ)(x)       M3 M3  
Gegenbauerd1 (#) (WGegenbauerd1) First derivative of Gegenbauer polynomial Cn(λ)(x)       M3 M3  
Gegenbauerd (#) (WGegenbauerd) K-th derivative of Gegenbauer polynomial Cn(λ)(x)       M3 M3  
Jacobi (#) (WJacobi) Jacobi polynomial Pn(α, β)(x)       M3 M3  
Jacobid1 (#) (WJacobid1) First derivative of Jacobi polynomial Pn(α, β)(x)       M3 M3  
Jacobid2 (#) (WJacobid2) Second derivative of Jacobi polynomial Pn(α, β)(x)       M3 M3  
Jacobid (#) (WJacobid) K-th derivative of Jacobi polynomial Pn(α, β)(x)       M3 M3  
C4. Elementary transcendental functions Expm1 (#) (WExpm1) exp(x)-1 B B   B B  
Exp2 (#) 2^x (base-2 exponent of x) B     B    
Log1p (#) (WLog1p) ln(1+x) B B   B B  
Log2 (#) log2(x) (base-2 logarithm of x) B     B    
Log10 log10(x) (base-10 logarithm of x) B     B    
Sqrt1pm1 (#) sqrt(1 + x) - 1 B     B    
Powm1 (#) x^y - 1 B     B    
Sinpi (#) sin(πx) B     B    
Cospi (#) cos(πx) B     B    
Acos arccos(x) B     B    
Asin arcsin(x) B     B    
Atan2 arctan2(y, x) B     B    
Cosh cosh(x) B     B    
Sinh sinh(x) B     B    
Tanh tanh(x) B     B    
Acosh arccosh(x) B     B    
Asinh arcsinh(x) B     B    
Atanh arctanh(x) B     B    
Cexp Complex exp(z) B     B    
Clog Complex ln(z) B     B    
Cexpm1 Complex exp(z)-1 B     B    
Clog1p Complex ln(1+z) B     B    
Ccos Complex cos(z) B     B    
Csin Complex sin(z) B     B    
Ctan Complex tan(z) B     B    
Cacos Complex arccos(z) B     B    
Casin Complex arcsin(z) B     B    
Catan Complex arctan(z) B     B    
Ccosh Complex cosh(z) B     B    
Csinh Complex sinh(z) B     B    
Ctanh Complex tanh(z) B     B    
Cacosh Complex arcosh(z) B     B    
Casinh Complex arsinh(z) B     B    
Catanh Complex artanh(z) B     B    
Ccot Complex cot(z) B     B    
C5. Exponential and logarithmic integrals Li (#) (WLi) Logarithmic integral li(x) B B   B B  
Ei (#) (WEi) Exponential integral Ei(x) B B   B B  
E1 (#) (WE_1) Exponential integral E1(x) B B   B B  
En (#) (WEn) Exponential integrals En(x) M3 M3   M3 M3  
Spence (#) (WSpence) Spence's function (dilogarithm function) Li2(x) M3 M3   M3 M3  
C6. Cosine and sine integrals Ci (#) (WCi) Cosine integral Ci(x) M3 M3   M3 M3  
Si (#) (WSi) Sine integral Si(x) M3 M3   M3 M3  
Chi (#) (WChi) Hyperbolic cosine integral Chi(x) M3 M3   M3 M3  
Shi (#) (WShi) Hyperbolic sine integral Shi(x) M3 M3   M3 M3  
C7a. Gamma functions Gamma Gamma function Γ(x) B     B    
Gamma1pm1 (#) Gamma function Γ(1+x) - 1 M3     M3    
Lngam (#) Logarithm of gamma function ln(Γ(x)) B     B    
Lngams (#) Logarithm of gamma function ln|Γ(x)| and sign of gamma function M3     M3    
Gamr (#) (WGamr) Reciprocal of gamma function 1/Γ(x) M3 M3   M3 M3  
Gamratio (#) Ratio of gamma functions Γ(a)/Γ(b) M3     M3    
Gamdratio (#) Ratio of gamma functions Γ(a)/Γ(a+δ) M3     M3    
Cgamma Gamma function Γ(z) (complex argument) M3     M3    
Clngam Logarithm of gamma function ln(Γ(z)) (complex argument) M3     M3    
Cgamr Reciprocal of gamma function 1/Γ(z) (complex argument) M3     M3    
Poch (#) (WPoch) Pochhammer's symbol (a)x M3 M3   M3 M3  
Poch1 (#) (WPoch1) Relative Pochhammer's symbol ((a)x - 1)/x M3 M3   M3 M3  
C7b. Beta functions Beta (#) (WBeta) Beta function B(a, b) M3 M3   M3 M3  
Lnbeta (#) (WLnbeta) Logarithm of beta function ln(B(a,b)) M3 M3   M3 M3  
Cbeta Beta function B(a, b) (complex argument) M3     M3    
Clnbeta Logarithm of beta function ln(B(a, b)) (complex argument) M3     M3    
C7c. Polygamma functions Digamma (#) (WDigamma) Digamma (or psi) function ψ(x) B B   B B  
Trigamma (#) (WTrigamma) Trigamma function ψ1(x) M3     M3 M3  
Polygamma (#) (WPolygamma) Polygamma function ψn(x) M3 M3   M3 M3  
Cdigamma Digamma (or psi) function ψ(z) (complex argument) M3     M3    
C7e. Incomplete Gamma functions Gami (#) (WGami) Incomplete gamma function γ(a, x) M3 M3   M3 M3  
Gamic (#) (WGamic) Complementary incomplete gamma function Γ(a, x) M3 M3   M3 M3  
Gamit (#) (WGamit) Tricomi's incomplete gamma function γ*(a, x) M3 M3   M3 M3  
Gammap (#) (WGammap) Normalized incomplete gamma function P(a, x) M3     M3 M3  
Gammaq (#) (WGammaq) Normalized complementary incomplete gamma function Q(a, x) M3     M3 M3  
Gammapi (#) (WGammapi) Inverse function of x for normalized incomplete gamma function P(a, x) M3     M3 M3  
Gammaqi (#) (WGammaqi) Inverse function of x for normalized complementary incomplete gamma function Q(a, x) M3     M3 M3  
Gammapia (#) (WGammapia) Inverse function of a for normalized incomplete gamma function P(a, x) M3     M3 M3  
Gammaqia (#) (WGammaqia) Inverse function of a for normalized complementary incomplete gamma function Q(a, x) M3     M3 M3  
Gammapd (#) (WGammapd) Derivative of normalized incomplete gamma function P(a, x) M3     M3 M3  
C7f. Incomplete Beta functions Betax (#) (WBetax) Incomplete beta function Bx(a, b) M3     M3 M3  
Betaxc (#) (WBetaxc) Compliment of incomplete beta function 1 - Bx(a, b) M3     M3 M3  
Ibeta (#) (WIbeta) Normalized incomplete beta function Ix(a, b) M3 M3   M3 M3  
Ibetac (#) (WIbetac) Normalized compliment of incomplete beta function 1 - Ix(a, b) M3     M3 M3  
Ibetai (#) (WIbetai) Normalized incomplete beta function Ix(a, b) inverse for x M3     M3 M3  
Ibetaci (#) (WIbetaci) Normalized compliment of incomplete beta function 1 - Ix(a, b) inverse for x M3     M3 M3  
Ibetaia (#) (WIbetaia) Normalized incomplete beta function Ix(a, b) inverse for a M3     M3 M3  
Ibetacia (#) (WIbetacia) Normalized compliment of incomplete beta function 1 - Ix(a, b) inverse for a M3     M3 M3  
Ibetaib (#) (WIbetaib) Normalized incomplete beta function Ix(a, b) inverse for b M3     M3 M3  
Ibetacib (#) (WIbetacib) Normalized compliment of incomplete beta function 1 - Ix(a, b) inverse for b M3     M3 M3  
Ibetad (#) (WIbetad) Derivative of normalized incomplete beta function Ix(a, b) M3     M3 M3  
C7g. Riemann zeta function Zeta (#) (WZeta) Riemann zeta function ζ(x) M3 M3   M3 M3  
C8. Error functions Erf Error function erf(x) B     B    
Erfc Complementary error function erfc(x) B     B    
Erfi (#) Error function erf(x) inverse M3     M3    
Erfci (#) Complementary error function erf(x) inverse M3     M3    
Dawson (#) (WDawson) Dawson's function F(x) M3 M3   M3 M3  
Fresc (#) (WFresc) Fresnel integral C(x) M3 M3   M3 M3  
Fress (#) (WFress) Fresnel integral S(x) M3 M3   M3 M3  
C10a. Bessel functions Besj0 (#) Bessel function of the first kind of order zero J0(x) B     B    
Besj1 (#) Bessel function of the first kind of order one J1(x) B     B    
Besjn (#) Bessel functions of the first kind of order n Jn(x) M3     M3    
Besjnu (#) (WBesj) Bessel function of the first kind of order ν Jν(x) (fractional order) B B   B B  
Besy0 (#) Bessel function of the second kind of order zero Y0(x) B     B    
Besy1 (#) Bessel function of the second kind of order one Y1(x) B     B    
Besyn (#) Bessel functions of the second kind of order n Yn(x) M3     M3    
Besynu (#) (WBesy) Bessel function of the second kind of order ν Yν(x) (fractional order) B B   B B  
Besjnd (#) Derivative J'n(x) of Bessel function of the first kind of order n Jn(x) M3     M3    
Besjnud (#) (WBesjd) Derivative J'ν(x) of Bessel function of the first kind of order ν Jν(x) (fractional order) M3 M3   M3 M3  
Besynd (#) Derivative Y'n(x) of modified Bessel functions of the second kind of order n Yn(x) M3     M3    
Besynud (#) (WBesyd) Derivative Y'ν(x) of Bessel function of the second kind of order ν Yν(x) (fractional order) M3 M3   M3 M3  
Sbesjn (#) Spherical Bessel function of the first kind jn(x) M3     M3    
Sbesjnu (#) (WSbesj) Spherical Bessel function of the first kind of order ν jν(x) (fractional order) M3 M3   M3 M3  
Sbesyn (#) Spherical Bessel function of the second kind yn(x) M3     M3    
Sbesynu (#) (WSbesy) Spherical Bessel function of the second kind of order ν yν(x) (fractional order) M3 M3   M3 M3  
Cbesh Sequence of Hankel functions Hν(m)(z) (fractional order) (complex argument) M3     M3    
Cbesj Sequence of Bessel functions of the first kind Jν(z) (fractional order) (complex argument) M3     M3    
Cbesy Sequence of Bessel functions of the second kind Yν(z) (fractional order) (complex argument) M3     M3    
C10b. Modified Bessel functions Besi0 (#) Modified Bessel function of the first kind of order zero I0(x) B     B    
Besi1 (#) Modified Bessel function of the first kind of order one I1(x) B     B    
Besin (#) Modified Bessel function of the first kind of order n In(x) M3     M3    
Besinu (#) (WBesi) Modified Bessel function of the first kind of order ν Iν(x) (fractional order) B B   B B  
Besk0 (#) Modified Bessel function of the second kind of order zero K0(x) B     B    
Besk1 (#) Modified Bessel function of the second kind of order one K1(x) B     B    
Beskn (#) Modified Bessel function of the second kind of order n Kn(x) M3     M3    
Besknu (#) (WBesk) Modified Bessel function of the second kind of order ν Kν(x) (fractional order) B B   B B  
Besind (#) Derivative I'n(x) of modified Bessel function of the first kind of order n In(x) M3     M3    
Besinud (#) (WBesid) Derivative I'ν(x) of modified Bessel function of the first kind of order ν Iν(x) (fractional order) M3 M3   M3 M3  
Besknd (#) Derivative K'n(x) of modified Bessel functions of the second kind of order n Kn(x) M3     M3    
Besknud (#) (WBeskd) Derivative K'ν(x) of modified Bessel function of the second kind of order ν Kν(x) (fractional order) M3 M3   M3 M3  
Sbesin (#) Modified spherical Bessel function of the first kind in(x) M3     M3    
Sbesinu (#) (WSbesi) Modified spherical Bessel function of the first kind of order ν iν(x) (fractional order) M3 M3   M3 M3  
Sbeskn (#) Modified spherical Bessel function of the second kind kn(x) M3     M3    
Sbesknu (#) (WSbesk) Modified spherical Bessel function of the second kind of order ν kν(x) (fractional order) M3 M3   M3 M3  
Cbesi Sequence of modified Bessel functions of the first kind Iν(z) (fractional order) (complex argument) M3     M3    
Cbesk Sequence of modified Bessel functions of the second kind Kν(z) (fractional order) (complex argument) M3     M3    
C10d. Airy functions Airyai (#) (WAiryai) Airy function Ai(x) M3 M3   M3 M3  
Airybi (#) (WAirybi) Airy function Bi(x) M3 M3   M3 M3  
Airyaid (#) (WAiryaid) Derivative Ai'(x) of Airy function Ai(x) M3 M3   M3 M3  
Airybid (#) (WAirybid) Derivative Bi'(x) of Airy function Bi(x) M3 M3   M3 M3  
Cairy Airy function Ai(x) or its derivative Ai'(z) (complex argument) M3     M3    
Cbiry Airy function Bi(x) or its derivative Bi'(z) (complex argument) M3     M3    
C11. Hypergeometric functions Hyp1f1 (#) (WHyp1f1) Hypergeometric function 1F1(a; b; z) (Kummer's function M(a, b, z)) M3     M3 M3  
Lhyp1f1 (#) (WLhyp1f1) Logarithm of hypergeometric function ln|1F1(a; b; z)| M3     M3 M3  
Hyp1f1r (#) (WHyp1f1r) Regularized hypergeometric functions 1F1(a; b; z)/Γ(b) M3     M3 M3  
Chu (#) (WChu) Confluent hypergeometric function U(a,b,x) M3     M3 M3  
Hyp2f1 (#) (WHyp2f1) Hypergeometric function 2F1(a1, a2; b; z) (Gaussian hypergeometric function) M3     M3 M3  
Hyp0f1 (#) (WHyp0f1) Hypergeometric function 0F1(; b; z) M3     M3 M3  
Hyp1f0 (#) (WHyp1f0) Hypergeometric function 1F0(a; z) M3     M3 M3  
Hyp2f0 (#) (WHyp2f0) Hypergeometric function 2F0(a1, a2; z) M3     M3 M3  
Hyppfq (WHyppfq) Hypergeometric function pFq(a1, ..., ap; b1, ..., bq; z)       M3    
C13. Jacobi elliptic functions Jelli Jacobi elliptic functions sn(u, k), cn(u,k), dn(u, k) M3     M3    
Jsn (#) (WJsn) Jacobi elliptic functions sn(u, k) M3 M3   M3 M3  
Jcn (#) (WJcn) Jacobi elliptic functions cn(u, k) M3 M3   M3 M3  
Jdn (#) (WJdn) Jacobi elliptic functions dn(u, k) M3 M3   M3 M3  
Jns (#) (WJns) Jacobi elliptic functions ns(u, k) M3     M3 M3  
Jnc (#) (WJnc) Jacobi elliptic functions nc(u, k) M3     M3 M3  
Jnd (#) (WJnd) Jacobi elliptic functions nd(u, k) M3     M3 M3  
Jsc (#) (WJsc) Jacobi elliptic functions sc(u, k) M3     M3 M3  
Jsd (#) (WJsd) Jacobi elliptic functions sd(u, k) M3     M3 M3  
Jdc (#) (WJdc) Jacobi elliptic functions dc(u, k) M3     M3 M3  
Jds (#) (WJds) Jacobi elliptic functions ds(u, k) M3     M3 M3  
Jcs (#) (WJcs) Jacobi elliptic functions cs(u, k) M3     M3 M3  
Jcd (#) (WJcd) Jacobi elliptic functions cd(u, k) M3     M3 M3  
Jtheta1 (#) (WJtheta1) Jacobi theta function θ1(x, q)       M3 M3  
Jheta1t (#) (WJtheta1t) Jacobi theta function θ1(x | τ)       M3 M3  
Jtheta2 (#) (WJtheta2) Jacobi theta function θ2(x, q)       M3 M3  
Jtheta2t (#) (WJtheta2t) Jacobi theta function θ2(x | τ)       M3 M3  
Jtheta3 (#) (WJtheta3) Jacobi theta function θ3(x, q)       M3 M3  
Jtheta3t (#) (WJtheta3t) Jacobi theta function θ3(x | τ)       M3 M3  
Jtheta3m1 (#) (WJtheta3m1) Jacobi theta function θ3(x, q) - 1       M3 M3  
Jtheta3m1t (#) (WJtheta3m1t) Jacobi theta function θ3(x | τ) - 1       M3 M3  
Jtheta4 (#) (WJtheta4) Jacobi theta function θ4(x, q)       M3 M3  
Jtheta4t (#) (WJtheta4t) Jacobi theta function θ4(x | τ)       M3 M3  
Jtheta4m1 (#) (WJtheta4m1) Jacobi theta function θ4(x, q) - 1       M3 M3  
Jtheta4m1t (#) (WJtheta4m1t) Jacobi theta function θ4(x | τ) - 1       M3 M3  
C14. Elliptic Integrals Celli1 (#) (WCelli1) Complete elliptic integral of the first kind K(k) B B   B B  
Celli2 (#) (WCelli2) Complete elliptic integral of the second kind E(k) B B   B B  
Celli3 (#) (WCelli3) Complete elliptic integral of the third kind P(n, k) B B   B B  
Elli1 (#) (WElli1) Incomplete elliptic integral of the first kind F(phi, k) M3 M3   M3 M3  
Elli2 (#) (WElli2) Incomplete elliptic integral of the second kind E(phi, k) M3 M3   M3 M3  
Elli3 (#) (WElli3) Incomplete elliptic integral of the third kind P(phi, n, k) M3 M3   M3 M3  
Rc (#) (WRc) Carlson form of elliptic integral RC(x, y) M3 M3   M3 M3  
Rd (#) (WRd) Carlson form of elliptic integral RD(x, y, z) M3 M3   M3 M3  
Rg (#) (WRg) Carlson form of elliptic integral RG(x, y, z) M3 M3   M3 M3  
Rf (#) (WRf) Carlson form of elliptic integral RF(x, y, z) M3 M3   M3 M3  
Rj (#) (WRj) Carlson form of elliptic integral RJ(x, y, z, p) M3 M3   M3 M3  
Jzeta (#) (WJzeta) Jacobi zeta function Z(φ, k) M3     M3 M3  
Hlambda (#) (WHlambda) Heuman lambda function Λ0(φ, k)       M3 M3  
C19. Other special functions Dconst (#) (WDconst) Numerical quantities B B   B B  
D. Linear algebra                  
D1. Elementary vector and matrix operations D1a. Elementary vector operations: BLAS1 Daxpy y <- ax + y M1     M1    
Dcopy y <- x M1     M1    
Ddot (#) x^T * y M1     M1    
Drotg Constructs Givens plane rotation M1     M1    
Drotmg Constructs modified Givens plane rotation M1     M1    
Drot Applies Givens plane rotation M1     M1    
Drotm Applies modified Givens plane rotation M1     M1    
Dscal x <- ax M1     M1    
Dswap y <-> x M1     M1    
Dasum (#) | X | (1-norm) M1     M1    
Dnrm2 (#) ||X||2 (2-norm of vector) M1     M1    
Zaxpy y <- ax + y (complex vector) M2     M2    
Zcopy y <- x (complex vector) M2     M2    
Zdotu x^T * y (complex vector) M2     M2    
Zdotc x^H * y (complex vector) M2     M2    
Zrotg Constructs Givens plane rotation (complex vector) M2     M2    
Zrot Applies Givens plane rotation (complex vector) M2     M2    
Zdrot Applies Givens plane rotation (complex vector) M2     M2    
Zdscal x <- ax (complex vector) M2     M2    
Zscal x <- ax (complex vector) (a is real number) M2     M2    
Zswap y <-> x (complex vector) M2     M2    
Dzasum |Re(x)|+|Im(x)| (1-norm) (complex vector) M2     M2    
Dznrm2 ||x||2 (2-norm) (complex vector) M2     M2    
D1a. Elementary vector operations: BLAS2 Dgemv y <- αAx+βy or y <- αA^Tx+βy M1     M1    
Dgbmv y <- αAx+βy or y <- αA^Tx+βy (band matrix) M1     M1    
Dsymv y <- αAx+βy (symmetric matrix) M1     M1    
Dsbmv y <- αAx+βy (symmetric band matrix) M1     M1    
Dspmv y <- αAx+βy (symmetric matrix) (packed form) M1     M1    
Dtrmv x <- Op(A)x (Op(A) = A or A^T) (triangular matrix) M1     M1    
Dtbmv x <- Op(A)x (Op(A) = A or A^T) (triangular band matrix) M1     M1    
Dtpmv x <- Op(A)x (Op(A) = A or A^T) (triangular matrix) (packed form) M1     M1    
Dtrsv Solution of Op(A)x = b (Op(A) = A or A^T) (triangular matrix) M1     M1    
Dtbsv Solution of Op(A)x = b (Op(A) = A or A^T) (triangular band matrix) M1     M1    
Dtpsv Solution of Op(A)x = b (Op(A) = A or A^T) (triangular matrix) (packed form) M1     M1    
Dger A <- αxy^T + A M1     M1    
Dsyr A <- αxx^T + A (symmetric matrix) M1     M1    
Dspr A <- αxx^T + A (symmetric matrix) (packed form) M1     M1    
Dsyr2 A <- αxy^T + αyx^T + A (symmetric matrix) M1     M1    
Dspr2 A <- αxy^T + αyx^T + A (symmetric matrix) (packed form) M1     M1    
Zgemv y <- αOp(A)x+βy (Op(A) = A, A^T or A^H) (complex matrix) M2     M2    
Zgbmv y <- αOp(A)x+βy (Op(A) = A, A^T or A^H) (complex band matrix) M2     M2    
Zhemv y <- αAx+βy (Hermitian matrix) M2     M2    
Zhbmv y <- αAx+βy (Hermitian band matrix) M2     M2    
Zhpmv y <- αAx+βy (Hermitian matrix) (packed form) M2     M2    
Zsymv y <- αAx+βy (complex symmetric matrix) M2     M2    
Zsbmv y <- αAx+βy (complex symmetric band matrix) M2     M2    
Zspmv y <- αAx+βy (complex symmetric matrix) (packed form) M2     M2    
Ztrmv x <- Op(A)x (Op(A) = A, A^T or A^H) (complex triangular matrix) M2     M2    
Ztbmv x <- Op(A)x (Op(A) = A, A^T or A^H) (complex triangular band matrix) M2     M2    
Ztpmv x <- Op(A)x (Op(A) = A, A^T or A^H) (complex triangular matrix) (packed form) M2     M2    
Ztrsv Solution of Op(A)x = b (Op(A) = A, A^T or A^H) (complex triangular matrix) M2     M2    
Ztbsv Solution of Op(A)x = b (Op(A) = A, A^T or A^H) (complex triangular band matrix) M2     M2    
Ztpsv Solution of Op(A)x = b  (Op(A) = A, A^T or A^H) (complex triangular matrix) (packed form) M2     M2    
Zgeru A <- αxy^T + A (complex matrix) M2     M2    
Zgerc A <- αxy^H + A (complex matrix) M2     M2    
Zher A <- αxx^H + A (Hermitian matrix) M2     M2    
Zhpr A <- αxx^H + A (Hermitian matrix) (packed form) M2     M2    
Zsyr A <- αxx^T + A (complex symmetric matrix) M2     M2    
Zspr A <- αxx^T + A (complex symmetric matrix) (packed form) M2     M2    
Zher2 A <- αxy^H + conjg(α)yx^H + A (Hermitian matrix) M2     M2    
Zhpr2 A <- αxy^H + conjg(α)yx^H + A (Hermitian matrix) (packed form) M2     M2    
Zsyr2 A <- αxy^T + αyx^T + A (complex symmetric matrix) M2     M2    
Zspr2 A <- αxy^T + αyx^T + A (complex symmetric matrix) (packed form) M2     M2    
D1b. Elementary matrix operations: BLAS3 Dgemm C <- αOp(A)Op(B) + βC (Op(X) = X, X^T) M1     M1    
Dsymm C <- αAB + βC or αBA + βC M1     M1    
Dtrmm B <- αOp(A)B or αBOp(A) (Op(A) = A or A^T) (triangular matrix) M1     M1    
Dtrsm Solution of Op(A)X = αB or XOp(A) = αB (Op(A) = A or A^T) M1     M1    
Dsyrk C <- αAA^T + βC or αA^TA + βC M1     M1    
Dsyr2k C <- αAB^T + αBA^T + βC or αA^TB + αB^TA + βC M1     M1    
Zgemm C <- αOp(A)Op(B) + βC (Op(X) = X, X^T or X^H) (complex matrix) M2     M2    
Zsymm C <- αAB + βC or αBA + βC (complex symmetric matrix) M2     M2    
Zhemm C <- αAB + βC or αBA + βC (Hermitian matrix) M2     M2    
Ztrmm B <- αOp(A)B or αBOp(A) (Op(A) = A, A^T or A^H) (complex triangular matrix) M2     M2    
Ztrsm Solution of Op(A)X = αB or XOp(A) = αB (Op(A) = A, A^T or A^H) (complex triangular matrix) M2     M2    
Zsyrk C <- αAA^T + βC or αA^TA + βC (complex symmetric matrix) M2     M2    
Zherk C <- αAA^H + βC or αA^HA + βC (Hermitian matrix) M2     M2    
Zsyr2k C <- αAB^T + αBA^T + βC or αA^TB + αB^TA + βC (complex symmetric matrix) M2     M2    
Zher2k C <- αAB^H + conjg(α)BA^H + βC or αA^HB + conjg(α)B^HA + βC (Hermitian matrix) M2     M2    
D1b. Elementary matrix operations: norm of matrix Dlange Norm of matrix (general matrix) B     B    
Dlangb Norm of matrix (band matrix) M1     M1    
Dlangt Norm of matrix (tridiagonal matrix) M1     M1    
Dlansy Norm of matrix (symmetric matrix) B     B    
Dlansb Norm of matrix (symmetric band matrix) M1     M1    
Dlansp Norm of matrix (symmetric matrix) (packed form) M1     M1    
Dlanst Norm of matrix (symmetric tridiagonal matrix) M1     M1    
Dlantr Norm of matrix (trapezoidal or triangular matrix) M1     M1    
Zlange Norm of matrix (complex matrix) M2     M2    
Zlangb Norm of matrix (complex band matrix) M2     M2    
Zlangt Norm of matrix (complex tridiagonal matrix) M2     M2    
Zlansy Norm of matrix (complex symmetric matrix) M2     M2    
Zlansb Norm of matrix (complex symmetric band matrix) M2     M2    
Zlansp Norm of matrix (complex symmetric matrix) (packed form) M2     M2    
Zlanhe Norm of matrix (Hermitian matrix) M2     M2    
Zlanhb Norm of matrix (Hermitian band matrix) M2     M2    
Zlanhp Norm of matrix (Hermitian matrix) (packed form) M2     M2    
Zlanht Norm of matrix (Hermitian tridiagonal matrix) M2     M2    
Zlantr Norm of matrix (complex trapezoidal or triangular matrix) M2     M2    
D2. Solution of systems of linear equations D2a. Solution of systems of linear equations (general matrices) Dgesv (WDgesv) (Simple driver) Solution of system of linear equations Ax = b B B   B B  
Dgetrf LU factorization of coefficient matrix M1     M1    
Dgetrs Solution of LU factorized system of linear equations M1     M1    
Dgetri Inverse matrix M1     M1    
Dgesvx (Expert driver)Solution of system of linear equations Ax = b M1     M1    
Dgecon Condition number of matrix B     B    
Dsgesv (Simple driver) Solution of system of linear equations Ax = b (mixed precision with iterative refinement) M1     M1    
Dgbsv (WDgbsv) (Simple driver) Solution of system of linear equations Ax = b (band matrix) M1 M1   M1 M1  
Dgbtrf LU factorization of coefficient matrix (band matrix) M1     M1    
Dgbtrs Solution of LU factorized system of linear equations (band matrix) M1     M1    
Dgbsvx (Expert driver) Solution of system of linear equations Ax = b (band matrix) M1     M1    
Dgbcon Condition number of matrix (band matrix) M1     M1    
Dgtsv (WDgtsv) (Simple driver) Solution of system of linear equations Ax = b (tridiagonal matrix) M1 M1   M1 M1  
Dgttrf LU factorization of coefficient matrix (tridiagonal matrix) M1     M1    
Dgttrs Solution of LU factorized system of linear equations (tridiagonal matrix) M1     M1    
Dgtsvx (Expert driver) Solution of system of linear equations Ax = b (tridiagonal matrix) M1     M1    
Dgtcon Condition number of matrix (tridiagonal matrix) M1     M1    
D2a3. Solution of systems of linear equations (triangular matrices) Dtrtrs (WDtrtrs) Solution of system of linear equations Ax = b (triangular matrix) M1 M1   M1 M1  
Dtrtri Inverse matrix (triangular matrix) M1     M1    
Dtrcon Condition number of matrix (triangular matrix) M1     M1    
Dtptrs Solution of system of linear equations Ax = b (triangular matrix) (packed form) M1     M1    
Dtptri Inverse matrix (triangular matrix) (packed form) M1     M1    
Dtpcon Condition number of matrix (triangular matrix) (packed form) M1     M1    
Dtbtrs Solution of system of linear equations Ax = b (triangular band matrix) M1     M1    
Dtbcon Condition number of matrix (triangular band matrix) M1     M1    
D2b1a. Solution of systems of linear equations (symmetric matrices) Dsysv (WDsysv) (Simple driver) Solution of system of linear equations Ax = b (symmetric matrix) M1 M1   M1 M1  
Dsytrf UDU^T or LDL^T factorization of coefficient matrix (symmetric matrix) M1     M1    
Dsytrs Solution of UDU^T or LDL^T factorized system of linear equations (symmetric matrix) M1     M1    
Dsytri Inverse matrix (symmetric matrix) M1     M1    
Dsysvx (Expert driver) Solution of system of linear equations Ax = b (symmetric matrix) M1     M1    
Dsycon Condition number of matrix (symmetric matrix) M1     M1    
Dspsv (Simple driver) Solution of system of linear equations Ax = b (symmetric matrix) (packed form) M1     M1    
Dsptrf UDU^T or LDL^T factorization of coefficient matrix (symmetric matrix) (packed form) M1     M1    
Dsptrs Solution of UDU^T or LDL^T factorized system of linear equations (symmetric matrix) (packed form) M1     M1    
Dsptri Inverse matrix (symmetric matrix) (packed form) M1     M1    
Dspsvx (Expert driver) Solution of system of linear equations Ax = b (symmetric matrix) (packed form) M1     M1    
Dspcon Condition number of matrix (symmetric matrix) (packed form) M1     M1    
D2b1b. Solution of systems of linear equations (symmetric positive definite matrices) Dposv (WDposv) (Simple driver) Solution of system of linear equations Ax = b (symmetric positive definite matrix) B B   B B  
Dpotrf Cholesky factorization of coefficient matrix (symmetric positive definite matrix) M1     M1    
Dpotrs Solution of Cholesky factorized system of linear equations (symmetric positive definite matrix) M1     M1    
Dpotri Inverse matrix (symmetric positive definite matrix) M1     M1    
Dposvx (Expert driver) Solution of system of linear equations Ax = b (symmetric positive definite matrix) M1     M1    
Dpocon Condition number of matrix (symmetric positive definite matrix) B     B    
Dsposv (Simple driver) Solution of system of linear equations Ax = b (symmetric positive definite matrix) (mixed precision with iterative refinement) M1     M1    
Dppsv (Simple driver) Solution of system of linear equations Ax = b (symmetric positive definite matrix) (packed form) M1     M1    
Dpptrf Cholesky factorization of coefficient matrix (symmetric positive definite matrix) (packed form) M1     M1    
Dpptrs Solution of Cholesky factorized system of linear equations (symmetric positive definite matrix) (packed form) M1     M1    
Dpptri Inverse matrix (symmetric positive definite matrix) (packed form) M1     M1    
Dppsvx (Expert driver) Solution of system of linear equations Ax = b (symmetric positive definite matrix) (packed form) M1     M1    
Dppcon Condition number of matrix (symmetric positive definite matrix) (packed form) M1     M1    
D2b2. Solution of systems of linear equations (symmetric positive definite banded matrices) Dpbsv (WDpbsv) (Simple driver) Solution of system of linear equations Ax = b (symmetric positive definite band matrix) M1 M1   M1 M1  
Dpbtrf Cholesky factorization of coefficient matrix (symmetric positive definite band matrix) M1     M1    
Dpbtrs Solution of Cholesky factorized system of linear equations (symmetric positive definite band matrix) M1     M1    
Dpbsvx (Expert driver) Solution of system of linear equations Ax = b (symmetric positive definite band matrix) M1     M1    
Dpbcon Condition number of matrix (symmetric positive definite band matrix) M1     M1    
Dptsv (WDptsv) (Simple driver) Solution of system of linear equations Ax = b (symmetric positive definite tridiagonal matrix) M1 M1   M1 M1  
Dpttrf LDL^T factorization of coefficient matrix (symmetric positive definite tridiagonal matrix) M1     M1    
Dpttrs Solution of LDL^T factorized system of linear equations (symmetric positive definite tridiagonal matrix) M1     M1    
Dptsvx (Expert driver) Solution of system of linear equations Ax = b (symmetric positive definite tridiagonal matrix) M1     M1    
Dptcon Condition number of matrix (symmetric positive definite tridiagonal matrix) M1     M1    
D2c. Solution of systems of linear equations (general complex matrices) Zgesv (WZgesv(2)) (Simple driver) Solution of system of linear equations Ax = b (complex matrix) M2 M2   M2 M2  
Zgetrf LU factorization of coefficient matrix (complex matrix) M2     M2    
Zgetrs Solution of LU factorized system of linear equations (complex matrix) M2     M2    
Zgetri Inverse matrix (complex matrix) M2     M2    
Zgesvx (Expert driver) Solution of system of linear equations Ax = b (complex matrix) M2     M2    
Zgecon Condition number of matrix (complex matrix) M2     M2    
Zcgesv (Simple driver) Solution of system of linear equations Ax = b (mixed precision with iterative refinement) (complex matrix) M2     M2    
Zgbsv (WZgbsv(2)) (Simple driver) Solution of system of linear equations Ax = b (complex band matrix) M2 M2   M2 M2  
Zgbtrf LU factorization of coefficient matrix (complex band matrix) M2     M2    
Zgbtrs Solution of LU factorized system of linear equations (complex band matrix) M2     M2    
Zgbsvx (Expert driver) Solution of system of linear equations Ax = b (complex band matrix) M2     M2    
Zgbcon Condition number of matrix (complex band matrix) M2     M2    
Zgtsv (WZgtsv(2)) (Simple driver) Solution of system of linear equations Ax = b (complex tridiagonal matrix) M2 M2   M2 M2  
Zgttrf LU factorization of coefficient matrix (complex tridiagonal matrix) M2     M2    
Zgttrs Solution of LU factorized system of linear equations (complex tridiagonal matrix) M2     M2    
Zgtsvx (Expert driver) Solution of system of linear equations Ax = b (complex tridiagonal matrix) M2     M2    
Zgtcon Condition number of matrix (complex tridiagonal matrix) M2     M2    
Zsysv (WZsysv(2)) (Simple driver) Solution of system of linear equations Ax = b (complex symmetric matrix) M2 M2   M2 M2  
Zsytrf UDU^H or LDL^H factorization of coefficient matrix (complex symmetric matrix) M2     M2    
Zsytrs Solution of UDU^H or LDL^H factorized system of linear equations (complex symmetric matrix) M2     M2    
Zsytri Inverse matrix (complex symmetric matrix) M2     M2    
Zsysvx (Expert driver) Solution of system of linear equations Ax = b (complex symmetric matrix) M2     M2    
Zsycon Condition number of matrix (complex symmetric matrix) M2     M2    
Zspsv (Simple driver) Solution of system of linear equations Ax = b (complex symmetric matrix) (packed form) M2     M2    
Zsptrf UDU^H or LDL^H factorization of coefficient matrix (complex symmetric matrix) (packed form) M2     M2    
Zsptrs Solution of UDU^H or LDL^H factorized system of linear equations (complex symmetric matrix) (packed form) M2     M2    
Zsptri Inverse matrix (complex symmetric matrix) (packed form) M2     M2    
Zspsvx (Expert driver) Solution of system of linear equations Ax = b (complex symmetric matrix) (packed form) M2     M2    
Zspcon Condition number of matrix (complex symmetric matrix) (packed form) M2     M2    
D2c3. Solution of systems of linear equations (triangular complex matrices) Ztrtrs (WZtrtrs(2)) Solution of system of linear equations Ax = b (complex triangular matrix) M2 M2   M2 M2  
Ztrtri Inverse matrix (complex triangular matrix) M2     M2    
Ztrcon Condition number of matrix (complex triangular matrix) M2     M2    
Ztptrs Solution of system of linear equations Ax = b (complex triangular matrix) (packed form) M2     M2    
Ztptri Inverse matrix (complex triangular matrix) (packed form) M2     M2    
Ztpcon Condition number of matrix (complex triangular matrix) (packed form) M2     M2    
Ztbtrs Solution of system of linear equations Ax = b (complex triangular band matrix) M2     M2    
Ztbcon Condition number of matrix (complex triangular band matrix) M2     M2    
D2d1a. Solution of systems of linear equations (Hermitian matrices) Zhesv (WZhesv(2)) (Simple driver) Solution of system of linear equations Ax = b (Hermitian matrix) M2 M2   M2 M2  
Zhetrf UDU^H or LDL^H factorization of coefficient matrix (Hermitian matrix) M2     M2    
Zhetrs Solution of UDU^H or LDL^H factorized system of linear equations (Hermitian matrix) M2     M2    
Zhetri Inverse matrix (Hermitian matrix) M2     M2    
Zhesvx (Expert driver) Solution of system of linear equations Ax = b (Hermitian matrix) M2     M2    
Zhecon Condition number of matrix (Hermitian matrix) M2     M2    
Zhpsv (Simple driver) Solution of system of linear equations Ax = b (Hermitian matrix) (packed form) M2     M2    
Zhptrf UDU^H or LDL^H factorization of coefficient matrix (Hermitian matrix) (packed form) M2     M2    
Zhptrs Solution of UDU^H or LDL^H factorized system of linear equations (Hermitian matrix) (packed form) M2     M2    
Zhptri Inverse matrix (Hermitian matrix) (packed form) M2     M2    
Zhpsvx (Expert driver) Solution of system of linear equations Ax = b (Hermitian matrix) (packed form) M2     M2    
Zhpcon Condition number of matrix (Hermitian matrix) (packed form) M2     M2    
D2d1b. Solution of systems of linear equations (positive definite Hermitian matrices) Zposv (WZposv(2)) (Simple driver) Solution of system of linear equations Ax = b (Hermitian positive definite matrix) M2 M2   M2 M2  
Zpotrf Cholesky factorization of coefficient matrix (Hermitian positive definite matrix) M2     M2    
Zpotrs Solution of Cholesky factorized system of linear equations (Hermitian positive definite matrix) M2     M2    
Zpotri Inverse matrix (Hermitian positive definite matrix) M2     M2    
Zposvx (Expert driver) Solution of system of linear equations Ax = b (Hermitian positive definite matrix) M2     M2    
Zpocon Condition number of matrix (Hermitian positive definite matrix) M2     M2    
Zcposv (Simple driver) Solution of system of linear equations Ax = b (Hermitian positive definite matrix) (mixed precision with iterative refinement) M2     M2    
Zppsv (Simple driver) Solution of system of linear equations Ax = b (Hermitian positive definite matrix) (packed form) M2     M2    
Zpptrf Cholesky factorization of coefficient matrix (Hermitian positive definite matrix) (packed form) M2     M2    
Zpptrs Solution of Cholesky factorized system of linear equations (Hermitian positive definite matrix) (packed form) M2     M2    
Zpptri Inverse matrix (Hermitian positive definite matrix) (packed form) M2     M2    
Zppsvx (Expert driver) Solution of system of linear equations Ax = b (Hermitian positive definite matrix) (packed form) M2     M2    
Zppcon Condition number of matrix (Hermitian positive definite matrix) (packed form) M2     M2    
D2d2. Solution of systems of linear equations (positive definite banded Hermitian matrices) Zpbsv (WZpbsv(2)) (Simple driver) Solution of system of linear equations Ax = b (Hermitian positive definite band matrix) M2 M2   M2 M2  
Zpbtrf Cholesky factorization of coefficient matrix (Hermitian positive definite band matrix) M2     M2    
Zpbtrs Solution of Cholesky factorized system of linear equations (Hermitian positive definite band matrix) M2     M2    
Zpbsvx (Expert driver) Solution of system of linear equations Ax = b (Hermitian positive definite band matrix) M2     M2    
Zpbcon Condition number of matrix (Hermitian positive definite band matrix) M2     M2    
Zptsv (WZptsv(2)) (Simple driver) Solution of system of linear equations Ax = b (Hermitian positive definite tridiagonal matrix) M2 M2   M2 M2  
Zpttrf LDL^H factorization of coefficient matrix (Hermitian positive definite tridiagonal matrix) M2     M2    
Zpttrs Solution of LDL^H factorized system of linear equations (Hermitian positive definite tridiagonal matrix) M2     M2    
Zptsvx (Expert driver) Solution of system of linear equations Ax = b (Hermitian positive definite tridiagonal matrix) M2     M2    
Zptcon Condition number of matrix (Hermitian positive definite tridiagonal matrix) M2     M2    
D4. Eigenvalues and eigenvectors D4a1. Ordinary eigenvalue problems (symmetric matrices) Dsyev (WDsyev) (Simple driver) Eigenvalues and eigenvectors (symmetric matrix) B B   B B  
Dsyevd (Divide and conquer driver) Eigenvalues and eigenvectors (symmetric matrix) M1     M1    
Dsyevr (MRRR driver) Eigenvalues and eigenvectors (symmetric matrix) M1     M1    
Dsyevx (Expert driver) Eigenvalues and eigenvectors (symmetric matrix) M1     M1    
Dsytrd Reduces a real symmetric matrix to tridiagonal form M1     M1    
Dorgtr Generates a transform matrix from a real symmetric matrix to tridiagonal form M1     M1    
Dormtr Multiplies by a transform matrix from a real symmetric matrix to tridiagonal form M1     M1    
Dsteqr Eigenvalues and eigenvectors of a symmetric tridiagonal matrix (QL or QR method) M1     M1    
Dsterf Eigenvalues of a symmetric tridiagonal matrix (QL or QR method) M1, M2     M1, M2    
Dstedc Eigenvalues and eigenvectors of a symmetric tridiagonal matrix (Divide and conquer method) M1     M1    
Dstemr Eigenvalues and eigenvectors of a symmetric tridiagonal matrix (MRRR method) M1     M1    
Dstebz Eigenvalues of a symmetric tridiagonal matrix (Bisection method) M1. M2     M1. M2    
Dstein Eigenvectors of a symmetric tridiagonal matrix (Inverse iteration method) M1     M1    
Dpteqr Eigenvalues and eigenvectors of a symmetric positive definite tridiagonal matrix M1     M1    
Dspev (Simple driver) Eigenvalues and eigenvectors (symmetric matrix) (packed form) M1     M1    
Dspevd (Divide and conquer driver) Eigenvalues and eigenvectors (symmetric matrix) (packed form) M1     M1    
Dspevx (Expert driver) Eigenvalues and eigenvectors (symmetric matrix) (packed form) M1     M1    
Dsptrd Reduces a real symmetric matrix to tridiagonal form (packed form) M1     M1    
Dopgtr Generates a transform matrix from a real symmetric matrix to tridiagonal form (packed form) M1     M1    
Dopmtr Multiplies by a transform matrix from a real symmetric matrix to tridiagonal form (packed form) M1     M1    
Dsbev (WDsbev) (Simple driver) Eigenvalues and eigenvectors (symmetric band matrix) M1 M1   M1 M1  
Dsbevd (Divide and conquer driver) Eigenvalues and eigenvectors (symmetric band matrix) M1     M1    
Dsbevx (Expert driver) Eigenvalues and eigenvectors (symmetric band matrix) M1     M1    
Dsbtrd Reduces a real symmetric band matrix to tridiagonal form M1     M1    
Dstev (WDstev) (Simple driver) Eigenvalues and eigenvectors (symmetric tridiagonal matrix) M1 M1   M1 M1  
Dstevd (Divide and conquer driver) Eigenvalues and eigenvectors (symmetric tridiagonal matrix) M1     M1    
Dstevr (MRRR driver) Eigenvalues and eigenvectors (symmetric tridiagonal matrix) M1     M1    
Dstevx (Expert driver) Eigenvalues and eigenvectors (symmetric tridiagonal matrix) M1     M1    
Ddisna Condition numbers for the eigenvectors M1. M2     M1. M2    
D4a2. Ordinary eigenvalue problems (general matrices) Dgeev (WDgeev) (Simple driver) Eigenvalues and eigenvectors M1 M1   M1 M1  
Dgeevx (Expert driver) Eigenvalues and eigenvectors M1     M1    
Dgehrd Reduces a real general matrix to upper Hessenberg form M1     M1    
Dgebal Balancing of a real general matrix M1     M1    
Dgebak Eigenvectors of original real general matrix by backward transformation on balanced matrix M1     M1    
Dorghr Generates a transform matrix to Hessenberg form M1     M1    
Dormhr Multiplies by a transform matrix to Hessenberg form M1     M1    
Dhseqr Eigenvalues and Schur factorization of Hessenberg matrix (QR method) M1     M1    
Dhsein Eigenvectors of Hessenberg matrix (Inverse iteration method) M1     M1    
Dtrevc3 Eigenvectors of quasi-triangular matrix of Schur factorization M1     M1    
Dtrexc Reordering of real Schur factorization of real matrix M1     M1    
Dtrsyl Solve real Sylvester matrix equation M1     M1    
Dtrsna Condition numbers for eigenvalues and/or eigenvectors of upper quasi-triangular matrix M1     M1    
Dtrsen Reordering of real Schur factorization of real matrix and condition numbers of cluster of eigenvalues and/or invariant subspace M1     M1    
Dgees (Simple driver) Schur decomposition M1     M1    
Dgees_r (Simple driver) Schur decomposition (reverse communication version) M1     M1    
Dgeesx (Expert driver) Schur decomposition M1     M1    
Dgeesx_r (Expert driver) Schur decomposition (reverse communication version) M1     M1    
D4a3. Ordinary eigenvalue problems (Hermitian matrices) Zheev (WZheev(2)) (Simple driver) Eigenvalues and eigenvectors (Hermitian matrix) M2 M2   M2 M2  
Zheevd (Divide and conquer driver) Eigenvalues and eigenvectors (Hermitian matrix) M2     M2    
Zheevr (MRRR driver) Eigenvalues and eigenvectors (Hermitian matrix) M2     M2    
Zheevx (Expert driver) Eigenvalues and eigenvectors (Hermitian matrix) M2     M2    
Zhetrd Reduces a real Hermitian matrix to tridiagonal form M2     M2    
Zungtr Generates a transform matrix from a complex Hermitian matrix to tridiagonal form M2     M2    
Zunmtr Multiplies by a transform matrix from a complex Hermitian matrix to tridiagonal form M2     M2    
Zsteqr Eigenvalues and eigenvectors of a symmetric tridiagonal matrix to which a Hermitian matrix was reduced (QL or QR method) M2     M2    
Zstedc Eigenvalues and eigenvectors of a symmetric tridiagonal matrix to which a Hermitian matrix was reduced (Divide and conquer method) M2     M2    
Zstemr Eigenvalues and eigenvectors of a symmetric tridiagonal matrix to which a Hermitian matrix was reduced (MRRR method) M2     M2    
Zstein Eigenvectors of a symmetric tridiagonal matrix to which a Hermitian matrix was reduced (Inverse iteration method) M2     M2    
Zpteqr Eigenvalues and eigenvectors of symmetric positive definite tridiagonal matrix to which a Hermitian matrix was reduced M2     M2    
Zhpev (Simple driver) Eigenvalues and eigenvectors (Hermitian matrix) (packed form) M2     M2    
Zhpevd (Divide and conquer driver) Eigenvalues and eigenvectors (Hermitian matrix) (packed form) M2     M2    
Zhpevx (Expert driver) Eigenvalues and eigenvectors (Hermitian matrix) (packed form) M2     M2    
Zhptrd Reduces a real Hermitian matrix to tridiagonal form (packed form) M2     M2    
Zupgtr Generates a transform matrix from a complex Hermitian matrix to tridiagonal form (packed form) M2     M2    
Zupmtr Multiplies by a transform matrix from a complex Hermitian matrix to tridiagonal form (packed form) M2     M2    
Zhbev (WZhbev(2)) (Simple driver) Eigenvalues and eigenvectors (Hermitian band matrix) M2 M2   M2 M2  
Zhbevd (Divide and conquer driver) Eigenvalues and eigenvectors (Hermitian band matrix) M2     M2    
Zhbevx (Expert driver) Eigenvalues and eigenvectors (Hermitian band matrix) M2     M2    
Zhbtrd Reduces a real Hermitian band matrix to tridiagonal form M2     M2    
D4a4. Ordinary eigenvalue problems (general complex matrices) Zgeev (WZgeev(2)) (Simple driver) Eigenvalues and eigenvectors (complex matrix) M2 M2   M2 M2  
Zgeevx (Expert driver) Eigenvalues and eigenvectors (complex matrix) M2     M2    
Zgehrd Reduces a general matrix to upper Hessenberg form (complex matrix) M2     M2    
Zgebal Balancing of a general matrix (complex matrix) M2     M2    
Zgebak Eigenvectors of original general matrix by backward transformation on balanced matrix (complex matrix) M2     M2    
Zunghr Generates a transform matrix to Hessenberg form (complex matrix) M2     M2    
Zunmhr Multiplies by a transform matrix to Hessenberg form (complex matrix) M2     M2    
Zhseqr Eigenvalues and Schur factorization of Hessenberg matrix (QR method) (complex matrix) M2     M2    
Zhsein Eigenvectors of Hessenberg matrix (Inverse iteration method) (complex matrix) M2     M2    
Ztrevc3 Eigenvectors of triangular matrix of Schur factorization (complex matrix) M2     M2    
Ztrexc Reordering of Schur factorization (complex matrix) M2     M2    
Ztrsyl Solve real Sylvester matrix equation (complex matrix) M2     M2    
Ztrsna Condition numbers for eigenvalues and/or eigenvectors of upper triangular matrix (complex matrix) M2     M2    
Ztrsen Reordering of real Schur factorization and condition numbers of cluster of eigenvalues and/or invariant subspace (complex matrix) M2     M2    
Zgees (Simple driver) Schur decomposition (complex matrix) M2     M2    
Zgees_r (Simple driver) Schur decomposition (complex matrix) (reverse communication version) M2     M2    
Zgeesx (Expert driver) Schur decomposition (complex matrix) M2     M2    
Zgeesx_r (Expert driver) Schur decomposition (complex matrix) (reverse communication version) M2     M2    
D4b1. Generalized eigenvalue problems (symmetric matrices) Dsygv (WDsygv) (Simple driver) Generalized eigenvalue problem (symmetric matrix) M1 M1   M1 M1  
Dsygvd (Divide and conquer driver) Generalized eigenvalue problem (symmetric matrix) M1     M1    
Dsygvx (Expert driver) Generalized eigenvalue problem (symmetric matrix) M1     M1    
Dspgv (Simple driver) Generalized eigenvalue problem (symmetric matrix) (packed form) M1     M1    
Dspgvd (Divide and conquer driver) Generalized eigenvalue problem (symmetric matrix) (packed form) M1     M1    
Dspgvx (Expert driver) Generalized eigenvalue problem (symmetric matrix) (packed form) M1     M1    
Dsbgv (WDsbgv) (Simple driver) Generalized eigenvalue problem (symmetric band matrix) M1 M1   M1 M1  
Dsbgvd (Divide and conquer driver) Generalized eigenvalue problem (symmetric band matrix) M1     M1    
Dsbgvx (Expert driver) Generalized eigenvalue problem (symmetric band matrix) M1     M1    
D4b2. Generalized eigenvalue problems (general matrices) Dggev (WDggev) (Simple driver) Generalized eigenvalue problem M1 M1   M1 M1  
Dggevx (Expert driver) Generalized eigenvalue problem M1     M1    
Dgges (Simple driver) Generalized Schur decomposition M1     M1    
Dgges_r (Simple driver) Generalized Schur decomposition (reverse communication version) M1     M1    
Dggesx (Expert driver) Generalized Schur decomposition M1     M1    
Dggesx_r (Expert driver) Generalized Schur decomposition (reverse communication version) M1     M1    
D4b3. Generalized eigenvalue problems (Hermitian matrices) Zhegv (WZhegv(2)) (Simple driver) Generalized eigenvalue problem (Hermitian matrix) M2 M2   M2 M2  
Zhegvd (Divide and conquer driver) Generalized eigenvalue problem (Hermitian matrix) M2     M2    
Zhegvx (Expert driver) Generalized eigenvalue problem (Hermitian matrix) M2     M2    
Zhpgv (Simple driver) Generalized eigenvalue problem (Hermitian matrix) (packed form) M2     M2    
Zhpgvd (Divide and conquer driver) Generalized eigenvalue problem (Hermitian matrix) (packed form) M2     M2    
Zhpgvx (Expert driver) Generalized eigenvalue problem (Hermitian matrix) (packed form) M2     M2    
Zhbgv (WZhbgv(2)) (Simple driver) Generalized eigenvalue problem (Hermitian band matrix) M2 M2   M2 M2  
Zhbgvd (Divide and conquer driver) Generalized eigenvalue problem (Hermitian band matrix) M2     M2    
Zhbgvx (Expert driver) Generalized eigenvalue problem (Hermitian band matrix) M2     M2    
D4b4. Generalized eigenvalue problems (general complex matrices) Zggev (WZggev(2)) (Simple driver) Generalized eigenvalue problem (complex matrix) M2 M2   M2 M2  
Zggevx (Expert driver) Generalized eigenvalue problem (complex matrix) M2     M2    
Zgges (Simple driver) Generalized Schur decomposition (complex matrix) M2     M2    
Zgges_r (Simple driver) Generalized Schur decomposition (complex matrix) (reverse communication version) M2     M2    
Zggesx (Expert driver) Generalized Schur decomposition (complex matrix) M2     M2    
Zggesx_r (Expert driver) Generalized Schur decomposition (complex matrix) (reverse communication version) M2     M2    
D5. QR factorization D5. QR factorization Dgeqp3 QR factorization with pivoting M1     M1    
Dgeqrf QR factorization M1     M1    
Dorgqr Generates matrix Q of QR factorization M1     M1    
Dormqr Multiplies matrix Q of QR factorization M1     M1    
Dgelqf LQ factorization M1     M1    
Dorglq Generates matrix Q of LQ factorization M1     M1    
Dormlq Multiplies matrix Q of LQ factorization M1     M1    
Zgeqp3 QR factorization with pivoting (complex matrix) M2     M2    
Zgeqrf QR factorization (complex matrix) M2     M2    
Zungqr Generates matrix Q of QR factorization (complex matrix) M2     M2    
Zunmqr Multiplies matrix Q of QR factorization (complex matrix) M2     M2    
Zgelqf LQ factorization (complex matrix) M2     M2    
Zunglq Generates matrix Q of LQ factorization (complex matrix) M2     M2    
Zunmlq Multiplies matrix Q of LQ factorization (complex matrix) M2     M2    
D6. Singular value decomposition D6. Singular value decomposition (SVD) Dgesvd (WDgesvd) (Simple driver) Singular value decomposition (SVD) M1 M1   M1 M1  
Dgesvdx (Expert driver) Singular value decomposition (SVD) M1     M1    
Dgesdd (Divide and conquer driver) Singular value decomposition (SVD) M1     M1    
Dgesvdq Singular value decomposition (SVD) (preconditioned QR method) M1     M1    
Dgejsv Singular value decomposition (SVD) (preconditioned Jacobi SVD algorithm) M1     M1    
Dggsvd3 (WDggsvd3) Generalized singular value decomposition (GSVD) M1 M1   M1 M1  
Zgesvd (WZgedvs(2)) (Simple driver) Singular value decomposition (SVD) (complex matrix) M2 M2   M2 M2  
Zgesvdx (Expert driver) Singular value decomposition (SVD) (complex matrix) M2     M2    
Zgesdd (Divide and conquer driver) Singular value decomposition (SVD) (complex matrix) M2     M2    
Zgesvdq Singular value decomposition (SVD) (preconditioned QR method) (complex matrix) M2     M2    
Zgejsv Singular value decomposition (SVD) (preconditioned Jacobi SVD algorithm) (complex matrix) M2     M2    
Zggsvd3(WZggsvd3(2)) Generalized singular value decomposition (GSVD) (complex matrix) M2 M2   M2 M2  
D9. Overdetermined or underdetermined systems of linear equations D9a. Overdetermined or underdetermined systems of linear equations (unconstrained) Dgels (WDgels) Full rank overdetermined or underdetermined linear systems B B   B B  
Dgetsls Full rank overdetermined or underdetermined linear systems (Tall skinny QR or Short wide LQ factorization) M1     M1    
Dgelsy (WDgelsy) Overdetermined or underdetermined linear systems (orthogonal factorization) M1 M1   M1 M1  
Dgelss (WDgelss) Overdetermined or underdetermined linear systems (SVD) M1 M1   M1 M1  
Dgelsd Overdetermined or underdetermined linear systems (SVD) (Divide and conquer method) M1     M1    
Zgels (WZgels(2)) Full rank overdetermined or underdetermined linear systems (complex matrix) M2 M2   M2 M2  
Zgetsls Full rank overdetermined or underdetermined linear systems (Tall skinny QR or Short wide LQ factorization) (complex matrix) M2     M2    
Zgelsy (WZgelsy(2)) Overdetermined or underdetermined linear systems (orthogonal factorization) (complex matrix) M2 M2   M2 M2  
Zgelss (WZgelss(2)) Overdetermined or underdetermined linear systems (SVD) (complex matrix) M2 M2   M2 M2  
Zgelsd Overdetermined or underdetermined linear systems (SVD) (Divide and conquer method) (complex matrix) M2     M2    
Dgecov Variance-covariance matrix of LLS factorized by Dgels B     B    
Dgecovy Variance-covariance matrix of LLS factorized by Dgelsy M1     M1    
Dgecovs Variance-covariance matrix of LLS factorized by Dgelss M1     M1    
Zgecov Variance-covariance matrix of LLS factorized by Zgels (complex matrix) M2     M2    
Zgecovy Variance-covariance matrix of LLS factorized by Zgelsy (complex matrix) M2     M2    
Zgecovs Variance-covariance matrix of LLS factorized by Zgelss (complex matrix) M2     M2    
D9b. Overdetermined or underdetermined systems of linear equations (constrained) Dgglse (WDgglse) Linear equality-constrained least squares (LSE) problem M1 M1   M1 M1  
Dggglm (WDggglm) General Gauss-Markov linear model (GLM) problem M1 M1   M1 M1  
Zgglse (WZgglse(2)) Linear equality-constrained least squares (LSE) problem (complex matrix) M2 M2   M2 M2  
Zggglm (WZggglm(2)) General Gauss-Markov linear model (GLM) problem (complex matrix) M2 M2   M2 M2  
E. Interpolation E. Interpolation (polynomial interpolation) Polint Polynomial interpolation M4     M4    
Polyvl Value of polynomial and derivatives M4     M4    
Polcof Coefficients of polynomial interpolation M4     M4    
Fitlag Iterative Lagrange interpolation M4     M4    
E. Interpolation (piecewise cubic Hermite interpolation / cubic spline interpolation) Pchim Piecewise cubic Hermite interpolation (default boundary conditions) M4     M4    
Pchic Piecewise cubic Hermite interpolation M4     M4    
Pchse (WPchse) Piecewise cubic spline interpolation ("not a not" condition) B B   B B  
Pchsp Piecewise cubic spline interpolation M4     M4    
Pchfe (WPchfe) Evaluation of function values for piecewise cubic Hermite (or cubic spline) interpolation B B   B B  
Pchfd Evaluation of function and derivative values for piecewise cubic Hermite (or cubic spline) interpolation M4     M4    
Chfev Cubic Hermite function values M4     M4    
Chfdv Cubic Hermite function and derivative values M4     M4    
Pchbs Piecewise cubic Hermite to B-spline conversion M4     M4    
Pchcm Monotonicity check for piecewise cubic Hermite function  M4     M4    
E. Interpolation (B-spline interpolation) Bint4 B-representation of cubic spline interpolation M4     M4    
Bintk B-representation of k-th order spline interpolation M4     M4    
Bvalue Evaluation of function or derivative value for B-representation of B-spline M4     M4    
Ppvalu Evaluation of function or derivative value for PP (piecewise polynomial) form of B-spline M4     M4    
Bsplpp B-representation to PP (piecewise polynomial) form of B-spline conversion M4     M4    
Bsplvn Compute the value of B-spline basis functions M4     M4    
Bsplvd Compute the value and the derivatives of B-spline basis functions M4     M4    
Bspldr Construct a divided difference table from B-representation for derivative calculation by Bsplev M4     M4    
Bsplev Evaluation of function and derivative values for B-representation of B-spline M4     M4    
Interv Compute Ileft for the input to Bsplvn and Bsplvd M4     M4    
Banfac LU factorization of banded coefficient matrix of system of linear equations (support routine for Bint4 and Bintk) M4     M4    
Banslv Solution of LU factorized system of linear equations  (support routine for Bint4 and Bintk) M4     M4    
E3a3. Quadrature involving fitted functions Pchia (WPchia) Integral of piecewise cubic Hermite / cubic spline function B B   B B  
Pchid Integral of piecewise cubic Hermite / cubic spline function (over an interval whoes endpoints are data points) M4     M4    
Bsqad Integral of B-representation of B-spline M4     M4    
Bfqad Integral of product of arbitrary function and B-representation of B-spline M4     M4    
Bfqad_r Integral of product of arbitrary function and B-representation of B-spline (reverse communication version) M4     M4    
Ppqad Integral of PP (piecewise polynomial) form of B-spline M4     M4    
Pfqad Integral of product of arbitrary function and PP (piecewise polynomial) form of B-spline M4     M4    
Pfqad_r Integral of product of arbitrary function and PP (piecewise polynomial) form of B-spline (reverse communication version) M4     M4    
F. Solution of nonlinear equations                  
F1a. Roots of polynomials F1a. Roots of polynomials Cpzero (WCpzero(2)) Roots of a polynomial (complex coefficients) (2nd order simultaneous iterative method) M3 M3   M3 M3  
Rpzero Roots of a polynomial (real coefficients) (2nd order simultaneous iterative method) M3     M3    
Rpzero2 (WRpzero2) Roots of a polynomial (real coefficients) (2nd order simultaneous iterative method) (Complex type is not used) B B   B B  
Cpqr79 (WCpqr79(2)) Roots of a polynomial (complex coefficients) (Companion matrix eigenvalues) M3 M3   M3 M3  
Rpqr79 Roots of a polynomial (real coefficients) (Companion matrix eigenvalues) M3     M3    
Dka (WDka(2)) Roots of a polynomial (complex coefficients) (3rd order Durand-Kerner-Aberth (DKA) method) M3 M3   M3 M3  
F1b. Solution of single general nonlinear equation F1b. Solution of single general nonlinear equation Dfzero Zero of the general nonlinear function B   B B   B
Dfzero_r Zero of the general nonlinear function (reverse communication version) B     B    
F2. Solution of a system of nonlinear equations F2. Solution of a system of nonlinear equations Hybrj Solution of a system of nonlinear equations by Powell hybrid method M3     M3    
Hybrj_r Solution of a system of nonlinear equations by Powell hybrid method (reverse communication version) M3     M3    
Hybrj1 Solution of a system of nonlinear equations by Powell hybrid method (simple driver) M3   M3 M3   M3
Hybrj1_r Solution of a system of nonlinear equations by Powell hybrid method (simple driver) (reverse communication version) M3     M3    
Hybrd Solution of a system of nonlinear equations by Powell hybrid method (Jacobian not required) M3     M3    
Hybrd_r Solution of a system of nonlinear equations by Powell hybrid method (Jacobian not required) (reverse communication version) M3     M3    
Hybrd1 Solution of a system of nonlinear equations by Powell hybrid method (Jacobian not required) (simple driver) B   B B   B
Hybrd1_r Solution of a system of nonlinear equations by Powell hybrid method (Jacobian not required) (simple driver) (reverse communication version) B     B    
Chkder Checks the gradient calculation (for Hybrj and Hybrj1) M3     M3    
Sos Solution of a system of nonlinear equations (Brown's method) M3   M3 M3   M3
Sos_r Solution of a system of nonlinear equations (Brown's method) (reverse communication version) M3     M3    
G. Optimization                  
G1a. Unconstrained optimization of a general univariate function G1a. Unconstrained optimization of a general univariate function Dfmin Minimum of a single variable general nonlinear function B   B B   B
Dfmin_r Minimum of a single variable general nonlinear function (reverse communication version) B     B    
G1b. Unconstrained optimization of a general multivariate function G1b. Unconstrained optimization of a general multivariate function Optif9 Minimum of a multivariable nonlinear function (quasi-Newton method or trust region method) M3     M3    
Optif9_r Minimum of a multivariable nonlinear function (quasi-Newton method or trust region method) (reverse communication version) M3     M3    
Optif0 Minimum of a multivariable nonlinear function (quasi-Newton method) (simple driver) B   B B   B
Optif0_r Minimum of a multivariable nonlinear function (quasi-Newton method) (simple driver) (reverse communication version) B     B    
Mng Minimum of a multivariable nonlinear function (trust region method) M3   M3 M3   M3
Mng_r Minimum of a multivariable nonlinear function (trust region method) (reverse communication version) M3     M3    
Mnf Minimum of a multivariable nonlinear function (trust region method) (gradient computed by finite differences) M3   M3 M3   M3
Mnf_r Minimum of a multivariable nonlinear function (trust region method) (gradient computed by finite differences) (reverse communication version) M3     M3    
Mnh Minimum of a multivariable nonlinear function (trust region method) (gradient and Hessian computed analytically) M3     M3    
Mnh_r Minimum of a multivariable nonlinear function (trust region method) (gradient and Hessian computed analytically) (reverse communication version) M3     M3    
Subplex Minimum of a multivariable nonlinear function (subspace-searching simplex method) M3   M3 M3   M3
Subplex_r Minimum of a multivariable nonlinear function (subspace-searching simplex method) (reverse communication version) M3     M3    
G2. Constrained optimization of a general multivariate function G2. Constrained optimization of a general multivariate function Mngb Minimization of multivariate function (trust region method) (simply bounded) M3     M3    
Mngb_r Minimization of multivariate function (trust region method) (simply bounded) (reverse communication version) M3     M3    
Mnfb Minimization of multivariate function (trust region method) (simply bounded) (gradient computed by finite differences) M3     M3    
Mnfb_r Minimization of multivariate function (trust region method) (simply bounded) (gradient computed by finite differences) (reverse communication version) M3     M3    
Mnhb Minimization of multivariate function (trust region method) (simply bounded) (gradient and Hessian computed analytically) M3     M3    
Mnhb_r Minimization of multivariate function (trust region method) (simply bounded) (gradient and Hessian computed analytically) (reverse communication version) M3     M3    
H. Differentiation, integration                  
H2. Quadrature                  
H2a1a. 1-D finite interval quadrature (user-defined integrand function) H2a1a. 1-D finite interval quadrature (fixed number of points) Qk15 Finite interval quadrature (15-point Gauss-Kronrod rule) B     B    
Qk15_r Finite interval quadrature (15-point Gauss-Kronrod rule) (reverse communication version) B     B    
Qk21 Finite interval quadrature (21-point Gauss-Kronrod rule) M4     M4    
Qk21_r Finite interval quadrature (21-point Gauss-Kronrod rule) (reverse communication version) M4     M4    
Qk31 Finite interval quadrature (31-point Gauss-Kronrod rule) M4     M4    
Qk31_r Finite interval quadrature (31-point Gauss-Kronrod rule) (reverse communication version) M4     M4    
Qk41 Finite interval quadrature (41-point Gauss-Kronrod rule) M4     M4    
Qk41_r Finite interval quadrature (41-point Gauss-Kronrod rule) (reverse communication version) M4     M4    
Qk51 Finite interval quadrature (51-point Gauss-Kronrod rule) M4     M4    
Qk51_r Finite interval quadrature (51-point Gauss-Kronrod rule) (reverse communication version) M4     M4    
Qk61 Finite interval quadrature (61-point Gauss-Kronrod rule) M4     M4    
Qk61_r Finite interval quadrature (61-point Gauss-Kronrod rule) (reverse communication version) M4     M4    
H2a1a. 1-D finite interval quadrature (automatic quadrature) Qng Finite interval automatic quadrature (21/43/87-point Gauss-Kronrod rule) M4     M4    
Qng_r Finite interval automatic quadrature (21/43/87-point Gauss-Kronrod rule) (reverse communication version) M4     M4    
Qag Finite interval adaptive quadrature (15/21/31/41/51/61-point Gauss-Kronrod rule) B   B B   B
Qag_r Finite interval adaptive quadrature (15/21/31/41/51/61-point Gauss-Kronrod rule) (reverse communication version) B     B    
Qags Finite interval adaptive quadrature with sigularities (21-point Gauss-Kronrod rule) M4   M4 M4   M4
Qags_r Finite interval adaptive quadrature with sigularities (21-point Gauss-Kronrod rule) (reverse communication version) M4     M4    
Defin Finite interval automatic quadrature (double exponential (DE) formula) M4   M4 M4   M4
Defin_r Finite interval automatic quadrature (double exponential (DE) formula) (reverse communication version) M4     M4    
H2a1b. 1-D finite interval quadrature (tabulated integrand) H2a1b. 1-D finite interval quadrature (tabulated integrand) Avint (WAvint) Finite interval quadrature for a function with tabulated data (approximation with overlapping parabolas) M4 M4   M4 M4  
H2a2a. 1-D finite interval quadrature (special integrand) (user-defined integrand function) H2a2a. 1-D finite interval quadrature (special integrand) Qagp Finite interval adaptive quadrature with known singular points (21-point Gauss-Kronrod rule) M4     M4    
Qagp_r Finite interval adaptive quadrature with known singular points (21-point Gauss-Kronrod rule) (reverse communication version) M4     M4    
Qawc Finite interval adaptive quadrature for Cauchy principal values (25-point Clenshaw-Curtis and 15-point Gauss-Kronrod rule) M4   M4 M4   M4
Qawc_r Finite interval adaptive quadrature for Cauchy principal values (25-point Clenshaw-Curtis and 15-point Gauss-Kronrod rule) (reverse communication version) M4     M4    
Qaws Finite interval adaptive quadrature for singular functions (25-point Clenshaw-Curtis and 15-point Gauss-Kronrod rule) M4   M4 M4   M4
Qaws_r Finite interval adaptive quadrature for singular functions (25-point Clenshaw-Curtis and 15-point Gauss-Kronrod rule) (reverse communication version) M4     M4    
Qawo Finite interval adaptive quadrature for oscillatory functions (25-point Clenshaw-Curtis and 15-point Gauss-Kronrod rule) M4   M4 M4   M4
Qawo_r Finite interval adaptive quadrature for oscillatory functions (25-point Clenshaw-Curtis and 15-point Gauss-Kronrod rule) (reverse communication version) M4     M4    
H2a3a. 1-D semi-infinite interval quadrature (user-defined integrand function) H2a3a. 1-D semi-infinite interval quadrature Qk15i Semi-infinite/infinite interval quadrature (15-point Gauss-Kronrod rule) M4     M4    
Qk15i_r Semi-infinite/infinite interval quadrature (15-point Gauss-Kronrod rule) (reverse communication version) M4     M4    
Qagi Semi-infinite/infinite interval adaptive quadrature (15-point Gauss-Kronrod rule) B   B B   B
Qagi_r Semi-infinite/infinite interval adaptive quadrature (15-point Gauss-Kronrod rule) (reverse communication version) B     B    
Qawf Semi-infinite interval adaptive quadrature for Fourier integrals (25-point Clenshaw-Curtis and 15-point Gauss-Kronrod rule) M4   M4 M4   M4
Qawf_r Semi-infinite interval adaptive quadrature for Fourier integrals (25-point Clenshaw-Curtis and 15-point Gauss-Kronrod rule) (reverse communication version) M4     M4    
Dehint Semi-infinite interval automatic quadrature (double exponential (DE) formula) M4   M4 M4   M4
Dehint_r Semi-infinite interval automatic quadrature (double exponential (DE) formula) (reverse communication version) M4     M4    
Deoint Semi-infinite interval automatic quadrature for Fourier integrals (double exponential (DE) formula) M4   M4 M4   M4
Deoint_r Semi-infinite interval automatic quadrature for Fourier integrals (double exponential (DE) formula) (reverse communication version) M4     M4    
H2a4. 1-D infinite interval quadrature (user-defined integrand function) H2a4. 1-D infinite interval quadrature Deiint Infinite interval automatic quadrature (double exponential (DE) formula) M4   M4 M4   M4
Deiint_r Infinite interval automatic quadrature (double exponential (DE) formula) (reverse communication version) M4   M4  
I. Differential and integral equations                  
I1. Ordinary differential equations I1a1. Initial value problem of ordinary differential equations (for non-stiff problem) Derkfa Initial value problem of ordinary differential equations (5(4)-th order Runge-Kutta-Fehlberg method)       B   B
Derkfa_r Initial value problem of ordinary differential equations (5(4)-th order Runge-Kutta-Fehlberg method) (reverse communication version)       B    
Dopri5a Initial value problem of ordinary differential equations (5(4)-th order Dorman-Prince method)       M4   M4
Dopri5a_r Initial value problem of ordinary differential equations (5(4)-th order Dorman-Prince method) (reverse communication version)       M4    
Dverka Initial value problem of ordinary differential equations (6(5)-th order Runge-Kutta-Verner method)       M4    
Dverka_r Initial value problem of ordinary differential equations (6(5)-th order Runge-Kutta-Verner method) (reverse communication version)       M4    
Dop853a Initial value problem of ordinary differential equations (8(5,3)-th order Dorman-Prince method)       M4   M4
Dop853a_r Initial value problem of ordinary differential equations (8(5,3)-th order Dorman-Prince method) (reverse communication version)       M4    
Deabm Initial value problem of ordinary differential equations (1~12-th order Adams-Bashforth-Moulton predictor-corrector method) M4   M4 M4   M4
Deabm_r Initial value problem of ordinary differential equations (1~12-th order Adams-Bashforth-Moulton predictor-corrector method) (reverse communication version) M4     M4    
Odexa Initial value problem of ordinary differential equations (extrapolation method (GBS algorithm))       M4    
Odexa_r Initial value problem of ordinary differential equations (extrapolation method (GBS algorithm)) (reverse communication version)       M4    
Dopn43 Initial value problem of ordinary differential equations (4(3)-th order Runge-Kutta-Nystrom method) (for second order differential equations)       B   B
Dopn43_r Initial value problem of ordinary differential equations (4(3)-th order Runge-Kutta-Nystrom method) (for second order differential equations) (reverse communication version)       B    
Dopn64 Initial value problem of ordinary differential equations (6(4)-th order Runge-Kutta-Nystrom method) (for second order differential equations)       M4   M4
Dopn64_r Initial value problem of ordinary differential equations (6(4)-th order Runge-Kutta-Nystrom method) (for second order differential equations) (reverse communication version)       M4    
Dopn86 Initial value problem of ordinary differential equations (8(6)-th order Runge-Kutta-Nystrom method) (for second order differential equations)       M4   M4
Dopn86_r Initial value problem of ordinary differential equations (8(6)-th order Runge-Kutta-Nystrom method) (for second order differential equations) (reverse communication version)       M4    
Dopn1210 Initial value problem of ordinary differential equations (12(10)-th order Runge-Kutta-Nystrom method) (for second order differential equations)       M4   M4
Dopn1210_r Initial value problem of ordinary differential equations (12(10)-th order Runge-Kutta-Nystrom method) (for second order differential equations) (reverse communication version)       M4    
Odex2a Initial value problem of ordinary differential equations (extrapolation method) (for second order differential equations)       M4    
Odex2a_r Initial value problem of ordinary differential equations (extrapolation method) (for second order differential equations) (reverse communication version)       M4    
Retarda Initial value problem of delay differential equations (5(4)-th order Dorman-Prince method)       M4    
Ylaga Initial value problem of delay differential equations (5(4)-th order Dorman-Prince method) (interpolation for back-values of solution)       M4    
Retarda_r Initial value problem of delay differential equations (5(4)-th order Dorman-Prince method) (reverse communication version)       M4    
Ylaga_r Initial value problem of delay differential equations (5(4)-th order Dorman-Prince method) (reverse communication version) (interpolation for back-values of solution)       M4    
Derkf Initial value problem of ordinary differential equations (5(4)-th order Runge-Kutta-Fehlberg method) B   B B    
Derkf_r Initial value problem of ordinary differential equations (5(4)-th order Runge-Kutta-Fehlberg method) (reverse communication version) B     B    
DerkfInt Initial value problem of ordinary differential equations (5(4)-th order Runge-Kutta-Fehlberg method) (interpolation for dense output) B     B    
Dopri5 Initial value problem of ordinary differential equations (5(4)-th order Dorman-Prince method) M4   M4 M4    
Contd5 Initial value problem of ordinary differential equations (5(4)-th order Dorman-Prince method) (interpolation for dense output) M4     M4    
Dopri5_r Initial value problem of ordinary differential equations (5(4)-th order Dorman-Prince method) (reverse communication version) M4     M4    
Contd5_r Initial value problem of ordinary differential equations (5(4)-th order Dorman-Prince method) (reverse communication version) (interpolation for dense output) M4     M4    
Dverk Initial value problem of ordinary differential equations (6(5)-th order Runge-Kutta-Verner method) M4     M4    
Dverk_r Initial value problem of ordinary differential equations (6(5)-th order Runge-Kutta-Verner method) (reverse communication version) M4     M4    
DverkInt Initial value problem of ordinary differential equations (6(5)-th order Runge-Kutta-Verner method) (interpolation for dense output) M4     M4    
Dop853 Initial value problem of ordinary differential equations (8(5,3)-th order Dorman-Prince method) M4   M4 M4    
Contd8 Initial value problem of ordinary differential equations (8(5,3)-th order Dorman-Prince method) (interpolation for dense output) M4     M4    
Dop853_r Initial value problem of ordinary differential equations (8(5,3)-th order Dorman-Prince method) (reverse communication version) M4     M4    
Contd8_r Initial value problem of ordinary differential equations (8(5,3)-th order Dorman-Prince method) (reverse communication version) (interpolation for dense output) M4     M4    
Odex Initial value problem of ordinary differential equations (extrapolation method (GBS algorithm)) M4     M4    
Contx1 Initial value problem of ordinary differential equations (extrapolation method (GBS algorithm)) (interpolation for dense output) M4     M4    
Odex_r Initial value problem of ordinary differential equations (extrapolation method (GBS algorithm)) (reverse communication version) M4     M4    
Contx1_r Initial value problem of ordinary differential equations (extrapolation method (GBS algorithm)) (reverse communication version) (interpolation for dense output) M4     M4    
Doprin Initial value problem of ordinary differential equations (7(6)-th order Runge-Kutta-Nystrom method) (for second order differential equations) M4     M4    
Doprin_r Initial value problem of ordinary differential equations (7(6)-th order Runge-Kutta-Nystrom method) (for second order differential equations) (reverse communication version) M4     M4    
Odex2 Initial value problem of ordinary differential equations (extrapolation method) (for second order differential equations) M4     M4    
Contx2 Initial value problem of ordinary differential equations (extrapolation method) (for second order differential equations) (interpolation for dense output) M4     M4    
Odex2_r Initial value problem of ordinary differential equations (extrapolation method) (for second order differential equations) (reverse communication version) M4     M4    
Contx2_r Initial value problem of ordinary differential equations (extrapolation method) (for second order differential equations) (reverse communication version) (interpolation for dense output) M4     M4    
Retard Initial value problem of delay differential equations (5(4)-th order Dorman-Prince method) M4     M4    
Ylag Initial value problem of delay differential equations (5(4)-th order Dorman-Prince method) (interpolation for back-values of solution) M4     M4    
Retard_r Initial value problem of delay differential equations (5(4)-th order Dorman-Prince method) (reverse communication version) M4     M4    
Ylag_r Initial value problem of delay differential equations (5(4)-th order Dorman-Prince method) (reverse communication version) (interpolation for back-values of solution) M4     M4    
I1a2. Initial value problem of ordinary differential equations (for stiff problem) Debdf Initial value problem of ordinary differential equations (1~5-th order backward differentiation formula (BDF)) M4   M4 M4   M4
Debdf_r Initial value problem of ordinary differential equations (1~5-th order backward differentiation formula (BDF)) (reverse communication version) M4     M4    
Radaua Initial value problem of ordinary differential equations (variable (5, 9, 13-th) order implicit Runge-Kutta method (Radau IIA))       M4   M4
Radaua_r Initial value problem of ordinary differential equations (variable (5, 9, 13-th) order implicit Runge-Kutta method (Radau IIA)) (reverse communication version)       M4    
Rodasa Initial value problem of ordinary differential equations (4(3)-th order Rosenbrock method)       M4    
Rodasa_r Initial value problem of ordinary differential equations (4(3)-th order Rosenbrock method) (reverse communication version)       M4    
Seulexa Initial value problem of ordinary differential equations (extrapolation method based on the linearly implicit Euler method)       M4    
Seulexa_r Initial value problem of ordinary differential equations (extrapolation method based on the linearly implicit Euler method) (reverse communication version)       M4    
Radau5 Initial value problem of ordinary differential equations (5-th order implicit Runge-Kutta method (Radau IIA)) M4     M4    
Contr5 Initial value problem of ordinary differential equations (5-th order implicit Runge-Kutta method (Radau IIA)) (interpolation for dense output) M4     M4    
Radau5_r Initial value problem of ordinary differential equations (5-th order implicit Runge-Kutta method (Radau IIA)) (reverse communication version) M4     M4    
Contr5_r Initial value problem of ordinary differential equations (5-th order implicit Runge-Kutta method (Radau IIA)) (reverse communication version) (interpolation for dense output) M4     M4    
Radaup Initial value problem of ordinary differential equations (5, 9, 13-th order implicit Runge-Kutta method (Radau IIA)) M4     M4    
Contrp Initial value problem of ordinary differential equations (5, 9, 13-th order implicit Runge-Kutta method (Radau IIA)) (interpolation for dense output) M4     M4    
Radaup_r Initial value problem of ordinary differential equations (5, 9, 13-th order implicit Runge-Kutta method (Radau IIA)) (reverse communication version) M4     M4    
Contrp_r Initial value problem of ordinary differential equations (5, 9, 13-th order implicit Runge-Kutta method (Radau IIA)) (reverse communication version) (interpolation for dense output) M4     M4    
Radau Initial value problem of ordinary differential equations (variable (5, 9, 13-th) order implicit Runge-Kutta method (Radau IIA)) M4   M4 M4    
Contra Initial value problem of ordinary differential equations (variable (5, 9, 13-th) order implicit Runge-Kutta method (Radau IIA)) (interpolation for dense output) M4     M4    
Radau_r Initial value problem of ordinary differential equations (variable (5, 9, 13-th) order implicit Runge-Kutta method (Radau IIA)) (reverse communication version) M4     M4    
Contra_r Initial value problem of ordinary differential equations (variable (5, 9, 13-th) order implicit Runge-Kutta method (Radau IIA)) (reverse communication version) (interpolation for dense output) M4     M4    
Rodas Initial value problem of ordinary differential equations (4(3)-th order Rosenbrock method) M4     M4    
Contro Initial value problem of ordinary differential equations (4(3)-th order Rosenbrock method) (interpolation for dense output) M4     M4    
Rodas_r Initial value problem of ordinary differential equations (4(3)-th order Rosenbrock method) (reverse communication version) M4     M4    
Contro_r Initial value problem of ordinary differential equations (4(3)-th order Rosenbrock method) (reverse communication version) (interpolation for dense output) M4     M4    
Seulex Initial value problem of ordinary differential equations (extrapolation method based on the linearly implicit Euler method) M4     M4    
Contex Initial value problem of ordinary differential equations (extrapolation method based on the linearly implicit Euler method) (interpolation for dense output) M4     M4    
Seulex_r Initial value problem of ordinary differential equations (extrapolation method based on the linearly implicit Euler method) (reverse communication version) M4     M4    
Contex_r Initial value problem of ordinary differential equations (extrapolation method based on the linearly implicit Euler method) (reverse communication version) (interpolation for dense output) M4     M4    
Dassl Solution of differential algebraic equation (DAE) (1~5-th order backward differentiation formula (BDF)) M4     M4    
Dassl_r Solution of differential algebraic equation (DAE) (1~5-th order backward differentiation formula (BDF)) (reverse communication version) M4     M4    
J. Integral transforms                  
J1. Fast Fourier transform (FFT) J1a1. One-dimensional real fast Fourier transforms Rfft1f (WRfft1f) One-dimensional real Fourier transform B B   B B  
Rfft1b (WRfft1b) One-dimensional real Fourier backward transform B B   B B  
Rfft1i Initialization of work data for Rfft1f and Rfft1b B     B    
Rfftmf One-dimensional real Fourier transform (multiple sequences) M3     M3    
Rfftmb One-dimensional real Fourier backward transform (multiple sequences) M3     M3    
Rfftmi Initialization of work data for Rfftmf and Rfftmb M3     M3    
J1a2. One-dimensional complex fast Fourier transforms Cfft1f (WCfft1f(2)) One-dimensional complex Fourier transform M3 M3   M3 M3  
Cfft1b (WCfft1b(2)) One-dimensional complex Fourier backward transform M3 M3   M3 M3  
Cfft1i Initialization of work data for Cfft1f and Cfft1b M3     M3    
Cfftmf One-dimensional complex Fourier transform (multiple sequences) M3     M3    
Cfftmb One-dimensional complex Fourier backward transform (multiple sequences) M3     M3    
Cfftmi Initialization of work data for Cfftmf and Cfftmb M3     M3    
J1a3. One-dimensional trigonometric fast Fourier transforms Sint1f (WSint1f) One-dimensional real sine transform M3 M3   M3 M3  
Sint1b (WSint1b) One-dimensional real sine backward transform M3 M3   M3 M3  
Sint1i Initialization of work data for Sint1f and Sint1b M3     M3    
Sintmf One-dimensional real sine transform (multiple sequences) M3     M3    
Sintmb One-dimensional real sine backward transform (multiple sequences) M3     M3    
Sintmi Initialization of work data for Sintmf and Sintmb M3     M3    
Cost1f (WCost1f) One-dimensional real cosine transform M3 M3   M3 M3  
Cost1b (WCost1b) One-dimensional real cosine backward transform M3 M3   M3 M3  
Cost1i Initialization of work data for Cost1f and Cost1b M3     M3    
Costmf One-dimensional real cosine transform (multiple sequences) M3     M3    
Costmb One-dimensional real cosine backward transform (multiple sequences) M3     M3    
Costmi Initialization of work data for Costmf and Costmb M3     M3    
J1a3. One-dimensional quarter trigonometric fast Fourier transforms Sinq1f One-dimensional real quarter sine transform M3     M3    
Sinq1b One-dimensional real quarter sine backward transform M3     M3    
Sinq1i Initialization of work data for Sinq1f and Sinq1b M3     M3    
Sinqmf One-dimensional real quarter sine transform (multiple sequences) M3     M3    
Sinqmb One-dimensional real quarter sine backward transform (multiple sequences) M3     M3    
Sinqmi Initialization of work data for Sinqmf and Sinqmb M3     M3    
Cosq1f One-dimensional real quarter cosine transform M3     M3    
Cosq1b One-dimensional real quarter cosine backward transform M3     M3    
Cosq1i Initialization of work data for Cosq1f and Cosq1b M3     M3    
Cosqmf One-dimensional real quarter cosine transform (multiple sequences) M3     M3    
Cosqmb One-dimensional real quarter cosine backward transform (multiple sequences) M3     M3    
Cosqmi Initialization of work data for Cosqmf and Cosqmb M3     M3    
J1b. Multidimensional fast Fourier transforms Rfft2f Two-dimensional real Fourier transform M3     M3    
Rfft2b Two-dimensional real Fourier backward transform M3     M3    
Rfft2i Initialization of work data for Rfft2f and Rfft2b M3     M3    
Rfft2c Full complex data of two-dimensional Fourier transform obtained by Rfft2f M3     M3    
Cfft2f Two-dimensional complex Fourier transform M3     M3    
Cfft2b Two-dimensional complex Fourier backward transform M3     M3    
Cfft2i Initialization of work data for Cfft2f and Cfft2b M3     M3    
K. Approximation                  
K1. Least squares approximation K1b1. Nonlinear least squares approximation Lmder Nonlinear least squares approximation (Levenberg-Marquardt method) M3     M3    
Lmder_r Nonlinear least squares approximation (Levenberg-Marquardt method) (reverse communication version) M3     M3    
Lmder1 Nonlinear least squares approximation (Levenberg-Marquardt method) (simple driver) M3   M3 M3   M3
Lmder1_r Nonlinear least squares approximation (Levenberg-Marquardt method) (simple driver) (reverse communication version) M3     M3    
Lmstr Nonlinear least squares approximation (Levenberg-Marquardt method) (limited storage version) M3     M3    
Lmstr_r Nonlinear least squares approximation (Levenberg-Marquardt method) (limited storage version) (reverse communication version) M3     M3    
Lmstr1 Nonlinear least squares approximation (Levenberg-Marquardt method) (limited storage version) (simple driver) M3     M3    
Lmstr1_r Nonlinear least squares approximation (Levenberg-Marquardt method) (limited storage version) (simple driver) (reverse communication version) M3     M3    
Lmdif Nonlinear least squares approximation (Levenberg-Marquardt method) (Jacobian not required) M3     M3    
Lmdif_r Nonlinear least squares approximation (Levenberg-Marquardt method) (Jacobian not required) (reverse communication version) M3     M3    
Lmdif1 Nonlinear least squares approximation (Levenberg-Marquardt method) (Jacobian not required) (simple driver) B   B B   B
Lmdif1_r Nonlinear least squares approximation (Levenberg-Marquardt method) (Jacobian not required) (simple driver) (reverse communication version) B     B    
Chkder Checks the gradient calculation (for Lmder, Lmder1, Lmstr and Lmstr1) (same as F2.) M3     M3    
Covar Variance covariance matrix calculation for Lmder, Lmder1, Lmstr, Lmstr1 and Lmdif) M3     M3    
N2g Nonlinear least squares approximation (Levenberg-Marquardt method) M3   M3 M3   M3
N2g_r Nonlinear least squares approximation (Levenberg-Marquardt method) (reverse communication version) M3     M3    
N2g1 Nonlinear least squares approximation (Levenberg-Marquardt method) (simple driver) M3     M3    
N2g1_r Nonlinear least squares approximation (Levenberg-Marquardt method) (simple driver) (reverse communication version) M3     M3    
N2f Nonlinear least squares approximation (adaptive algorithm) (Jacobian not required) M3   M3 M3   M3
N2f_r Nonlinear least squares approximation (adaptive algorithm) (Jacobian not required) (reverse communication version) M3     M3    
N2f1 Nonlinear least squares approximation (adaptive algorithm) (Jacobian not required) (simple driver) M3     M3    
N2f1_r Nonlinear least squares approximation (adaptive algorithm) (Jacobian not required) (simple driver) (reverse communication version) M3     M3    
N2p Nonlinear least squares approximation (adaptive algorithm) (limited storage version) M3     M3    
N2p_r Nonlinear least squares approximation (adaptive algorithm) (limited storage version) (reverse communication version) M3     M3    
K1b2. Constrained nonlinear least squares approximation N2gb Nonlinear least squares approximation (adaptive algorithm) (simply bounded) M3     M3    
N2gb_r Nonlinear least squares approximation (adaptive algorithm) (simply bounded) (reverse communication version) M3     M3    
N2fb Nonlinear least squares approximation (adaptive algorithm) (simply bounded) (Jacobian not required) M3     M3    
N2fb_r Nonlinear least squares approximation (adaptive algorithm) (simply bounded) (Jacobian not required) (reverse communication version) M3     M3    
N2pb Nonlinear least squares approximation (adaptive algorithm) (simply bounded) (limited storage version) M3     M3    
N2pb_r Nonlinear least squares approximation (adaptive algorithm) (simply bounded) (limited storage version) (reverse communication version) M3     M3    
L. Statistics, probability                  
L6. Random number generation                  
L6a21. Uniform random numbers L6a21. Uniform random numbers (Mersenne-Twister) InitGenrand Initialization with seed for random number generator (Mersenne Twister) B     B    
InitByArray Initialization with array of integers for random number generator (Mersenne Twister) M4     M4    
GenrandInt32 (#) Unsigned 32 bit integer random number (Mersenne Twister) B     B    
GenrandInt31 (#) Unsigned 31 bit integer random number (Mersenne Twister) B     B    
GenrandReal1 (#) 32 bit real random number in [0,1] (Mersenne Twister) M4     M4    
GenrandReal2 (#) 32 bit real random number in [0,1) (Mersenne Twister) M4     M4    
GenrandReal3 (#) 32 bit real random number in (0,1) (Mersenne Twister) M4     M4    
GenrandReal53 (#) 53 bit real random number in [0,1) (Mersenne Twister) B     B    
InitGenrand64 Initialization of random number generator (64 bit Mersenne Twister) M4@     M4@    
InitByArray64 Initialization with array of integers for random number generator (64 bit Mersenne Twister) M4@     M4@    
Genrand64Int64 (#) Unsigned 64 bit integer random number (64 bit Mersenne Twister) M4@     M4@    
Genrand64Int63 (#) Unsigned 63 bit integer random number (64 bit Mersenne Twister) M4@     M4@    
Genrand64Real1 (#) Double precision real random number in [0, 1] (64 bit Mersenne Twister) M4@     M4@    
Genrand64Real2 (#) Double precision real random number in [0, 1) (64 bit Mersenne Twister) M4@     M4@    
Genrand64Real3 (#) Double precision real random number in (0, 1) (64 bit Mersenne Twister) M4@     M4@    
L6a21. Uniform random numbers (Lagged Fibonacci method) RanStart Initialization for integer random number generator (Lagged Fibonacci method) M4     M4    
RanArray Unsigned 30 bit integer random numbers (Lagged Fibonacci method) M4     M4    
RanArrNext (#) Unsigned 30 bit integer random number (Lagged Fibonacci method) M4     M4    
RanfStart Initialization for real random number generator (Lagged Fibonacci method) M4     M4    
RanfArray 53 bit real random numbers in [0,1) (Lagged Fibonacci method) M4     M4    
RanfArrNext (#) 53 bit real random number in [0,1) (Lagged Fibonacci method) M4     M4    
L6a21. Uniform random numbers (Linear congruential method) Srand48 Initialization with 32-bit seed for Drand48, Lrand48 and Mrand48 (Linear congruential method) M4     M4    
Seed48 Initialization with 48-bit seed for Drand48, Lrand48 and Mrand48 (Linear congruential method) M4     M4    
Lcong48 Set up parameters for random number generators (Linear congruential method) M4     M4    
Drand48 (#) 48 bit real random number in [0,1) (Linear congruential method) M4     M4    
Erand48 48 bit real random number in [0,1) (Linear congruential method) M4     M4    
Lrand48 (#) Unsigned 31 bit integer random number (Linear congruential method) M4     M4    
Nrand48 Unsigned 31 bit integer random number (Linear congruential method) M4     M4    
Mrand48 (#) Signed 32 bit integer random number (Linear congruential method) M4     M4    
Jrand48 Signed 32 bit integer random number (Linear congruential method) M4     M4    
L6a14. Normal random numbers L6a14. Normal random numbers InitNorm Initialization of normal random number generator (Ziggurat method) M4     M4    
GenrandNorm (#) 53 bit real normal random number (Ziggurat method) (Mersenne Twister) M4   M4  
RanfArrNextNorm (#) 53 bit real normal random number (Ziggurat method) (Lagged Fibonacci method) M4     M4    
Drand48Norm (#) 48 bit real normal random number (Ziggurat method) (Linear congruential method) M4     M4    
L6a5. Exponential random numbers L6a5. Exponential random numbers InitExp Initialization of exponential random number generator (Ziggurat method) M4     M4    
GenrandExp (#) 53 bit real exponential random number (Ziggurat method) (Mersenne Twister) M4   M4  
RanfArrNextExp (#) 53 bit real exponential random number (Ziggurat method) (Lagged Fibonacci method) M4     M4    
Drand48Exp (#) 48 bit real exponential random number (Ziggurat method) (Linear congruential method) M4     M4    
L6a7. Gamma random numbers L6a7. Gamma random numbers.  GenrandGam (#) 53 bit real gamma random number (Squeeze method of Marsaglia and Tsang) (Mersenne Twister) M4     M4    
RanfArrNextGam (#) 53 bit real gamma random number (Squeeze method of Marsaglia and Tsang) (Lagged Fibonacci method) M4     M4    
Drand48Gam (#) 48 bit real gamma random number (Squeeze method of Marsaglia and Tsang) (Linear congruential method) M4     M4    
R. Service routines                  
R1. Machine-dependent constants R1. Machine-dependent constants Dlamch (#) Machine parameters (double precision floating-point arithmetic) B     B    
D1mach (#) Machine parameters (double precision floating-point arithmetic) B     B    
Slamch (#) Machine parameters (single precision floating-point arithmetic) B     B    
R1mach (#) Machine parameters (single precision floating-point arithmetic) B     B    
I1mach (#) Machine parameters (integer machine dependent constants) B     B    
Z. Others                  
Z1. Test matrix generation Z1. Test matrix generation Dlatms Generates random matrices with specified singular values or symmetric random matrices with specified eigenvalues M1     M1    
Dlatmt Generates random matrices with specified singular values or symmetric random matrices with specified eigenvalues (with specified rank of matrix) M1     M1    
Dlatme Generates random non-symmetric square matrices with specified eigenvalues M1     M1    
Dlatmr Generates random matrices with specified diagonal elements M1     M1    
Zlatms Generates random matrices with specified singular values or symmetric random matrices with specified eigenvalues (complex matrix) M2     M2    
Zlatmt Generates random matrices with specified singular values or symmetric random matrices with specified eigenvalues (with specified rank of matrix) (complex matrix) M2     M2    
Zlatme Generates random non-symmetric square matrices with specified eigenvalues (complex matrix) M2     M2    
Zlatmr Generates random matrices with specified diagonal elements (complex matrix) M2     M2    
Category (Sparse matrix calculation) VBA routine name Functions V6.1 V7.0
VBA WS Solver VBA WS Solver
Sparse BLAS D1a. Sparse BLAS CsrDusmv y <- αAx + βy or y <- αA^Tx + βy (CSR)       B    
CsrDussv Solution of Ax = b or A^Tx = b (Triangular matrices) (CSR)       B    
CsrDusmm C <- αAB + βC or C <- αA^TB + βC (CSR)       B    
CsrDussm Solution of AX = B or A^TX = B (Triangular matrices) (CSR)       B    
SsrDusmv y <- αAx + βy (Symmetric matrix) (CSR)       B    
CscDusmv y <- αAx + βy or y <- αA^Tx + βy (CSC)       M5    
CscDussv Solution of Ax = b or A^Tx = b (Triangular matrices) (CSC)       M5    
CscDusmm C <- αAB + βC or C <- αA^TB + βC (CSC)       M5    
CscDussm Solution of AX = B or A^TX = B (Triangular matrices) (CSC)       M5    
SscDusmv y <- αAx + βy (Symmetric matrix) (CSC)       M5    
D1a. Sparse BLAS (Complex) CsrZusmv y <- αAx + βy, y <- αA^Tx + βy or y <- αA^Hx + βy (Complex matrices) (CSR)       M5    
CsrZussv Solution of Ax = b, A^Tx = b or A^Hx = b (Complex triangular matrices) (CSR)       M5    
CsrZusmm C <- αAB + βC, C <- αA^TB + βC or C <- αA^HB + βC (Complex matrices) (CSR)       M5    
CsrZussm Solution of AX = B, A^TX = B or A^HX = B (Complex triangular matrices) (CSR)       M5    
HsrZusmv y <- αAx + βy (Hermitian matrix) (CSR)       M5    
SsrZusmv y <- αAx + βy (Complex symmetric matrix) (CSR)       M5    
CscZusmv y <- αAx + βy, y <- αA^Tx + βy or y <- αA^Hx + βy (Complex matrices) (CSC)       M5    
CscZussv Solution of Ax = b, A^Tx = b or A^Hx = b (Complex triangular matrices) (CSC)       M5    
CscZusmm C <- αAB + βC, C <- αA^TB + βC or C <- αA^HB + βC (Complex matrices) (CSC)       M5    
CscZussm Solution of AX = B, A^TX = B or A^HX = B (Complex triangular matrices) (CSC)       M5    
HscZusmv y <- αAx + βy (Hermitian matrix) (CSC)       M5    
SscZusmv y <- αAx + βy (Complex symmetric matrix) (CSC)       M5    
Elementary vector and matrix operations D1b9. Matrix storage mode conversion CooCsr COO -> CSR       B    
CooCsc COO -> CSC       M5    
CsrCoo CSR -> COO       B    
CscCoo CSC -> COO       M5    
CsrCsc CSR -> CSC       M5    
CscCsr CSC -> CSR       M5    
SsrCsr SSR (CSR sparse matrix packed form) -> CSR (symmetric full matrix)       B    
SscCsc SSC (CSC sparse matrix packed form) -> CSC (symmetric full matrix)       M5    
CsrSsr CSR (symmetric full matrix) -> SSR (CSR sparse matrix packed form)       B    
CscSsc CSC (symmetric full matrix) -> SSC (CSC sparse matrix packed form)       M5    
CsrDense CSR -> dense matrix       B    
CscDense CSC -> dense matrix       M5    
CooDense COO -> dense matrix       M5    
DenseCsr Dense matrix -> CSR       B    
DenseCsc Dense matrix -> CSC       M5    
DenseCoo Dense matrix -> COO       M5    
D1b. Other matrix operations CsxDiag Diagonal elements of sparse matrix (CSC/CSR)       M5    
CsxDiagInd Indices to diagonal elements of sparse matrix (CSC/CSR)       M5    
CsrTrans Transpose of sparse matrix (CSR)       B    
CscTrans Transpose of sparse matrix (CSC)       M5    
CsxSort Sort elements of sparse matrix (CSC/CSR)       M5    
CsrDusadd C <- αA + βB (CSR)       M5    
CscDusadd C <- αA + βB (CSC)       M5    
DenseNnz Number of nonzero elements in dense matrix       M5    
D1b9. Matrix storage mode conversion (Complex matrices) ZCooCsr COO -> CSR (Complex matrices)       M5    
ZCooCsc COO -> CSC (Complex matrices)       M5    
ZCsrCoo CSR -> COO (Complex matrices)       M5    
ZCscCoo CSC -> COO (Complex matrices)       M5    
ZCsrCsc CSR -> CSC (Complex matrices)       M5    
ZCscCsr CSC -> CSR (Complex matrices)       M5    
ZHsrCsr HSR (CSR Hermitian sparse matrix packed form) -> CSR (Hermitian full matrix) (Complex matrices)       M5    
ZHscCsc HSC (CSC Hermitian sparse matrix packed form) -> CSC (Hermitian full matrix) (Complex matrices)       M5    
ZSsrCsr SSR (CSR sparse matrix packed form) -> CSR (symmetric full matrix) (Complex matrices)       M5    
ZSscCsc SSC (CSC sparse matrix packed form) -> CSC (symmetric full matrix) (Complex matrices)       M5    
ZCsrSsr CSR (symmetric full matrix) -> SSR (CSR sparse matrix packed form) (Complex matrices)       M5    
ZCscSsc CSC (symmetric full matrix) -> SSC (CSC sparse matrix packed form) (Complex matrices)       M5    
ZCsrDense CSR -> dense matrix (Complex matrices)       M5    
ZCscDense CSC -> dense matrix (Complex matrices)       M5    
ZCooDense COO -> dense matrix (Complex matrices)       M5    
ZDenseCsr Dense matrix -> CSR (Complex matrices)       M5    
ZDenseCsc Dense matrix -> CSC (Complex matrices)       M5    
ZDenseCoo Dense matrix -> COO (Complex matrices)       M5    
D1b. Other matrix operations (Complex matrices) ZCsxDiag Diagonal elements of sparse matrix (Complex matrices) (CSC/CSR)       M5    
ZCsrTrans Transpose of sparse matrix (Complex matrices) (CSR)       M5    
ZCscTrans Transpose of sparse matrix (Complex matrices) (CSC)       M5    
ZCsxSort Sort elements of sparse matrix (Complex matrices) (CSC/CSR)       M5    
CsrZusadd C <- αA + βB (Complex matrices) (CSR)       M5    
CscZusadd C <- αA + βB (Complex matrices) (CSC)       M5    
CsrDzusadd C <- αA + βB (α and C are complex numbers) (CSR)       M5    
CscDzusadd C <- αA + βB (α and C are complex numbers) (CSC)       M5    
ZDenseNnz Number of nonzero elements in dense matrix (Complex matrices)       M5    
Solution of systems of linear equations (iterative solvers) D2a. General matrices Bicg1 BICG method (simplex driver)       B    
Bicg BICG method (driver)       M5    
Bicg_r BICG method (Reverse communication version)       M5    
Bicg_s BICG method (Subroutine version)       M5    
Cgs CGS method (driver)       M5    
Cgs_r CGS method (Reverse communication version)       M5    
Cgs_s CGS method (Subroutine version)       M5    
Diom DIOM (driver)       M5    
Diom_r DIOM (Reverse communication version)       M5    
Diom_s DIOM (Subroutine version)       M5    
Dqgmres DQGMRES method (driver)       M5    
Dqgmres_r DQGMRES method (Reverse communication version)       M5    
Dqgmres_s DQGMRES method (Subroutine version)       M5    
Fgmres FGMRES method (driver)       M5    
Fgmres_r FGMRES method (Reverse communication version)       M5    
Fgmres_s FGMRES method (Subroutine version)       M5    
Fom FOM (driver)       M5    
Fom_r FOM (Reverse communication version)       M5    
Fom_s FOM (Subroutine version)       M5    
Gcr GCR method (driver)       M5    
Gcr_r GCR method (Reverse communication version)       M5    
Gcr_s GCR method (Subroutine version)       M5    
Gpbicg GPBICG method (driver)       M5    
Gpbicg_r GPBICG method (Reverse communication version)       M5    
Gpbicg_s GPBICG method (Subroutine version)       M5    
Orthomin ORTHOMIN method (driver)       M5    
Orthomin_r ORTHOMIN method (Reverse communication version)       M5    
Orthomin_s ORTHOMIN method (Subroutine version)       M5    
Qmr QMR method (driver)       M5    
Qmr_r QMR method (Reverse communication version)       M5    
Qmr_s QMR method (Subroutine version)       M5    
Tfqmr TFQMR method (driver)       M5    
Tfqmr_r TFQMR method (Reverse communication version)       M5    
Tfqmr_s TFQMR method (Subroutine version)       M5    
Sor SOR method (driver)       M5    
Sor_r SOR method (Reverse communication version)       M5    
Sor_s SOR method (Subroutine version)       M5    
D2a. General matrices (Preconditioners) CsxDs Initialize diagonal scaling preconditioner (CSC/CSR)       M5    
CsxDsSolve Diagonal scaling preconditioner (CSC/CSR)       M5    
CsxSsor Initialize symmetric successive over-relaxation (SSOR) preconditioner (CSC/CSR)       M5    
CsrSsorSolve Symmetric successive over-relaxation (SSOR) preconditioner (CSR)       M5    
CscSsorSolve Symmetric successive over-relaxation (SSOR) preconditioner (CSC)       M5    
CsrIlu0 Initialize incomplete LU decomposition without fill-in (ILU0) preconditioner (CSR)       M5    
CscIlu0 Initialize incomplete LU decomposition without fill-in (ILU0) preconditioner (CSC)       M5    
CsrIlu Initialize incomplete LU decomposition with level (ILU(p)) preconditioner (CSR)       M5    
CscIlu Initialize incomplete LU decomposition with level (ILU(p)) preconditioner (CSC)       M5    
CsrIluSolve Incomplete LU decomposition preconditioner (ILU) (CSR)       M5    
CscIluSolve Incomplete LU decomposition preconditioner (ILU) (CSC)       M5    
D2b. Symmetric matrices Cg1 CG method (symmetric positive definite matrix) (simple driver)       B    
Cg CG method (symmetric positive definite matrix) (driver)       M5    
Cg_r CG method (symmetric positive definite matrix) (Reverse communication version)       M5    
Cg_s CG method (symmetric positive definite matrix) (Subroutine version)       M5    
Cr CR method (symmetric matrix) (driver)       M5    
Cr_r CR method (symmetric matrix) (Reverse communication version)       M5    
Cr_s CR method (symmetric matrix) (Subroutine version)       M5    
D2a. Symmetric matrices (Preconditioners) SsrIc0 Initialize incomplete Cholesky decomposition without fill-in (IC0) preconditioner (symmetric positive definite matrix) (CSR)       M5    
SscIc0 Initialize incomplete Cholesky decomposition without fill-in (IC0) preconditioner (symmetric positive definite matrix) (CSC)       M5    
SsrIcSolve Incomplete Cholesky decomposition preconditioner (IC) (symmetric positive definite matrix) (CSR)       M5    
SscIcSolve Incomplete Cholesky decomposition preconditioner (IC) (symmetric positive definite matrix) (CSC)       M5    
SsrSsorSolve Symmetric successive over-relaxation (SSOR) preconditioner (Symmetric matrix) (CSR)       M5    
SscSsorSolve Symmetric successive over-relaxation (SSOR) preconditioner (Symmetric matrix) (CSC)       M5    
D2c. Complex general matrices ZBicg BICG method (Complex matrices) (driver)       M5    
ZBicg_r BICG method (Complex matrices) (Reverse communication version)       M5    
ZBicg_s BICG method (Complex matrices) (Subroutine version)       M5    
ZCgs CGS method (Complex matrices) (driver)       M5    
ZCgs_r CGS method (Complex matrices) (Reverse communication version)       M5    
ZCgs_s CGS method (Complex matrices) (Subroutine version)       M5    
ZDiom DIOM (Complex matrices) (driver)       M5    
ZDiom_r DIOM (Complex matrices) (Reverse communication version)       M5    
ZDiom_s DIOM (Complex matrices) (Subroutine version)       M5    
ZDqgmres DQGMRES method (Complex matrices) (driver)       M5    
ZDqgmres_r DQGMRES method (Complex matrices) (Reverse communication version)       M5    
ZDqgmres_s DQGMRES method (Complex matrices) (Subroutine version)       M5    
ZFgmres FGMRES method (Complex matrices) (driver)       M5    
ZFgmres_r FGMRES method (Complex matrices) (Reverse communication version)       M5    
ZFgmres_s FGMRES method (Complex matrices) (Subroutine version)       M5    
ZFom FOM (Complex matrices) (driver)       M5    
ZFom_r FOM (Complex matrices) (Reverse communication version)       M5    
ZFom_s FOM (Complex matrices) (Subroutine version)       M5    
ZGcr GCR method (Complex matrices) (driver)       M5    
ZGcr_r GCR method (Complex matrices) (Reverse communication version)       M5    
ZGcr_s GCR method (Complex matrices) (Subroutine version)       M5    
ZGpbicg GPBICG method (Complex matrices) (driver)       M5    
ZGpbicg_r GPBICG method (Complex matrices) (Reverse communication version)       M5    
ZGpbicg_s GPBICG method (Complex matrices) (Subroutine version)       M5    
ZOrthomin ORTHOMIN method (Complex matrices) (driver)       M5    
ZOrthomin_r ORTHOMIN method (Complex matrices) (Reverse communication version)       M5    
ZOrthomin_s ORTHOMIN method (Complex matrices) (Subroutine version)       M5    
ZQmr QMR method (Complex matrices) (driver)       M5    
ZQmr_r QMR method (Complex matrices) (Reverse communication version)       M5    
ZQmr_s QMR method (Complex matrices) (Subroutine version)       M5    
ZTfqmr TFQMR method (Complex matrices) (driver)       M5    
ZTfqmr_r TFQMR method (Complex matrices) (Reverse communication version)       M5    
ZTfqmr_s TFQMR method (Complex matrices) (Subroutine version)       M5    
ZSor SOR method (Complex matrices) (driver)       M5    
ZSor_r SOR method (Complex matrices) (Reverse communication version)       M5    
ZSor_s SOR method (Complex matrices) (Subroutine version)       M5    
D2c. Complex general matrices (Preconditioners) ZCsxDs Initialize diagonal scaling preconditioner (Complex matrices) (CSC/CSR)       M5    
ZCsxDsSolve Diagonal scaling preconditioner (Complex matrices) (CSC/CSR)       M5    
ZCsxSsor Initialize symmetric successive over-relaxation (SSOR) preconditioner (Complex matrices) (CSC/CSR)       M5    
ZCsrSsorSolve Symmetric successive over-relaxation (SSOR) preconditioner (Complex symmetric matrices) (CSR)       M5    
ZCscSsorSolve Symmetric successive over-relaxation (SSOR) preconditioner (Complex symmetric matrices) (CSC)       M5    
ZCsrIlu0 Initialize incomplete LU decomposition without fill-in (ILU0) preconditioner (Complex matrices) (CSR)       M5    
ZCscIlu0 Initialize incomplete LU decomposition without fill-in (ILU0) preconditioner (Complex matrices) (CSC)       M5    
ZCsrIlu Initialize incomplete LU decomposition with level (ILU(p)) preconditioner (Complex matrices) (CSR)       M5    
ZCscIlu Initialize incomplete LU decomposition with level (ILU(p)) preconditioner (Complex matrices) (CSC)       M5    
ZCsrIluSolve Incomplete LU decomposition (ILU) preconditioner (Complex matrices) (CSR)       M5    
ZCscIluSolve Incomplete LU decomposition (ILU) preconditioner (Complex matrices) (CSC)       M5    
D2c. Complex symmetric matrices ZCocg COCG method (Complex symmetric matrices) (driver)       M5    
ZCocg_r COCG method (Complex symmetric matrices) (Reverse communication version)       M5    
ZCocg_s COCG method (Complex symmetric matrices) (Subroutine version)       M5    
Zcocr COCR method (Complex symmetric matrices) (driver)       M5    
ZCocr_r COCR method (Complex symmetric matrices) (Reverse communication version)       M5    
ZCocr_s COCR method (Complex symmetric matrices) (Subroutine version)       M5    
D2c. Complex symmetric matrices (Preconditioners) ZSsrSsorSolve Symmetric successive over-relaxation (SSOR) preconditioner (Complex symmetric matrices) (CSR)       M5    
ZSscSsorSolve Symmetric successive over-relaxation (SSOR) preconditioner (Complex symmetric matrices) (CSC)       M5    
D2d. Hermitian matrices ZCg CG method (Hermitian positive definite matrices) (driver)       M5    
ZCg_r CG method (Hermitian positive definite matrices) (Reverse communication version)       M5    
ZCg_s CG method (Hermitian positive definite matrices) (Subroutine version)       M5    
ZCr CR method (Hermitian matrices) (driver)       M5    
ZCr_r CR method (Hermitian matrices) (Reverse communication version)       M5    
ZCr_s CR method (Hermitian matrices) (Subroutine version)       M5    
D2d. Hermitian matrices (Preconditioners) ZHsrIc0 Initialize incomplete Cholesky decomposition without fill-in (IC0) preconditioner (Hermitian positive definite matrices) (CSR)       M5    
ZHscIc0 Initialize incomplete Cholesky decomposition without fill-in (IC0) preconditioner (Hermitian positive definite matrices) (CSC)       M5    
ZHsrIcSolve Incomplete Cholesky preconditioner (IC) (Hermitian positive definite matrix) (CSR)       M5    
ZHscIcSolve Incomplete Cholesky preconditioner (IC) (Hermitian positive definite matrix) (CSC)       M5    
ZHsrSsorSolve Symmetric successive over-relaxation (SSOR) preconditioner (Hermitian matrices) (CSR)       M5    
ZHscSsorSolve Symmetric successive over-relaxation (SSOR) preconditioner (Hermitian matrices) (CSC)       M5    
Solution of systems of linear equations (direct methods) D2. SuperLU (General matrices) Dgssv Solves the system of linear equations (direct method) (sparse matrix) (SuperLU) (simple driver)       M5    
Dgssvx Solves the system of linear equations (direct method) (sparse matrix) (SuperLU) (expert driver)       M5    
D2. SuperLU (Complex matrices) Zgssv Solves the system of linear equations (direct method) (complex sparse matrix) (SuperLU) (simple driver)       M5    
Zgssvx Solves the system of linear equations (direct method) (complex sparse matrix) (SuperLU) (expert driver)       M5    
Eigenvalues and eigenvectors (sparse matrices) D4. Arpack (General sparse matrices) Dgsev Eigenvalues and eigenvectors of a general sparse matrix (implicitly restarted Arnoldi method (IRAM)) (Arpack) (driver)       M5    
Dgsgv Generalized eigenvalue problem of a general sparse matrix (implicitly restarted Arnoldi method (IRAM)) (Arpack) (driver)       M5    
Dnaupd Arnoldi factorization of a general sparse matrix (Arpack)       M5    
Dneupd Approximate eigenvalues and eigenvectors of a general sparse matrix from Arnoldi factorization (Arpack)       M5    
D4. Arpack (Symmetric sparse matrices) Dssev Eigenvalues and eigenvectors of a symmetric sparse matrix (implicitly restarted Lanczos method (IRLM)) (Arpack) (driver)       M5    
Dssgv Generalized eigenvalue problem of a symmetric sparse matrix (implicitly restarted Lanczos method (IRLM)) (Arpack) (driver)       M5    
Dsaupd Lanczos factorization of a symmetric sparse matrix (Arpack)       M5    
Dseupd Approximate eigenvalues and eigenvectors of a symmetric sparse matrix from Lanczos factorization (Arpack)       M5    
D4. Arpack (Complex sparse matrices) Zgsev Eigenvalues and eigenvectors of a complex sparse matrix (implicitly restarted Arnoldi method (IRAM)) (Arpack) (driver)       M5    
Zgsgv Generalized eigenvalue problem of a complex sparse matrix (implicitly restarted Arnoldi method (IRAM)) (Arpack) (simple driver)       M5    
Znaupd Arnoldi factorization of a complex sparse matrix (Arpack)       M5    
Zneupd Approximate eigenvalues and eigenvectors of a general sparse matrix from Arnoldi factorization (complex matrix) (Arpack)       M5    
Singular value decomposition (sparse matrices) D6. Arpack (General matrices) Dgssvd Singular value decomposition (SVD) of a general sparse matrix (implicitly restarted Lanczos method (IRLM)) (Arpack) (driver)       M5    
Zgssvd Singular value decomposition (SVD) of a complex sparse matrix (implicitly restarted Arnoldi method (IRAM)) (Arpack) (driver)       M5    
Differential and integral equations I2. Partial differential equations Fem2p (*) Assemble finite element matrix of Poisson equation (2D) in CSR sparse matrix format       B    
Fem3p (*) Assemble finite element matrix of Poisson equation (3D) in CSR sparse matrix format       M5    
Mesh23 (*) Generates simple rectangular mesh for FEM (2D) (Triangular element)       B    
Mesh24 (*) Generates simple rectangular mesh for FEM (2D) (4-node quadrangle element)       M5    
Mesh29 (*) Generates simple rectangular mesh for FEM (2D) (9-node quadrangle element)       M5    
Mesh34 (*) Generates simple rectangular mesh for FEM (3D) (4-node tetrahedral element)       M5    
Mesh35 (*) Generates simple rectangular mesh for FEM (3D) (5-node pentahedral (pyramid) element)       M5    
Mesh36 (*) Generates simple rectangular mesh for FEM (3D) (6-node pentahedral (prism) element)       M5    
Mesh38 (*) Generates simple rectangular mesh for FEM (3D) (8-node tetrahedral element)       M5    
Data handling N1. Input, output of data MMRead Read a matrix from the Matrix Market file       B    
MMReadInfo Read matrix information from Matrix Market file       B    
MMWrite Write a matrix to the Matrix Market file       B    
HBRead Read a matrix from the Harwell-Boeing file       M5    
HBRead1 Read a matrix from the Harwell-Boeing file (simple driver)       M5    
HBReadInfo Read matrix information from Harwell-Boeing file       M5    
HBReadInfo1 Read matrix information from Harwell-Boeing file (simple driver)       M5    
HBWrite Write a matrix to the file in Harwell-Boeing format       M5    
HBWrite1 Write a matrix to the Harwell-Boeing file (Simple driver)       M5    
ZMMRead Read a matrix from the Matrix Market file (complex matrix)       M5    
ZMMWrite Write a matrix to the Matrix Market file (complex matrix)       M5    
ZHBRead Read a matrix from the Harwell-Boeing file (complex matrix)       M5    
ZHBRead1 Read a matrix from the Harwell-Boeing file (complex matrix) (simple driver)       M5    
ZHBWrite Write a matrix to the file in Harwell-Boeing format (complex matrix)       M5    
ZHBWrite1 Write a matrix to the Harwell-Boeing file (complex matrix) (simple driver)       M5    
ReadGmsh22 (*) Read FEM mesh information from Gmsh file (Version 2.2)       B    
ReadMsh2 (*) Read FEM mesh information from FreeFEM++ msh (2D) file       M5    
ReadMsh3 (*) Read FEM mesh information from FreeFEM++ msh (3D) file       M5    
ReadMesh (*) Read FEM mesh information from FreeFEM++ mesh file       M5    
WriteGmsh22 (*) Write FEM mesh information to Gmsh file (Version 2.2)       B    
WriteCsv2 (*) Write coordinates and their values to CSV file (2D)       M5    
WriteCsv3 (*) Write coordinates and their values to CSV file (3D)       M5    
WriteVtkug (*) VTK file output (Solution of PDE by FEM)       B    
Service routines R2. Error checking CsrCheck Check sparse matrix (CSR)       B    
CscCheck Check sparse matrix (CSC)       M5    
ZCsrCheck Check complex sparse matrix (CSR)       M5    
ZCscCheck Check complex sparse matrix (CSC)       M5