sdrvrfp.f3

Usage

Functions/Subroutines

subroutine sdrvrfp (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)
SDRVRFP

Function/Subroutine Documentation

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)

SDRVRFP

Purpose:

!>
!> 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.
!>

Parameters

NOUT

!>          NOUT is INTEGER
!>                The unit number for output.
!>

!>          NN is INTEGER
!>                The number of values of N contained in the vector NVAL.
!>

NVAL

!>          NVAL is INTEGER array, dimension (NN)
!>                The values of the matrix dimension N.
!>

NNS

!>          NNS is INTEGER
!>                The number of values of NRHS contained in the vector NSVAL.
!>

NSVAL

!>          NSVAL is INTEGER array, dimension (NNS)
!>                The values of the number of right-hand sides NRHS.
!>

NNT

!>          NNT is INTEGER
!>                The number of values of MATRIX TYPE contained in the vector NTVAL.
!>

NTVAL

!>          NTVAL is INTEGER array, dimension (NNT)
!>                The values of matrix type (between 0 and 9 for PO/PP/PF matrices).
!>

THRESH

!>          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.
!>

!>          A is REAL array, dimension (NMAX*NMAX)
!>

ASAV

!>          ASAV is REAL array, dimension (NMAX*NMAX)
!>

AFAC

!>          AFAC is REAL array, dimension (NMAX*NMAX)
!>

AINV

!>          AINV is REAL array, dimension (NMAX*NMAX)
!>

!>          B is REAL array, dimension (NMAX*MAXRHS)
!>

BSAV

!>          BSAV is REAL array, dimension (NMAX*MAXRHS)
!>

XACT

!>          XACT is REAL array, dimension (NMAX*MAXRHS)
!>

!>          X is REAL array, dimension (NMAX*MAXRHS)
!>

ARF

!>          ARF is REAL array, dimension ((NMAX*(NMAX+1))/2)
!>

ARFINV

!>          ARFINV is REAL array, dimension ((NMAX*(NMAX+1))/2)
!>

S_WORK_SLATMS

!>          S_WORK_SLATMS is REAL array, dimension ( 3*NMAX )
!>

S_WORK_SPOT01

!>          S_WORK_SPOT01 is REAL array, dimension ( NMAX )
!>

S_TEMP_SPOT02

!>          S_TEMP_SPOT02 is REAL array, dimension ( NMAX*MAXRHS )
!>

S_TEMP_SPOT03

!>          S_TEMP_SPOT03 is REAL array, dimension ( NMAX*NMAX )
!>

S_WORK_SLANSY

!>          S_WORK_SLANSY is REAL array, dimension ( NMAX )
!>

S_WORK_SPOT02

!>          S_WORK_SPOT02 is REAL array, dimension ( NMAX )
!>

S_WORK_SPOT03

!>          S_WORK_SPOT03 is REAL array, dimension ( NMAX )
!>

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.

Definition at line 232 of file sdrvrfp.f.