ScaLAPACK 2.1  2.1
ScaLAPACK: Scalable Linear Algebra PACKage
pdbrddriver.f
Go to the documentation of this file.
1  PROGRAM pdbrddriver
2 *
3 * -- ScaLAPACK testing driver (version 1.7) --
4 * University of Tennessee, Knoxville, Oak Ridge National Laboratory,
5 * and University of California, Berkeley.
6 * March 13, 2000
7 *
8 * Purpose
9 * =======
10 *
11 * PDBRDDRIVER is the main test program for the DOUBLE PRECISION
12 * ScaLAPACK BRD (bidiagonal reduction) routines.
13 *
14 * The program must be driven by a short data file. An annotated
15 * example of a data file can be obtained by deleting the first 3
16 * characters from the following 13 lines:
17 * 'ScaLAPACK BRD computation input file'
18 * 'PVM machine'
19 * 'BRD.out' output file name
20 * 6 device out
21 * 3 number of problems sizes
22 * 16 20 18 values of M
23 * 16 18 20 values of N
24 * 3 number of NB's
25 * 2 3 5 values of NB
26 * 7 number of process grids (ordered pairs of P & Q)
27 * 1 2 1 4 2 3 8 values of P
28 * 1 2 4 1 3 2 1 values of Q
29 * 1.0 threshold
30 *
31 *
32 * Internal Parameters
33 * ===================
34 *
35 * TOTMEM INTEGER, default = 2000000
36 * TOTMEM is a machine-specific parameter indicating the
37 * maximum amount of available memory in bytes.
38 * The user should customize TOTMEM to his platform. Remember
39 * to leave room in memory for the operating system, the BLACS
40 * buffer, etc. For example, on a system with 8 MB of memory
41 * per process (e.g., one processor on an Intel iPSC/860), the
42 * parameters we use are TOTMEM=6200000 (leaving 1.8 MB for OS,
43 * code, BLACS buffer, etc). However, for PVM, we usually set
44 * TOTMEM = 2000000. Some experimenting with the maximum value
45 * of TOTMEM may be required.
46 *
47 * INTGSZ INTEGER, default = 4 bytes.
48 * DBLESZ INTEGER, default = 8 bytes.
49 * INTGSZ and DBLESZ indicate the length in bytes on the
50 * given platform for an integer and a double precision real.
51 * MEM DOUBLE PRECISION array, dimension ( TOTMEM / DBLESZ )
52 *
53 * All arrays used by SCALAPACK routines are allocated from
54 * this array and referenced by pointers. The integer IPA,
55 * for example, is a pointer to the starting element of MEM for
56 * the matrix A.
57 *
58 * =====================================================================
59 *
60 * .. Parameters ..
61  INTEGER block_cyclic_2d, csrc_, ctxt_, dlen_, dtype_,
62  $ lld_, mb_, m_, nb_, n_, rsrc_
63  parameter( block_cyclic_2d = 1, dlen_ = 9, dtype_ = 1,
64  $ ctxt_ = 2, m_ = 3, n_ = 4, mb_ = 5, nb_ = 6,
65  $ rsrc_ = 7, csrc_ = 8, lld_ = 9 )
66  INTEGER dblesz, memsiz, ntests, totmem
67  DOUBLE PRECISION padval
68  parameter( dblesz = 8, totmem = 2000000,
69  $ memsiz = totmem / dblesz, ntests = 20,
70  $ padval = -9923.0d+0 )
71 * ..
72 * .. Local Scalars ..
73  LOGICAL check
74  CHARACTER*6 passed
75  CHARACTER*80 outfile
76  INTEGER i, iam, iaseed, ictxt, imidpad, info, ipa, ipd,
77  $ ipe, ipostpad, iprepad, iptp, iptq, ipw, j, k,
78  $ kfail, kpass, kskip, ktests, lwork, m, maxmn,
79  $ minmn, mnp, mnq, mp, mycol, myrow, n, nb,
80  $ ndiag, ngrids, nmat, nnb, noffd, nout, npcol,
81  $ nprocs, nprow, nq, workbrd, worksiz
82  REAL thresh
83  DOUBLE PRECISION anorm, fresid, nops, tmflops
84 * ..
85 * .. Local Arrays ..
86  INTEGER desca( dlen_ ), ierr( 1 ), nbval( ntests ),
87  $ mval( ntests ), nval( ntests ),
88  $ pval( ntests ), qval( ntests )
89  DOUBLE PRECISION ctime( 1 ), mem( memsiz ), wtime( 1 )
90 * ..
91 * .. External Subroutines ..
92  EXTERNAL blacs_barrier, blacs_exit, blacs_get,
93  $ blacs_gridexit, blacs_gridinfo, blacs_gridinit,
94  $ blacs_pinfo, descinit, igsum2d, pdchekpad,
95  $ pdbrdinfo, pdfillpad, pdlafchk,
96  $ pdmatgen, pdgebdrv, pdgebrd, slboot,
97  $ slcombine, sltimer
98 * ..
99 * .. External Functions ..
100  INTEGER iceil, numroc
101  DOUBLE PRECISION pdlange
102  EXTERNAL iceil, numroc, pdlange
103 * ..
104 * .. Intrinsic Functions ..
105  INTRINSIC dble, max, min
106 * ..
107 * .. Data statements ..
108  DATA ktests, kpass, kfail, kskip / 4*0 /
109 * ..
110 * .. Executable Statements ..
111 *
112 * Get starting information
113 *
114  CALL blacs_pinfo( iam, nprocs )
115  iaseed = 100
116  CALL pdbrdinfo( outfile, nout, nmat, mval, ntests, nval, ntests,
117  $ nnb, nbval, ntests, ngrids, pval, ntests, qval,
118  $ ntests, thresh, mem, iam, nprocs )
119  check = ( thresh.GE.0.0e+0 )
120 *
121 * Print headings
122 *
123  IF( iam.EQ.0 ) THEN
124  WRITE( nout, fmt = * )
125  WRITE( nout, fmt = 9995 )
126  WRITE( nout, fmt = 9994 )
127  WRITE( nout, fmt = * )
128  END IF
129 *
130 * Loop over different process grids
131 *
132  DO 30 i = 1, ngrids
133 *
134  nprow = pval( i )
135  npcol = qval( i )
136 *
137 * Make sure grid information is correct
138 *
139  ierr( 1 ) = 0
140  IF( nprow.LT.1 ) THEN
141  IF( iam.EQ.0 )
142  $ WRITE( nout, fmt = 9999 ) 'GRID', 'nprow', nprow
143  ierr( 1 ) = 1
144  ELSE IF( npcol.LT.1 ) THEN
145  IF( iam.EQ.0 )
146  $ WRITE( nout, fmt = 9999 ) 'GRID', 'npcol', npcol
147  ierr( 1 ) = 1
148  ELSE IF( nprow*npcol.GT.nprocs ) THEN
149  IF( iam.EQ.0 )
150  $ WRITE( nout, fmt = 9998 ) nprow*npcol, nprocs
151  ierr( 1 ) = 1
152  END IF
153 *
154  IF( ierr( 1 ).GT.0 ) THEN
155  IF( iam.EQ.0 )
156  $ WRITE( nout, fmt = 9997 ) 'grid'
157  kskip = kskip + 1
158  GO TO 30
159  END IF
160 *
161 * Define process grid
162 *
163  CALL blacs_get( -1, 0, ictxt )
164  CALL blacs_gridinit( ictxt, 'Row-major', nprow, npcol )
165  CALL blacs_gridinfo( ictxt, nprow, npcol, myrow, mycol )
166 *
167  IF( myrow.GE.nprow .OR. mycol.GE.npcol )
168  $ GO TO 30
169 *
170 * Go to bottom of loop if this case doesn't use my process
171 *
172  DO 20 j = 1, nmat
173 *
174  m = mval( j )
175  n = nval( j )
176 *
177 * Make sure matrix information is correct
178 *
179  ierr( 1 ) = 0
180  IF( m.LT.1 ) THEN
181  IF( iam.EQ.0 )
182  $ WRITE( nout, fmt = 9999 ) 'MATRIX', 'M', m
183  ierr( 1 ) = 1
184  ELSE IF( n.LT.1 ) THEN
185  IF( iam.EQ.0 )
186  $ WRITE( nout, fmt = 9999 ) 'MATRIX', 'N', n
187  ierr( 1 ) = 1
188  END IF
189 *
190 * Make sure no one had error
191 *
192  CALL igsum2d( ictxt, 'All', ' ', 1, 1, ierr, 1, -1, 0 )
193 *
194  IF( ierr( 1 ).GT.0 ) THEN
195  IF( iam.EQ.0 )
196  $ WRITE( nout, fmt = 9997 ) 'matrix'
197  kskip = kskip + 1
198  GO TO 20
199  END IF
200 *
201 * Loop over different blocking sizes
202 *
203  DO 10 k = 1, nnb
204 *
205  nb = nbval( k )
206 *
207 * Make sure nb is legal
208 *
209  ierr( 1 ) = 0
210  IF( nb.LT.1 ) THEN
211  ierr( 1 ) = 1
212  IF( iam.EQ.0 )
213  $ WRITE( nout, fmt = 9999 ) 'NB', 'NB', nb
214  END IF
215 *
216 * Check all processes for an error
217 *
218  CALL igsum2d( ictxt, 'All', ' ', 1, 1, ierr, 1, -1, 0 )
219 *
220  IF( ierr( 1 ).GT.0 ) THEN
221  IF( iam.EQ.0 )
222  $ WRITE( nout, fmt = 9997 ) 'NB'
223  kskip = kskip + 1
224  GO TO 10
225  END IF
226 *
227 * Padding constants
228 *
229  mp = numroc( m, nb, myrow, 0, nprow )
230  nq = numroc( n, nb, mycol, 0, npcol )
231  mnp = numroc( min( m, n ), nb, myrow, 0, nprow )
232  mnq = numroc( min( m, n ), nb, mycol, 0, npcol )
233  IF( check ) THEN
234  iprepad = max( nb, mp )
235  imidpad = nb
236  ipostpad = max( nb, nq )
237  ELSE
238  iprepad = 0
239  imidpad = 0
240  ipostpad = 0
241  END IF
242 *
243 * Initialize the array descriptor for the matrix A
244 *
245  CALL descinit( desca, m, n, nb, nb, 0, 0, ictxt,
246  $ max( 1, mp )+imidpad, ierr( 1 ) )
247 *
248  CALL igsum2d( ictxt, 'All', ' ', 1, 1, ierr, 1, -1, 0 )
249 *
250  IF( ierr( 1 ).LT.0 ) THEN
251  IF( iam.EQ.0 )
252  $ WRITE( nout, fmt = 9997 ) 'descriptor'
253  kskip = kskip + 1
254  GO TO 10
255  END IF
256 *
257 * Assign pointers into MEM for SCALAPACK arrays, A is
258 * allocated starting at position MEM( IPREPAD+1 )
259 *
260  IF( m.GE.n ) THEN
261  ndiag = mnq
262  noffd = mnp
263  ELSE
264  ndiag = mnp
265  noffd = numroc( min( m, n )-1, nb, mycol, 0, npcol )
266  END IF
267 *
268  ipa = iprepad + 1
269  ipd = ipa + desca( lld_ )*nq + ipostpad + iprepad
270  ipe = ipd + ndiag + ipostpad + iprepad
271  iptq = ipe + noffd + ipostpad + iprepad
272  iptp = iptq + mnq + ipostpad + iprepad
273  ipw = iptp + mnp + ipostpad + iprepad
274 *
275 * Calculate the amount of workspace required for the
276 * reduction
277 *
278  lwork = nb*( mp+nq+1 ) + nq
279  workbrd = lwork + ipostpad
280  worksiz = workbrd
281 *
282 * Figure the amount of workspace required by the check
283 *
284  IF( check ) THEN
285  worksiz = max( lwork, 2*nb*( mp+nq+nb ) ) + ipostpad
286  END IF
287 *
288 * Check for adequate memory for problem size
289 *
290  ierr( 1 ) = 0
291  IF( ipw+worksiz.GT.memsiz ) THEN
292  IF( iam.EQ.0 )
293  $ WRITE( nout, fmt = 9996 ) 'Bidiagonal reduction',
294  $ ( ipw+worksiz )*dblesz
295  ierr( 1 ) = 1
296  END IF
297 *
298 * Check all processes for an error
299 *
300  CALL igsum2d( ictxt, 'All', ' ', 1, 1, ierr, 1, -1, 0 )
301 *
302  IF( ierr( 1 ).GT.0 ) THEN
303  IF( iam.EQ.0 )
304  $ WRITE( nout, fmt = 9997 ) 'MEMORY'
305  kskip = kskip + 1
306  GO TO 10
307  END IF
308 *
309 * Generate the matrix A
310 *
311  CALL pdmatgen( ictxt, 'No', 'No', desca( m_ ),
312  $ desca( n_ ), desca( mb_ ), desca( nb_ ),
313  $ mem( ipa ), desca( lld_ ), desca( rsrc_ ),
314  $ desca( csrc_ ), iaseed, 0, mp, 0, nq,
315  $ myrow, mycol, nprow, npcol )
316 *
317 * Need Infinity-norm of A for checking
318 *
319  IF( check ) THEN
320  CALL pdfillpad( ictxt, mp, nq, mem( ipa-iprepad ),
321  $ desca( lld_ ), iprepad, ipostpad,
322  $ padval )
323  CALL pdfillpad( ictxt, ndiag, 1, mem( ipd-iprepad ),
324  $ ndiag, iprepad, ipostpad, padval )
325  CALL pdfillpad( ictxt, noffd, 1, mem( ipe-iprepad ),
326  $ noffd, iprepad, ipostpad, padval )
327  CALL pdfillpad( ictxt, mnq, 1, mem( iptq-iprepad ),
328  $ mnq, iprepad, ipostpad, padval )
329  CALL pdfillpad( ictxt, mnp, 1, mem( iptp-iprepad ),
330  $ mnp, iprepad, ipostpad, padval )
331  CALL pdfillpad( ictxt, worksiz-ipostpad, 1,
332  $ mem( ipw-iprepad ), worksiz-ipostpad,
333  $ iprepad, ipostpad, padval )
334  anorm = pdlange( 'I', m, n, mem( ipa ), 1, 1, desca,
335  $ mem( ipw ) )
336  CALL pdchekpad( ictxt, 'PDLANGE', mp, nq,
337  $ mem( ipa-iprepad ), desca( lld_ ),
338  $ iprepad, ipostpad, padval )
339  CALL pdchekpad( ictxt, 'PDLANGE', worksiz-ipostpad,
340  $ 1, mem( ipw-iprepad ),
341  $ worksiz-ipostpad, iprepad, ipostpad,
342  $ padval )
343  CALL pdfillpad( ictxt, workbrd-ipostpad, 1,
344  $ mem( ipw-iprepad ), workbrd-ipostpad,
345  $ iprepad, ipostpad, padval )
346  END IF
347 *
348  CALL slboot()
349  CALL blacs_barrier( ictxt, 'All' )
350  CALL sltimer( 1 )
351 *
352 * Reduce to bidiagonal form
353 *
354  CALL pdgebrd( m, n, mem( ipa ), 1, 1, desca, mem( ipd ),
355  $ mem( ipe ), mem( iptq ), mem( iptp ),
356  $ mem( ipw ), lwork, info )
357 *
358  CALL sltimer( 1 )
359 *
360  IF( check ) THEN
361 *
362 * Check for memory overwrite
363 *
364  CALL pdchekpad( ictxt, 'PDGEBRD', mp, nq,
365  $ mem( ipa-iprepad ), desca( lld_ ),
366  $ iprepad, ipostpad, padval )
367  CALL pdchekpad( ictxt, 'PDGEBRD', ndiag, 1,
368  $ mem( ipd-iprepad ), ndiag, iprepad,
369  $ ipostpad, padval )
370  CALL pdchekpad( ictxt, 'PDGEBRD', noffd, 1,
371  $ mem( ipe-iprepad ), noffd, iprepad,
372  $ ipostpad, padval )
373  CALL pdchekpad( ictxt, 'PDGEBRD', mnq, 1,
374  $ mem( iptq-iprepad ), mnq, iprepad,
375  $ ipostpad, padval )
376  CALL pdchekpad( ictxt, 'PDGEBRD', mnp, 1,
377  $ mem( iptp-iprepad ), mnp, iprepad,
378  $ ipostpad, padval )
379  CALL pdchekpad( ictxt, 'PDGEBRD', workbrd-ipostpad,
380  $ 1, mem( ipw-iprepad ),
381  $ workbrd-ipostpad, iprepad,
382  $ ipostpad, padval )
383  CALL pdfillpad( ictxt, worksiz-ipostpad, 1,
384  $ mem( ipw-iprepad ), worksiz-ipostpad,
385  $ iprepad, ipostpad, padval )
386 *
387 * Compute fctres = ||A-Q*B*P|| / (||A|| * N * eps)
388 *
389  CALL pdgebdrv( m, n, mem( ipa ), 1, 1, desca,
390  $ mem( ipd ), mem( ipe ), mem( iptq ),
391  $ mem( iptp ), mem( ipw ), ierr( 1 ) )
392  CALL pdlafchk( 'No', 'No', m, n, mem( ipa ), 1, 1,
393  $ desca, iaseed, anorm, fresid,
394  $ mem( ipw ) )
395 *
396 * Check for memory overwrite
397 *
398  CALL pdchekpad( ictxt, 'PDGEBDRV', mp, nq,
399  $ mem( ipa-iprepad ), desca( lld_ ),
400  $ iprepad, ipostpad, padval )
401  CALL pdchekpad( ictxt, 'PDGEBDRV', ndiag, 1,
402  $ mem( ipd-iprepad ), ndiag, iprepad,
403  $ ipostpad, padval )
404  CALL pdchekpad( ictxt, 'PDGEBDRV', noffd, 1,
405  $ mem( ipe-iprepad ), noffd, iprepad,
406  $ ipostpad, padval )
407  CALL pdchekpad( ictxt, 'PDGEBDRV', worksiz-ipostpad,
408  $ 1, mem( ipw-iprepad ),
409  $ worksiz-ipostpad, iprepad,
410  $ ipostpad, padval )
411 *
412 * Test residual and detect NaN result
413 *
414  IF( fresid.LE.thresh .AND. fresid-fresid.EQ.0.0d+0
415  $ .AND. ierr( 1 ).EQ.0 ) THEN
416  kpass = kpass + 1
417  passed = 'PASSED'
418  ELSE
419  IF( myrow.EQ.0 .AND. mycol.EQ.0 )
420  $ WRITE( nout, fmt = 9986 ) fresid
421 *
422  kfail = kfail + 1
423  passed = 'FAILED'
424  END IF
425 *
426  IF( myrow.EQ.0 .AND. mycol.EQ.0 .AND. ierr( 1 ).NE.0 )
427  $ WRITE( nout, fmt = * )
428  $ 'D or E copies incorrect ...'
429  ELSE
430 *
431 * Don't perform the checking, only the timing operation
432 *
433  kpass = kpass + 1
434  fresid = fresid - fresid
435  passed = 'BYPASS'
436 *
437  END IF
438 *
439 * Gather maximum of all CPU and WALL clock timings
440 *
441  CALL slcombine( ictxt, 'All', '>', 'W', 1, 1, wtime )
442  CALL slcombine( ictxt, 'All', '>', 'C', 1, 1, ctime )
443 *
444 * Print results
445 *
446  IF( myrow.EQ.0 .AND. mycol.EQ.0 ) THEN
447 *
448 * BRD requires 8/3 N^3 floating point operations
449 *
450  maxmn = max( m, n )
451  minmn = min( m, n )
452  nops = 4.0d+0 * dble( minmn ) * dble( minmn ) *
453  $ ( dble( maxmn ) - dble( minmn ) / 3.0d+0 )
454  nops = nops / 1.0d+6
455 *
456 * Print WALL time
457 *
458  IF( wtime( 1 ).GT.0.0d+0 ) THEN
459  tmflops = nops / wtime( 1 )
460  ELSE
461  tmflops = 0.0d+0
462  END IF
463  IF( wtime( 1 ).GE.0.0d+0 )
464  $ WRITE( nout, fmt = 9993 ) 'WALL', m, n, nb, nprow,
465  $ npcol, wtime( 1 ), tmflops, fresid, passed
466 *
467 * Print CPU time
468 *
469  IF( ctime( 1 ).GT.0.0d+0 ) THEN
470  tmflops = nops / ctime( 1 )
471  ELSE
472  tmflops = 0.0d+0
473  END IF
474  IF( ctime( 1 ).GE.0.0d+0 )
475  $ WRITE( nout, fmt = 9993 ) 'CPU ', m, n, nb, nprow,
476  $ npcol, ctime( 1 ), tmflops, fresid, passed
477  END IF
478  10 CONTINUE
479  20 CONTINUE
480 *
481  CALL blacs_gridexit( ictxt )
482  30 CONTINUE
483 *
484 * Print ending messages and close output file
485 *
486  IF( iam.EQ.0 ) THEN
487  ktests = kpass + kfail + kskip
488  WRITE( nout, fmt = * )
489  WRITE( nout, fmt = 9992 ) ktests
490  IF( check ) THEN
491  WRITE( nout, fmt = 9991 ) kpass
492  WRITE( nout, fmt = 9989 ) kfail
493  ELSE
494  WRITE( nout, fmt = 9990 ) kpass
495  END IF
496  WRITE( nout, fmt = 9988 ) kskip
497  WRITE( nout, fmt = * )
498  WRITE( nout, fmt = * )
499  WRITE( nout, fmt = 9987 )
500  IF( nout.NE.6 .AND. nout.NE.0 ) CLOSE ( nout )
501  END IF
502 *
503  CALL blacs_exit( 0 )
504 *
505  9999 FORMAT( 'ILLEGAL ', a6, ': ', a5, ' = ', i3,
506  $ '; It should be at least 1' )
507  9998 FORMAT( 'ILLEGAL GRID: nprow*npcol = ', i4, '. It can be at most',
508  $ i4 )
509  9997 FORMAT( 'Bad ', a6, ' parameters: going on to next test case.' )
510  9996 FORMAT( 'Unable to perform ', a, ': need TOTMEM of at least',
511  $ i11 )
512  9995 FORMAT( 'TIME M N NB P Q BRD Time ',
513  $ ' MFLOPS Residual CHECK' )
514  9994 FORMAT( '---- ------ ------ --- ----- ----- --------- ',
515  $ '----------- -------- ------' )
516  9993 FORMAT( a4, 1x, i6, 1x, i6, 1x, i3, 1x, i5, 1x, i5, 1x, f9.2, 1x,
517  $ f11.2, 1x, f8.2, 1x, a6 )
518  9992 FORMAT( 'Finished', i4, ' tests, with the following results:' )
519  9991 FORMAT( i5, ' tests completed and passed residual checks.' )
520  9990 FORMAT( i5, ' tests completed without checking.' )
521  9989 FORMAT( i5, ' tests completed and failed residual checks.' )
522  9988 FORMAT( i5, ' tests skipped because of illegal input values.' )
523  9987 FORMAT( 'END OF TESTS.' )
524  9986 FORMAT( '||A - Q*B*P|| / (||A|| * N * eps) = ', g25.7 )
525 *
526  stop
527 *
528 * End of PDBRDDRIVER
529 *
530  END
max
#define max(A, B)
Definition: pcgemr.c:180
pdgebdrv
subroutine pdgebdrv(M, N, A, IA, JA, DESCA, D, E, TAUQ, TAUP, WORK, INFO)
Definition: pdgebdrv.f:3
pdchekpad
subroutine pdchekpad(ICTXT, MESS, M, N, A, LDA, IPRE, IPOST, CHKVAL)
Definition: pdchekpad.f:3
sltimer
subroutine sltimer(I)
Definition: sltimer.f:47
pdbrdinfo
subroutine pdbrdinfo(SUMMRY, NOUT, NMAT, MVAL, LDMVAL, NVAL, LDNVAL, NNB, NBVAL, LDNBVAL, NGRIDS, PVAL, LDPVAL, QVAL, LDQVAL, THRESH, WORK, IAM, NPROCS)
Definition: pdbrdinfo.f:5
pdgebrd
subroutine pdgebrd(M, N, A, IA, JA, DESCA, D, E, TAUQ, TAUP, WORK, LWORK, INFO)
Definition: pdgebrd.f:3
descinit
subroutine descinit(DESC, M, N, MB, NB, IRSRC, ICSRC, ICTXT, LLD, INFO)
Definition: descinit.f:3
slboot
subroutine slboot()
Definition: sltimer.f:2
pdmatgen
subroutine pdmatgen(ICTXT, AFORM, DIAG, M, N, MB, NB, A, LDA, IAROW, IACOL, ISEED, IROFF, IRNUM, ICOFF, ICNUM, MYROW, MYCOL, NPROW, NPCOL)
Definition: pdmatgen.f:4
pdlange
double precision function pdlange(NORM, M, N, A, IA, JA, DESCA, WORK)
Definition: pdlange.f:3
pdlafchk
subroutine pdlafchk(AFORM, DIAG, M, N, A, IA, JA, DESCA, IASEED, ANORM, FRESID, WORK)
Definition: pdlafchk.f:3
numroc
integer function numroc(N, NB, IPROC, ISRCPROC, NPROCS)
Definition: numroc.f:2
pdfillpad
subroutine pdfillpad(ICTXT, M, N, A, LDA, IPRE, IPOST, CHKVAL)
Definition: pdfillpad.f:2
pdbrddriver
program pdbrddriver
Definition: pdbrddriver.f:1
slcombine
subroutine slcombine(ICTXT, SCOPE, OP, TIMETYPE, N, IBEG, TIMES)
Definition: sltimer.f:267
min
#define min(A, B)
Definition: pcgemr.c:181
iceil
integer function iceil(INUM, IDENOM)
Definition: iceil.f:2