Routine PCLARFT() Documented By Universal Report

Routine: PCLARFT()

# lines:	543
# code:	543
# comment:	0
# blank:	0
# Variables:	47
# Callers:	12
# Callings:	0
# Words:	254
# Keywords:	168

..
     .. Array Arguments ..
     ..
  Purpose
  =======
  PCLARFT 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 distributed matrix 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 distributed matrix V, and
     H  =  I - V' * T * V
  Notes
  =====
  Each global data object is described by an associated description
  vector.  This vector stores the information required to establish
  the mapping between an object element and its corresponding process
  and memory location.
  Let A be a generic term for any 2D block cyclicly distributed array.
  Such a global array has an associated description vector DESCA.
  In the following comments, the character _ should be read as
  "of the global array".
  NOTATION        STORED IN      EXPLANATION
  --------------- -------------- --------------------------------------
  DTYPE_A(global) DESCA( DTYPE_ )The descriptor type.  In this case,
                                 DTYPE_A = 1.
  CTXT_A (global) DESCA( CTXT_ ) The BLACS context handle, indicating
                                 the BLACS process grid A is distribu-
                                 ted over. The context itself is glo-
                                 bal, but the handle (the integer
                                 value) may vary.
  M_A    (global) DESCA( M_ )    The number of rows in the global
                                 array A.
  N_A    (global) DESCA( N_ )    The number of columns in the global
                                 array A.
  MB_A   (global) DESCA( MB_ )   The blocking factor used to distribute
                                 the rows of the array.
  NB_A   (global) DESCA( NB_ )   The blocking factor used to distribute
                                 the columns of the array.
  RSRC_A (global) DESCA( RSRC_ ) The process row over which the first
                                 row of the array A is distributed.
  CSRC_A (global) DESCA( CSRC_ ) The process column over which the
                                 first column of the array A is
                                 distributed.
  LLD_A  (local)  DESCA( LLD_ )  The leading dimension of the local
                                 array.  LLD_A >= MAX(1,LOCr(M_A)).
  Let K be the number of rows or columns of a distributed matrix,
  and assume that its process grid has dimension p x q.
  LOCr( K ) denotes the number of elements of K that a process
  would receive if K were distributed over the p processes of its
  process column.
  Similarly, LOCc( K ) denotes the number of elements of K that a
  process would receive if K were distributed over the q processes of
  its process row.
  The values of LOCr() and LOCc() may be determined via a call to the
  ScaLAPACK tool function, NUMROC:
          LOCr( M ) = NUMROC( M, MB_A, MYROW, RSRC_A, NPROW ),
          LOCc( N ) = NUMROC( N, NB_A, MYCOL, CSRC_A, NPCOL ).
  An upper bound for these quantities may be computed by:
          LOCr( M ) <= ceil( ceil(M/MB_A)/NPROW )*MB_A
          LOCc( N ) <= ceil( ceil(N/NB_A)/NPCOL )*NB_A
  Arguments
  =========
  DIRECT  (global 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)
          = 'B': H = H(k) . . . H(2) H(1) (Backward)
  STOREV  (global input) CHARACTER*1
          Specifies how the vectors which define the elementary
          reflectors are stored (see also Further Details):
          = 'C': columnwise
          = 'R': rowwise
  N       (global input) INTEGER
          The order of the block reflector H. N >= 0.
  K       (global input) INTEGER
          The order of the triangular factor T (= the number of
          elementary reflectors). 1 <= K <= MB_V (= NB_V).
  V       (input/output) COMPLEX pointer into the local memory
          to an array of local dimension (LOCr(IV+N-1),LOCc(JV+K-1))
          if STOREV = 'C', and (LOCr(IV+K-1),LOCc(JV+N-1)) if
          STOREV = 'R'. The distributed matrix V contains the
          Householder vectors. See further details.
  IV      (global input) INTEGER
          The row index in the global array V indicating the first
          row of sub( V ).
  JV      (global input) INTEGER
          The column index in the global array V indicating the
          first column of sub( V ).
  DESCV   (global and local input) INTEGER array of dimension DLEN_.
          The array descriptor for the distributed matrix V.
  TAU     (local input) COMPLEX, array, dimension LOCr(IV+K-1)
          if INCV = M_V, and LOCc(JV+K-1) otherwise. This array
          contains the Householder scalars related to the Householder
          vectors.  TAU is tied to the distributed matrix V.
  T       (local output) COMPLEX array, dimension (NB_V,NB_V)
          if STOREV = 'Col', and (MB_V,MB_V) otherwise. It contains
          the k-by-k triangular factor of the block reflector asso-
          ciated with V. If DIRECT = 'F', T is upper triangular;
          if DIRECT = 'B', T is lower triangular.
  WORK    (local workspace) COMPLEX array,
                                          dimension (K*(K-1)/2)
  Further Details
  ===============
  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( IV:IV+N-1,    (  1       )    V( IV:IV+K-1,    (  1 v1 v1 v1 v1 )
     JV:JV+K-1 ) = ( v1  1    )       JV:JV+N-1 ) = (     1 v2 v2 v2 )
                   ( v1 v2  1 )                     (        1 v3 v3 )
                   ( v1 v2 v3 )
                   ( v1 v2 v3 )
  DIRECT = 'B' and STOREV = 'C':   DIRECT = 'B' and STOREV = 'R':
  V( IV:IV+N-1,    ( v1 v2 v3 )    V( IV:IV+K-1,    ( v1 v1  1       )
     JV:JV+K-1 ) = ( v1 v2 v3 )       JV:JV+N-1 ) = ( v2 v2 v2  1    )
                   (  1 v2 v3 )                     ( v3 v3 v3 v3  1 )
                   (     1 v3 )
                   (        1 )
  =====================================================================
     .. Parameters ..

Display dynamic version

FontSize: - +

001        SUBROUTINE PCLARFT( DIRECT , STOREV , N , K , V , IV , JV , DESCV , TAU ,
002       $T , WORK )
003  
004  *     -- ScaLAPACK auxiliary routine(version 1.7) --
005  *     University of Tennessee , Knoxville , Oak Ridge National Laboratory ,
006  *     and University of California , Berkeley.
007  *     May 1 , 1997
008  
009  *     .. Scalar Arguments ..
010        CHARACTER DIRECT , STOREV
011        INTEGER IV , JV , K , N
012        INTEGER BLOCK_CYCLIC_2D , CSRC_ , CTXT_ , DLEN_ , DTYPE_ ,
013       $LLD_ , MB_ , M_ , NB_ , N_ , RSRC_
014        PARAMETER( BLOCK_CYCLIC_2D = 1 , DLEN_ = 9 , DTYPE_ = 1 ,
015       $CTXT_ = 2 , M_ = 3 , N_ = 4 , MB_ = 5 , NB_ = 6 ,
016       $RSRC_ = 7 , CSRC_ = 8 , LLD_ = 9 )
017        COMPLEX ONE , ZERO
018        PARAMETER( ONE =( 1.0E + 0 , 0.0E + 0 ) ,
019       $ZERO =( 0.0E + 0 , 0.0E + 0 ) )
020  *     ..
021  *     .. Local Scalars ..
022        LOGICAL FORWARD
023        INTEGER ICOFF , ICTXT , II , IIV , IROFF , IVCOL , IVROW ,
024       $ITMP0 , ITMP1 , IW , JJ , JJV , LDV , MICOL , MIROW ,
025       $MYCOL , MYROW , NP , NPCOL , NPROW , NQ
026        COMPLEX VII
027  *     ..
028  *     .. External Subroutines ..
029        EXTERNAL BLACS_GRIDINFO , CCOPY , CGEMV , CGSUM2D ,
030       $CLACGV , CLASET , CTRMV , INFOG2L
031  *     ..
032  *     .. External Functions ..
033        LOGICAL LSAME
034        INTEGER INDXG2P , NUMROC
035        EXTERNAL INDXG2P , LSAME , NUMROC
036  *     ..
037  *     .. Intrinsic Functions ..
038        INTRINSIC MOD
039  *     ..
040  *     .. Executable Statements ..
041  
042  *     Quick return if possible
043  
044        IF( N.LE.0 .OR. K.LE.0 )
044               
045       $    RETURN
046  
047            ICTXT = DESCV( CTXT_ )
048            CALL BLACS_GRIDINFO( ICTXT , NPROW , NPCOL , MYROW , MYCOL )
049  
050            FORWARD = LSAME( DIRECT , 'F' )
051            CALL INFOG2L( IV , JV , DESCV , NPROW , NPCOL , MYROW , MYCOL ,
052       $    IIV , JJV , IVROW , IVCOL )
053  
054            IF( LSAME( STOREV , 'C' ) .AND. MYCOL.EQ.IVCOL ) THEN
055  
055                   
056                IW = 1
057                LDV = DESCV( LLD_ )
058                IROFF = MOD( IV - 1 , DESCV( MB_ ) )
059  
060                IF( FORWARD ) THEN
061  
062  *                 DIRECT = 'Forward' , STOREV = 'Columnwise'
063  
063                       
064                    NP = NUMROC( N + IROFF , DESCV( MB_ ) , MYROW , IVROW , NPROW )
065                    IF( MYROW.EQ.IVROW ) THEN
065                           
066                        NP = NP - IROFF
067                        II = IIV + 1
068                    ELSE
068                           
069                        II = IIV
070                    END IF
071                    IF( IROFF + 1.EQ.DESCV( MB_ ) ) THEN
071                           
072                        MIROW = MOD( IVROW + 1 , NPROW )
073                    ELSE
073                           
074                        MIROW = IVROW
075                    END IF
076                    ITMP0 = 0
077  
078                    DO 10 JJ = JJV + 1 , JJV + K - 1
079  
079                           
080                        IF( MYROW.EQ.MIROW ) THEN
080                               
081                            VII = V( II + (JJ - 1)*LDV )
082                            V( II + (JJ - 1)*LDV ) = ONE
083                        END IF
084  
085  *                     T(1 : i - 1 , i) = - tau( jv + i - 1 ) *
086  *                     V(iv + i - 1 : iv + n - 1 , jv : jv + i - 2)' * V(iv + i - 1 : iv + n - 1 , jv + i - 1)
087  
088                        ITMP0 = ITMP0 + 1
089                        IF( NP - II + IIV.GT.0 ) THEN
089                               
090                            CALL CGEMV( 'Conjugate transpose' , NP - II + IIV , ITMP0 ,
091       $                    - TAU( JJ ) , V( II + (JJV - 1)*LDV ) , LDV ,
092       $                    V( II + (JJ - 1)*LDV ) , 1 , ZERO ,
093       $                    WORK( IW ) , 1 )
094                        ELSE
094                               
095                            CALL CLASET( 'All' , ITMP0 , 1 , ZERO , ZERO , WORK( IW ) ,
096       $                    ITMP0 )
097                        END IF
098  
099                        IW = IW + ITMP0
100                        IF( MYROW.EQ.MIROW ) THEN
100                               
101                            V( II + (JJ - 1)*LDV ) = VII
102                            II = II + 1
103                        END IF
104  
105                        IF( MOD( IV + ITMP0 , DESCV( MB_ ) ).EQ.0 )
105                               
106       $                    MIROW = MOD( MIROW + 1 , NPROW )
107  
108     10             CONTINUE
109  
109                       
110                    CALL CGSUM2D( ICTXT , 'Columnwise' , ' ' , IW - 1 , 1 , WORK , IW - 1 ,
111       $            IVROW , MYCOL )
112  
113                    IF( MYROW.EQ.IVROW ) THEN
114  
114                           
115                        IW = 1
116                        ITMP0 = 0
117                        ITMP1 = 1
118  
119                        T( ITMP1 ) = TAU( JJV )
120  
121                        DO 20 JJ = JJV + 1 , JJV + K - 1
122  
123  *                         T(1 : j - 1 , j) = T(1 : j - 1 , 1 : j - 1) * T(1 : j - 1 , j)
124  
124                               
125                            ITMP0 = ITMP0 + 1
126                            ITMP1 = ITMP1 + DESCV( NB_ )
127                            CALL CCOPY( ITMP0 , WORK( IW ) , 1 , T( ITMP1 ) , 1 )
128                            IW = IW + ITMP0
129  
130                            CALL CTRMV( 'Upper' , 'No transpose' , 'Non - unit' ,
131       $                    ITMP0 , T , DESCV( NB_ ) , T( ITMP1 ) , 1 )
132                            T(ITMP1 + ITMP0) = TAU( JJ )
133  
134     20                 CONTINUE
135  
135                       
136                    END IF
137  
138                ELSE
139  
140  *                 DIRECT = 'Backward' , STOREV = 'Columnwise'
141  
141                       
142                    NP = NUMROC( N + IROFF - 1 , DESCV( MB_ ) , MYROW , IVROW , NPROW )
143                    IF( MYROW.EQ.IVROW )
143                           
144       $                NP = NP - IROFF
145                        MIROW = INDXG2P( IV + N - 2 , DESCV( MB_ ) , MYROW ,
146       $                DESCV( RSRC_ ) , NPROW )
147                        II = IIV + NP - 1
148                        ITMP0 = 0
149  
150                        DO 30 JJ = JJV + K - 2 , JJV , - 1
151  
151                               
152                            IF( MYROW.EQ.MIROW ) THEN
152                                   
153                                VII = V( II + (JJ - 1)*LDV )
154                                V( II + (JJ - 1)*LDV ) = ONE
155                            END IF
156  
157  *                         T(1 : i - 1 , i) = - tau( jv + i - 1 ) *
158  *                         V(iv : iv + n - k + i - 1 , jv + i : jv + k - 1)' * V(iv : iv + n - k + i - 1 , jv + i - 1)
159  
160                            ITMP0 = ITMP0 + 1
161                            IF( II - IIV + 1.GT.0 ) THEN
161                                   
162                                CALL CGEMV( 'Conjugate transpose' , II - IIV + 1 , ITMP0 ,
163       $                        - TAU( JJ ) , V( IIV + JJ*LDV ) , LDV ,
164       $                        V( IIV + (JJ - 1)*LDV ) , 1 , ZERO ,
165       $                        WORK( IW ) , 1 )
166                            ELSE
166                                   
167                                CALL CLASET( 'All' , ITMP0 , 1 , ZERO , ZERO , WORK( IW ) ,
168       $                        ITMP0 )
169                            END IF
170  
171                            IW = IW + ITMP0
172                            IF( MYROW.EQ.MIROW ) THEN
172                                   
173                                V( II + (JJ - 1)*LDV ) = VII
174                                II = II - 1
175                            END IF
176  
177                            IF( MOD( IV + N - ITMP0 - 2 , DESCV(MB_) ).EQ.0 )
177                                   
178       $                        MIROW = MOD( MIROW + NPROW - 1 , NPROW )
179  
180     30                 CONTINUE
181  
181                           
182                        CALL CGSUM2D( ICTXT , 'Columnwise' , ' ' , IW - 1 , 1 , WORK , IW - 1 ,
183       $                IVROW , MYCOL )
184  
185                        IF( MYROW.EQ.IVROW ) THEN
186  
186                               
187                            IW = 1
188                            ITMP0 = 0
189                            ITMP1 = K + 1 + (K - 1) * DESCV( NB_ )
190  
191                            T( ITMP1 - 1 ) = TAU( JJV + K - 1 )
192  
193                            DO 40 JJ = JJV + K - 2 , JJV , - 1
194  
195  *                             T(j + 1 : k , j) = T(j + 1 : k , j + 1 : k) * T(j + 1 : k , j)
196  
196                                   
197                                ITMP0 = ITMP0 + 1
198                                ITMP1 = ITMP1 - DESCV( NB_ ) - 1
199                                CALL CCOPY( ITMP0 , WORK( IW ) , 1 , T( ITMP1 ) , 1 )
200                                IW = IW + ITMP0
201  
202                                CALL CTRMV( 'Lower' , 'No transpose' , 'Non - unit' ,
203       $                        ITMP0 , T( ITMP1 + DESCV( NB_ ) ) ,
204       $                        DESCV( NB_ ) , T( ITMP1 ) , 1 )
205                                T( ITMP1 - 1 ) = TAU( JJ )
206  
207     40                     CONTINUE
208  
208                           
209                        END IF
210  
211                    END IF
212  
213                ELSE IF( LSAME( STOREV , 'R' ) .AND. MYROW.EQ.IVROW ) THEN
214  
214                       
215                    IW = 1
216                    LDV = DESCV( LLD_ )
217                    ICOFF = MOD( JV - 1 , DESCV( NB_ ) )
218  
219                    IF( FORWARD ) THEN
220  
221  *                     DIRECT = 'Forward' , STOREV = 'Rowwise'
222  
222                           
223                        NQ = NUMROC( N + ICOFF , DESCV( NB_ ) , MYCOL , IVCOL , NPCOL )
224                        IF( MYCOL.EQ.IVCOL ) THEN
224                               
225                            NQ = NQ - ICOFF
226                            JJ = JJV + 1
227                        ELSE
227                               
228                            JJ = JJV
229                        END IF
230                        IF( ICOFF + 1.EQ.DESCV( NB_ ) ) THEN
230                               
231                            MICOL = MOD( IVCOL + 1 , NPCOL )
232                        ELSE
232                               
233                            MICOL = IVCOL
234                        END IF
235                        ITMP0 = 0
236  
237                        DO 50 II = IIV + 1 , IIV + K - 1
238  
238                               
239                            IF( MYCOL.EQ.MICOL ) THEN
239                                   
240                                VII = V( II + (JJ - 1)*LDV )
241                                V( II + (JJ - 1)*LDV ) = ONE
242                            END IF
243  
244  *                         T(1 : i - 1 , i) = - tau( iv + i - 1 ) *
245  *                         V(iv + i - 1 , jv + i - 1 : jv + n - 1) * V(iv : iv + i - 2 , jv + i - 1 : jv + n - 1)'
246  
247                            ITMP0 = ITMP0 + 1
248                            IF( NQ - JJ + JJV.GT.0 ) THEN
248                                   
249                                CALL CLACGV( NQ - JJ + JJV , V( II + (JJ - 1)*LDV ) , LDV )
250                                CALL CGEMV( 'No transpose' , ITMP0 , NQ - JJ + JJV ,
251       $                        - TAU(II) , V( IIV + (JJ - 1)*LDV ) , LDV ,
252       $                        V( II + (JJ - 1)*LDV ) , LDV , ZERO ,
253       $                        WORK( IW ) , 1 )
254                                CALL CLACGV( NQ - JJ + JJV , V( II + (JJ - 1)*LDV ) , LDV )
255                            ELSE
255                                   
256                                CALL CLASET( 'All' , ITMP0 , 1 , ZERO , ZERO ,
257       $                        WORK( IW ) , ITMP0 )
258                            END IF
259  
260                            IW = IW + ITMP0
261                            IF( MYCOL.EQ.MICOL ) THEN
261                                   
262                                V( II + (JJ - 1)*LDV ) = VII
263                                JJ = JJ + 1
264                            END IF
265  
266                            IF( MOD( JV + ITMP0 , DESCV( NB_ ) ).EQ.0 )
266                                   
267       $                        MICOL = MOD( MICOL + 1 , NPCOL )
268  
269     50                 CONTINUE
270  
270                           
271                        CALL CGSUM2D( ICTXT , 'Rowwise' , ' ' , IW - 1 , 1 , WORK , IW - 1 ,
272       $                MYROW , IVCOL )
273  
274                        IF( MYCOL.EQ.IVCOL ) THEN
275  
275                               
276                            IW = 1
277                            ITMP0 = 0
278                            ITMP1 = 1
279  
280                            T( ITMP1 ) = TAU( IIV )
281  
282                            DO 60 II = IIV + 1 , IIV + K - 1
283  
284  *                             T(1 : i - 1 , i) = T(1 : i - 1 , 1 : i - 1) * T(1 : i - 1 , i)
285  
285                                   
286                                ITMP0 = ITMP0 + 1
287                                ITMP1 = ITMP1 + DESCV( MB_ )
288                                CALL CCOPY( ITMP0 , WORK( IW ) , 1 , T( ITMP1 ) , 1 )
289                                IW = IW + ITMP0
290  
291                                CALL CTRMV( 'Upper' , 'No transpose' , 'Non - unit' ,
292       $                        ITMP0 , T , DESCV( MB_ ) , T( ITMP1 ) , 1 )
293                                T( ITMP1 + ITMP0 ) = TAU( II )
294  
295     60                     CONTINUE
296  
296                           
297                        END IF
298  
299                    ELSE
300  
301  *                     DIRECT = 'Backward' , STOREV = 'Rowwise'
302  
302                           
303                        NQ = NUMROC( N + ICOFF - 1 , DESCV( NB_ ) , MYCOL , IVCOL , NPCOL )
304                        IF( MYCOL.EQ.IVCOL )
304                               
305       $                    NQ = NQ - ICOFF
306                            MICOL = INDXG2P( JV + N - 2 , DESCV( NB_ ) , MYCOL ,
307       $                    DESCV( CSRC_ ) , NPCOL )
308                            JJ = JJV + NQ - 1
309                            ITMP0 = 0
310  
311                            DO 70 II = IIV + K - 2 , IIV , - 1
312  
312                                   
313                                IF( MYCOL.EQ.MICOL ) THEN
313                                       
314                                    VII = V( II + (JJ - 1)*LDV )
315                                    V( II + (JJ - 1)*LDV ) = ONE
316                                END IF
317  
318  *                             T(i + 1 : k , i) = - tau( iv + i - 1 ) *
319  *                             V(iv + i : iv + k - 1 , jv : jv + n - k + i - 1)' * V(iv + i - 1 , jv : jv + n - k + i - 1)'
320  
321                                ITMP0 = ITMP0 + 1
322                                IF( JJ - JJV + 1.GT.0 ) THEN
322                                       
323                                    CALL CLACGV( JJ - JJV + 1 , V( II + (JJV - 1)*LDV ) , LDV )
324                                    CALL CGEMV( 'No transpose' , ITMP0 , JJ - JJV + 1 ,
325       $                            - TAU( II ) , V( II + 1 + (JJV - 1)*LDV ) , LDV ,
326       $                            V( II + (JJV - 1)*LDV ) , LDV , ZERO ,
327       $                            WORK( IW ) , 1 )
328                                    CALL CLACGV( JJ - JJV + 1 , V( II + (JJV - 1)*LDV ) , LDV )
329                                ELSE
329                                       
330                                    CALL CLASET( 'All' , ITMP0 , 1 , ZERO , ZERO ,
331       $                            WORK( IW ) , ITMP0 )
332                                END IF
333  
334                                IW = IW + ITMP0
335                                IF( MYCOL.EQ.MICOL ) THEN
335                                       
336                                    V( II + (JJ - 1)*LDV ) = VII
337                                    JJ = JJ - 1
338                                END IF
339  
340                                IF( MOD( JV + N - ITMP0 - 2 , DESCV( NB_ ) ).EQ.0 )
340                                       
341       $                            MICOL = MOD( MICOL + NPCOL - 1 , NPCOL )
342  
343     70                     CONTINUE
344  
344                               
345                            CALL CGSUM2D( ICTXT , 'Rowwise' , ' ' , IW - 1 , 1 , WORK , IW - 1 ,
346       $                    MYROW , IVCOL )
347  
348                            IF( MYCOL.EQ.IVCOL ) THEN
349  
349                                   
350                                IW = 1
351                                ITMP0 = 0
352                                ITMP1 = K + 1 + (K - 1) * DESCV( MB_ )
353  
354                                T( ITMP1 - 1 ) = TAU( IIV + K - 1 )
355  
356                                DO 80 II = IIV + K - 2 , IIV , - 1
357  
358  *                                 T(i + 1 : k , i) = T(i + 1 : k , i + 1 : k) * T(i + 1 : k , i)
359  
359                                       
360                                    ITMP0 = ITMP0 + 1
361                                    ITMP1 = ITMP1 - DESCV( MB_ ) - 1
362                                    CALL CCOPY( ITMP0 , WORK( IW ) , 1 , T( ITMP1 ) , 1 )
363                                    IW = IW + ITMP0
364  
365                                    CALL CTRMV( 'Lower' , 'No transpose' , 'Non - unit' ,
366       $                            ITMP0 , T( ITMP1 + DESCV( MB_ ) ) ,
367       $                            DESCV( MB_ ) , T( ITMP1 ) , 1 )
368                                    T( ITMP1 - 1 ) = TAU( II )
369  
370     80                         CONTINUE
371  
371                               
372                            END IF
373  
374                        END IF
375  
376                    END IF
377  
378                    RETURN
379  
380  *                 End of PCLARFT
381  
382                END3396

Variables in Routine PCLARFT()

Summary Report
Data Type	Quantity	Size(byte)
CHARACTER	2	2
COMPLEX	3	12
INTEGER	38	152
LOGICAL	2	2
REAL	2	8
TOTAL	47	176

List of Variables

CHARACTER

DIRECT STOREV

COMPLEX

ONE VII ZERO

INTEGER

BLOCK_CYCLIC_2D CSRC_ CTXT_ DLEN_ DTYPE_
ICOFF ICTXT II IIV INDXG2P
IROFF ITMP0 ITMP1 IV IVCOL
IVROW IW JJ JJV JV
K LDV LLD_ M_ MB_
MICOL MIROW MYCOL MYROW N
N_ NB_ NP NPCOL NPROW
NQ NUMROC RSRC_

LOGICAL

FORWARD LSAME

REAL

T V

Variables Dependence Graph
Put the mouse over a right hand side variable to display the corresponding line of the dependence
-
- -
FORWARD <--- LSAMEFORWARD = LSAME( DIRECT, 'F' ), DIRECTFORWARD = LSAME( DIRECT, 'F' )

ICOFF <--- JVICOFF = MOD( JV-1, DESCV( NB_ ) ), NB_ICOFF = MOD( JV-1, DESCV( NB_ ) )

ICTXT <--- CTXT_ICTXT = DESCV( CTXT_ )

II <--- IIII = II + 1{2II = II - 1}, IIVII = IIV + NP - 1{2DO 50 II = IIV+1, IIV+K-1, 3DO 60 II = IIV+1, IIV+K-1, 4DO 70 II = IIV+K-2, IIV, -1, 5DO 80 II = IIV+K-2, IIV, -1, 6II = IIV + 1, 7II = IIV}, KDO 50 II = IIV+1, IIV+K-1{2DO 60 II = IIV+1, IIV+K-1, 3DO 70 II = IIV+K-2, IIV, -1, 4DO 80 II = IIV+K-2, IIV, -1}, NPII = IIV + NP - 1

IROFF <--- IVIROFF = MOD( IV-1, DESCV( MB_ ) ), MB_IROFF = MOD( IV-1, DESCV( MB_ ) )

ITMP0 <--- ITMP0ITMP0 = ITMP0 + 1{2ITMP0 = ITMP0 + 1, 3ITMP0 = ITMP0 + 1, 4ITMP0 = ITMP0 + 1, 5ITMP0 = ITMP0 + 1, 6ITMP0 = ITMP0 + 1, 7ITMP0 = ITMP0 + 1, 8ITMP0 = ITMP0 + 1}

ITMP1 <--- ITMP1ITMP1 = ITMP1 + DESCV( NB_ ){2ITMP1 = ITMP1 - DESCV( NB_ ) - 1, 3ITMP1 = ITMP1 + DESCV( MB_ ), 4ITMP1 = ITMP1 - DESCV( MB_ ) - 1}, KITMP1 = K + 1 + (K-1) * DESCV( NB_ ){2ITMP1 = K + 1 + (K-1) * DESCV( MB_ )}, MB_ITMP1 = ITMP1 + DESCV( MB_ ){2ITMP1 = K + 1 + (K-1) * DESCV( MB_ ), 3ITMP1 = ITMP1 - DESCV( MB_ ) - 1}, NB_ITMP1 = ITMP1 + DESCV( NB_ ){2ITMP1 = K + 1 + (K-1) * DESCV( NB_ ), 3ITMP1 = ITMP1 - DESCV( NB_ ) - 1}

IW <--- ITMP0IW = IW + ITMP0{2IW = IW + ITMP0, 3IW = IW + ITMP0, 4IW = IW + ITMP0, 5IW = IW + ITMP0, 6IW = IW + ITMP0, 7IW = IW + ITMP0, 8IW = IW + ITMP0}, IWIW = IW + ITMP0{2IW = IW + ITMP0, 3IW = IW + ITMP0, 4IW = IW + ITMP0, 5IW = IW + ITMP0, 6IW = IW + ITMP0, 7IW = IW + ITMP0, 8IW = IW + ITMP0}

JJ <--- JJJJ = JJ + 1{2JJ = JJ - 1}, JJVDO 20 JJ = JJV+1, JJV+K-1{2DO 30 JJ = JJV+K-2, JJV, -1, 3DO 40 JJ = JJV+K-2, JJV, -1, 4JJ = JJV + 1, 5JJ = JJV, 6JJ = JJV + NQ - 1, 7DO 10 JJ = JJV+1, JJV+K-1}, KDO 20 JJ = JJV+1, JJV+K-1{2DO 30 JJ = JJV+K-2, JJV, -1, 3DO 40 JJ = JJV+K-2, JJV, -1, 4DO 10 JJ = JJV+1, JJV+K-1}, NQJJ = JJV + NQ - 1

LDV <--- LLD_LDV = DESCV( LLD_ ){2LDV = DESCV( LLD_ )}

MICOL <--- INDXG2P, IVCOLMICOL = MOD( IVCOL+1, NPCOL ){2MICOL = IVCOL}, CSRC_, JV, MYCOL, N, NB_, NPCOLMICOL = MOD( IVCOL+1, NPCOL ){2}

MIROW <--- INDXG2P, IV, IVROWMIROW = MOD( IVROW+1, NPROW ){2MIROW = IVROW}, MB_, MYROW, N, NPROW{2MIROW = MOD( IVROW+1, NPROW )}, RSRC_

NP <--- IROFF{2, 3NP = NP - IROFF}, IVROW{2}, MB_{2}, MYROW{2}, N{2}, NPNP = NP - IROFF, NPROW{2}, NUMROC{2}

NQ <--- IVCOL{2}, MYCOL{2}, N{2}, NB_{2}, NPCOL{2}, NQNQ = NQ - ICOFF, NUMROC{2}, ICOFF{2NQ = NQ - ICOFF, 3}

T <--- IIT( ITMP1+ITMP0 ) = TAU( II ){2T( ITMP1-1 ) = TAU( II )}, IIVT( ITMP1 ) = TAU( IIV ){2T( ITMP1-1 ) = TAU( IIV+K-1 )}, JJT(ITMP1+ITMP0) = TAU( JJ ){2T( ITMP1-1 ) = TAU( JJ )}, JJVT( ITMP1 ) = TAU( JJV ){2T( ITMP1-1 ) = TAU( JJV+K-1 )}, KT( ITMP1-1 ) = TAU( JJV+K-1 ){2T( ITMP1-1 ) = TAU( IIV+K-1 )}

V <--- ONEV( II+(JJ-1)*LDV ) = ONE{2V( II+(JJ-1)*LDV ) = ONE, 3V( II+(JJ-1)*LDV ) = ONE, 4V( II+(JJ-1)*LDV ) = ONE}, VIIV( II+(JJ-1)*LDV ) = VII{2V( II+(JJ-1)*LDV ) = VII, 3V( II+(JJ-1)*LDV ) = VII, 4V( II+(JJ-1)*LDV ) = VII}

VII <--- IIVII = V( II+(JJ-1)*LDV ){2VII = V( II+(JJ-1)*LDV ), 3VII = V( II+(JJ-1)*LDV ), 4VII = V( II+(JJ-1)*LDV )}, JJVII = V( II+(JJ-1)*LDV ){2VII = V( II+(JJ-1)*LDV ), 3VII = V( II+(JJ-1)*LDV ), 4VII = V( II+(JJ-1)*LDV )}, LDVVII = V( II+(JJ-1)*LDV ){2VII = V( II+(JJ-1)*LDV ), 3VII = V( II+(JJ-1)*LDV ), 4VII = V( II+(JJ-1)*LDV )}, VVII = V( II+(JJ-1)*LDV ){2VII = V( II+(JJ-1)*LDV ), 3VII = V( II+(JJ-1)*LDV ), 4VII = V( II+(JJ-1)*LDV )}

Analysis elements of the routine PCLARFT()
Put the mouse over each element to display detailed matching information

Assigned variables
		BLOCK_CYCLIC_2D , CSRC_ , CTXT_ , DIRECT , DLEN_ , DTYPE_ , FORWARD , ICOFF , ICTXT , II , IROFF , ITMP0 , ITMP1 , IW , JJ , LDV , LLD_ , M_ , MB_ , MICOL , MIROW , N_ , NB_ , NP , NQ , ONE , RSRC_ , STOREV , VII , ZERO

Active variables
		BLOCK_CYCLIC_2D , CSRC_ , CTXT_ , DESCV , DIRECT , DLEN_ , DTYPE_ , FORWARD , ICOFF , ICTXT , II , IIV , INDXG2P , IROFF , ITMP0 , ITMP1 , IV , IVCOL , IVROW , IW , JJ , JJV , JV , K , LDV , LLD_ , LSAME , M_ , MB_ , MICOL , MIROW , MYCOL , MYROW , N , N_ , NB_ , NP , NPCOL , NPROW , NQ , NUMROC , ONE , RSRC_ , STOREV , T , TAU , V , VII , WORK , ZERO

Accessed arrays [ array name : associated index ]
	DESCV	: CSRC_ , CTXT_ , LLD_ , LLD_ , MB_ , MB_ , MB_ , MB_ , MB_ , MB_ , MB_ , MB_ , MB_ , MB_ , MB_ , MB_ , MB_ , NB_ , NB_ , NB_ , NB_ , NB_ , NB_ , NB_ , NB_ , NB_ , NB_ , NB_ , NB_ , NB_ , RSRC_
	LSAME	: DIRECT, 'F' , STOREV, 'C' , STOREV, 'R'
	NUMROC	: N+ICOFF, DESCV( NB_ ), MYCOL, IVCOL, NPCOL , N+ICOFF-1, DESCV( NB_ ), MYCOL, IVCOL, NPCOL , N+IROFF, DESCV( MB_ ), MYROW, IVROW, NPROW , N+IROFF-1, DESCV( MB_ ), MYROW, IVROW, NPROW
	T	: 1:i-1,1:i-1 , 1:i-1,i , 1:i-1,i , 1:i-1,i , 1:i-1,i , 1:j-1,1:j-1 , 1:j-1,j , i+1:k,i , i+1:k,i , i+1:k,i+1:k , ITMP1 , ITMP1 , ITMP1 , ITMP1 , ITMP1 , ITMP1 , ITMP1 , ITMP1 , ITMP1 , ITMP1 , ITMP1+DESCV( MB_ ) , ITMP1+DESCV( NB_ ) , ITMP1+ITMP0 , ITMP1+ITMP0 , ITMP1-1 , ITMP1-1 , ITMP1-1 , ITMP1-1 , j+1:k,j , j+1:k,j+1:k
	TAU	: II , II , II , II , IIV , IIV+K-1 , iv+i-1 , iv+i-1 , JJ , JJ , JJ , JJ , JJV , JJV+K-1 , jv+i-1 , jv+i-1
	V	: II+(JJ-1)LDV , II+(JJ-1)LDV , II+(JJ-1)LDV , II+(JJ-1)LDV , II+(JJ-1)LDV , II+(JJ-1)LDV , II+(JJ-1)LDV , II+(JJ-1)LDV , II+(JJ-1)LDV , II+(JJ-1)LDV , II+(JJ-1)LDV , II+(JJ-1)LDV , II+(JJ-1)LDV , II+(JJ-1)LDV , II+(JJ-1)LDV , II+(JJ-1)LDV , II+(JJV-1)LDV , II+(JJV-1)LDV , II+(JJV-1)LDV , II+(JJV-1)LDV , II+1+(JJV-1)LDV , IIV+(JJ-1)LDV , IIV+(JJ-1)LDV , IIV+JJLDV , iv:iv+i-2,jv+i-1:jv+n-1 , iv:iv+n-k+i-1,jv+i:jv+k-1 , iv:iv+n-k+i-1,jv+i-1 , iv+i:iv+k-1,jv:jv+n-k+i-1 , iv+i-1,jv:jv+n-k+i-1 , iv+i-1,jv+i-1:jv+n-1 , iv+i-1:iv+n-1,jv:jv+i-2 , iv+i-1:iv+n-1,jv+i-1
	WORK	: IW , IW , IW , IW , IW , IW , IW , IW , IW , IW , IW , IW

Conditional statements [ statement : associated predicate ]
	do	: `(` 10 JJ = JJV + 1 , JJV + K - 1 `)` , `(` 20 JJ = JJV + 1 , JJV + K - 1 `)` , `(` 30 JJ = JJV + K - 2 , JJV , - 1 `)` , `(` 40 JJ = JJV + K - 2 , JJV , - 1 `)` , `(` 50 II = IIV + 1 , IIV + K - 1 `)` , `(` 60 II = IIV + 1 , IIV + K - 1 `)` , `(` 70 II = IIV + K - 2 , IIV , - 1 `)` , `(` 80 II = IIV + K - 2 , IIV , - 1 `)`
	if	: `(` possible `)` , `(` N.LE.0 .OR. K.LE.0 `)` , `(` (LSAME( STOREV , 'C' ) .AND. MYCOL.EQ.IVCOL ) `)` , `(` FORWARD `)` , `(` MYROW.EQ.IVROW `)` , `(` (IROFF + 1.EQ.DESCV( MB_ ) ) `)` , `(` MYROW.EQ.MIROW `)` , `(` NP-II+IIV.GT.0 `)` , `(` MYROW.EQ.MIROW `)` , `(` (MOD( IV + ITMP0 , DESCV( MB_ ) ).EQ.0 ) `)` , `(` MYROW.EQ.IVROW `)` , `(` MYROW.EQ.IVROW `)` , `(` MYROW.EQ.MIROW `)` , `(` II-IIV + 1.GT.0 `)` , `(` MYROW.EQ.MIROW `)` , `(` (MOD( IV + N - ITMP0 - 2 , DESCV(MB_) ).EQ.0 ) `)` , `(` MYROW.EQ.IVROW `)` , `(` (LSAME( STOREV , 'R' ) .AND. MYROW.EQ.IVROW ) `)` , `(` FORWARD `)` , `(` MYCOL.EQ.IVCOL `)` , `(` (ICOFF + 1.EQ.DESCV( NB_ ) ) `)` , `(` MYCOL.EQ.MICOL `)` , `(` NQ-JJ+JJV.GT.0 `)` , `(` MYCOL.EQ.MICOL `)` , `(` (MOD( JV + ITMP0 , DESCV( NB_ ) ).EQ.0 ) `)` , `(` MYCOL.EQ.IVCOL `)` , `(` MYCOL.EQ.IVCOL `)` , `(` MYCOL.EQ.MICOL `)` , `(` JJ-JJV + 1.GT.0 `)` , `(` MYCOL.EQ.MICOL `)` , `(` (MOD( JV + N - ITMP0 - 2 , DESCV( NB_ ) ).EQ.0 ) `)` , `(` MYCOL.EQ.IVCOL `)`