Note: the dimension of the array
AP
must be at least
max1,N×N+1/2.
On entry: the
n by
n triangular matrix
A, packed by columns.
More precisely,
- if UPLO='U', the upper triangle of A must be stored with element Aij in APi+jj-1/2 for i≤j;
- if UPLO='L', the lower triangle of A must be stored with element Aij in APi+2n-jj-1/2 for i≥j.
If DIAG='U', the diagonal elements of A are assumed to be 1, and are not referenced; the same storage scheme is used whether DIAG='N' or ‘U’.