TESTING/LIN/dqrt03.f(3)

πŸ“¦
version 3.12.1 Β· 2 versions
, in arch rolling

πŸ’– NAME πŸ’–

okay bestie, welcome to the page for TESTING/LIN/dqrt03.f β€” and yes, that's the whole name, no cap. This one lives in the Library Functions Manual and it's giving "LAPACK testing routine realness" 🫢✨


πŸ¦‹ SYNOPSIS πŸ¦‹

Here's the cast list, fr fr:

   Functions/Subroutines
       subroutine dqrt03 (m, n, k, af, c, cc, q, lda, tau, work, lwork, rwork,
           result)
           DQRT03

🌈 Function/Subroutine Documentation 🌈

This is the full signature, character-for-character, because we do not mess with signatures πŸ’…:

subroutine dqrt03 (integer m, integer n, integer k, double precision,
    dimension( lda, * ) af, double precision, dimension( lda, * ) c, double
    precision, dimension( lda, * ) cc, double precision, dimension( lda, * )
    q, integer lda, double precision, dimension( * ) tau, double precision,
    dimension( lwork ) work, integer lwork, double precision, dimension( * )
    rwork, double precision, dimension( * ) result)

DQRT03

🎯 Purpose 🎯

DQRT03 is the vibe check for DORMQR, which computes Q*C, Q'*C, C*Q or C*Q'.

It compares the results of a call to DORMQR with the results of forming Q explicitly by a call to DORGQR and then performing matrix multiplication by a call to DGEMM. Two techniques, one showdown 😭πŸ”₯

🧸 Parameters 🧸

  • M β€” M is INTEGER. The order of the orthogonal matrix Q. M >= 0.
  • N β€” N is INTEGER. The number of rows or columns of the matrix C; C is m-by-n if Q is applied from the left, or n-by-m if Q is applied from the right. N >= 0.
  • K β€” K is INTEGER. The number of elementary reflectors whose product defines the orthogonal matrix Q. M >= K >= 0.
  • AF β€” AF is DOUBLE PRECISION array, dimension (LDA,N). Details of the QR factorization of an m-by-n matrix, as returned by DGEQRF. See DGEQRF for further details.
  • C β€” C is DOUBLE PRECISION array, dimension (LDA,N).
  • CC β€” CC is DOUBLE PRECISION array, dimension (LDA,N).
  • Q β€” Q is DOUBLE PRECISION array, dimension (LDA,M).
  • LDA β€” LDA is INTEGER. The leading dimension of the arrays AF, C, CC, and Q.
  • TAU β€” TAU is DOUBLE PRECISION array, dimension (min(M,N)). The scalar factors of the elementary reflectors corresponding to the QR factorization in AF.
  • WORK β€” WORK is DOUBLE PRECISION array, dimension (LWORK).
  • LWORK β€” LWORK is INTEGER. The length of WORK. LWORK must be at least M, and should be M*NB, where NB is the blocksize for this environment.
  • RWORK β€” RWORK is DOUBLE PRECISION array, dimension (M).
  • RESULT β€” RESULT is DOUBLE PRECISION array, dimension (4). The test ratios compare two techniques for multiplying a random matrix C by an m-by-m orthogonal matrix Q.
RESULT(1) = norm( Q*C - Q*C )  / ( M * norm(C) * EPS )
RESULT(2) = norm( C*Q - C*Q )  / ( M * norm(C) * EPS )
RESULT(3) = norm( Q'*C - Q'*C )/ ( M * norm(C) * EPS )
RESULT(4) = norm( C*Q' - C*Q' )/ ( M * norm(C) * EPS )

πŸ‘‘ Author πŸ‘‘

The academic credits, as listed:

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

Definition at line 134 of file dqrt03.f.


πŸ“œ Author πŸ“œ

Generated automatically by Doxygen for LAPACK from the source code.


okay that's the whole page, bestie πŸ«ΆπŸ’– This cutie comes from TESTING/LIN/dqrt03.f(3) with footer energy from LAPACK Version 3.12.0 β€” generated automatically by Doxygen for LAPACK from the source code, periodt πŸŽ€βœ¨

πŸ“¦ where this page comes from