asa007

asa007, a Fortran77 code which computes the inverse of a symmetric positive definite (SPD) matrix, by Michael Healy.

This is a version of Applied Statistics Algorithm 7.

The algorithm implemented here uses a compressed storage for both the matrix A and the factor U. This saves some storage, but can make computations a little awkward.

Licensing:

The information on this web page is distributed under the MIT license.

Languages:

asa007 is available in a C version and a C++ version and a Fortran90 version and a MATLAB version and an Octave version and a Python version.

Related Data and Programs:

asa007_test

asa006, a Fortran77 library which computes the Cholesky factorization of a symmetric positive definite matrix. A version of this routine is used by ASA007.

ASA047, a Fortran77 library which implements the Nelder-Mead minimization algorithm, and uses a version of ASA007 for the computation.

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

LINPACK, a Fortran77 library which includes routines to compute the inverse of symmetric positive matrix.

Author:

Michael Healy

Reference:

1. Michael Healy,
   Algorithm AS 7: Inversion of a Positive Semi-Definite Symmetric Matrix,
   Applied Statistics,
   Volume 17, Number 2, 1968, pages 198-199.

Source Code:

• asa007.f, the source code.
• asa007.sh, commands to compile the source code.

Examples and Tests:

• asa007_prb.f, a sample calling program.
• asa007_prb.sh, commands to compile and run the sample program.
• asa007_prb_output.txt, the output file.

List of Routines:

• CHOL computes the Cholesky factorization of a PDS matrix.
• SYMINV computes the inverse of a symmetric matrix.
• TIMESTAMP prints out the current YMDHMS date as a timestamp.

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