Example (from Program LA_PBSVX_EXAMPLE)

The results below are computed with $\epsilon = 1.19209 \times 10^{-7}$.
$A$ and $B$ are the matrices from Example 2 for LA_PBSV.

The call:

  CALL LA_PBSVX( A, B, X, FACT='E', EQUED=EQUED, S=S, & 


FERR=FERR, BERR=BERR, RCOND=RCOND ) 

EQUED, S, FERR, BERR and RCOND on exit:

$${\bf EQUED} = \mbox{'N'} \ \ \ \ \begin{array}{c} {\bf S} \\... ... \end{array} \ \ \ \ \ \ {\bf RCOND} = 3.15913 \times 10^{-1}$$

$$\begin{array}{c} {\bf FERR} \\ \begin{array}{\vert rrr\vert}... ... 4.61455 \times 10^{-8} & 0 \\ \hline \end{array} \end{array}$$

The forward and backward errors of the three solution vectors are:

$$\left( \begin{array}{ccc} 2.05706 \times 10^{-6} & 2.05706 ... ...imes 10^{-8} & 4.61455 \times 10^{-8} & 0 \end{array} \right).$$

The solution of the system $A\,X = B$ is:

$$X = \left( \begin{array}{rrr} 1.00000 & 2.00000 & 3.00000 ... ...3.00000 \\ 1.00000 & 2.00000 & 3.00000 \end{array} \right).$$