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
|
/* lapack.h -- C/C++ prototypes for (some) BLAS+LAPACK+wrapper routines */
/* $Header$ */
/*
* prerequisites:
* "cctk.h"
* "config.h" // for "integer" = Fortran integer
*/
#ifdef __cplusplus
extern "C"
{
#endif
/*
* ***** BLAS *****
*/
integer CCTK_FCALL
CCTK_FNAME(isamax)(const integer* N, const float SX[], const integer* incx);
integer CCTK_FCALL
CCTK_FNAME(idamax)(const integer* N, const double DX[], const integer* incx);
/*
* ***** LAPACK *****
*/
void CCTK_FCALL
CCTK_FNAME(sgesv)(const integer* N, const integer* NRHS,
float A[], const int* LDA,
integer IPIV[],
float B[], const integer* LDB, integer* info);
void CCTK_FCALL
CCTK_FNAME(dgesv)(const integer* N, const integer* NRHS,
double A[], const int* LDA,
integer IPIV[],
double B[], const integer* LDB, integer* info);
/*
* ***** wrappers (for passing character-string args) *****
*/
/* norm_int = 0 for infinity-norm, 1 for 1-norm */
void CCTK_FCALL
CCTK_FNAME(sgecon_wrapper)(const integer* norm_int,
const integer* N,
float A[], const integer* LDA,
const float* anorm, float* rcond,
float WORK[], integer IWORK[],
integer* info);
void CCTK_FCALL
CCTK_FNAME(dgecon_wrapper)(const integer* norm_int,
const integer* N,
double A[], const integer* LDA,
const double* anorm, double* rcond,
double WORK[], integer IWORK[],
integer* info);
#ifdef __cplusplus
} /* extern "C" */
#endif
|