1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
|
/* ---------------------------------------------------------------------
*
* -- PBLAS routine (version 2.0) --
* University of Tennessee, Knoxville, Oak Ridge National Laboratory,
* and University of California, Berkeley.
* April 1, 1998
*
* ---------------------------------------------------------------------
*/
/*
* This file includes the standard C libraries, as well as system depen-
* dent include files. All PBLAS routines include this file.
*
* ---------------------------------------------------------------------
* Machine Specific PBLAS macros
* ---------------------------------------------------------------------
*/
#define _HAL_ 0
#define _T3D_ 1
#define _T3E_ 2
#ifdef T3D
#define _MACH_ _T3D_
#endif
#ifdef T3E
#define _MACH_ _T3E_
#endif
#ifndef _MACH_
#define _MACH_ _HAL_
#endif
/*
* CBRATIO is the ratio of the transfer cost per element for the combine
* sum to one process and the broadcast operation. This value is used
* within the Level 3 PBLAS routines to decide on which algorithm to se-
* lect.
*/
#define CBRATIO 1.3
/*
* ---------------------------------------------------------------------
* Include files
* ---------------------------------------------------------------------
*/
#include <stdio.h>
#include <stdlib.h>
#ifdef __STDC__
#include <stdarg.h>
#else
#include <varargs.h>
#endif
#if( ( _MACH_ == _T3D_ ) || ( _MACH_ == _T3E_ ) )
#include <fortran.h>
#endif
/*
* ---------------------------------------------------------------------
* FORTRAN <-> C interface
* ---------------------------------------------------------------------
*
* These macros identifies how the PBLAS will be called as follows:
*
* _F2C_ADD_: the FORTRAN compiler expects the name of C functions to be
* in all lower case and to have an underscore postfixed it (Suns, Intel
* compilers expect this).
*
* _F2C_NOCHANGE: the FORTRAN compiler expects the name of C functions
* to be in all lower case (IBM RS6K compilers do this).
*
* _F2C_UPCASE: the FORTRAN compiler expects the name of C functions
* to be in all upcase. (Cray compilers expect this).
*
* _F2C_F77ISF2C: the FORTRAN compiler in use is f2c, a FORTRAN to C
* converter.
*/
#define _F2C_ADD_ 0
#define _F2C_NOCHANGE 1
#define _F2C_UPCASE 2
#define _F2C_F77ISF2C 3
#ifdef UpCase
#define _F2C_CALL_ _F2C_UPCASE
#endif
#ifdef NoChange
#define _F2C_CALL_ _F2C_NOCHANGE
#endif
#ifdef Add_
#define _F2C_CALL_ _F2C_ADD_
#endif
#ifdef f77IsF2C
#define _F2C_CALL_ _F2C_F77ISF2C
#endif
#ifndef _F2C_CALL_
#define _F2C_CALL_ _F2C_ADD_
#endif
/*
* ---------------------------------------------------------------------
* TYPE DEFINITIONS AND CONVERSION UTILITIES
* ---------------------------------------------------------------------
*/
#ifndef Int
#define Int int
#endif
#if( ( _MACH_ == _T3D_ ) || ( _MACH_ == _T3E_ ) )
#define float double
/* Type of character argument in a FORTRAN call */
#define F_CHAR_T _fcd
/* Character conversion utilities */
#define F2C_CHAR(a) ( _fcdtocp( (a) ) )
#define C2F_CHAR(a) ( _cptofcd( (a), 1 ) )
/* Type of FORTRAN functions */
#define F_VOID_FCT void fortran /* Subroutine */
#define F_INTG_FCT Int fortran /* INTEGER function */
#else /* Type of character argument in a FORTRAN call */
typedef char * F_CHAR_T;
/* Character conversion utilities */
#define F2C_CHAR(a) (a)
#define C2F_CHAR(a) (a)
/* Type of FORTRAN functions */
#define F_VOID_FCT void /* Subroutine */
#define F_INTG_FCT Int /* INTEGER function */
#endif
/*
* ----------------------------------------------------------------------
* #typedef definitions
* ---------------------------------------------------------------------
*/
typedef float cmplx [2];
typedef double cmplx16[2];
#define REAL_PART 0
#define IMAG_PART 1
#ifdef __STDC__
typedef void (*GESD2D_T) ( Int, Int, Int,
char *, Int, Int,
Int );
typedef void (*GERV2D_T) ( Int, Int, Int,
char *, Int, Int,
Int );
typedef void (*GEBS2D_T) ( Int, char *, char *,
Int, Int, char *,
Int );
typedef void (*GEBR2D_T) ( Int, char *, char *,
Int, Int, char *,
Int, Int, Int );
typedef void (*GSUM2D_T) ( Int, char *, char *,
Int, Int, char *,
Int, Int, Int );
typedef F_VOID_FCT (*MMADD_T) ( Int *, Int *, char *,
char *, Int *, char *,
char *, Int * );
typedef F_VOID_FCT (*MMSHFT_T) ( Int *, Int *, Int *,
char *, Int * );
typedef F_VOID_FCT (*VVDOT_T) ( Int *, char *, char *,
Int *, char *, Int * );
typedef F_VOID_FCT (*VVSET_T) ( Int *, char *, char *,
Int * );
typedef F_VOID_FCT (*TZPAD_T) ( F_CHAR_T, F_CHAR_T, Int *,
Int *, Int *, char *,
char *, char *, Int * );
typedef F_VOID_FCT (*TZPADCPY_T) ( F_CHAR_T, F_CHAR_T, Int *,
Int *, Int *, char *,
Int *, char *, Int * );
typedef F_VOID_FCT (*TZSET_T) ( F_CHAR_T, Int *, Int *,
Int *, char *, char *,
char *, Int * );
typedef F_VOID_FCT (*TZSCAL_T) ( F_CHAR_T, Int *, Int *,
Int *, char *, char *,
Int * );
typedef F_VOID_FCT (*AXPY_T) ( Int *, char *, char *,
Int *, char *, Int * );
typedef F_VOID_FCT (*COPY_T) ( Int *, char *, Int *,
char *, Int * );
typedef F_VOID_FCT (*SWAP_T) ( Int *, char *, Int *,
char *, Int * );
typedef F_VOID_FCT (*GEMV_T) ( F_CHAR_T, Int *, Int *,
char *, char *, Int *,
char *, Int *, char *,
char *, Int * );
typedef F_VOID_FCT (*AGEMV_T) ( F_CHAR_T, Int *, Int *,
char *, char *, Int *,
char *, Int *, char *,
char *, Int * );
typedef F_VOID_FCT (*SYMV_T) ( F_CHAR_T, Int *, char *,
char *, Int *, char *,
Int *, char *, char *,
Int * );
typedef F_VOID_FCT (*ASYMV_T) ( F_CHAR_T, Int *, char *,
char *, Int *, char *,
Int *, char *, char *,
Int * );
typedef F_VOID_FCT (*HEMV_T) ( F_CHAR_T, Int *, char *,
char *, Int *, char *,
Int *, char *, char *,
Int * );
typedef F_VOID_FCT (*AHEMV_T) ( F_CHAR_T, Int *, char *,
char *, Int *, char *,
Int *, char *, char *,
Int * );
typedef F_VOID_FCT (*TRMV_T) ( F_CHAR_T, F_CHAR_T, F_CHAR_T,
Int *, char *, Int *,
char *, Int * );
typedef F_VOID_FCT (*ATRMV_T) ( F_CHAR_T, F_CHAR_T, F_CHAR_T,
Int *, char *, char *,
Int *, char *, Int *,
char *, char *, Int * );
typedef F_VOID_FCT (*TRSV_T) ( F_CHAR_T, F_CHAR_T, F_CHAR_T,
Int *, char *, Int *,
char *, Int * );
typedef F_VOID_FCT (*GERC_T) ( Int *, Int *, char *,
char *, Int *, char *,
Int *, char *, Int * );
typedef F_VOID_FCT (*GERU_T) ( Int *, Int *, char *,
char *, Int *, char *,
Int *, char *, Int * );
typedef F_VOID_FCT (*SYR_T) ( F_CHAR_T, Int *, char *,
char *, Int *, char *,
Int * );
typedef F_VOID_FCT (*HER_T) ( F_CHAR_T, Int *, char *,
char *, Int *, char *,
Int * );
typedef F_VOID_FCT (*SYR2_T) ( F_CHAR_T, Int *, char *,
char *, Int *, char *,
Int *, char *, Int * );
typedef F_VOID_FCT (*HER2_T) ( F_CHAR_T, Int *, char *,
char *, Int *, char *,
Int *, char *, Int * );
typedef F_VOID_FCT (*GEMM_T) ( F_CHAR_T, F_CHAR_T, Int *,
Int *, Int *, char *,
char *, Int *, char *,
Int *, char *, char *,
Int * );
typedef F_VOID_FCT (*SYMM_T) ( F_CHAR_T, F_CHAR_T, Int *,
Int *, char *, char *,
Int *, char *, Int *,
char *, char *, Int * );
typedef F_VOID_FCT (*HEMM_T) ( F_CHAR_T, F_CHAR_T, Int *,
Int *, char *, char *,
Int *, char *, Int *,
char *, char *, Int * );
typedef F_VOID_FCT (*SYRK_T) ( F_CHAR_T, F_CHAR_T, Int *,
Int *, char *, char *,
Int *, char *, char *,
Int * );
typedef F_VOID_FCT (*HERK_T) ( F_CHAR_T, F_CHAR_T, Int *,
Int *, char *, char *,
Int *, char *, char *,
Int * );
typedef F_VOID_FCT (*SYR2K_T) ( F_CHAR_T, F_CHAR_T, Int *,
Int *, char *, char *,
Int *, char *, Int *,
char *, char *, Int * );
typedef F_VOID_FCT (*HER2K_T) ( F_CHAR_T, F_CHAR_T, Int *,
Int *, char *, char *,
Int *, char *, Int *,
char *, char *, Int * );
typedef F_VOID_FCT (*TRMM_T) ( F_CHAR_T, F_CHAR_T, F_CHAR_T,
F_CHAR_T, Int *, Int *,
char *, char *, Int *,
char *, Int * );
typedef F_VOID_FCT (*TRSM_T) ( F_CHAR_T, F_CHAR_T, F_CHAR_T,
F_CHAR_T, Int *, Int *,
char *, char *, Int *,
char *, Int * );
#else
typedef void (*GESD2D_T) ();
typedef void (*GERV2D_T) ();
typedef void (*GEBS2D_T) ();
typedef void (*GEBR2D_T) ();
typedef void (*GSUM2D_T) ();
typedef F_VOID_FCT (*MMADD_T) ();
typedef F_VOID_FCT (*MMSHFT_T) ();
typedef F_VOID_FCT (*VVDOT_T) ();
typedef F_VOID_FCT (*VVSET_T) ();
typedef F_VOID_FCT (*TZPAD_T) ();
typedef F_VOID_FCT (*TZPADCPY_T) ();
typedef F_VOID_FCT (*TZSET_T) ();
typedef F_VOID_FCT (*TZSCAL_T) ();
typedef F_VOID_FCT (*AXPY_T) ();
typedef F_VOID_FCT (*COPY_T) ();
typedef F_VOID_FCT (*SWAP_T) ();
typedef F_VOID_FCT (*GEMV_T) ();
typedef F_VOID_FCT (*AGEMV_T) ();
typedef F_VOID_FCT (*SYMV_T) ();
typedef F_VOID_FCT (*ASYMV_T) ();
typedef F_VOID_FCT (*HEMV_T) ();
typedef F_VOID_FCT (*AHEMV_T) ();
typedef F_VOID_FCT (*TRMV_T) ();
typedef F_VOID_FCT (*ATRMV_T) ();
typedef F_VOID_FCT (*TRSV_T) ();
typedef F_VOID_FCT (*GERC_T) ();
typedef F_VOID_FCT (*GERU_T) ();
typedef F_VOID_FCT (*SYR_T) ();
typedef F_VOID_FCT (*HER_T) ();
typedef F_VOID_FCT (*SYR2_T) ();
typedef F_VOID_FCT (*HER2_T) ();
typedef F_VOID_FCT (*GEMM_T) ();
typedef F_VOID_FCT (*SYMM_T) ();
typedef F_VOID_FCT (*HEMM_T) ();
typedef F_VOID_FCT (*SYRK_T) ();
typedef F_VOID_FCT (*HERK_T) ();
typedef F_VOID_FCT (*SYR2K_T) ();
typedef F_VOID_FCT (*HER2K_T) ();
typedef F_VOID_FCT (*TRMM_T) ();
typedef F_VOID_FCT (*TRSM_T) ();
#endif
typedef struct
{
char type; /* Encoding of the data type */
Int usiz; /* length in bytes of elementary data type */
Int size; /* length in bytes of data type */
char * zero,
* one,
* negone; /* pointers to contants of correct type */
GESD2D_T Cgesd2d; /* BLACS functions */
GERV2D_T Cgerv2d;
GEBS2D_T Cgebs2d;
GEBR2D_T Cgebr2d;
GSUM2D_T Cgsum2d;
MMADD_T Fmmadd; /* Addition functions */
MMADD_T Fmmcadd;
MMADD_T Fmmtadd;
MMADD_T Fmmtcadd;
MMADD_T Fmmdda;
MMADD_T Fmmddac;
MMADD_T Fmmddat;
MMADD_T Fmmddact;
MMSHFT_T Fcshft; /* Shift functions */
MMSHFT_T Frshft;
VVDOT_T Fvvdotu; /* Dot functions */
VVDOT_T Fvvdotc;
TZPAD_T Ftzpad; /* Array pad function */
TZPADCPY_T Ftzpadcpy;
VVSET_T Fset;
TZSCAL_T Ftzscal; /* Scaling functions */
TZSCAL_T Fhescal;
TZSCAL_T Ftzcnjg;
AXPY_T Faxpy; /* Level 1 BLAS */
COPY_T Fcopy;
SWAP_T Fswap;
GEMV_T Fgemv; /* Level 2 BLAS */
SYMV_T Fsymv;
HEMV_T Fhemv;
TRMV_T Ftrmv;
TRSV_T Ftrsv;
AGEMV_T Fagemv;
ASYMV_T Fasymv;
AHEMV_T Fahemv;
ATRMV_T Fatrmv;
GERC_T Fgerc;
GERU_T Fgeru;
SYR_T Fsyr;
HER_T Fher;
SYR2_T Fsyr2;
HER2_T Fher2;
GEMM_T Fgemm; /* Level 3 BLAS */
SYMM_T Fsymm;
HEMM_T Fhemm;
SYRK_T Fsyrk;
HERK_T Fherk;
SYR2K_T Fsyr2k;
HER2K_T Fher2k;
TRMM_T Ftrmm;
TRSM_T Ftrsm;
} PBTYP_T;
#ifdef __STDC__
typedef void (*TZSYR_T) ( PBTYP_T *, char *, Int,
Int, Int, Int,
char *, char *, Int,
char *, Int, char *,
Int );
typedef void (*TZSYR2_T) ( PBTYP_T *, char *, Int,
Int, Int, Int,
char *, char *, Int,
char *, Int, char *,
Int, char *, Int,
char *, Int );
typedef void (*TZTRM_T) ( PBTYP_T *, char *, char *,
char *, char *, Int,
Int, Int, Int,
char *, char *, Int,
char *, Int, char *,
Int );
typedef void (*TZSYM_T) ( PBTYP_T *, char *, char *,
Int, Int, Int,
Int, char *, char *,
Int, char *, Int,
char *, Int, char *,
Int, char *, Int );
#else
typedef void (*TZSYR_T) ();
typedef void (*TZSYR2_T) ();
typedef void (*TZTRM_T) ();
typedef void (*TZSYM_T) ();
#endif
typedef struct
{
Int offd; /* Global diagonal offset */
Int lcmt00; /* LCM value of first block */
Int mp; /* Local number of rows */
Int imb1; /* Size of first row block (global) */
Int imbloc; /* Size of first local row block */
Int mb; /* Row block size */
Int lmbloc; /* Size of last local row block */
Int mblks; /* Number of local row blocks */
Int iupp; /* LCM row bound for first diagonal block */
Int upp; /* LCM row bound for diagonal block */
Int prow; /* Relative row process coordinate */
Int nprow; /* Number of process rows */
Int nq; /* Local number of columns */
Int inb1; /* Size of first column block (global) */
Int inbloc; /* Size of first local column block */
Int nb; /* Column block size */
Int lnbloc; /* Size of last local column block */
Int nblks; /* Number of local column blocks */
Int ilow; /* LCM column bound for first diagonal block */
Int low; /* LCM column bound for diagonal block */
Int pcol; /* Relative column process coordinate */
Int npcol; /* Number of process columns */
Int lcmb; /* Least common multiple of nprow * mb and npcol * nb */
} PB_VM_T;
/*
* ---------------------------------------------------------------------
* #define macro constants
* ---------------------------------------------------------------------
*/
#define INT 'I' /* type identifiers */
#define SREAL 'S'
#define DREAL 'D'
#define SCPLX 'C'
#define DCPLX 'Z'
#define crot_ CROT
|
