next up previous contents index
Next: Example 2 (from Program Up: Examples Previous: Examples   Contents   Index

Example 1

The packed storage scheme is illustrated below for the matrix


\begin{displaymath}
A = \left( \begin{array}{rrrr}
a_{11} & a_{12} & a_{13} & a...
...{34} \\
a_{41} & a_{42} & a_{43} & a_{44} \end{array} \right)
\end{displaymath}

where $a_{ij} = \mbox{conjg}(a_{ji}) \;\; \mbox{for} \;\; i,j = 1,
\cdots, 4.$

If ${\bf UPLO} =$ 'U' then:


\begin{displaymath}
\begin{array}{c} {\bf AP} \;\; \mbox{on} \;\; \mbox{entry} ...
...} &
a_{24} & a_{34} & a_{44} \\ \hline \end{array} \end{array}\end{displaymath}

If ${\bf UPLO} =$ 'L' then:


\begin{displaymath}
\begin{array}{c} {\bf AP} \;\; \mbox{on} \;\; \mbox{entry} ...
...} &
a_{33} & a_{43} & a_{44} \\ \hline \end{array} \end{array}\end{displaymath}



Susan Blackford 2001-08-19