TESTING/LIN/sdrvrfp.f(3)
πΈ TESTING/LIN/sdrvrfp.f(3) β SDRVRFP π
okay so bestie, this is the Library Functions Manual page for the LAPACK test driver sdrvrfp β welcome to the group chat β¨
πΈ NAME πΈ
TESTING/LIN/sdrvrfp.f β the LAPACK test driver known as SDRVRFP, and it is literally here to test the RFP routines π
β¨ SYNOPSIS β¨
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
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
SDRVRFP
π§ Purpose
okay so SDRVRFP tests the LAPACK RFP routines: SPFTRF, SPFTRS, and SPFTRI π
this testing routine follows the same tests as DDRVPO (the test for the full format Symmetric Positive Definite solver) β no cap, that's the vibe π
the tests are performed in Full Format, and conversion back and forth from full format to RFP format is performed using the routines STRTTF and STFTTR β¨
first up, a specific matrix A of size N is created. there are nine types of different matrixes possible, and they are literally:
1.Diagonal2.Random,CNDNUM = 2*3.First row and column zero*4.Last row and column zero*5.Middle row and column zero6.Random,CNDNUM = sqrt(0.1/EPS)7.Random,CNDNUM = 0.1/EPS8.Scaled near underflow9.Scaled near overflow
(* - 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βNN is INTEGER; the number of values ofNcontained in the vectorNVALπ’NVALβNVAL is INTEGER array, dimension (NN); the values of the matrix dimensionNπNNSβNNS is INTEGER; the number of values ofNRHScontained in the vectorNSVALβ¨NSVALβNSVAL is INTEGER array, dimension (NNS); the values of the number of right-hand sidesNRHSπNNTβNNT is INTEGER; the number of values ofMATRIX TYPEcontained in the vectorNTVALπ¦NTVALβNTVAL is INTEGER array, dimension (NNT); the values of matrix type (between0and9forPO/PP/PFmatrices) πΈTHRESHβTHRESH is REAL; the threshold value for the test ratios. A result is included in the output file ifRESULT >= THRESH. To have every test ratio printed, useTHRESH = 0π₯Aβ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βB is REAL array, dimension (NMAX*MAXRHS)π¦BSAVβBSAV is REAL array, dimension (NMAX*MAXRHS)πΈXACTβXACT is REAL array, dimension (NMAX*MAXRHS)π₯Xβ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
Definition at line 232 of file sdrvrfp.f.
π«Ά Author (auto-generated)
Generated automatically by Doxygen for LAPACK from the source code.
and that's the whole page, bestie! the footer's version/source energy is LAPACK Version 3.12.0 for page TESTING/LIN/sdrvrfp.f(3) β no extra date shown in this footer, just that version/source stamp. byeee π
β¨
π¦ where this page comes from
- arch rolling: lapack-doc 3.12.1-2 Β· custom
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms