Next: Band Storage Up: Matrix Storage Schemes Previous: Conventional Storage Contents Index

##
Packed Storage

Symmetric, Hermitian or triangular matrices may be stored more compactly, if the relevant triangle (again as specified by UPLO) is packed by columns in a one-dimensional array. In LAPACK, arrays that hold matrices in packed storage, have names ending in `P'. So:

- if UPLO = `U', aij is stored in AP(i+j(j-1)/2) for ;

- if UPLO = `L', aij is stored in AP( i+(2n-j)(j-1)/2) for .

For example:

UPLO | Triangular matrix A | Packed storage in array AP
`U' |  |
`L' |  |

Note that for real or complex symmetric matrices, packing the upper triangle by columns is equivalent to packing the lower triangle by rows; packing the lower triangle by columns is equivalent to packing the upper triangle by rows. For complex Hermitian matrices, packing the upper triangle by columns is equivalent to packing the conjugate of the lower triangle by rows; packing the lower triangle by columns is equivalent to packing the conjugate of the upper triangle by rows.

Next: Band Storage Up: Matrix Storage Schemes Previous: Conventional Storage Contents Index