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

Example 3 (from Program LA_HESV_EXAMPLE)


\begin{displaymath}
\hspace{-0.50 cm}
A = \left( \begin{array}{ccccc}
1 & 4 +...
... \\ 17 + 43i \\ 47 + 23i \\ 39 + 31i \\ 56 \end{array} \right)
\end{displaymath}

Arrays ${\bf A}$ and ${\bf B}$ on entry:

\begin{displaymath}
\begin{array}{c} {\bf A} \\
\begin{array}{\vert ccccc\ver...
...(47,23) \\ (39,31) \\ (56,0)
\\ \hline \end{array} \end{array}\end{displaymath}

Elements marked $\ast$ are not used by the routine.

The call:
CALL LA_HESV( A, B, IPIV
=IPIV )

${\bf A}$, ${\bf B}$ and ${\bf IPIV}$ on exit:

\begin{displaymath}
\begin{array}{c} {\bf A} \\
\begin{array}{\vert ccc} \hli...
...ast \\
\ast & \ast & \ast \\
\hline \end{array} \end{array}\end{displaymath}


\begin{displaymath}
\begin{array}{c} {\bf A (continued)} \\
\begin{array}{cc\...
... \\
\ast & ( 3.00000, 0 ) \\
\hline \end{array} \end{array}\end{displaymath}


\begin{displaymath}
\begin{array}{cc} {\bf B} \\
\begin{array}{\vert r\vert} ...
...ne 1 \\ -2 \\ -2 \\ -2 \\ -2 \\ \hline
\end{array} \end{array}\end{displaymath}

The solution of the system $A\,x = b$ is:

\begin{displaymath}x = \left( \begin{array}{rrr}
1.00000 + 1.00000i \\
1.0000...
...00000 + 1.00000i \\
1.00000 + 1.00000i
\end{array} \right). \end{displaymath}



Susan Blackford 2001-08-19