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
|
#include "./pblas.h"
#ifdef __STDC__
typedef void (*CPYPTR)(Int, Int, float *, Int, float *, Int);
#define SLVOID void
#else
typedef void (*CPYPTR)();
#define SLVOID char
#endif
#define ErrPrnt fprintf(stderr, "line %d of file %s\n",__LINE__, __FILE__);
#define Mdescset(desc, m, n, mb, nb, rsrc, csrc, ictxt, lld) \
{ \
(desc)[DT_] = BLOCK_CYCLIC_2D; \
(desc)[CTXT_] = (ictxt); \
(desc)[M_] = (m); \
(desc)[N_] = (n); \
(desc)[MB_] = (mb); \
(desc)[NB_] = (nb); \
(desc)[RSRC_] = (rsrc); \
(desc)[CSRC_] = (csrc); \
(desc)[LLD_] = (lld); \
}
#define MCindxg2p(IG, nb, srcproc, nprocs) \
( ((srcproc) + (IG)/(nb)) % nprocs )
typedef struct {double r, i;} DCOMPLEX;
typedef struct {float r, i;} SCOMPLEX;
#define Mmalloc(M_ptr, M_type, M_elt, M_i, M_ctxt) \
{ \
void pberror_(); \
(M_ptr) = ( M_type * ) malloc((M_elt)*(sizeof(M_type))); \
if (!(M_ptr)) \
{ \
if ((M_elt) > 0) \
{ \
(M_i) = 1; \
fprintf(stderr, "Not enough memory on line %d of file %s!!\n", \
__LINE__, __FILE__); \
pberror_(&(M_ctxt), __FILE__, &(M_i)); \
} \
} \
}
|
