LINPACK_C
Linear Algebra Library
Single Precision Complex

LINPACK_C is a library of MATLAB routines, using single precision complex arithmetic, which can solve systems of linear equations for a variety of matrix types and storage modes.

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.

Related Data and Programs:

BLAS1 is a library of MATLAB routines for vector operations, used by LINPACK.

LAPACK is a FORTRAN90 library of routines which can solve linear systems of equations or problems involving eigenvalues.

LINPACK_C is also available in a C++ version and a FORTRAN77 version and a FORTRAN90 version.

LINPACK_D is a version of LINPACK for double precision real arithmetic.

LINPACK_S is a version of LINPACK for single precision real arithmetic.

LINPACK_Z is a version of LINPACK for double precision complex arithmetic.

LINPLUS is a library of MATLAB routines similar to LINPACK but covering some other storage formats.

TEMPLATES is a library of MATLAB routines that implement various iterative methods for solving a linear system.

TEST_MAT is a FORTRAN90 collection of test matrices with known inverses, determinants, eigenvalues and so on.

Reference:

Jack Dongarra, Jim Bunch, Cleve Moler, Pete Stewart,
LINPACK User's Guide,
SIAM, 1979,
ISBN13: 978-0-898711-72-1,
LC: QA214.L56.
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.

Tar File:

A GZIP'ed TAR file of the contents of this directory is available. This is only done as a convenience for users who want ALL the files, and don't want to download them individually. This is not a convenience for me, so don't be surprised if the tar file is somewhat out of date.

linpack_c_m_src.tar.csh the commands used to create the GZIP'ed TAR file.
linpack_c_m_src.tar.gz the GZIP'ed TAR file.

Source Code:

cabs1.m, returns the L1 norm of a complex number.
cchdc.m, computes the Cholesky decomposition of a positive definite matrix;
cchdd.m, downdates a Cholesky decomposition of a positive definite matrix;
cchex.m, updates a Cholesky decomposition of a positive definite matrix;
cchud.m, updates a Cholesky decomposition of a positive definite matrix;
cdotc.m, computes the conjugated dot product of two vectors.
cgbco.m, factors a band matrix and estimates its condition.
cgbdi.m, computes the determinant of a band matrix factored by CGBCO or CGBFA.
cgbfa.m, factors a complex band matrix by elimination.
cgbsl.m, solves a complex band system factored by CGBCO or CGBFA.
cgeco.m, factors a complex matrix and estimates its condition.
cgedi.m, computes the determinant and inverse of a matrix.
cgefa.m, factors a complex matrix by Gaussian elimination.
cgesl.m, solves a complex system factored by CGECO or CGEFA.
cgtsl.m, solves a complex tridiagonal system.
chico.m, factors a Hermitian matrix and estimates its condition.
chidi.m, computes inverse, determinant and interia of a Hermitian matrix.
chifa.m, factors a Hermitian matrix by Gaussian elimination.
chisl.m, solves a linear system factored by CHIFA or CHICO.
chpco.m, factors a Hermitian packed matrix and estimates its condition.
chpdi.m, computes inverse, determinant and interia of a Hermitian packed matrix.
chpfa.m, factors a Hermitian packed matrix by Gaussian elimination.
chpsl.m, solves a linear system factored by CHPFA or CHPCO.
cpbco.m, factors a Hermitian positive definite band matrix and estimates its condition.
cpbdi.m, computes inverse, determinant and interia of a Hermitian positive definite band matrix.
cpbfa.m, factors a Hermitian positive definite band matrix by Gaussian elimination.
cpbsl.m, solves a linear system factored by CPBFA or CPBCO.
cpoco.m, factors a Hermitian positive definite matrix and estimates its condition.
cpodi.m, computes inverse, determinant and interia of a Hermitian positive definite matrix.
cpofa.m, factors a Hermitian positive definite matrix by Gaussian elimination.
cposl.m, solves a linear system factored by CPOFA or CPOCO.
cppco.m, factors a Hermitian positive definite packed matrix and estimates its condition.
cppdi.m, computes inverse, determinant and interia of a Hermitian positive definite packed matrix.
cppfa.m, factors a Hermitian positive definite packed matrix by Gaussian elimination.
cppsl.m, solves a linear system factored by CPPFA or CPPCO.
cptsl.m, solves a Hermitian positive definite tridiagonal system.
cqrdc.m, computes the QR decomposition of a general rectangular complex matrix.
cqrsl.m, uses the QR decomposition to solve linear systems and other tasks.
crotg.m, computes a complex Givens rotation.
csico.m, factors a complex symmetric indefinite matrix and estimates its condition.
csidi.m, computes inverse and determinant of a complex symmetric indefinite matrix.
csign1.m, is a complex transfer-of-sign function.
csign2.m, is a complex transfer-of-sign function.
csrot.m, applies a complex Givens rotation.
csvdc.m, computes the singular value decomposition of a general complex matrix.
ctrco.m, returns the condition of a complex triangular matrix.
ctrdi.m, computes the inverse or determinant of a complex triangular matrix.
ctrsl.m, solves a complex triangular linear system.
i4_huge.m, returns a "huge" I4.
icamax.m, indexes the complex vector element of maximum absolute value.
r4_sign.m, returns the sign of an R4.
scasum.m, takes the sum of the absolute values of a complex vector.
scnrm2.m, computes the Euclidean norm of a complex vector.
srotg.m, computes a real Givens rotation.

Examples and Tests:

linpack_c_test.m, calls all the tests;
linpack_c_test.out, output from the tests.
linpack_c_test01.m, tests CCHDC;
linpack_c_test02.m, tests CCHEX;
linpack_c_test03.m, tests CCHUD;
linpack_c_test04.m, tests CGBCO;
linpack_c_test05.m, tests CGBFA and CGBSL;
linpack_c_test06.m, tests CGBFA and CGBDI;
linpack_c_test07.m, tests CGECO;
linpack_c_test08.m, tests CGEFA and CGESL;
linpack_c_test09.m, tests CGEFA and CGEDI;
linpack_c_test10.m, tests CGTSL;
linpack_c_test11.m, tests CHICO;
linpack_c_test12.m, tests CHIFA and CHISL;
linpack_c_test13.m, tests CHIFA and CHIDI;
linpack_c_test14.m, tests CHPCO;
linpack_c_test15.m, tests CHPFA and CHPSL;
linpack_c_test16.m, tests CHPFA and CHPDI;
linpack_c_test17.m, tests CPBCO;
linpack_c_test18.m, tests CPBFA and CPBSL;
linpack_c_test19.m, tests CPBFA and CPBDI;
linpack_c_test20.m, tests CPOCO;
linpack_c_test21.m, tests CPOFA and CPOSL;
linpack_c_test22.m, tests CPOFA and CPODI;
linpack_c_test23.m, tests CPPCO;
linpack_c_test24.m, tests CPPFA and CPPSL;
linpack_c_test25.m, tests CPPFA and CPPDI;
linpack_c_test26.m, tests CPTSL;
linpack_c_test27.m, tests CQRDC and CQRSL;
linpack_c_test28.m, tests CSICO;
linpack_c_test29.m, tests CSIFA and CSISL;
linpack_c_test30.m, tests CSIFA and CSIDI;
linpack_c_test31.m, tests CSPCO;
linpack_c_test32.m, tests CSPFA and CSPSL;
linpack_c_test33.m, tests CSPFA and CSPDI;
linpack_c_test34.m, tests CSVDC;
linpack_c_test35.m, tests CTRCO;
linpack_c_test36.m, tests CTRDI;
linpack_c_test37.m, tests CTRSL;
c4_uniform_01.m, returns a unit pseudorandom C4.
c4mat_uniform_01.m, returns a unit pseudorandom C4MAT.
c4vec_uniform_01.m, returns a unit pseudorandom C4VEC.
r4_uniform_01.m, returns a unit pseudorandom R4.
timestamp.m, prints the current YMDHMS date as a timestamp.

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

Last revised on 13 May 2007.