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. Diagonal
• 2. Random, CNDNUM = 2
• *3. First row and column zero
• *4. Last row and column zero
• *5. Middle row and column zero
• 6. Random, CNDNUM = sqrt(0.1/EPS)
• 7. Random, CNDNUM = 0.1/EPS
• 8. Scaled near underflow
• 9. 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 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 — 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 💅✨