next up previous contents index
Next: Examples Up: Symmetric/Hermitian Positive Definite Linear Previous: Purpose   Contents   Index

Arguments

AB
(input/output) REAL or COMPLEX array, shape $(:,:)$ with $size({\bf AB},1) = kd+1$ and $size({\bf AB},2) = n $, where $kd$ is the number of superdiagonals or subdiagonals in the band and $n$ is the order of $A$.
On entry, the upper (if UPLO = 'U') or lower (if UPLO = 'L') triangle of matrix $A$ in band storage. The $(kd+1)$ diagonals of $A$ are stored in the rows of AB so that the $j^{th}$ column of $A$ is stored in the $j^{th}$ column of ${\bf AB}$ as follows:

\begin{displaymath}
\begin{array}{c\vert c\vert c}
A_{i,j} & i,j & {\bf UPLO} ...
...1 \leq j \leq n \end{array} & \mbox{'L'} \\ \hline
\end{array}\end{displaymath}

On exit, the factor $U$ or $L$ from the Cholesky factorization $A = U^HU = LL^H$ in the same storage format as $A$.

B
(input/output) REAL or COMPLEX array, shape $(:,:)$ with $size({\bf B},1) = n$ or shape $(:)$ with $size({\bf B}) = n$.
On entry, the matrix $B$.
On exit, the solution matrix $X$.

UPLO
Optional (input) CHARACTER(LEN=1)

\begin{optionarg}
\item[{= 'U':}] Upper triangle of A is stored;
\item[{= 'L':}] Lower triangle of A is stored.
\end{optionarg}
Default value: 'U'.

INFO
Optional (output) INTEGER.

\begin{infoarg}
\item[{$=$\ 0:}] successful exit.
\item[{$<$\ 0:}] if ${\bf IN...
...n could not be
completed and the solution could not be computed.
\end{infoarg}
If INFO is not present and an error occurs, then the program is terminated with an error message.
References: [1] and [17,9,20].
next up previous contents index
Next: Examples Up: Symmetric/Hermitian Positive Definite Linear Previous: Purpose   Contents   Index
Susan Blackford 2001-08-19