|
|
| |
| # lines: |
411 | | # code: |
411 | | # comment: | 0 | |
# blank: | 0 |
| # Variables: | 53 |
| # Callers: | 1 |
| # Callings: | 0 |
| # Words: | 182 |
| # Keywords: | 104 |
|
|
|
|
|
..
.. Array Arguments ..
..
Purpose
=======
PCLASE2 initializes an M-by-N distributed matrix sub( A ) denoting
A(IA:IA+M-1,JA:JA+N-1) to BETA on the diagonal and ALPHA on the
offdiagonals. PCLASE2 requires that only dimension of the matrix
operand is distributed.
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
=========
UPLO (global input) CHARACTER
Specifies the part of the distributed matrix sub( A ) to be
set:
= 'U': Upper triangular part is set; the strictly lower
triangular part of sub( A ) is not changed;
= 'L': Lower triangular part is set; the strictly upper
triangular part of sub( A ) is not changed;
Otherwise: All of the matrix sub( A ) is set.
M (global input) INTEGER
The number of rows to be operated on i.e the number of rows
of the distributed submatrix sub( A ). M >= 0.
N (global input) INTEGER
The number of columns to be operated on i.e the number of
columns of the distributed submatrix sub( A ). N >= 0.
ALPHA (global input) COMPLEX
The constant to which the offdiagonal elements are to be
set.
BETA (global input) COMPLEX
The constant to which the diagonal elements are to be set.
A (local output) COMPLEX pointer into the local memory
to an array of dimension (LLD_A,LOCc(JA+N-1)). This array
contains the local pieces of the distributed matrix sub( A )
to be set. On exit, the leading M-by-N submatrix sub( A )
is set as follows:
if UPLO = 'U', A(IA+i-1,JA+j-1) = ALPHA, 1<=i<=j-1, 1<=j<=N,
if UPLO = 'L', A(IA+i-1,JA+j-1) = ALPHA, j+1<=i<=M, 1<=j<=N,
otherwise, A(IA+i-1,JA+j-1) = ALPHA, 1<=i<=M, 1<=j<=N,
IA+i.NE.JA+j,
and, for all UPLO, A(IA+i-1,JA+i-1) = BETA, 1<=i<=min(M,N).
IA (global input) INTEGER
The row index in the global array A indicating the first
row of sub( A ).
JA (global input) INTEGER
The column index in the global array A indicating the
first column of sub( A ).
DESCA (global and local input) INTEGER array of dimension DLEN_.
The array descriptor for the distributed matrix A.
=====================================================================
.. Parameters ..
|
|
|
|
001 SUBROUTINE PCLASE2( UPLO , M , N , ALPHA , BETA , A , IA , JA , DESCA )
002
003 * -- ScaLAPACK auxiliary routine(version 1.7) --
004 * University of Tennessee , Knoxville , Oak Ridge National Laboratory ,
005 * and University of California , Berkeley.
006 * May 1 , 1997
007
008 * .. Scalar Arguments ..
009 CHARACTER UPLO
010 INTEGER IA , JA , M , N
011 COMPLEX ALPHA , BETA
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 * ..
018 * .. Local Scalars ..
019 INTEGER HEIGHT , IACOL , IAROW , IBASE , ICOFFA , II , IIA ,
020 $IIBEG , IIEND , IINXT , ILEFT , IRIGHT , IROFFA ,
021 $ITOP , JJ , JJA , JJBEG , JJEND , JJNXT , LDA , MBA ,
022 $MP , MPA , MYCOL , MYDIST , MYROW , NBA , NPCOL ,
023 $NPROW , NQ , NQA , WIDE
024 * ..
025 * .. External Subroutines ..
026 EXTERNAL BLACS_GRIDINFO , CLASET , INFOG2L
027 * ..
028 * .. External Functions ..
029 LOGICAL LSAME
030 INTEGER ICEIL , NUMROC
031 EXTERNAL ICEIL , LSAME , NUMROC
032 * ..
033 * .. Intrinsic Functions ..
034 INTRINSIC MAX , MIN , MOD
035 * ..
036 * .. Executable Statements ..
037
038 IF( M.EQ.0 .OR. N.EQ.0 )
038  
039 $ RETURN
040
041 * Get grid parameters
042
043 CALL BLACS_GRIDINFO( DESCA( CTXT_ ) , NPROW , NPCOL , MYROW , MYCOL )
044
045 CALL INFOG2L( IA , JA , DESCA , NPROW , NPCOL , MYROW , MYCOL , IIA , JJA ,
046 $ IAROW , IACOL )
047 MBA = DESCA( MB_ )
048 NBA = DESCA( NB_ )
049 LDA = DESCA( LLD_ )
050 IROFFA = MOD( IA - 1 , MBA )
051 ICOFFA = MOD( JA - 1 , NBA )
052
053 IF( N.LE.( NBA - ICOFFA ) ) THEN
054
055 * It is assumed that the local columns JJA : JJA + N - 1 of the matrix
056 * A are in the same process column(IACOL).
057
058 * N
059 * JJA JJA + N - 1
060 * / --------------------- \
061 * IROFFA| | | |
062 * \ |...................| |( IAROW )
063 * IIA |x | | MB_A
064 * | x | |
065 * | -- x ---------------- | /
066 * | x |
067 * | x | ITOP
068 * | x | |
069 * | x | /------- \
070 * | ------- x ----------- | | ------- x ----------- |
071 * | x | | x |
072 * | x | | x |
073 * | x | | x |
074 * | x | | x |
075 * | ------------ x ------ | | ------------ x ------ |
076 * | x | \____________ /
077 * | x | |
078 * | x | IBASE
079 * | x |
080 * | ----------------- x - | Local picture
081 * | x|
082 * | |
083 * | |
084 * | |
085 * | ------------------- |
086 * | |
087 * . .
088 * . .
089 * .(IACOL) .
090
090
091 IF( MYCOL.EQ.IACOL ) THEN
092
092
093 MPA = NUMROC( M + IROFFA , MBA , MYROW , IAROW , NPROW )
094 IF( MPA.LE.0 )
094
095 $ RETURN
096 IF( MYROW.EQ.IAROW )
096
097 $ MPA = MPA - IROFFA
098 MYDIST = MOD( MYROW - IAROW + NPROW , NPROW )
099 ITOP = MYDIST * MBA - IROFFA
100
101 IF( LSAME( UPLO , 'U' ) ) THEN
102
102
103 ITOP = MAX( 0 , ITOP )
104 IIBEG = IIA
105 IIEND = IIA + MPA - 1
106 IINXT = MIN( ICEIL( IIBEG , MBA ) * MBA , IIEND )
107
108 10 CONTINUE
109 IF(( N - ITOP ).GT.0 ) THEN
109
110 CALL CLASET( UPLO , IINXT - IIBEG + 1 , N - ITOP , ALPHA , BETA ,
111 $ A( IIBEG + (JJA + ITOP - 1)*LDA ) , LDA )
112 MYDIST = MYDIST + NPROW
113 ITOP = MYDIST * MBA - IROFFA
114 IIBEG = IINXT + 1
115 IINXT = MIN( IINXT + MBA , IIEND )
116 GO TO 10
117 END IF
118
119 ELSE IF( LSAME( UPLO , 'L' ) ) THEN
120
120
121 II = IIA
122 JJ = JJA
123 MP = MPA
124 IBASE = MIN( ITOP + MBA , N )
125 ITOP = MIN( MAX( 0 , ITOP ) , N )
126
127 20 CONTINUE
128 IF( JJ.LE.( JJA + N - 1 ) ) THEN
128
129 HEIGHT = IBASE - ITOP
130 CALL CLASET( 'All' , MP , ITOP - JJ + JJA , ALPHA , ALPHA ,
131 $ A( II + (JJ - 1)*LDA ) , LDA )
132 CALL CLASET( UPLO , MP , HEIGHT , ALPHA , BETA ,
133 $ A( II + (JJA + ITOP - 1)*LDA ) , LDA )
134 MP = MAX( 0 , MP - HEIGHT )
135 II = II + HEIGHT
136 JJ = JJA + IBASE
137 MYDIST = MYDIST + NPROW
138 ITOP = MYDIST * MBA - IROFFA
139 IBASE = MIN( ITOP + MBA , N )
140 ITOP = MIN( ITOP , N )
141 GO TO 20
142 END IF
143
144 ELSE
145
145
146 II = IIA
147 JJ = JJA
148 MP = MPA
149 IBASE = MIN( ITOP + MBA , N )
150 ITOP = MIN( MAX( 0 , ITOP ) , N )
151
152 30 CONTINUE
153 IF( JJ.LE.( JJA + N - 1 ) ) THEN
153
154 HEIGHT = IBASE - ITOP
155 CALL CLASET( 'All' , MPA , ITOP - JJ + JJA , ALPHA , ALPHA ,
156 $ A( IIA + (JJ - 1)*LDA ) , LDA )
157 CALL CLASET( 'All' , MPA - MP , HEIGHT , ALPHA , ALPHA ,
158 $ A( IIA + (JJA + ITOP - 1)*LDA ) , LDA )
159 CALL CLASET( 'All' , MP , HEIGHT , ALPHA , BETA ,
160 $ A( II + (JJA + ITOP - 1)*LDA ) , LDA )
161 MP = MAX( 0 , MP - HEIGHT )
162 II = II + HEIGHT
163 JJ = JJA + IBASE
164 MYDIST = MYDIST + NPROW
165 ITOP = MYDIST * MBA - IROFFA
166 IBASE = MIN( ITOP + MBA , N )
167 ITOP = MIN( ITOP , N )
168 GO TO 30
169 END IF
170
171 END IF
172
173 END IF
174
175 ELSE IF( M.LE.( MBA - IROFFA ) ) THEN
176
177 * It is assumed that the local rows IIA : IIA + M - 1 of the matrix A
178 * are in the same process row(IAROW).
179
180 * ICOFFA
181 * / \JJA
182 * IIA ------------------ .... --------
183 * | .x | | | / | | \
184 * | . x | | | ILEFT| | | |
185 * | . x | | | | | |
186 * | . x | | \ x | |
187 * | . |x | | |x | | IRIGHT
188 * | . | x | | | x | |
189 * (IAROW) | . | x | | | x | |
190 * | . | x| | | x| |
191 * | . | x | | x /
192 * | . | |x | | |
193 * | . | | x | | |
194 * | . | | x | | |
195 * | . | | x| | |
196 * IIA + M - 1 ------------------ .... -------
197 * NB_A
198 * (IACOL) Local picture
199
199
200 IF( MYROW.EQ.IAROW ) THEN
201
201
202 NQA = NUMROC( N + ICOFFA , NBA , MYCOL , IACOL , NPCOL )
203 IF( NQA.LE.0 )
203
204 $ RETURN
205 IF( MYCOL.EQ.IACOL )
205
206 $ NQA = NQA - ICOFFA
207 MYDIST = MOD( MYCOL - IACOL + NPCOL , NPCOL )
208 ILEFT = MYDIST * NBA - ICOFFA
209
210 IF( LSAME( UPLO , 'L' ) ) THEN
211
211
212 ILEFT = MAX( 0 , ILEFT )
213 JJBEG = JJA
214 JJEND = JJA + NQA - 1
215 JJNXT = MIN( ICEIL( JJBEG , NBA ) * NBA , JJEND )
216
217 40 CONTINUE
218 IF(( M - ILEFT ).GT.0 ) THEN
218
219 CALL CLASET( UPLO , M - ILEFT , JJNXT - JJBEG + 1 , ALPHA ,
220 $ BETA , A( IIA + ILEFT + (JJBEG - 1)*LDA ) , LDA )
221 MYDIST = MYDIST + NPCOL
222 ILEFT = MYDIST * NBA - ICOFFA
223 JJBEG = JJNXT + 1
224 JJNXT = MIN( JJNXT + NBA , JJEND )
225 GO TO 40
226 END IF
227
228 ELSE IF( LSAME( UPLO , 'U' ) ) THEN
229
229
230 II = IIA
231 JJ = JJA
232 NQ = NQA
233 IRIGHT = MIN( ILEFT + NBA , M )
234 ILEFT = MIN( MAX( 0 , ILEFT ) , M )
235
236 50 CONTINUE
237 IF( II.LE.( IIA + M - 1 ) ) THEN
237
238 WIDE = IRIGHT - ILEFT
239 CALL CLASET( 'All' , ILEFT - II + IIA , NQ , ALPHA , ALPHA ,
240 $ A( II + (JJ - 1)*LDA ) , LDA )
241 CALL CLASET( UPLO , WIDE , NQ , ALPHA , BETA ,
242 $ A( IIA + ILEFT + (JJ - 1)*LDA ) , LDA )
243 NQ = MAX( 0 , NQ - WIDE )
244 II = IIA + IRIGHT
245 JJ = JJ + WIDE
246 MYDIST = MYDIST + NPCOL
247 ILEFT = MYDIST * NBA - ICOFFA
248 IRIGHT = MIN( ILEFT + NBA , M )
249 ILEFT = MIN( ILEFT , M )
250 GO TO 50
251 END IF
252
253 ELSE
254
254
255 II = IIA
256 JJ = JJA
257 NQ = NQA
258 IRIGHT = MIN( ILEFT + NBA , M )
259 ILEFT = MIN( MAX( 0 , ILEFT ) , M )
260
261 60 CONTINUE
262 IF( II.LE.( IIA + M - 1 ) ) THEN
262
263 WIDE = IRIGHT - ILEFT
264 CALL CLASET( 'All' , ILEFT - II + IIA , NQA , ALPHA , ALPHA ,
265 $ A( II + (JJA - 1)*LDA ) , LDA )
266 CALL CLASET( 'All' , WIDE , NQA - NQ , ALPHA , ALPHA ,
267 $ A( IIA + ILEFT + (JJA - 1)*LDA ) , LDA )
268 CALL CLASET( 'All' , WIDE , NQ , ALPHA , BETA ,
269 $ A( IIA + ILEFT + (JJ - 1)*LDA ) , LDA )
270 NQ = MAX( 0 , NQ - WIDE )
271 II = IIA + IRIGHT
272 JJ = JJ + WIDE
273 MYDIST = MYDIST + NPCOL
274 ILEFT = MYDIST * NBA - ICOFFA
275 IRIGHT = MIN( ILEFT + NBA , M )
276 ILEFT = MIN( ILEFT , M )
277 GO TO 60
278 END IF
279
280 END IF
281
282 END IF
283
284 END IF
285
286 RETURN
287
288 * End of PCLASE2
289
290 END42
21
|
|
Variables in Routine PCLASE2()
| Summary Report |
| Data Type | Quantity | Size(byte) |
| CHARACTER | 1 | 1 |
| COMPLEX | 2 | 8 |
| INTEGER | 49 | 196 |
| LOGICAL | 1 | 1 |
| TOTAL | 53 | 206 |
List of Variables
CHARACTER
COMPLEX
INTEGER
| BLOCK_CYCLIC_2D | CSRC_ | CTXT_ | DLEN_ | DTYPE_ |
| HEIGHT | IA | IACOL | IAROW | IBASE |
| ICEIL | ICOFFA | II | IIA | IIBEG |
| IIEND | IINXT | ILEFT | IRIGHT | IROFFA |
| ITOP | JA | JJ | JJA | JJBEG |
| JJEND | JJNXT | LDA | LLD_ | M |
| M_ | MB_ | MBA | MP | MPA |
| MYCOL | MYDIST | MYROW | N | N_ |
| NB_ | NBA | NPCOL | NPROW | NQ |
| NQA | NUMROC | RSRC_ | WIDE | |
LOGICAL
Variables Dependence Graph
Put the mouse over a right hand side variable to display the corresponding line of the dependence | | - | | - | - | | HEIGHT | <--- | IBASEHEIGHT = IBASE - ITOP{2HEIGHT = IBASE - ITOP}, ITOPHEIGHT = IBASE - ITOP{2HEIGHT = IBASE - ITOP} |
| IBASE | <--- | ITOPIBASE = MIN( ITOP+MBA, N ){2IBASE = MIN( ITOP + MBA, N ), 3IBASE = MIN( ITOP+MBA, N ), 4IBASE = MIN( ITOP + MBA, N )}, MBAIBASE = MIN( ITOP+MBA, N ){2IBASE = MIN( ITOP + MBA, N ), 3IBASE = MIN( ITOP+MBA, N ), 4IBASE = MIN( ITOP + MBA, N )}, NIBASE = MIN( ITOP+MBA, N ){2IBASE = MIN( ITOP + MBA, N ), 3IBASE = MIN( ITOP+MBA, N ), 4IBASE = MIN( ITOP + MBA, N )} |
| ICOFFA | <--- | JAICOFFA = MOD( JA-1, NBA ), NBAICOFFA = MOD( JA-1, NBA ) |
| II | <--- | IIII = II + HEIGHT{2II = II + HEIGHT}, IIAII = IIA{2II = IIA, 3II = IIA, 4II = IIA + IRIGHT, 5II = IIA, 6II = IIA + IRIGHT}, IRIGHTII = IIA + IRIGHT{2II = IIA + IRIGHT}, HEIGHTII = II + HEIGHT{2II = II + HEIGHT} |
| IIBEG | <--- | IIAIIBEG = IIA, IINXTIIBEG = IINXT +1 |
| IIEND | <--- | IIAIIEND = IIA + MPA - 1, MPAIIEND = IIA + MPA - 1 |
| IINXT | <--- | ICEILIINXT = MIN( ICEIL( IIBEG, MBA ) * MBA, IIEND ), IIBEGIINXT = MIN( ICEIL( IIBEG, MBA ) * MBA, IIEND ), IIENDIINXT = MIN( ICEIL( IIBEG, MBA ) * MBA, IIEND ){2IINXT = MIN( IINXT+MBA, IIEND )}, IINXTIINXT = MIN( IINXT+MBA, IIEND ), MBAIINXT = MIN( ICEIL( IIBEG, MBA ) * MBA, IIEND ){2IINXT = MIN( IINXT+MBA, IIEND )} |
| ILEFT | <--- | ICOFFAILEFT = MYDIST * NBA - ICOFFA{2ILEFT = MYDIST * NBA - ICOFFA, 3ILEFT = MYDIST * NBA - ICOFFA, 4ILEFT = MYDIST * NBA - ICOFFA}, ILEFTILEFT = MAX( 0, ILEFT ){2ILEFT = MIN( MAX( 0, ILEFT ), M ), 3ILEFT = MIN( ILEFT, M ), 4ILEFT = MIN( MAX( 0, ILEFT ), M ), 5ILEFT = MIN( ILEFT, M )}, MILEFT = MIN( MAX( 0, ILEFT ), M ){2ILEFT = MIN( ILEFT, M ), 3ILEFT = MIN( MAX( 0, ILEFT ), M ), 4ILEFT = MIN( ILEFT, M )}, MYDISTILEFT = MYDIST * NBA - ICOFFA{2ILEFT = MYDIST * NBA - ICOFFA, 3ILEFT = MYDIST * NBA - ICOFFA, 4ILEFT = MYDIST * NBA - ICOFFA}, NBAILEFT = MYDIST * NBA - ICOFFA{2ILEFT = MYDIST * NBA - ICOFFA, 3ILEFT = MYDIST * NBA - ICOFFA, 4ILEFT = MYDIST * NBA - ICOFFA} |
| IRIGHT | <--- | ILEFTIRIGHT = MIN( ILEFT+NBA, M ){2IRIGHT = MIN( ILEFT + NBA, M ), 3IRIGHT = MIN( ILEFT+NBA, M ), 4IRIGHT = MIN( ILEFT + NBA, M )}, MIRIGHT = MIN( ILEFT+NBA, M ){2IRIGHT = MIN( ILEFT + NBA, M ), 3IRIGHT = MIN( ILEFT+NBA, M ), 4IRIGHT = MIN( ILEFT + NBA, M )}, NBAIRIGHT = MIN( ILEFT+NBA, M ){2IRIGHT = MIN( ILEFT + NBA, M ), 3IRIGHT = MIN( ILEFT+NBA, M ), 4IRIGHT = MIN( ILEFT + NBA, M )} |
| IROFFA | <--- | MBAIROFFA = MOD( IA-1, MBA ), IAIROFFA = MOD( IA-1, MBA ) |
| ITOP | <--- | IROFFAITOP = MYDIST * MBA - IROFFA{2ITOP = MYDIST * MBA - IROFFA, 3ITOP = MYDIST * MBA - IROFFA, 4ITOP = MYDIST * MBA - IROFFA}, ITOPITOP = MAX( 0, ITOP ){2ITOP = MIN( MAX( 0, ITOP ), N ), 3ITOP = MIN( ITOP, N ), 4ITOP = MIN( MAX( 0, ITOP ), N ), 5ITOP = MIN( ITOP, N )}, MBAITOP = MYDIST * MBA - IROFFA{2ITOP = MYDIST * MBA - IROFFA, 3ITOP = MYDIST * MBA - IROFFA, 4ITOP = MYDIST * MBA - IROFFA}, MYDISTITOP = MYDIST * MBA - IROFFA{2ITOP = MYDIST * MBA - IROFFA, 3ITOP = MYDIST * MBA - IROFFA, 4ITOP = MYDIST * MBA - IROFFA}, NITOP = MIN( MAX( 0, ITOP ), N ){2ITOP = MIN( ITOP, N ), 3ITOP = MIN( MAX( 0, ITOP ), N ), 4ITOP = MIN( ITOP, N )} |
| JJ | <--- | IBASEJJ = JJA + IBASE{2JJ = JJA + IBASE}, JJJJ = JJ + WIDE{2JJ = JJ + WIDE}, JJAJJ = JJA{2JJ = JJA + IBASE, 3JJ = JJA, 4JJ = JJA + IBASE, 5JJ = JJA, 6JJ = JJA}, WIDEJJ = JJ + WIDE{2JJ = JJ + WIDE} |
| JJBEG | <--- | JJAJJBEG = JJA, JJNXTJJBEG = JJNXT +1 |
| JJEND | <--- | JJAJJEND = JJA + NQA - 1, NQAJJEND = JJA + NQA - 1 |
| JJNXT | <--- | ICEILJJNXT = MIN( ICEIL( JJBEG, NBA ) * NBA, JJEND ), JJBEGJJNXT = MIN( ICEIL( JJBEG, NBA ) * NBA, JJEND ), JJENDJJNXT = MIN( ICEIL( JJBEG, NBA ) * NBA, JJEND ){2JJNXT = MIN( JJNXT+NBA, JJEND )}, JJNXTJJNXT = MIN( JJNXT+NBA, JJEND ), NBAJJNXT = MIN( ICEIL( JJBEG, NBA ) * NBA, JJEND ){2JJNXT = MIN( JJNXT+NBA, JJEND )} |
| LDA | <--- | LLD_LDA = DESCA( LLD_ ) |
| MBA | <--- | MB_MBA = DESCA( MB_ ) |
| MP | <--- | MPMP = MAX( 0, MP - HEIGHT ){2MP = MAX( 0, MP - HEIGHT )}, MPAMP = MPA{2MP = MPA}, HEIGHTMP = MAX( 0, MP - HEIGHT ){2MP = MAX( 0, MP - HEIGHT )} |
| MPA | <--- | IAROWMPA = NUMROC( M+IROFFA, MBA, MYROW, IAROW, NPROW ), IROFFAMPA = NUMROC( M+IROFFA, MBA, MYROW, IAROW, NPROW ), MMPA = NUMROC( M+IROFFA, MBA, MYROW, IAROW, NPROW ), MBAMPA = NUMROC( M+IROFFA, MBA, MYROW, IAROW, NPROW ), MYROWMPA = NUMROC( M+IROFFA, MBA, MYROW, IAROW, NPROW ), NPROWMPA = NUMROC( M+IROFFA, MBA, MYROW, IAROW, NPROW ), NUMROCMPA = NUMROC( M+IROFFA, MBA, MYROW, IAROW, NPROW ) |
| MYDIST | <--- | IACOLMYDIST = MOD( MYCOL-IACOL+NPCOL, NPCOL ), IAROWMYDIST = MOD( MYROW-IAROW+NPROW, NPROW ), MYCOLMYDIST = MOD( MYCOL-IACOL+NPCOL, NPCOL ), MYDISTMYDIST = MYDIST + NPROW{2MYDIST = MYDIST + NPROW, 3MYDIST = MYDIST + NPROW, 4MYDIST = MYDIST + NPCOL, 5MYDIST = MYDIST + NPCOL, 6MYDIST = MYDIST + NPCOL}, MYROWMYDIST = MOD( MYROW-IAROW+NPROW, NPROW ), NPCOLMYDIST = MOD( MYCOL-IACOL+NPCOL, NPCOL ){2MYDIST = MYDIST + NPCOL, 3MYDIST = MYDIST + NPCOL, 4MYDIST = MYDIST + NPCOL}, NPROWMYDIST = MYDIST + NPROW{2MYDIST = MYDIST + NPROW, 3MYDIST = MYDIST + NPROW, 4MYDIST = MOD( MYROW-IAROW+NPROW, NPROW )} |
| NBA | <--- | NB_NBA = DESCA( NB_ ) |
| NQ | <--- | NQNQ = MAX( 0, NQ - WIDE ){2NQ = MAX( 0, NQ - WIDE )}, NQANQ = NQA{2NQ = NQA}, WIDENQ = MAX( 0, NQ - WIDE ){2NQ = MAX( 0, NQ - WIDE )} |
| NQA | <--- | IACOLNQA = NUMROC( N+ICOFFA, NBA, MYCOL, IACOL, NPCOL ), ICOFFANQA = NUMROC( N+ICOFFA, NBA, MYCOL, IACOL, NPCOL ), MYCOLNQA = NUMROC( N+ICOFFA, NBA, MYCOL, IACOL, NPCOL ), NNQA = NUMROC( N+ICOFFA, NBA, MYCOL, IACOL, NPCOL ), NBANQA = NUMROC( N+ICOFFA, NBA, MYCOL, IACOL, NPCOL ), NPCOLNQA = NUMROC( N+ICOFFA, NBA, MYCOL, IACOL, NPCOL ), NUMROCNQA = NUMROC( N+ICOFFA, NBA, MYCOL, IACOL, NPCOL ) |
| WIDE | <--- | ILEFTWIDE = IRIGHT - ILEFT{2WIDE = IRIGHT - ILEFT}, IRIGHTWIDE = IRIGHT - ILEFT{2WIDE = IRIGHT - ILEFT} |
|
|
Analysis elements of the routine PCLASE2()
Put the mouse over each element to display detailed matching information
 Assigned variables |
| | | BLOCK_CYCLIC_2D , CSRC_ , CTXT_ , DLEN_ , DTYPE_ , HEIGHT , IBASE , ICOFFA , II , IIBEG , IIEND , IINXT , ILEFT , IRIGHT , IROFFA , ITOP , JJ , JJBEG , JJEND , JJNXT , LDA , LLD_ , M_ , MB_ , MBA , MP , MPA , MYDIST , N_ , NB_ , NBA , NQ , NQA , RSRC_ , WIDE |
|
| Active variables |
| | | A , ALPHA , BETA , BLOCK_CYCLIC_2D , CSRC_ , CTXT_ , DESCA , DLEN_ , DTYPE_ , HEIGHT , IA , IACOL , IAROW , IBASE , ICEIL , ICOFFA , II , IIA , IIBEG , IIEND , IINXT , ILEFT , IRIGHT , IROFFA , ITOP , JA , JJ , JJA , JJBEG , JJEND , JJNXT , LDA , LLD_ , LSAME , M , M_ , MB_ , MBA , MP , MPA , MYCOL , MYDIST , MYROW , N , N_ , NB_ , NBA , NPCOL , NPROW , NQ , NQA , NUMROC , RSRC_ , UPLO , WIDE |
|
| Accessed arrays [ array name : associated index ] |
| |  A | : II+(JJ-1)*LDA , II+(JJ-1)*LDA , II+(JJA+ITOP-1)*LDA , II+(JJA+ITOP-1)*LDA , II+(JJA-1)*LDA , IIA+(JJ-1)*LDA , IIA+(JJA+ITOP-1)*LDA , IIA+ILEFT+(JJ-1)*LDA , IIA+ILEFT+(JJ-1)*LDA , IIA+ILEFT+(JJA-1)*LDA , IIA+ILEFT+(JJBEG-1)*LDA , IIBEG+(JJA+ITOP-1)*LDA |
| | DESCA | : CTXT_ , LLD_ , MB_ , NB_ |
| | ICEIL | : IIBEG, MBA , JJBEG, NBA |
| | LSAME | : UPLO, 'L' , UPLO, 'L' , UPLO, 'U' , UPLO, 'U' |
| | NUMROC | : M+IROFFA, MBA, MYROW, IAROW, NPROW , N+ICOFFA, NBA, MYCOL, IACOL, NPCOL |
|
| Conditional statements [ statement : associated predicate ] |
| | if | : ( M.EQ.0 .OR. N.EQ.0 ) , ( (N.LE.( NBA - ICOFFA ) ) ) , ( MYCOL.EQ.IACOL ) , ( MPA.LE.0 ) , ( MYROW.EQ.IAROW ) , ( (LSAME( UPLO , 'U' ) ) ) , ( (( N-ITOP ).GT.0 ) ) , ( (LSAME( UPLO , 'L' ) ) ) , ( (JJ.LE.( JJA + N - 1 ) ) ) , ( (JJ.LE.( JJA + N - 1 ) ) ) , ( (M.LE.( MBA - IROFFA ) ) ) , ( MYROW.EQ.IAROW ) , ( NQA.LE.0 ) , ( MYCOL.EQ.IACOL ) , ( (LSAME( UPLO , 'L' ) ) ) , ( (( M-ILEFT ).GT.0 ) ) , ( (LSAME( UPLO , 'U' ) ) ) , ( (II.LE.( IIA + M - 1 ) ) ) , ( (II.LE.( IIA + M - 1 ) ) ) |
|
| List of variables |
ALPHA
BETA
BLOCK_CYCLIC_2D
CSRC_
CTXT_
DLEN_
DTYPE_
| HEIGHT
IA
IACOL
IAROW
IBASE
ICEIL
ICOFFA
II
| IIA
IIBEG
IIEND
IINXT
ILEFT
IRIGHT
IROFFA
ITOP
| JA
JJ
JJA
JJBEG
JJEND
JJNXT
LDA
LLD_
| LSAME
M
M_
MB_
MBA
MP
MPA
MYCOL
| MYDIST
MYROW
N
N_
NB_
NBA
NPCOL
NPROW
| NQ
NQA
NUMROC
RSRC_
UPLO
WIDE | | close
| |
ALPHA
BETA
BLOCK_CYCLIC_2D
CSRC_
CTXT_
DLEN_
DTYPE_
HEIGHT
IA
IACOL
IAROW
IBASE
ICEIL
ICOFFA
II
IIA
IIBEG
IIEND
IINXT
ILEFT
IRIGHT
IROFFA
ITOP
JA
JJ
JJA
JJBEG
JJEND
JJNXT
LDA
LLD_
LSAME
M
M_
MB_
MBA
MP
MPA
MYCOL
MYDIST
MYROW
N
N_
NB_
NBA
NPCOL
NPROW
NQ
NQA
NUMROC
RSRC_
UPLO
WIDE
| |
