ddrvrfp.f3

Usage

Functions/Subroutines

subroutine ddrvrfp (nout, nn, nval, nns, nsval, nnt, ntval, thresh, a, asav, afac, ainv, b, bsav, xact, x, arf, arfinv, d_work_dlatms, d_work_dpot01, d_temp_dpot02, d_temp_dpot03, d_work_dlansy, d_work_dpot02, d_work_dpot03)
DDRVRFP

Function/Subroutine Documentation

subroutine ddrvrfp (integer nout, integer nn, integer, dimension( nn ) nval, integer nns, integer, dimension( nns ) nsval, integer nnt, integer, dimension( nnt ) ntval, double precision thresh, double precision, dimension( * ) a, double precision, dimension( * ) asav, double precision, dimension( * ) afac, double precision, dimension( * ) ainv, double precision, dimension( * ) b, double precision, dimension( * ) bsav, double precision, dimension( * ) xact, double precision, dimension( * ) x, double precision, dimension( * ) arf, double precision, dimension( * ) arfinv, double precision, dimension( * ) d_work_dlatms, double precision, dimension( * ) d_work_dpot01, double precision, dimension( * ) d_temp_dpot02, double precision, dimension( * ) d_temp_dpot03, double precision, dimension( * ) d_work_dlansy, double precision, dimension( * ) d_work_dpot02, double precision, dimension( * ) d_work_dpot03)

DDRVRFP

Purpose:

!>
!> DDRVRFP tests the LAPACK RFP routines:
!>     DPFTRF, DPFTRS, and DPFTRI.
!>
!> 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 DTRTTF and
!> DTFTTR.
!>
!> 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 DPFTRF, 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 DPFTRF 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 DOUBLE PRECISION
!>                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 DOUBLE PRECISION array, dimension (NMAX*NMAX)
!>

ASAV

!>          ASAV is DOUBLE PRECISION array, dimension (NMAX*NMAX)
!>

AFAC

!>          AFAC is DOUBLE PRECISION array, dimension (NMAX*NMAX)
!>

AINV

!>          AINV is DOUBLE PRECISION array, dimension (NMAX*NMAX)
!>

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

BSAV

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

XACT

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

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

ARF

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

ARFINV

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

D_WORK_DLATMS

!>          D_WORK_DLATMS is DOUBLE PRECISION array, dimension ( 3*NMAX )
!>

D_WORK_DPOT01

!>          D_WORK_DPOT01 is DOUBLE PRECISION array, dimension ( NMAX )
!>

D_TEMP_DPOT02

!>          D_TEMP_DPOT02 is DOUBLE PRECISION array, dimension ( NMAX*MAXRHS )
!>

D_TEMP_DPOT03

!>          D_TEMP_DPOT03 is DOUBLE PRECISION array, dimension ( NMAX*NMAX )
!>

D_WORK_DLANSY

!>          D_WORK_DLANSY is DOUBLE PRECISION array, dimension ( NMAX )
!>

D_WORK_DPOT02

!>          D_WORK_DPOT02 is DOUBLE PRECISION array, dimension ( NMAX )
!>

D_WORK_DPOT03

!>          D_WORK_DPOT03 is DOUBLE PRECISION array, dimension ( NMAX )
!>

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.

Definition at line 232 of file ddrvrfp.f.