next up previous contents index
Next: Examples Up: Symmetric Indefinite Linear Systems Previous: Purpose   Contents   Index

Arguments

AP
(input/output) REAL or COMPLEX array, shape $(:)$ with $size(AP) = n(n+1)/2$, where $n$ is the order of $A$.
On entry, the upper or lower triangle of matrix $A$ in packed storage. The elements are stored columnwise as follows:

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

On exit, the block diagonal matrix $D$ and the multipliers used to obtain $U$ or $L$ from the factorization of $A$, stored as a packed triangular matrix 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'.

IPIV
Optional (output) INTEGER array, shape $(:)$ with $size({\bf IPIV})=n$.
Details of the row and column interchanges and the block structure of $D$.
If ${\bf IPIV}_k > 0$, then rows and columns $k$ and ${\bf IPIV}_k$ were interchanged, and $D_{k,k}$ is a $1 \times 1$ diagonal block.
If ${\bf IPIV}_k < 0$, then there are two cases:

\begin{numbersec}
\item If ${\bf UPLO} =$\ 'U' and ${\bf IPIV}_k = {\bf IPIV}_{...
...nd $D_{k:k+1,k:k+1}$\ is a
\hbox{$2 \times 2$} diagonal block.
\end{numbersec}

INFO
Optional (output) INTEGER.

\begin{infoarg}
\item[{$=$\ 0:}] successful exit
\item[{$<$\ 0:}] if ${\bf INF...
...l matrix $D$\ is
singular, so 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,21].
next up previous contents index
Next: Examples Up: Symmetric Indefinite Linear Systems Previous: Purpose   Contents   Index
Susan Blackford 2001-08-19