/home/abui.../zdrvrfp.f(3)

🌸 NAME 🌸

/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/TESTING/LIN/zdrvrfp.f

okay bestie, this is the zdrvrfp.f(3) page from the Library Functions Manual — and yes, we're gonna get through all of it together. 💖

---

✨ SYNOPSIS ✨

🌈 Functions/Subroutines 🌈

       subroutine ZDRVRFP (nout, nn, nval, nns, nsval, nnt, ntval, thresh, a,
           asav, afac, ainv, b, bsav, xact, x, arf, arfinv, z_work_zlatms,
           z_work_zpot02, z_work_zpot03, d_work_zlatms, d_work_zlanhe,
           d_work_zpot01, d_work_zpot02, d_work_zpot03)
           ZDRVRFP

---

💖 Function/Subroutine Documentation 💖

   subroutine ZDRVRFP (integer nout, integer nn, integer, dimension( nn ) nval,
       integer nns, integer, dimension( nns ) nsval, integer nnt, integer,
       dimension( nnt ) ntval, double precision thresh, complex*16, dimension(
       * ) a, complex*16, dimension( * ) asav, complex*16, dimension( * ) afac,
       complex*16, dimension( * ) ainv, complex*16, dimension( * ) b,
       complex*16, dimension( * ) bsav, complex*16, dimension( * ) xact,
       complex*16, dimension( * ) x, complex*16, dimension( * ) arf,
       complex*16, dimension( * ) arfinv, complex*16, dimension( * )
       z_work_zlatms, complex*16, dimension( * ) z_work_zpot02, complex*16,
       dimension( * ) z_work_zpot03, double precision, dimension( * )
       d_work_zlatms, double precision, dimension( * ) d_work_zlanhe, double
       precision, dimension( * ) d_work_zpot01, double precision, dimension( *
       ) d_work_zpot02, double precision, dimension( * ) d_work_zpot03)
       ZDRVRFP

🎯 Purpose 🎯

okay so ZDRVRFP is the vibe-check routine for the LAPACK RFP routines: ZPFTRF, ZPFTRS, and ZPFTRI. It runs the same tests as ZDRVPO — which is the test for the full-format Symmetric Positive Definite solver. slay. 💅

The tests are performed in Full Format, and converting back and forth between full format and RFP format is done using ZTRTTF and ZTFTTR.

First, a specific matrix A of size N is created. There are nine types of different matrixes possible:

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 ZPFTRF, 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 ZPFTRF 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 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 — A is COMPLEX*16 array, dimension (NMAX*NMAX)
• ASAV — ASAV is COMPLEX*16 array, dimension (NMAX*NMAX)
• AFAC — AFAC is COMPLEX*16 array, dimension (NMAX*NMAX)
• AINV — AINV is COMPLEX*16 array, dimension (NMAX*NMAX)
• B — B is COMPLEX*16 array, dimension (NMAX*MAXRHS)
• BSAV — BSAV is COMPLEX*16 array, dimension (NMAX*MAXRHS)
• XACT — XACT is COMPLEX*16 array, dimension (NMAX*MAXRHS)
• X — X is COMPLEX*16 array, dimension (NMAX*MAXRHS)
• ARF — ARF is COMPLEX*16 array, dimension ((NMAX*(NMAX+1))/2)
• ARFINV — ARFINV is COMPLEX*16 array, dimension ((NMAX*(NMAX+1))/2)
• Z_WORK_ZLATMS — Z_WORK_ZLATMS is COMPLEX*16 array, dimension ( 3*NMAX )
• Z_WORK_ZPOT02 — Z_WORK_ZPOT02 is COMPLEX*16 array, dimension ( NMAX*MAXRHS )
• Z_WORK_ZPOT03 — Z_WORK_ZPOT03 is COMPLEX*16 array, dimension ( NMAX*NMAX )
• D_WORK_ZLATMS — D_WORK_ZLATMS is DOUBLE PRECISION array, dimension ( NMAX )
• D_WORK_ZLANHE — D_WORK_ZLANHE is DOUBLE PRECISION array, dimension ( NMAX )
• D_WORK_ZPOT01 — D_WORK_ZPOT01 is DOUBLE PRECISION array, dimension ( NMAX )
• D_WORK_ZPOT02 — D_WORK_ZPOT02 is DOUBLE PRECISION array, dimension ( NMAX )
• D_WORK_ZPOT03 — D_WORK_ZPOT03 is DOUBLE PRECISION array, dimension ( NMAX )

🌟 Author 🌟

• Univ. of Tennessee
• Univ. of California Berkeley
• Univ. of Colorado Denver
• NAG Ltd.

Definition at line 238 of file zdrvrfp.f.

---

🫶 Author 🫶

Generated automatically by Doxygen for LAPACK from the source code.

---

mwah 💋 this page is part of LAPACK Version 3.12.0, and its footer path is /home/abui.../zdrvrfp.f(3). iykyk 💅