LINPACK_Z
Linear Algebra Library
Double Precision Complex

LINPACK_Z is a MATLAB library which solves systems of linear equations for a variety of matrix types and storage modes, by Jack Dongarra, Jim Bunch, Cleve Moler, Pete Stewart.

MATLAB already provides a wide set of linear equation solvers. This (partial) set of LINPACK routines is provided just for testing and comparison.

LINPACK has officially been superseded by the LAPACK library. The LAPACK library uses more modern algorithms and code structure. However, the LAPACK library can be extraordinarily complex; what is done in a single LINPACK routine may correspond to 10 or 20 utility routines in LAPACK. This is fine if you treat LAPACK as a black box. But if you wish to learn how the algorithm works, or to adapt it, or to convert the code to another language, this is a real drawback. This is one reason I still keep a copy of LINPACK around.

Versions of LINPACK in various arithmetic precisions are available through the NETLIB web site.

Licensing:

The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.

Languages:

LINPACK_Z is available in a C++ version and a FORTRAN77 version and a FORTRAN90 version and a MATLAB version.

Related Data and Programs:

BLAS1_Z, a MATLAB library which contains basic linear algebra routines for vector-vector operations, using double precision complex arithmetic.

LAPACK_EXAMPLES, a FORTRAN90 program which demonstrates the use of the LAPACK linear algebra library.

LINPACK_C, a MATLAB library which solves linear systems using sihgle precision complex arithmetic;

LINPACK_D, a MATLAB library which solves linear systems using double precision real arithmetic;

LINPACK_S, a MATLAB library which solves linear systems using single precision real arithmetic;

LINPLUS, a MATLAB library which carries out simple manipulations of matrices in a variety of formats.

TEST_MAT, a MATLAB library which defines test matrices.

TEST_MATRIX, a MATLAB library which contains a collection of test matrices by Nick Higham.

Author:

Original FORTRAN77 version by Jack Dongarra, Jim Bunch, Cleve Moler, Pete Stewart. MATLAB version by John Burkardt.

Reference:

1. Jack Dongarra, Jim Bunch, Cleve Moler, Pete Stewart,
   LINPACK User's Guide,
   SIAM, 1979,
   ISBN13: 978-0-898711-72-1,
   LC: QA214.L56.
2. Charles Lawson, Richard Hanson, David Kincaid, Fred Krogh,
   Algorithm 539, Basic Linear Algebra Subprograms for Fortran Usage,
   ACM Transactions on Mathematical Software,
   Volume 5, Number 3, September 1979, pages 308-323.

Source Code:

• drotg.m, computes a real Givens rotation.
• dzasum.m, takes the sum of the absolute values of a complex vector.
• dznrm2.m, computes the Euclidean norm of a complex vector.
\
• i4_huge.m, returns a "huge" I4.
• izamax.m, indexes the complex vector element of maximum absolute value.
• r8_sign.m, returns the sign of an R8.
• zabs1.m, returns the L1 norm of a complex number.
• zchdc.m, computes the Cholesky decomposition of a positive definite matrix;
• zchdd.m, downdates a Cholesky decomposition of a positive definite matrix;
• zchex.m, updates a Cholesky decomposition of a positive definite matrix;
• zchud.m, updates a Cholesky decomposition of a positive definite matrix;
• zdotc.m, computes the conjugated dot product of two vectors.
• zdrot.m, applies a complex Givens rotation.
• zgbco.m, factors a band matrix and estimates its condition.
• zgbdi.m, computes the determinant of a band matrix factored by ZGBCO or ZGBFA.
• zgbfa.m, factors a complex band matrix by elimination.
• zgbsl.m, solves a complex band system factored by ZGBCO or ZGBFA.
• zgeco.m, factors a complex matrix and estimates its condition.
• zgedi.m, computes the determinant and inverse of a matrix.
• zgefa.m, factors a complex matrix by Gaussian elimination.
• zgesl.m, solves a complex system factored by ZGECO or ZGEFA.
• zgtsl.m, solves a complex tridiagonal system.
• zhico.m, factors a Hermitian matrix and estimates its condition.
• zhidi.m, computes inverse, determinant and interia of a Hermitian matrix.
• zhifa.m, factors a Hermitian matrix by Gaussian elimination.
• zhisl.m, solves a linear system factored by ZHIFA or ZHICO.
• zhpco.m, factors a Hermitian packed matrix and estimates its condition.
• zhpdi.m, computes inverse, determinant and interia of a Hermitian packed matrix.
• zhpfa.m, factors a Hermitian packed matrix by Gaussian elimination.
• zhpsl.m, solves a linear system factored by ZHPFA or ZHPCO.
• zpbco.m, factors a Hermitian positive definite band matrix and estimates its condition.
• zpbdi.m, computes inverse, determinant and interia of a Hermitian positive definite band matrix.
• zpbfa.m, factors a Hermitian positive definite band matrix by Gaussian elimination.
• zpbsl.m, solves a linear system factored by ZPBFA or ZPBCO.
• zpoco.m, factors a Hermitian positive definite matrix and estimates its condition.
• zpodi.m, computes inverse, determinant and interia of a Hermitian positive definite matrix.
• zpofa.m, factors a Hermitian positive definite matrix by Gaussian elimination.
• zposl.m, solves a linear system factored by ZPOFA or ZPOCO.
• zppco.m, factors a Hermitian positive definite packed matrix and estimates its condition.
• zppdi.m, computes inverse, determinant and interia of a Hermitian positive definite packed matrix.
• zppfa.m, factors a Hermitian positive definite packed matrix by Gaussian elimination.
• zppsl.m, solves a linear system factored by ZPPFA or ZPPCO.
• zptsl.m, solves a Hermitian positive definite tridiagonal system.
• zqrdc.m, computes the QR decomposition of a general rectangular complex matrix.
• zqrsl.m, uses the QR decomposition to solve linear systems and other tasks.
• zrotg.m, computes a complex Givens rotation.
• zsico.m, factors a complex symmetric indefinite matrix and estimates its condition.
• zsidi.m, computes inverse and determinant of a complex symmetric indefinite matrix.
• zsign1.m, is a complex transfer-of-sign function.
• zsign2.m, is a complex transfer-of-sign function.
• zsvdc.m, computes the singular value decomposition of a general complex matrix.
• ztrco.m, returns the condition of a complex triangular matrix.
• ztrdi.m, computes the inverse or determinant of a complex triangular matrix.
• ztrsl.m, solves a complex triangular linear system.

Examples and Tests:

• linpack_z_test.m, calls all the tests;
• linpack_z_test_output.txt, output from the tests.
• linpack_z_test01.m, tests ZCHDC;
• linpack_z_test02.m, tests ZCHEX;
• linpack_z_test03.m, tests ZCHUD;
• linpack_z_test04.m, tests ZGBCO;
• linpack_z_test05.m, tests ZGBFA and ZGBSL;
• linpack_z_test06.m, tests ZGBFA and ZGBDI;
• linpack_z_test07.m, tests ZGECO;
• linpack_z_test08.m, tests ZGEFA and ZGESL;
• linpack_z_test09.m, tests ZGEFA and ZGEDI;
• linpack_z_test10.m, tests ZGTSL;
• linpack_z_test11.m, tests ZHICO;
• linpack_z_test12.m, tests ZHIFA and XZISL;
• linpack_z_test13.m, tests ZHIFA and ZHIDI;
• linpack_z_test14.m, tests ZHPCO;
• linpack_z_test15.m, tests ZHPFA and ZHPSL;
• linpack_z_test16.m, tests ZHPFA and ZHPDI;
• linpack_z_test17.m, tests ZPBCO;
• linpack_z_test18.m, tests ZPBFA and ZPBSL;
• linpack_z_test19.m, tests ZPBFA and ZPBDI;
• linpack_z_test20.m, tests ZPOCO;
• linpack_z_test21.m, tests ZPOFA and ZPOSL;
• linpack_z_test22.m, tests ZPOFA and ZPODI;
• linpack_z_test23.m, tests ZPPCO;
• linpack_z_test24.m, tests ZPPFA and ZPPSL;
• linpack_z_test25.m, tests ZPPFA and ZPPDI;
• linpack_z_test26.m, tests ZPTSL;
• linpack_z_test27.m, tests ZQRDC and ZQRSL;
• linpack_z_test28.m, tests ZSICO;
• linpack_z_test29.m, tests ZSIFA and ZSISL;
• linpack_z_test30.m, tests ZSIFA and ZSIDI;
• linpack_z_test31.m, tests ZSPCO;
• linpack_z_test32.m, tests ZSPFA and ZSPSL;
• linpack_z_test33.m, tests ZSPFA and ZSPDI;
• linpack_z_test34.m, tests ZSVDC;
• linpack_z_test345.m, tests ZSVDC;
• linpack_z_test35.m, tests ZTRCO;
• linpack_z_test36.m, tests ZTRDI;
• linpack_z_test37.m, tests ZTRSL;
• c8_uniform_01.m, returns a unit pseudorandom C8.
• c4mat_uniform_01.m, returns a unit pseudorandom C8MAT.
• c4vec_uniform_01.m, returns a unit pseudorandom C8VEC.
• r8_uniform_01.m, returns a unit pseudorandom R8.
• timestamp.m, prints the current YMDHMS date as a timestamp.

You can go up one level to the MATLAB source codes.