1SPPRFS(1)                LAPACK routine (version 3.2)                SPPRFS(1)
2
3
4

NAME

6       SPPRFS - improves the computed solution to a system of linear equations
7       when the coefficient matrix is symmetric positive definite and  packed,
8       and provides error bounds and backward error estimates for the solution
9

SYNOPSIS

11       SUBROUTINE SPPRFS( UPLO,  N, NRHS, AP, AFP, B, LDB, X, LDX, FERR, BERR,
12                          WORK, IWORK, INFO )
13
14           CHARACTER      UPLO
15
16           INTEGER        INFO, LDB, LDX, N, NRHS
17
18           INTEGER        IWORK( * )
19
20           REAL           AFP( * ), AP( * ), B( LDB, * ), BERR( * ),  FERR(  *
21                          ), WORK( * ), X( LDX, * )
22

PURPOSE

24       SPPRFS  improves  the computed solution to a system of linear equations
25       when the coefficient matrix is symmetric positive definite and  packed,
26       and  provides  error  bounds and backward error estimates for the solu‐
27       tion.
28

ARGUMENTS

30       UPLO    (input) CHARACTER*1
31               = 'U':  Upper triangle of A is stored;
32               = 'L':  Lower triangle of A is stored.
33
34       N       (input) INTEGER
35               The order of the matrix A.  N >= 0.
36
37       NRHS    (input) INTEGER
38               The number of right hand sides, i.e., the number of columns  of
39               the matrices B and X.  NRHS >= 0.
40
41       AP      (input) REAL array, dimension (N*(N+1)/2)
42               The  upper  or lower triangle of the symmetric matrix A, packed
43               columnwise in a linear array.  The j-th column of A  is  stored
44               in  the array AP as follows: if UPLO = 'U', AP(i + (j-1)*j/2) =
45               A(i,j) for 1<=i<=j; if UPLO = 'L',  AP(i  +  (j-1)*(2n-j)/2)  =
46               A(i,j) for j<=i<=n.
47
48       AFP     (input) REAL array, dimension (N*(N+1)/2)
49               The  triangular factor U or L from the Cholesky factorization A
50               = U**T*U or A = L*L**T, as computed  by  SPPTRF/CPPTRF,  packed
51               columnwise in a linear array in the same format as A (see AP).
52
53       B       (input) REAL array, dimension (LDB,NRHS)
54               The right hand side matrix B.
55
56       LDB     (input) INTEGER
57               The leading dimension of the array B.  LDB >= max(1,N).
58
59       X       (input/output) REAL array, dimension (LDX,NRHS)
60               On  entry,  the  solution  matrix X, as computed by SPPTRS.  On
61               exit, the improved solution matrix X.
62
63       LDX     (input) INTEGER
64               The leading dimension of the array X.  LDX >= max(1,N).
65
66       FERR    (output) REAL array, dimension (NRHS)
67               The estimated forward error bound for each solution vector X(j)
68               (the  j-th  column  of the solution matrix X).  If XTRUE is the
69               true solution corresponding to X(j), FERR(j)  is  an  estimated
70               upper bound for the magnitude of the largest element in (X(j) -
71               XTRUE) divided by the magnitude of the largest element in X(j).
72               The  estimate  is as reliable as the estimate for RCOND, and is
73               almost always a slight overestimate of the true error.
74
75       BERR    (output) REAL array, dimension (NRHS)
76               The componentwise relative backward error of each solution vec‐
77               tor  X(j) (i.e., the smallest relative change in any element of
78               A or B that makes X(j) an exact solution).
79
80       WORK    (workspace) REAL array, dimension (3*N)
81
82       IWORK   (workspace) INTEGER array, dimension (N)
83
84       INFO    (output) INTEGER
85               = 0:  successful exit
86               < 0:  if INFO = -i, the i-th argument had an illegal value
87

PARAMETERS

89       ITMAX is the maximum number of steps of iterative refinement.
90
91
92
93 LAPACK routine (version 3.2)    November 2008                       SPPRFS(1)
Impressum