.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\&.