ZLARZT - Online Linux Manual Page

Section : 1
Updated : November 2008
Source : LAPACK routine (version 3.2)
Note : LAPACK routine (version 3.2)

NAMEZLARZT - forms the triangular factor T of a complex block reflector H of order > n, which is defined as a product of k elementary reflectors

SYNOPSISSUBROUTINE ZLARZT(  DIRECT, STOREV, N, K, V, LDV, TAU, T, LDT )  CHARACTER DIRECT, STOREV  INTEGER K, LDT, LDV, N  COMPLEX*16 T( LDT, * ), TAU( * ), V( LDV, * )

PURPOSEZLARZT forms the triangular factor T of a complex block reflector H of order > n, which is defined as a product of k elementary reflectors. If DIRECT = 'F', H = H(1) H(2) . . . H(k) and T is upper triangular; If DIRECT = 'B', H = H(k) . . . H(2) H(1) and T is lower triangular. If STOREV = 'C', the vector which defines the elementary reflector H(i) is stored in the i-th column of the array V, and

   H = I - V * T * V'
If STOREV = 'R', the vector which defines the elementary reflector H(i) is stored in the i-th row of the array V, and

   H = I - V' * T * V
Currently, only STOREV = 'R' and DIRECT = 'B' are supported.

ARGUMENTSDIRECT (input) CHARACTER*1  Specifies the order in which the elementary reflectors are multiplied to form the block reflector:
= 'F': H = H(1) H(2) . . . H(k) (Forward, not supported yet)
= 'B': H = H(k) . . . H(2) H(1) (Backward)
STOREV (input) CHARACTER*1  Specifies how the vectors which define the elementary reflectors are stored (see also Further Details):
= 'R': rowwise
N (input) INTEGER  The order of the block reflector H. N >= 0. K (input) INTEGER  The order of the triangular factor T (= the number of elementary reflectors). K >= 1. V (input/output) COMPLEX*16 array, dimension  (LDV,K) if STOREV = 'C' (LDV,N) if STOREV = 'R' The matrix V. See further details. LDV (input) INTEGER  The leading dimension of the array V. If STOREV = 'C', LDV >= max(1,N); if STOREV = 'R', LDV >= K. TAU (input) COMPLEX*16 array, dimension (K)  TAU(i) must contain the scalar factor of the elementary reflector H(i). T (output) COMPLEX*16 array, dimension (LDT,K)  The k by k triangular factor T of the block reflector. If DIRECT = 'F', T is upper triangular; if DIRECT = 'B', T is lower triangular. The rest of the array is not used. LDT (input) INTEGER  The leading dimension of the array T. LDT >= K.

FURTHER DETAILSBased on contributions by

  A. Petitet, Computer Science Dept., Univ. of Tenn., Knoxville, USA The shape of the matrix V and the storage of the vectors which define the H(i) is best illustrated by the following example with n = 5 and k = 3. The elements equal to 1 are not stored; the corresponding array elements are modified but restored on exit. The rest of the array is not used.
DIRECT = 'F' and STOREV = 'C': DIRECT = 'F' and STOREV = 'R':
                                            ______V_____

       ( v1 v2 v3 ) / ( v1 v2 v3 ) ( v1 v1 v1 v1 v1 . . . . 1 )
   V = ( v1 v2 v3 ) ( v2 v2 v2 v2 v2 . . . 1 )
       ( v1 v2 v3 ) ( v3 v3 v3 v3 v3 . . 1 )
       ( v1 v2 v3 )

          . . .

          . . .

          1 . .

             1 .

                1
DIRECT = 'B' and STOREV = 'C': DIRECT = 'B' and STOREV = 'R':
                                                      ______V_____
          1 / . 1 ( 1 . . . . v1 v1 v1 v1 v1 )
          . . 1 ( . 1 . . . v2 v2 v2 v2 v2 )
          . . . ( . . 1 . . v3 v3 v3 v3 v3 )
          . . .

       ( v1 v2 v3 )

       ( v1 v2 v3 )

   V = ( v1 v2 v3 )

       ( v1 v2 v3 )

       ( v1 v2 v3 )
0
Johanes Gumabo
Data Size   :   13,382 byte
man-zlarzt.lBuild   :   2025-03-22, 13:26   :  
Visitor Screen   :   x
Visitor Counter ( page / site )   :   3 / 1,209,981
Visitor ID   :     :  
Visitor IP   :   18.119.28.119   :  
Visitor Provider   :   AMAZON-02   :  
Provider Position ( lat x lon )   :   39.962500 x -83.006100   :   x
Provider Accuracy Radius ( km )   :   1000   :  
Provider City   :   Columbus   :  
Provider Province   :   Ohio ,   :   ,
Provider Country   :   United States   :  
Provider Continent   :   North America   :  
Visitor Recorder   :   Version   :  
Visitor Recorder   :   Library   :  
Online Linux Manual Page   :   Version   :   Online Linux Manual Page - Fedora.40 - march=x86-64 - mtune=generic - 25.03.22
Online Linux Manual Page   :   Library   :   lib_c - 24.10.03 - march=x86-64 - mtune=generic - Fedora.40
Online Linux Manual Page   :   Library   :   lib_m - 24.10.03 - march=x86-64 - mtune=generic - Fedora.40
Data Base   :   Version   :   Online Linux Manual Page Database - 24.04.13 - march=x86-64 - mtune=generic - fedora-38
Data Base   :   Library   :   lib_c - 23.02.07 - march=x86-64 - mtune=generic - fedora.36


Very long time ago, I have the best tutor, Wenzel Svojanovsky. If someone knows the email address of Wenzel Svojanovsky, please send an email to johanesgumabo@gmail.com.
Help me, linux0001.com will expire on July 16, 2025. I have no money to renew it. View detail

If error, please print screen and send to johanes_gumabo@yahoo.co.id
Under development. Support me via PayPal.