LINPACK_C
Linear Algebra Library
Single Precision Complex

LINPACK_C 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_C is available in a C++ version and a FORTRAN77 version and a FORTRAN90 version and a MATLAB version.

Related Data and Programs:

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

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

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;

LINPACK_Z, a MATLAB library which solves linear systems using double precision complex arithmetic;

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

TEMPLATES, a MATLAB library which carries out simple versions of various iterative solvers.

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:

• 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_output.txt, 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_test345.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.