hankel_cholesky

hankel_cholesky, a Fortran90 code which can compute the upper Cholesky factor of a symmetric positive definite (SPD) Hankel matrix H, that is, H = R' * R.

A Hankel matrix is a matrix which is constant along all antidiagonals. A schematic of a 5x5 symmetric Hankel matrix would be:

        a  b  c  d  e
        b  c  d  e  f
        c  d  e  f  g
        d  e  f  g  h
        e  f  g  h  i
      

Licensing:

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

Languages:

hankel_cholesky 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:

hankel_cholesky_test

asa006, a Fortran90 code which computes the Cholesky factorization of a symmetric positive definite (SPD) matrix, by Michael Healy. This is a version of Applied Statistics Algorithm 6;

hankel_spd, a Fortran90 code which can compute a lower triangular matrix L which is the Cholesky factor of a symmetric positive definite (SPD) Hankel matrix H, that is, H = L * L'.

toeplitz_cholesky, a Fortran90 code which computes the Cholesky factorization of a symmetric positive definite (SPD) Toeplitz matrix.

Reference:

1. James Phillips,
   The triangular decomposition of Hankel matrices,
   Mathematics of Computation,
   Volume 25, Number 115, July 1971, pages 599-602.

Source Code:

• hankel_cholesky.f90, the source file.
• hankel_cholesky.sh, compiles the source file.