.TH "TESTING/LIN/sdrvrfp.f" 3 "Version 3.12.0" "LAPACK" \" -*- nroff -*-
.ad l
.nh
.SH NAME
TESTING/LIN/sdrvrfp.f
.SH SYNOPSIS
.br
.PP
.SS "Functions/Subroutines"

.in +1c
.ti -1c
.RI "subroutine \fBsdrvrfp\fP (nout, nn, nval, nns, nsval, nnt, ntval, thresh, a, asav, afac, ainv, b, bsav, xact, x, arf, arfinv, s_work_slatms, s_work_spot01, s_temp_spot02, s_temp_spot03, s_work_slansy, s_work_spot02, s_work_spot03)"
.br
.RI "\fBSDRVRFP\fP "
.in -1c
.SH "Function/Subroutine Documentation"
.PP 
.SS "subroutine sdrvrfp (integer nout, integer nn, integer, dimension( nn ) nval, integer nns, integer, dimension( nns ) nsval, integer nnt, integer, dimension( nnt ) ntval, real thresh, real, dimension( * ) a, real, dimension( * ) asav, real, dimension( * ) afac, real, dimension( * ) ainv, real, dimension( * ) b, real, dimension( * ) bsav, real, dimension( * ) xact, real, dimension( * ) x, real, dimension( * ) arf, real, dimension( * ) arfinv, real, dimension( * ) s_work_slatms, real, dimension( * ) s_work_spot01, real, dimension( * ) s_temp_spot02, real, dimension( * ) s_temp_spot03, real, dimension( * ) s_work_slansy, real, dimension( * ) s_work_spot02, real, dimension( * ) s_work_spot03)"

.PP
\fBSDRVRFP\fP 
.PP
\fBPurpose:\fP
.RS 4

.PP
.nf
 SDRVRFP tests the LAPACK RFP routines:
     SPFTRF, SPFTRS, and SPFTRI\&.

 This testing routine follow the same tests as DDRVPO (test for the full
 format Symmetric Positive Definite solver)\&.

 The tests are performed in Full Format, conversion back and forth from
 full format to RFP format are performed using the routines STRTTF and
 STFTTR\&.

 First, a specific matrix A of size N is created\&. There is nine types of
 different matrixes possible\&.
  1\&. Diagonal                        6\&. Random, CNDNUM = sqrt(0\&.1/EPS)
  2\&. Random, CNDNUM = 2              7\&. Random, CNDNUM = 0\&.1/EPS
 *3\&. First row and column zero       8\&. Scaled near underflow
 *4\&. Last row and column zero        9\&. Scaled near overflow
 *5\&. Middle row and column zero
 (* - tests error exits from SPFTRF, no test ratios are computed)
 A solution XACT of size N-by-NRHS is created and the associated right
 hand side B as well\&. Then SPFTRF is called to compute L (or U), the
 Cholesky factor of A\&. Then L (or U) is used to solve the linear system
 of equations AX = B\&. This gives X\&. Then L (or U) is used to compute the
 inverse of A, AINV\&. The following four tests are then performed:
 (1) norm( L*L' - A ) / ( N * norm(A) * EPS ) or
     norm( U'*U - A ) / ( N * norm(A) * EPS ),
 (2) norm(B - A*X) / ( norm(A) * norm(X) * EPS ),
 (3) norm( I - A*AINV ) / ( N * norm(A) * norm(AINV) * EPS ),
 (4) ( norm(X-XACT) * RCOND ) / ( norm(XACT) * EPS ),
 where EPS is the machine precision, RCOND the condition number of A, and
 norm( \&. ) the 1-norm for (1,2,3) and the inf-norm for (4)\&.
 Errors occur when INFO parameter is not as expected\&. Failures occur when
 a test ratios is greater than THRES\&.
.fi
.PP
 
.RE
.PP
\fBParameters\fP
.RS 4
\fINOUT\fP 
.PP
.nf
          NOUT is INTEGER
                The unit number for output\&.
.fi
.PP
.br
\fINN\fP 
.PP
.nf
          NN is INTEGER
                The number of values of N contained in the vector NVAL\&.
.fi
.PP
.br
\fINVAL\fP 
.PP
.nf
          NVAL is INTEGER array, dimension (NN)
                The values of the matrix dimension N\&.
.fi
.PP
.br
\fINNS\fP 
.PP
.nf
          NNS is INTEGER
                The number of values of NRHS contained in the vector NSVAL\&.
.fi
.PP
.br
\fINSVAL\fP 
.PP
.nf
          NSVAL is INTEGER array, dimension (NNS)
                The values of the number of right-hand sides NRHS\&.
.fi
.PP
.br
\fINNT\fP 
.PP
.nf
          NNT is INTEGER
                The number of values of MATRIX TYPE contained in the vector NTVAL\&.
.fi
.PP
.br
\fINTVAL\fP 
.PP
.nf
          NTVAL is INTEGER array, dimension (NNT)
                The values of matrix type (between 0 and 9 for PO/PP/PF matrices)\&.
.fi
.PP
.br
\fITHRESH\fP 
.PP
.nf
          THRESH is REAL
                The threshold value for the test ratios\&.  A result is
                included in the output file if RESULT >= THRESH\&.  To have
                every test ratio printed, use THRESH = 0\&.
.fi
.PP
.br
\fIA\fP 
.PP
.nf
          A is REAL array, dimension (NMAX*NMAX)
.fi
.PP
.br
\fIASAV\fP 
.PP
.nf
          ASAV is REAL array, dimension (NMAX*NMAX)
.fi
.PP
.br
\fIAFAC\fP 
.PP
.nf
          AFAC is REAL array, dimension (NMAX*NMAX)
.fi
.PP
.br
\fIAINV\fP 
.PP
.nf
          AINV is REAL array, dimension (NMAX*NMAX)
.fi
.PP
.br
\fIB\fP 
.PP
.nf
          B is REAL array, dimension (NMAX*MAXRHS)
.fi
.PP
.br
\fIBSAV\fP 
.PP
.nf
          BSAV is REAL array, dimension (NMAX*MAXRHS)
.fi
.PP
.br
\fIXACT\fP 
.PP
.nf
          XACT is REAL array, dimension (NMAX*MAXRHS)
.fi
.PP
.br
\fIX\fP 
.PP
.nf
          X is REAL array, dimension (NMAX*MAXRHS)
.fi
.PP
.br
\fIARF\fP 
.PP
.nf
          ARF is REAL array, dimension ((NMAX*(NMAX+1))/2)
.fi
.PP
.br
\fIARFINV\fP 
.PP
.nf
          ARFINV is REAL array, dimension ((NMAX*(NMAX+1))/2)
.fi
.PP
.br
\fIS_WORK_SLATMS\fP 
.PP
.nf
          S_WORK_SLATMS is REAL array, dimension ( 3*NMAX )
.fi
.PP
.br
\fIS_WORK_SPOT01\fP 
.PP
.nf
          S_WORK_SPOT01 is REAL array, dimension ( NMAX )
.fi
.PP
.br
\fIS_TEMP_SPOT02\fP 
.PP
.nf
          S_TEMP_SPOT02 is REAL array, dimension ( NMAX*MAXRHS )
.fi
.PP
.br
\fIS_TEMP_SPOT03\fP 
.PP
.nf
          S_TEMP_SPOT03 is REAL array, dimension ( NMAX*NMAX )
.fi
.PP
.br
\fIS_WORK_SLANSY\fP 
.PP
.nf
          S_WORK_SLANSY is REAL array, dimension ( NMAX )
.fi
.PP
.br
\fIS_WORK_SPOT02\fP 
.PP
.nf
          S_WORK_SPOT02 is REAL array, dimension ( NMAX )
.fi
.PP
.br
\fIS_WORK_SPOT03\fP 
.PP
.nf
          S_WORK_SPOT03 is REAL array, dimension ( NMAX )
.fi
.PP
 
.RE
.PP
\fBAuthor\fP
.RS 4
Univ\&. of Tennessee 
.PP
Univ\&. of California Berkeley 
.PP
Univ\&. of Colorado Denver 
.PP
NAG Ltd\&. 
.RE
.PP

.PP
Definition at line \fB232\fP of file \fBsdrvrfp\&.f\fP\&.
.SH "Author"
.PP 
Generated automatically by Doxygen for LAPACK from the source code\&.