|
|
| |
| # lines: |
331 | | # code: |
331 | | # comment: | 0 | |
# blank: | 0 |
| # Variables: | 54 |
| # Callers: | 1 |
| # Callings: | 0 |
| # Words: | 121 |
| # Keywords: | 56 |
|
|
|
|
|
..
.. Array Arguments ..
..
Purpose
=======
ZLAREF applies one or several Householder reflectors of size 3
to one or two matrices (if column is specified) on either their
rows or columns.
Arguments
=========
TYPE (global input) CHARACTER*1
If 'R': Apply reflectors to the rows of the matrix
(apply from left)
Otherwise: Apply reflectors to the columns of the matrix
Unchanged on exit.
A (global input/output) COMPLEX*16 array, (LDA,*)
On entry, the matrix to receive the reflections.
The updated matrix on exit.
LDA (local input) INTEGER
On entry, the leading dimension of A. Unchanged on exit.
WANTZ (global input) LOGICAL
If .TRUE., then apply any column reflections to Z as well.
If .FALSE., then do no additional work on Z.
Z (global input/output) COMPLEX*16 array, (LDZ,*)
On entry, the second matrix to receive column reflections.
This is changed only if WANTZ is set.
LDZ (local input) INTEGER
On entry, the leading dimension of Z. Unchanged on exit.
BLOCK (global input) LOGICAL
If .TRUE., then apply several reflectors at once and read
their data from the VECS array.
If .FALSE., apply the single reflector given by V2, V3,
T1, T2, and T3.
IROW1 (local input/output) INTEGER
On entry, the local row element of A.
Undefined on output.
ICOL1 (local input/output) INTEGER
On entry, the local column element of A.
Undefined on output.
ISTART (global input) INTEGER
Specifies the "number" of the first reflector. This is
used as an index into VECS if BLOCK is set.
ISTART is ignored if BLOCK is .FALSE..
ISTOP (global input) INTEGER
Specifies the "number" of the last reflector. This is
used as an index into VECS if BLOCK is set.
ISTOP is ignored if BLOCK is .FALSE..
ITMP1 (local input) INTEGER
Starting range into A. For rows, this is the local
first column. For columns, this is the local first row.
ITMP2 (local input) INTEGER
Ending range into A. For rows, this is the local last
column. For columns, this is the local last row.
LILOZ
LIHIZ (local input) INTEGER
These serve the same purpose as ITMP1,ITMP2 but for Z
when WANTZ is set.
VECS (global input) COMPLEX*16 array of size 3*N (matrix size)
This holds the size 3 reflectors one after another and this
is only accessed when BLOCK is .TRUE.
V2
V3
T1
T2
T3 (global input/output) COMPLEX*16
This holds information on a single size 3 Householder
reflector and is read when BLOCK is .FALSE., and
overwritten when BLOCK is .TRUE.
Further Details
===============
Implemented by: M. Fahey, May 28, 1999
=====================================================================
.. Local Scalars ..
|
|
|
|
001 SUBROUTINE ZLAREF( TYPE , A , LDA , WANTZ , Z , LDZ , BLOCK , IROW1 ,
002 $ICOL1 , ISTART , ISTOP , ITMP1 , ITMP2 , LILOZ ,
003 $LIHIZ , VECS , V2 , V3 , T1 , T2 , T3 )
004
005 * -- ScaLAPACK routine(version 1.7) --
006 * Univ. of Tennessee , Univ. of California Berkeley , NAG Ltd. ,
007 * Courant Institute , Argonne National Lab , and Rice University
008 * May 28 , 1999
009
010 * .. Scalar Arguments ..
011 LOGICAL BLOCK , WANTZ
012 CHARACTER TYPE
013 INTEGER ICOL1 , IROW1 , ISTART , ISTOP , ITMP1 , ITMP2 , LDA ,
014 $LDZ , LIHIZ , LILOZ
015 COMPLEX*16 T1 , T2 , T3 , V2 , V3
016 INTEGER J , K
017 COMPLEX*16 A1 , A11 , A2 , A22 , A3 , A4 , A5 , B1 , B2 , B3 , B4 ,
018 $B5 , H11 , H22 , SUM , SUM1 , SUM2 , SUM3 , T12 , T13 ,
019 $T22 , T23 , T32 , T33 , TMP1 , TMP2 , TMP3 , V22 , V23 ,
020 $V32 , V33
021 * ..
022 * .. External Functions ..
023 LOGICAL LSAME
024 EXTERNAL LSAME
025 * ..
026 * .. Intrinsic Functions ..
027 INTRINSIC DCONJG , MOD
028 * ..
029 * .. Executable Statements ..
030
031 IF( LSAME( TYPE , 'R' ) ) THEN
031  
032 IF( BLOCK ) THEN
032
033 DO 30 K = ISTART , ISTOP - MOD( ISTOP - ISTART + 1 , 3 ) , 3
033
034 V2 = VECS(( K - 1 )*3 + 1 )
035 V3 = VECS(( K - 1 )*3 + 2 )
036 T1 = VECS(( K - 1 )*3 + 3 )
037 V22 = VECS(( K - 1 )*3 + 4 )
038 V32 = VECS(( K - 1 )*3 + 5 )
039 T12 = VECS(( K - 1 )*3 + 6 )
040 V23 = VECS(( K - 1 )*3 + 7 )
041 V33 = VECS(( K - 1 )*3 + 8 )
042 T13 = VECS(( K - 1 )*3 + 9 )
043 T2 = T1*V2
044 T3 = T1*V3
045 T22 = T12*V22
046 T32 = T12*V32
047 T23 = T13*V23
048 T33 = T13*V33
049 DO 10 J = ITMP1 , ITMP2 - MOD( ITMP2 - ITMP1 + 1 , 2 ) , 2
049
050 A1 = A( IROW1 , J )
051 A2 = A( IROW1 + 1 , J )
052 A3 = A( IROW1 + 2 , J )
053 A4 = A( IROW1 + 3 , J )
054 A5 = A( IROW1 + 4 , J )
055 B1 = A( IROW1 , J + 1 )
056 B2 = A( IROW1 + 1 , J + 1 )
057 B3 = A( IROW1 + 2 , J + 1 )
058 B4 = A( IROW1 + 3 , J + 1 )
059 B5 = A( IROW1 + 4 , J + 1 )
060 SUM1 = DCONJG( T1 )*A1 + DCONJG( T2 )*A2 +
061 $ DCONJG( T3 )*A3
062 A( IROW1 , J ) = A1 - SUM1
063 H11 = A2 - SUM1*V2
064 H22 = A3 - SUM1*V3
065 TMP1 = DCONJG( T1 )*B1 + DCONJG( T2 )*B2 +
066 $ DCONJG( T3 )*B3
067 A( IROW1 , J + 1 ) = B1 - TMP1
068 A11 = B2 - TMP1*V2
069 A22 = B3 - TMP1*V3
070 SUM2 = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 +
071 $ DCONJG( T32 )*A4
072 A( IROW1 + 1 , J ) = H11 - SUM2
073 H11 = H22 - SUM2*V22
074 H22 = A4 - SUM2*V32
075 TMP2 = DCONJG( T12 )*A11 + DCONJG( T22 )*A22 +
076 $ DCONJG( T32 )*B4
077 A( IROW1 + 1 , J + 1 ) = A11 - TMP2
078 A11 = A22 - TMP2*V22
079 A22 = B4 - TMP2*V32
080 SUM3 = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 +
081 $ DCONJG( T33 )*A5
082 A( IROW1 + 2 , J ) = H11 - SUM3
083 A( IROW1 + 3 , J ) = H22 - SUM3*V23
084 A( IROW1 + 4 , J ) = A5 - SUM3*V33
085 TMP3 = DCONJG( T13 )*A11 + DCONJG( T23 )*A22 +
086 $ DCONJG( T33 )*B5
087 A( IROW1 + 2 , J + 1 ) = A11 - TMP3
088 A( IROW1 + 3 , J + 1 ) = A22 - TMP3*V23
089 A( IROW1 + 4 , J + 1 ) = B5 - TMP3*V33
090 10 CONTINUE
090
091 DO 20 J = ITMP2 - MOD( ITMP2 - ITMP1 + 1 , 2 ) + 1 , ITMP2
091
092 SUM = DCONJG( T1 )*A( IROW1 , J ) +
093 $ DCONJG( T2 )*A( IROW1 + 1 , J ) +
094 $ DCONJG( T3 )*A( IROW1 + 2 , J )
095 A( IROW1 , J ) = A( IROW1 , J ) - SUM
096 H11 = A( IROW1 + 1 , J ) - SUM*V2
097 H22 = A( IROW1 + 2 , J ) - SUM*V3
098 SUM = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 +
099 $ DCONJG( T32 )*A( IROW1 + 3 , J )
100 A( IROW1 + 1 , J ) = H11 - SUM
101 H11 = H22 - SUM*V22
102 H22 = A( IROW1 + 3 , J ) - SUM*V32
103 SUM = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 +
104 $ DCONJG( T33 )*A( IROW1 + 4 , J )
105 A( IROW1 + 2 , J ) = H11 - SUM
106 A( IROW1 + 3 , J ) = H22 - SUM*V23
107 A( IROW1 + 4 , J ) = A( IROW1 + 4 , J ) - SUM*V33
108 20 CONTINUE
108
109 IROW1 = IROW1 + 3
110 30 CONTINUE
110
111 DO 50 K = ISTOP - MOD( ISTOP - ISTART + 1 , 3 ) + 1 , ISTOP
111
112 V2 = VECS(( K - 1 )*3 + 1 )
113 V3 = VECS(( K - 1 )*3 + 2 )
114 T1 = VECS(( K - 1 )*3 + 3 )
115 T2 = T1*V2
116 T3 = T1*V3
117 DO 40 J = ITMP1 , ITMP2
117
118 SUM = DCONJG( T1 )*A( IROW1 , J ) +
119 $ DCONJG( T2 )*A( IROW1 + 1 , J ) +
120 $ DCONJG( T3 )*A( IROW1 + 2 , J )
121 A( IROW1 , J ) = A( IROW1 , J ) - SUM
122 A( IROW1 + 1 , J ) = A( IROW1 + 1 , J ) - SUM*V2
123 A( IROW1 + 2 , J ) = A( IROW1 + 2 , J ) - SUM*V3
124 40 CONTINUE
124
125 IROW1 = IROW1 + 1
126 50 CONTINUE
126
127 ELSE
127
128 DO 60 J = ITMP1 , ITMP2
128
129 SUM = DCONJG( T1 )*A( IROW1 , J ) +
130 $ DCONJG( T2 )*A( IROW1 + 1 , J ) +
131 $ DCONJG( T3 )*A( IROW1 + 2 , J )
132 A( IROW1 , J ) = A( IROW1 , J ) - SUM
133 A( IROW1 + 1 , J ) = A( IROW1 + 1 , J ) - SUM*V2
134 A( IROW1 + 2 , J ) = A( IROW1 + 2 , J ) - SUM*V3
135 60 CONTINUE
135
136 END IF
137 ELSE
138
139 * Do column transforms
140
140
141 IF( BLOCK ) THEN
141
142 DO 90 K = ISTART , ISTOP - MOD( ISTOP - ISTART + 1 , 3 ) , 3
142
143 V2 = VECS(( K - 1 )*3 + 1 )
144 V3 = VECS(( K - 1 )*3 + 2 )
145 T1 = VECS(( K - 1 )*3 + 3 )
146 V22 = VECS(( K - 1 )*3 + 4 )
147 V32 = VECS(( K - 1 )*3 + 5 )
148 T12 = VECS(( K - 1 )*3 + 6 )
149 V23 = VECS(( K - 1 )*3 + 7 )
150 V33 = VECS(( K - 1 )*3 + 8 )
151 T13 = VECS(( K - 1 )*3 + 9 )
152 T2 = T1*V2
153 T3 = T1*V3
154 T22 = T12*V22
155 T32 = T12*V32
156 T23 = T13*V23
157 T33 = T13*V33
158 DO 70 J = ITMP1 , ITMP2
158
159 SUM = T1*A( J , ICOL1 ) + T2*A( J , ICOL1 + 1 ) +
160 $ T3*A( J , ICOL1 + 2 )
161 A( J , ICOL1 ) = A( J , ICOL1 ) - SUM
162 H11 = A( J , ICOL1 + 1 ) - SUM*DCONJG( V2 )
163 H22 = A( J , ICOL1 + 2 ) - SUM*DCONJG( V3 )
164 SUM = T12*H11 + T22*H22 + T32*A( J , ICOL1 + 3 )
165 A( J , ICOL1 + 1 ) = H11 - SUM
166 H11 = H22 - SUM*DCONJG( V22 )
167 H22 = A( J , ICOL1 + 3 ) - SUM*DCONJG( V32 )
168 SUM = T13*H11 + T23*H22 + T33*A( J , ICOL1 + 4 )
169 A( J , ICOL1 + 2 ) = H11 - SUM
170 A( J , ICOL1 + 3 ) = H22 - SUM*DCONJG( V23 )
171 A( J , ICOL1 + 4 ) = A( J , ICOL1 + 4 ) - SUM*DCONJG( V33 )
172 70 CONTINUE
172
173 IF( WANTZ ) THEN
173
174 DO 80 J = LILOZ , LIHIZ
174
175 SUM = T1*Z( J , ICOL1 ) + T2*Z( J , ICOL1 + 1 ) +
176 $ T3*Z( J , ICOL1 + 2 )
177 Z( J , ICOL1 ) = Z( J , ICOL1 ) - SUM
178 H11 = Z( J , ICOL1 + 1 ) - SUM*DCONJG( V2 )
179 H22 = Z( J , ICOL1 + 2 ) - SUM*DCONJG( V3 )
180 SUM = T12*H11 + T22*H22 + T32*Z( J , ICOL1 + 3 )
181 Z( J , ICOL1 + 1 ) = H11 - SUM
182 H11 = H22 - SUM*DCONJG( V22 )
183 H22 = Z( J , ICOL1 + 3 ) - SUM*DCONJG( V32 )
184 SUM = T13*H11 + T23*H22 + T33*Z( J , ICOL1 + 4 )
185 Z( J , ICOL1 + 2 ) = H11 - SUM
186 Z( J , ICOL1 + 3 ) = H22 - SUM*DCONJG( V23 )
187 Z( J , ICOL1 + 4 ) = Z( J , ICOL1 + 4 ) -
188 $ SUM*DCONJG( V33 )
189 80 CONTINUE
189
190 END IF
191 ICOL1 = ICOL1 + 3
192 90 CONTINUE
192
193 DO 120 K = ISTOP - MOD( ISTOP - ISTART + 1 , 3 ) + 1 , ISTOP
193
194 V2 = VECS(( K - 1 )*3 + 1 )
195 V3 = VECS(( K - 1 )*3 + 2 )
196 T1 = VECS(( K - 1 )*3 + 3 )
197 T2 = T1*V2
198 T3 = T1*V3
199 DO 100 J = ITMP1 , ITMP2
199
200 SUM = T1*A( J , ICOL1 ) + T2*A( J , ICOL1 + 1 ) +
201 $ T3*A( J , ICOL1 + 2 )
202 A( J , ICOL1 ) = A( J , ICOL1 ) - SUM
203 A( J , ICOL1 + 1 ) = A( J , ICOL1 + 1 ) - SUM*DCONJG( V2 )
204 A( J , ICOL1 + 2 ) = A( J , ICOL1 + 2 ) - SUM*DCONJG( V3 )
205 100 CONTINUE
205
206 IF( WANTZ ) THEN
206
207 DO 110 J = LILOZ , LIHIZ
207
208 SUM = T1*Z( J , ICOL1 ) + T2*Z( J , ICOL1 + 1 ) +
209 $ T3*Z( J , ICOL1 + 2 )
210 Z( J , ICOL1 ) = Z( J , ICOL1 ) - SUM
211 Z( J , ICOL1 + 1 ) = Z( J , ICOL1 + 1 ) -
212 $ SUM*DCONJG( V2 )
213 Z( J , ICOL1 + 2 ) = Z( J , ICOL1 + 2 ) -
214 $ SUM*DCONJG( V3 )
215 110 CONTINUE
215
216 END IF
217 ICOL1 = ICOL1 + 1
218 120 CONTINUE
218
219 ELSE
219
220 DO 130 J = ITMP1 , ITMP2
220
221 SUM = T1*A( J , ICOL1 ) + T2*A( J , ICOL1 + 1 ) +
222 $ T3*A( J , ICOL1 + 2 )
223 A( J , ICOL1 ) = A( J , ICOL1 ) - SUM
224 A( J , ICOL1 + 1 ) = A( J , ICOL1 + 1 ) - SUM*DCONJG( V2 )
225 A( J , ICOL1 + 2 ) = A( J , ICOL1 + 2 ) - SUM*DCONJG( V3 )
226 130 CONTINUE
226
227 END IF
228 END IF
229 RETURN
230
231 * End of ZLAREF
232
233 END7
34
|
|
Variables in Routine ZLAREF()
| Summary Report |
| Data Type | Quantity | Size(byte) |
| CHARACTER | 1 | 1 |
| COMPLEX*16 | 36 | ? |
| INTEGER | 12 | 48 |
| LOGICAL | 3 | 3 |
| REAL | 2 | 8 |
| TOTAL | 54 | 60 |
List of Variables
CHARACTER
COMPLEX*16
| A1 | A11 | A2 | A22 | A3 |
| A4 | A5 | B1 | B2 | B3 |
| B4 | B5 | H11 | H22 | SUM |
| SUM1 | SUM2 | SUM3 | T1 | T12 |
| T13 | T2 | T22 | T23 | T3 |
| T32 | T33 | TMP1 | TMP2 | TMP3 |
| V2 | V22 | V23 | V3 | V32 |
| V33 | | | | |
INTEGER
| ICOL1 | IROW1 | ISTART | ISTOP | ITMP1 |
| ITMP2 | J | K | LDA | LDZ |
| LIHIZ | LILOZ | | | |
LOGICAL
REAL
Variables Dependence Graph
Put the mouse over a right hand side variable to display the corresponding line of the dependence | | - | | - | - | | A | <--- | AA( IROW1+4, J ) = A( IROW1+4, J ) - SUM*V33{2A( IROW1, J ) = A( IROW1, J ) - SUM, 3A( IROW1+1, J ) = A( IROW1+1, J ) - SUM*V2, 4A( IROW1+2, J ) = A( IROW1+2, J ) - SUM*V3, 5A( IROW1, J ) = A( IROW1, J ) - SUM, 6A( IROW1+1, J ) = A( IROW1+1, J ) - SUM*V2, 7A( IROW1+2, J ) = A( IROW1+2, J ) - SUM*V3, 8A( J, ICOL1 ) = A( J, ICOL1 ) - SUM, 9A( J, ICOL1+4 ) = A( J, ICOL1+4 ) - SUM*DCONJG( V33 ), 10A( J, ICOL1 ) = A( J, ICOL1 ) - SUM, 11A( J, ICOL1+1 ) = A( J, ICOL1+1 ) - SUM*DCONJG( V2 ), 12A( J, ICOL1+2 ) = A( J, ICOL1+2 ) - SUM*DCONJG( V3 ), 13A( J, ICOL1 ) = A( J, ICOL1 ) - SUM, 14A( J, ICOL1+1 ) = A( J, ICOL1+1 ) - SUM*DCONJG( V2 ), 15A( J, ICOL1+2 ) = A( J, ICOL1+2 ) - SUM*DCONJG( V3 ), 16A( IROW1, J ) = A( IROW1, J ) - SUM}, B5A( IROW1+4, J+1 ) = B5 - TMP3*V33, H11A( IROW1+1, J ) = H11 - SUM{2A( IROW1+2, J ) = H11 - SUM, 3A( J, ICOL1+1 ) = H11 - SUM, 4A( J, ICOL1+2 ) = H11 - SUM, 5A( IROW1+1, J ) = H11 - SUM2, 6A( IROW1+2, J ) = H11 - SUM3}, H22A( IROW1+3, J ) = H22 - SUM*V23{2A( J, ICOL1+3 ) = H22 - SUM*DCONJG( V23 ), 3A( IROW1+3, J ) = H22 - SUM3*V23}, ICOL1A( J, ICOL1 ) = A( J, ICOL1 ) - SUM{2A( J, ICOL1+4 ) = A( J, ICOL1+4 ) - SUM*DCONJG( V33 ), 3A( J, ICOL1 ) = A( J, ICOL1 ) - SUM, 4A( J, ICOL1+1 ) = A( J, ICOL1+1 ) - SUM*DCONJG( V2 ), 5A( J, ICOL1+2 ) = A( J, ICOL1+2 ) - SUM*DCONJG( V3 ), 6A( J, ICOL1 ) = A( J, ICOL1 ) - SUM, 7A( J, ICOL1+1 ) = A( J, ICOL1+1 ) - SUM*DCONJG( V2 ), 8A( J, ICOL1+2 ) = A( J, ICOL1+2 ) - SUM*DCONJG( V3 )}, IROW1A( IROW1+4, J ) = A( IROW1+4, J ) - SUM*V33{2A( IROW1, J ) = A( IROW1, J ) - SUM, 3A( IROW1+1, J ) = A( IROW1+1, J ) - SUM*V2, 4A( IROW1+2, J ) = A( IROW1+2, J ) - SUM*V3, 5A( IROW1, J ) = A( IROW1, J ) - SUM, 6A( IROW1+1, J ) = A( IROW1+1, J ) - SUM*V2, 7A( IROW1+2, J ) = A( IROW1+2, J ) - SUM*V3, 8A( IROW1, J ) = A( IROW1, J ) - SUM}, A1A( IROW1, J ) = A1 - SUM1, JA( IROW1+4, J ) = A( IROW1+4, J ) - SUM*V33{2A( IROW1, J ) = A( IROW1, J ) - SUM, 3A( IROW1+1, J ) = A( IROW1+1, J ) - SUM*V2, 4A( IROW1+2, J ) = A( IROW1+2, J ) - SUM*V3, 5A( IROW1, J ) = A( IROW1, J ) - SUM, 6A( IROW1+1, J ) = A( IROW1+1, J ) - SUM*V2, 7A( IROW1+2, J ) = A( IROW1+2, J ) - SUM*V3, 8A( J, ICOL1 ) = A( J, ICOL1 ) - SUM, 9A( J, ICOL1+4 ) = A( J, ICOL1+4 ) - SUM*DCONJG( V33 ), 10A( J, ICOL1 ) = A( J, ICOL1 ) - SUM, 11A( J, ICOL1+1 ) = A( J, ICOL1+1 ) - SUM*DCONJG( V2 ), 12A( J, ICOL1+2 ) = A( J, ICOL1+2 ) - SUM*DCONJG( V3 ), 13A( J, ICOL1 ) = A( J, ICOL1 ) - SUM, 14A( J, ICOL1+1 ) = A( J, ICOL1+1 ) - SUM*DCONJG( V2 ), 15A( J, ICOL1+2 ) = A( J, ICOL1+2 ) - SUM*DCONJG( V3 ), 16A( IROW1, J ) = A( IROW1, J ) - SUM}, A11A( IROW1+1, J+1 ) = A11 - TMP2{2A( IROW1+2, J+1 ) = A11 - TMP3}, SUMA( IROW1+1, J ) = H11 - SUM{2A( IROW1+2, J ) = H11 - SUM, 3A( IROW1+3, J ) = H22 - SUM*V23, 4A( IROW1+4, J ) = A( IROW1+4, J ) - SUM*V33, 5A( IROW1, J ) = A( IROW1, J ) - SUM, 6A( IROW1+1, J ) = A( IROW1+1, J ) - SUM*V2, 7A( IROW1+2, J ) = A( IROW1+2, J ) - SUM*V3, 8A( IROW1, J ) = A( IROW1, J ) - SUM, 9A( IROW1+1, J ) = A( IROW1+1, J ) - SUM*V2, 10A( IROW1+2, J ) = A( IROW1+2, J ) - SUM*V3, 11A( J, ICOL1 ) = A( J, ICOL1 ) - SUM, 12A( J, ICOL1+1 ) = H11 - SUM, 13A( J, ICOL1+2 ) = H11 - SUM, 14A( J, ICOL1+3 ) = H22 - SUM*DCONJG( V23 ), 15A( J, ICOL1+4 ) = A( J, ICOL1+4 ) - SUM*DCONJG( V33 ), 16A( J, ICOL1 ) = A( J, ICOL1 ) - SUM, 17A( J, ICOL1+1 ) = A( J, ICOL1+1 ) - SUM*DCONJG( V2 ), 18A( J, ICOL1+2 ) = A( J, ICOL1+2 ) - SUM*DCONJG( V3 ), 19A( J, ICOL1 ) = A( J, ICOL1 ) - SUM, 20A( J, ICOL1+1 ) = A( J, ICOL1+1 ) - SUM*DCONJG( V2 ), 21A( J, ICOL1+2 ) = A( J, ICOL1+2 ) - SUM*DCONJG( V3 ), 22A( IROW1, J ) = A( IROW1, J ) - SUM}, SUM1A( IROW1, J ) = A1 - SUM1, SUM2A( IROW1+1, J ) = H11 - SUM2, SUM3A( IROW1+2, J ) = H11 - SUM3{2A( IROW1+3, J ) = H22 - SUM3*V23, 3A( IROW1+4, J ) = A5 - SUM3*V33}, TMP1A( IROW1, J+1 ) = B1 - TMP1, TMP2A( IROW1+1, J+1 ) = A11 - TMP2, TMP3A( IROW1+2, J+1 ) = A11 - TMP3{2A( IROW1+3, J+1 ) = A22 - TMP3*V23, 3A( IROW1+4, J+1 ) = B5 - TMP3*V33}, V2A( IROW1+1, J ) = A( IROW1+1, J ) - SUM*V2{2A( IROW1+1, J ) = A( IROW1+1, J ) - SUM*V2, 3A( J, ICOL1+1 ) = A( J, ICOL1+1 ) - SUM*DCONJG( V2 ), 4A( J, ICOL1+1 ) = A( J, ICOL1+1 ) - SUM*DCONJG( V2 )}, V23A( IROW1+3, J ) = H22 - SUM*V23{2A( J, ICOL1+3 ) = H22 - SUM*DCONJG( V23 ), 3A( IROW1+3, J ) = H22 - SUM3*V23, 4A( IROW1+3, J+1 ) = A22 - TMP3*V23}, A22A( IROW1+3, J+1 ) = A22 - TMP3*V23, V3A( IROW1+2, J ) = A( IROW1+2, J ) - SUM*V3{2A( IROW1+2, J ) = A( IROW1+2, J ) - SUM*V3, 3A( J, ICOL1+2 ) = A( J, ICOL1+2 ) - SUM*DCONJG( V3 ), 4A( J, ICOL1+2 ) = A( J, ICOL1+2 ) - SUM*DCONJG( V3 )}, V33A( IROW1+4, J ) = A( IROW1+4, J ) - SUM*V33{2A( J, ICOL1+4 ) = A( J, ICOL1+4 ) - SUM*DCONJG( V33 ), 3A( IROW1+4, J ) = A5 - SUM3*V33, 4A( IROW1+4, J+1 ) = B5 - TMP3*V33}, A5A( IROW1+4, J ) = A5 - SUM3*V33, B1A( IROW1, J+1 ) = B1 - TMP1 |
| A1 | <--- | AA1 = A( IROW1, J ), IROW1A1 = A( IROW1, J ), JA1 = A( IROW1, J ) |
| A11 | <--- | B2A11 = B2 - TMP1*V2, TMP1A11 = B2 - TMP1*V2, TMP2A11 = A22 - TMP2*V22, V2A11 = B2 - TMP1*V2, V22A11 = A22 - TMP2*V22, A22A11 = A22 - TMP2*V22 |
| A2 | <--- | AA2 = A( IROW1+1, J ), IROW1A2 = A( IROW1+1, J ), JA2 = A( IROW1+1, J ) |
| A22 | <--- | B3A22 = B3 - TMP1*V3, B4A22 = B4 - TMP2*V32, TMP1A22 = B3 - TMP1*V3, TMP2A22 = B4 - TMP2*V32, V3A22 = B3 - TMP1*V3, V32A22 = B4 - TMP2*V32 |
| A3 | <--- | AA3 = A( IROW1+2, J ), IROW1A3 = A( IROW1+2, J ), JA3 = A( IROW1+2, J ) |
| A4 | <--- | AA4 = A( IROW1+3, J ), IROW1A4 = A( IROW1+3, J ), JA4 = A( IROW1+3, J ) |
| A5 | <--- | AA5 = A( IROW1+4, J ), IROW1A5 = A( IROW1+4, J ), JA5 = A( IROW1+4, J ) |
| B1 | <--- | AB1 = A( IROW1, J+1 ), IROW1B1 = A( IROW1, J+1 ), JB1 = A( IROW1, J+1 ) |
| B2 | <--- | AB2 = A( IROW1+1, J+1 ), IROW1B2 = A( IROW1+1, J+1 ), JB2 = A( IROW1+1, J+1 ) |
| B3 | <--- | AB3 = A( IROW1+2, J+1 ), IROW1B3 = A( IROW1+2, J+1 ), JB3 = A( IROW1+2, J+1 ) |
| B4 | <--- | AB4 = A( IROW1+3, J+1 ), IROW1B4 = A( IROW1+3, J+1 ), JB4 = A( IROW1+3, J+1 ) |
| B5 | <--- | AB5 = A( IROW1+4, J+1 ), IROW1B5 = A( IROW1+4, J+1 ), JB5 = A( IROW1+4, J+1 ) |
| H11 | <--- | AH11 = A( J, ICOL1+1 ) - SUM*DCONJG( V2 ){2H11 = A( IROW1+1, J ) - SUM*V2}, H22H11 = H22 - SUM*V22{2H11 = H22 - SUM*DCONJG( V22 ), 3H11 = H22 - SUM*DCONJG( V22 ), 4H11 = H22 - SUM2*V22}, ICOL1H11 = A( J, ICOL1+1 ) - SUM*DCONJG( V2 ){2H11 = Z( J, ICOL1+1 ) - SUM*DCONJG( V2 )}, IROW1H11 = A( IROW1+1, J ) - SUM*V2, JH11 = A( J, ICOL1+1 ) - SUM*DCONJG( V2 ){2H11 = Z( J, ICOL1+1 ) - SUM*DCONJG( V2 ), 3H11 = A( IROW1+1, J ) - SUM*V2}, SUMH11 = H22 - SUM*V22{2H11 = A( J, ICOL1+1 ) - SUM*DCONJG( V2 ), 3H11 = H22 - SUM*DCONJG( V22 ), 4H11 = Z( J, ICOL1+1 ) - SUM*DCONJG( V2 ), 5H11 = H22 - SUM*DCONJG( V22 ), 6H11 = A( IROW1+1, J ) - SUM*V2}, SUM1H11 = A2 - SUM1*V2, SUM2H11 = H22 - SUM2*V22, A2H11 = A2 - SUM1*V2, V2H11 = A( J, ICOL1+1 ) - SUM*DCONJG( V2 ){2H11 = Z( J, ICOL1+1 ) - SUM*DCONJG( V2 ), 3H11 = A2 - SUM1*V2, 4H11 = A( IROW1+1, J ) - SUM*V2}, V22H11 = H22 - SUM*V22{2H11 = H22 - SUM*DCONJG( V22 ), 3H11 = H22 - SUM*DCONJG( V22 ), 4H11 = H22 - SUM2*V22}, ZH11 = Z( J, ICOL1+1 ) - SUM*DCONJG( V2 ) |
| H22 | <--- | AH22 = A( IROW1+3, J ) - SUM*V32{2H22 = A( J, ICOL1+2 ) - SUM*DCONJG( V3 ), 3H22 = A( J, ICOL1+3 ) - SUM*DCONJG( V32 ), 4H22 = A( IROW1+2, J ) - SUM*V3}, ICOL1H22 = A( J, ICOL1+2 ) - SUM*DCONJG( V3 ){2H22 = A( J, ICOL1+3 ) - SUM*DCONJG( V32 ), 3H22 = Z( J, ICOL1+2 ) - SUM*DCONJG( V3 ), 4H22 = Z( J, ICOL1+3 ) - SUM*DCONJG( V32 )}, IROW1H22 = A( IROW1+3, J ) - SUM*V32{2H22 = A( IROW1+2, J ) - SUM*V3}, JH22 = A( IROW1+3, J ) - SUM*V32{2H22 = A( J, ICOL1+2 ) - SUM*DCONJG( V3 ), 3H22 = A( J, ICOL1+3 ) - SUM*DCONJG( V32 ), 4H22 = Z( J, ICOL1+2 ) - SUM*DCONJG( V3 ), 5H22 = Z( J, ICOL1+3 ) - SUM*DCONJG( V32 ), 6H22 = A( IROW1+2, J ) - SUM*V3}, SUMH22 = A( IROW1+3, J ) - SUM*V32{2H22 = A( J, ICOL1+2 ) - SUM*DCONJG( V3 ), 3H22 = A( J, ICOL1+3 ) - SUM*DCONJG( V32 ), 4H22 = Z( J, ICOL1+2 ) - SUM*DCONJG( V3 ), 5H22 = Z( J, ICOL1+3 ) - SUM*DCONJG( V32 ), 6H22 = A( IROW1+2, J ) - SUM*V3}, SUM1H22 = A3 - SUM1*V3, SUM2H22 = A4 - SUM2*V32, V3H22 = A( J, ICOL1+2 ) - SUM*DCONJG( V3 ){2H22 = Z( J, ICOL1+2 ) - SUM*DCONJG( V3 ), 3H22 = A3 - SUM1*V3, 4H22 = A( IROW1+2, J ) - SUM*V3}, V32H22 = A( IROW1+3, J ) - SUM*V32{2H22 = A( J, ICOL1+3 ) - SUM*DCONJG( V32 ), 3H22 = Z( J, ICOL1+3 ) - SUM*DCONJG( V32 ), 4H22 = A4 - SUM2*V32}, ZH22 = Z( J, ICOL1+2 ) - SUM*DCONJG( V3 ){2H22 = Z( J, ICOL1+3 ) - SUM*DCONJG( V32 )}, A3H22 = A3 - SUM1*V3, A4H22 = A4 - SUM2*V32 |
| ICOL1 | <--- | ICOL1ICOL1 = ICOL1 + 3{2ICOL1 = ICOL1 + 1} |
| IROW1 | <--- | IROW1IROW1 = IROW1 + 3{2IROW1 = IROW1 + 1} |
| J | <--- | ITMP1DO 40 J = ITMP1, ITMP2{2DO 60 J = ITMP1, ITMP2, 3DO 70 J = ITMP1, ITMP2, 4DO 100 J = ITMP1, ITMP2, 5DO 130 J = ITMP1, ITMP2, 6DO 10 J = ITMP1, ITMP2 - MOD( ITMP2-ITMP1+1, 2 ), 2, 7DO 20 J = ITMP2 - MOD( ITMP2-ITMP1+1, 2 ) + 1, ITMP2}, ITMP2DO 40 J = ITMP1, ITMP2{2DO 60 J = ITMP1, ITMP2, 3DO 70 J = ITMP1, ITMP2, 4DO 100 J = ITMP1, ITMP2, 5DO 130 J = ITMP1, ITMP2, 6DO 10 J = ITMP1, ITMP2 - MOD( ITMP2-ITMP1+1, 2 ), 2, 7DO 20 J = ITMP2 - MOD( ITMP2-ITMP1+1, 2 ) + 1, ITMP2}, LIHIZDO 80 J = LILOZ, LIHIZ{2DO 110 J = LILOZ, LIHIZ}, LILOZDO 80 J = LILOZ, LIHIZ{2DO 110 J = LILOZ, LIHIZ} |
| K | <--- | ISTARTDO 50 K = ISTOP - MOD( ISTOP-ISTART+1, 3 ) + 1, ISTOP{2DO 90 K = ISTART, ISTOP - MOD( ISTOP-ISTART+1, 3 ), 3, 3DO 120 K = ISTOP - MOD( ISTOP-ISTART+1, 3 ) + 1, ISTOP, 4DO 30 K = ISTART, ISTOP - MOD( ISTOP-ISTART+1, 3 ), 3}, ISTOPDO 50 K = ISTOP - MOD( ISTOP-ISTART+1, 3 ) + 1, ISTOP{2DO 90 K = ISTART, ISTOP - MOD( ISTOP-ISTART+1, 3 ), 3, 3DO 120 K = ISTOP - MOD( ISTOP-ISTART+1, 3 ) + 1, ISTOP, 4DO 30 K = ISTART, ISTOP - MOD( ISTOP-ISTART+1, 3 ), 3} |
| SUM | <--- | ASUM = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 +{2SUM = DCONJG( T1 )*A( IROW1, J ) +, 3SUM = DCONJG( T1 )*A( IROW1, J ) +, 4SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 5SUM = T12*H11 + T22*H22 + T32*A( J, ICOL1+3 ), 6SUM = T13*H11 + T23*H22 + T33*A( J, ICOL1+4 ), 7SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 8SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 9SUM = DCONJG( T1 )*A( IROW1, J ) +, 10SUM = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 +}, H11SUM = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 +{2SUM = T12*H11 + T22*H22 + T32*A( J, ICOL1+3 ), 3SUM = T13*H11 + T23*H22 + T33*A( J, ICOL1+4 ), 4SUM = T12*H11 + T22*H22 + T32*Z( J, ICOL1+3 ), 5SUM = T13*H11 + T23*H22 + T33*Z( J, ICOL1+4 ), 6SUM = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 +}, H22SUM = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 +{2SUM = T12*H11 + T22*H22 + T32*A( J, ICOL1+3 ), 3SUM = T13*H11 + T23*H22 + T33*A( J, ICOL1+4 ), 4SUM = T12*H11 + T22*H22 + T32*Z( J, ICOL1+3 ), 5SUM = T13*H11 + T23*H22 + T33*Z( J, ICOL1+4 ), 6SUM = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 +}, ICOL1SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +{2SUM = T12*H11 + T22*H22 + T32*A( J, ICOL1+3 ), 3SUM = T13*H11 + T23*H22 + T33*A( J, ICOL1+4 ), 4SUM = T1*Z( J, ICOL1 ) + T2*Z( J, ICOL1+1 ) +, 5SUM = T12*H11 + T22*H22 + T32*Z( J, ICOL1+3 ), 6SUM = T13*H11 + T23*H22 + T33*Z( J, ICOL1+4 ), 7SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 8SUM = T1*Z( J, ICOL1 ) + T2*Z( J, ICOL1+1 ) +, 9SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +}, IROW1SUM = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 +{2SUM = DCONJG( T1 )*A( IROW1, J ) +, 3SUM = DCONJG( T1 )*A( IROW1, J ) +, 4SUM = DCONJG( T1 )*A( IROW1, J ) +, 5SUM = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 +}, JSUM = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 +{2SUM = DCONJG( T1 )*A( IROW1, J ) +, 3SUM = DCONJG( T1 )*A( IROW1, J ) +, 4SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 5SUM = T12*H11 + T22*H22 + T32*A( J, ICOL1+3 ), 6SUM = T13*H11 + T23*H22 + T33*A( J, ICOL1+4 ), 7SUM = T1*Z( J, ICOL1 ) + T2*Z( J, ICOL1+1 ) +, 8SUM = T12*H11 + T22*H22 + T32*Z( J, ICOL1+3 ), 9SUM = T13*H11 + T23*H22 + T33*Z( J, ICOL1+4 ), 10SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 11SUM = T1*Z( J, ICOL1 ) + T2*Z( J, ICOL1+1 ) +, 12SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 13SUM = DCONJG( T1 )*A( IROW1, J ) +, 14SUM = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 +}, T1SUM = DCONJG( T1 )*A( IROW1, J ) +{2SUM = DCONJG( T1 )*A( IROW1, J ) +, 3SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 4SUM = T1*Z( J, ICOL1 ) + T2*Z( J, ICOL1+1 ) +, 5SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 6SUM = T1*Z( J, ICOL1 ) + T2*Z( J, ICOL1+1 ) +, 7SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 8SUM = DCONJG( T1 )*A( IROW1, J ) +}, T12SUM = T12*H11 + T22*H22 + T32*A( J, ICOL1+3 ){2SUM = T12*H11 + T22*H22 + T32*Z( J, ICOL1+3 ), 3SUM = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 +}, T13SUM = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 +{2SUM = T13*H11 + T23*H22 + T33*A( J, ICOL1+4 ), 3SUM = T13*H11 + T23*H22 + T33*Z( J, ICOL1+4 )}, T2SUM = DCONJG( T1 )*A( IROW1, J ) +{2SUM = DCONJG( T1 )*A( IROW1, J ) +, 3SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 4SUM = T1*Z( J, ICOL1 ) + T2*Z( J, ICOL1+1 ) +, 5SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 6SUM = T1*Z( J, ICOL1 ) + T2*Z( J, ICOL1+1 ) +, 7SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 8SUM = DCONJG( T1 )*A( IROW1, J ) +}, T22SUM = T12*H11 + T22*H22 + T32*A( J, ICOL1+3 ){2SUM = T12*H11 + T22*H22 + T32*Z( J, ICOL1+3 ), 3SUM = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 +}, T23SUM = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 +{2SUM = T13*H11 + T23*H22 + T33*A( J, ICOL1+4 ), 3SUM = T13*H11 + T23*H22 + T33*Z( J, ICOL1+4 )}, T3SUM = DCONJG( T1 )*A( IROW1, J ) +{2SUM = DCONJG( T1 )*A( IROW1, J ) +, 3SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 4SUM = T1*Z( J, ICOL1 ) + T2*Z( J, ICOL1+1 ) +, 5SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 6SUM = T1*Z( J, ICOL1 ) + T2*Z( J, ICOL1+1 ) +, 7SUM = T1*A( J, ICOL1 ) + T2*A( J, ICOL1+1 ) +, 8SUM = DCONJG( T1 )*A( IROW1, J ) +}, T32SUM = T12*H11 + T22*H22 + T32*A( J, ICOL1+3 ){2SUM = T12*H11 + T22*H22 + T32*Z( J, ICOL1+3 ), 3SUM = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 +}, T33SUM = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 +{2SUM = T13*H11 + T23*H22 + T33*A( J, ICOL1+4 ), 3SUM = T13*H11 + T23*H22 + T33*Z( J, ICOL1+4 )}, ZSUM = T1*Z( J, ICOL1 ) + T2*Z( J, ICOL1+1 ) +{2SUM = T12*H11 + T22*H22 + T32*Z( J, ICOL1+3 ), 3SUM = T13*H11 + T23*H22 + T33*Z( J, ICOL1+4 ), 4SUM = T1*Z( J, ICOL1 ) + T2*Z( J, ICOL1+1 ) +} |
| SUM1 | <--- | A1SUM1 = DCONJG( T1 )*A1 + DCONJG( T2 )*A2 +, T1SUM1 = DCONJG( T1 )*A1 + DCONJG( T2 )*A2 +, T2SUM1 = DCONJG( T1 )*A1 + DCONJG( T2 )*A2 +, A2SUM1 = DCONJG( T1 )*A1 + DCONJG( T2 )*A2 +, T3SUM1 = DCONJG( T1 )*A1 + DCONJG( T2 )*A2 +, A3SUM1 = DCONJG( T1 )*A1 + DCONJG( T2 )*A2 + |
| SUM2 | <--- | H11SUM2 = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 +, H22SUM2 = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 +, T12SUM2 = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 +, T22SUM2 = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 +, T32SUM2 = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 +, A4SUM2 = DCONJG( T12 )*H11 + DCONJG( T22 )*H22 + |
| SUM3 | <--- | H11SUM3 = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 +, H22SUM3 = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 +, T13SUM3 = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 +, T23SUM3 = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 +, T33SUM3 = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 +, A5SUM3 = DCONJG( T13 )*H11 + DCONJG( T23 )*H22 + |
| T1 | <--- | KT1 = VECS( ( K-1 )*3+3 ){2T1 = VECS( ( K-1 )*3+3 ), 3T1 = VECS( ( K-1 )*3+3 ), 4T1 = VECS( ( K-1 )*3+3 )} |
| T12 | <--- | KT12 = VECS( ( K-1 )*3+6 ){2T12 = VECS( ( K-1 )*3+6 )} |
| T13 | <--- | KT13 = VECS( ( K-1 )*3+9 ){2T13 = VECS( ( K-1 )*3+9 )} |
| T2 | <--- | T1T2 = T1*V2{2T2 = T1*V2, 3T2 = T1*V2, 4T2 = T1*V2}, V2T2 = T1*V2{2T2 = T1*V2, 3T2 = T1*V2, 4T2 = T1*V2} |
| T22 | <--- | T12T22 = T12*V22{2T22 = T12*V22}, V22T22 = T12*V22{2T22 = T12*V22} |
| T23 | <--- | T13T23 = T13*V23{2T23 = T13*V23}, V23T23 = T13*V23{2T23 = T13*V23} |
| T3 | <--- | T1T3 = T1*V3{2T3 = T1*V3, 3T3 = T1*V3, 4T3 = T1*V3}, V3T3 = T1*V3{2T3 = T1*V3, 3T3 = T1*V3, 4T3 = T1*V3} |
| T32 | <--- | T12T32 = T12*V32{2T32 = T12*V32}, V32T32 = T12*V32{2T32 = T12*V32} |
| T33 | <--- | T13T33 = T13*V33{2T33 = T13*V33}, V33T33 = T13*V33{2T33 = T13*V33} |
| TMP1 | <--- | B2TMP1 = DCONJG( T1 )*B1 + DCONJG( T2 )*B2 +, B3TMP1 = DCONJG( T1 )*B1 + DCONJG( T2 )*B2 +, T1TMP1 = DCONJG( T1 )*B1 + DCONJG( T2 )*B2 +, T2TMP1 = DCONJG( T1 )*B1 + DCONJG( T2 )*B2 +, T3TMP1 = DCONJG( T1 )*B1 + DCONJG( T2 )*B2 +, B1TMP1 = DCONJG( T1 )*B1 + DCONJG( T2 )*B2 + |
| TMP2 | <--- | B4TMP2 = DCONJG( T12 )*A11 + DCONJG( T22 )*A22 +, A11TMP2 = DCONJG( T12 )*A11 + DCONJG( T22 )*A22 +, T12TMP2 = DCONJG( T12 )*A11 + DCONJG( T22 )*A22 +, T22TMP2 = DCONJG( T12 )*A11 + DCONJG( T22 )*A22 +, T32TMP2 = DCONJG( T12 )*A11 + DCONJG( T22 )*A22 +, A22TMP2 = DCONJG( T12 )*A11 + DCONJG( T22 )*A22 + |
| TMP3 | <--- | B5TMP3 = DCONJG( T13 )*A11 + DCONJG( T23 )*A22 +, A11TMP3 = DCONJG( T13 )*A11 + DCONJG( T23 )*A22 +, T13TMP3 = DCONJG( T13 )*A11 + DCONJG( T23 )*A22 +, T23TMP3 = DCONJG( T13 )*A11 + DCONJG( T23 )*A22 +, T33TMP3 = DCONJG( T13 )*A11 + DCONJG( T23 )*A22 +, A22TMP3 = DCONJG( T13 )*A11 + DCONJG( T23 )*A22 + |
| V2 | <--- | KV2 = VECS( ( K-1 )*3+1 ){2V2 = VECS( ( K-1 )*3+1 ), 3V2 = VECS( ( K-1 )*3+1 ), 4V2 = VECS( ( K-1 )*3+1 )} |
| V22 | <--- | KV22 = VECS( ( K-1 )*3+4 ){2V22 = VECS( ( K-1 )*3+4 )} |
| V23 | <--- | KV23 = VECS( ( K-1 )*3+7 ){2V23 = VECS( ( K-1 )*3+7 )} |
| V3 | <--- | KV3 = VECS( ( K-1 )*3+2 ){2V3 = VECS( ( K-1 )*3+2 ), 3V3 = VECS( ( K-1 )*3+2 ), 4V3 = VECS( ( K-1 )*3+2 )} |
| V32 | <--- | KV32 = VECS( ( K-1 )*3+5 ){2V32 = VECS( ( K-1 )*3+5 )} |
| V33 | <--- | KV33 = VECS( ( K-1 )*3+8 ){2V33 = VECS( ( K-1 )*3+8 )} |
| Z | <--- | H11Z( J, ICOL1+1 ) = H11 - SUM{2Z( J, ICOL1+2 ) = H11 - SUM}, H22Z( J, ICOL1+3 ) = H22 - SUM*DCONJG( V23 ), ICOL1Z( J, ICOL1 ) = Z( J, ICOL1 ) - SUM{2Z( J, ICOL1+4 ) = Z( J, ICOL1+4 ) -, 3Z( J, ICOL1 ) = Z( J, ICOL1 ) - SUM, 4Z( J, ICOL1+1 ) = Z( J, ICOL1+1 ) -, 5Z( J, ICOL1+2 ) = Z( J, ICOL1+2 ) -}, JZ( J, ICOL1 ) = Z( J, ICOL1 ) - SUM{2Z( J, ICOL1+4 ) = Z( J, ICOL1+4 ) -, 3Z( J, ICOL1 ) = Z( J, ICOL1 ) - SUM, 4Z( J, ICOL1+1 ) = Z( J, ICOL1+1 ) -, 5Z( J, ICOL1+2 ) = Z( J, ICOL1+2 ) -}, SUMZ( J, ICOL1 ) = Z( J, ICOL1 ) - SUM{2Z( J, ICOL1+1 ) = H11 - SUM, 3Z( J, ICOL1+2 ) = H11 - SUM, 4Z( J, ICOL1+3 ) = H22 - SUM*DCONJG( V23 ), 5Z( J, ICOL1+4 ) = Z( J, ICOL1+4 ) -, 6Z( J, ICOL1 ) = Z( J, ICOL1 ) - SUM, 7Z( J, ICOL1+1 ) = Z( J, ICOL1+1 ) -, 8Z( J, ICOL1+2 ) = Z( J, ICOL1+2 ) -}, V2Z( J, ICOL1+1 ) = Z( J, ICOL1+1 ) -, V23Z( J, ICOL1+3 ) = H22 - SUM*DCONJG( V23 ), V3Z( J, ICOL1+2 ) = Z( J, ICOL1+2 ) -, V33Z( J, ICOL1+4 ) = Z( J, ICOL1+4 ) -, ZZ( J, ICOL1 ) = Z( J, ICOL1 ) - SUM{2Z( J, ICOL1+4 ) = Z( J, ICOL1+4 ) -, 3Z( J, ICOL1 ) = Z( J, ICOL1 ) - SUM, 4Z( J, ICOL1+1 ) = Z( J, ICOL1+1 ) -, 5Z( J, ICOL1+2 ) = Z( J, ICOL1+2 ) -} |
|
|
Analysis elements of the routine ZLAREF()
Put the mouse over each element to display detailed matching information
 Assigned variables |
| | | A1 , A11 , A2 , A22 , A3 , A4 , A5 , B1 , B2 , B3 , B4 , B5 , H11 , H22 , ICOL1 , IROW1 , J , K , SUM , SUM1 , SUM2 , SUM3 , T1 , T12 , T13 , T2 , T22 , T23 , T3 , T32 , T33 , TMP1 , TMP2 , TMP3 , V2 , V22 , V23 , V3 , V32 , V33 |
|
| Active variables |
| | | A , A1 , A11 , A2 , A22 , A3 , A4 , A5 , B1 , B2 , B3 , B4 , B5 , BLOCK , H11 , H22 , ICOL1 , IROW1 , ISTART , ISTOP , ITMP1 , ITMP2 , J , K , LDA , LDZ , LIHIZ , LILOZ , LSAME , SUM , SUM1 , SUM2 , SUM3 , T1 , T12 , T13 , T2 , T22 , T23 , T3 , T32 , T33 , TMP1 , TMP2 , TMP3 , TYPE , V2 , V22 , V23 , V3 , V32 , V33 , VECS , WANTZ , Z |
|
| Accessed arrays [ array name : associated index ] |
| |  A | : IROW1, J , IROW1, J , IROW1, J , IROW1, J , IROW1, J , IROW1, J , IROW1, J , IROW1, J , IROW1, J+1 , IROW1, J+1 , IROW1+1, J , IROW1+1, J , IROW1+1, J , IROW1+1, J , IROW1+1, J , IROW1+1, J , IROW1+1, J , IROW1+1, J , IROW1+1, J , IROW1+1, J+1 , IROW1+1, J+1 , IROW1+2, J , IROW1+2, J , IROW1+2, J , IROW1+2, J , IROW1+2, J , IROW1+2, J , IROW1+2, J , IROW1+2, J , IROW1+2, J , IROW1+2, J+1 , IROW1+2, J+1 , IROW1+3, J , IROW1+3, J , IROW1+3, J , IROW1+3, J , IROW1+3, J , IROW1+3, J+1 , IROW1+3, J+1 , IROW1+4, J , IROW1+4, J , IROW1+4, J , IROW1+4, J , IROW1+4, J+1 , IROW1+4, J+1 , J, ICOL1 , J, ICOL1 , J, ICOL1 , J, ICOL1 , J, ICOL1 , J, ICOL1 , J, ICOL1+1 , J, ICOL1+1 , J, ICOL1+1 , J, ICOL1+1 , J, ICOL1+1 , J, ICOL1+1 , J, ICOL1+1 , J, ICOL1+2 , J, ICOL1+2 , J, ICOL1+2 , J, ICOL1+2 , J, ICOL1+2 , J, ICOL1+2 , J, ICOL1+2 , J, ICOL1+3 , J, ICOL1+3 , J, ICOL1+3 , J, ICOL1+4 , J, ICOL1+4 |
| | LSAME | : TYPE, 'R' |
| | VECS | : ( K-1 )*3+1 , ( K-1 )*3+1 , ( K-1 )*3+1 , ( K-1 )*3+1 , ( K-1 )*3+2 , ( K-1 )*3+2 , ( K-1 )*3+2 , ( K-1 )*3+2 , ( K-1 )*3+3 , ( K-1 )*3+3 , ( K-1 )*3+3 , ( K-1 )*3+3 , ( K-1 )*3+4 , ( K-1 )*3+4 , ( K-1 )*3+5 , ( K-1 )*3+5 , ( K-1 )*3+6 , ( K-1 )*3+6 , ( K-1 )*3+7 , ( K-1 )*3+7 , ( K-1 )*3+8 , ( K-1 )*3+8 , ( K-1 )*3+9 , ( K-1 )*3+9 |
| | Z | : J, ICOL1 , J, ICOL1 , J, ICOL1 , J, ICOL1 , J, ICOL1+1 , J, ICOL1+1 , J, ICOL1+1 , J, ICOL1+1 , J, ICOL1+1 , J, ICOL1+2 , J, ICOL1+2 , J, ICOL1+2 , J, ICOL1+2 , J, ICOL1+2 , J, ICOL1+3 , J, ICOL1+3 , J, ICOL1+3 , J, ICOL1+4 , J, ICOL1+4 |
|
| Conditional statements [ statement : associated predicate ] |
| | do | : ( 30 K = ISTART , ISTOP - MOD( ISTOP - ISTART + 1 , 3 ) , 3 ) , ( 10 J = ITMP1 , ITMP2 - MOD( ITMP2 - ITMP1 + 1 , 2 ) , 2 ) , ( 20 J = ITMP2 - MOD( ITMP2 - ITMP1 + 1 , 2 ) + 1 , ITMP2 ) , ( 50 K = ISTOP - MOD( ISTOP - ISTART + 1 , 3 ) + 1 , ISTOP ) , ( 40 J = ITMP1 , ITMP2 ) , ( 60 J = ITMP1 , ITMP2 ) , ( column transforms ) , ( 90 K = ISTART , ISTOP - MOD( ISTOP - ISTART + 1 , 3 ) , 3 ) , ( 70 J = ITMP1 , ITMP2 ) , ( 80 J = LILOZ , LIHIZ ) , ( 120 K = ISTOP - MOD( ISTOP - ISTART + 1 , 3 ) + 1 , ISTOP ) , ( 100 J = ITMP1 , ITMP2 ) , ( 110 J = LILOZ , LIHIZ ) , ( 130 J = ITMP1 , ITMP2 ) |
| | if | : ( (LSAME( TYPE , 'R' ) ) ) , ( BLOCK ) , ( BLOCK ) , ( WANTZ ) , ( WANTZ ) |
|
| List of variables |
A
A1
A11
A2
A22
A3
A4
| A5
B1
B2
B3
B4
B5
BLOCK
H11
| H22
ICOL1
IROW1
ISTART
ISTOP
ITMP1
ITMP2
J
| K
LDA
LDZ
LIHIZ
LILOZ
LSAME
SUM
SUM1
| SUM2
SUM3
T1
T12
T13
T2
T22
T23
| T3
T32
T33
TMP1
TMP2
TMP3
TYPE
V2
| V22
V23
V3
V32
V33
WANTZ
Z | | close
| |
A
A1
A11
A2
A22
A3
A4
A5
B1
B2
B3
B4
B5
BLOCK
H11
H22
ICOL1
IROW1
ISTART
ISTOP
ITMP1
ITMP2
J
K
LDA
LDZ
LIHIZ
LILOZ
LSAME
SUM
SUM1
SUM2
SUM3
T1
T12
T13
T2
T22
T23
T3
T32
T33
TMP1
TMP2
TMP3
TYPE
V2
V22
V23
V3
V32
V33
WANTZ
Z
| |
