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
|
// Copyright (c) 2017-2023, University of Tennessee. All rights reserved.
// SPDX-License-Identifier: BSD-3-Clause
// This program is free software: you can redistribute it and/or modify it under
// the terms of the BSD 3-Clause license. See the accompanying LICENSE file.
#ifndef LAPACK_CONFIG_H
#define LAPACK_CONFIG_H
#include "blas/defines.h"
#include "lapack/defines.h"
#include <stdlib.h>
#if defined(BLAS_ILP64) && ! defined(LAPACK_ILP64)
#define LAPACK_ILP64
#endif
#ifndef lapack_int
#ifdef LAPACK_ILP64
typedef int64_t lapack_int;
#else
typedef int lapack_int;
#endif
/* #define so that #ifdef works. */
#define lapack_int lapack_int
#endif
#ifndef lapack_logical
#define lapack_logical lapack_int
#endif
/* f2c, hence MacOS Accelerate, returns double instead of float
* for sdot, slange, clange, etc. */
#if defined(BLAS_HAVE_ACCELERATE) || defined(BLAS_HAVE_F2C)
typedef double lapack_float_return;
#else
typedef float lapack_float_return;
#endif
/* --------------------------------------------------------------------------
* Complex types */
#if defined(LAPACK_COMPLEX_CUSTOM)
/* If user defines LAPACK_COMPLEX_CUSTOM, then the user must define:
* lapack_complex_float
* lapack_complex_double
* lapack_complex_float_real(z)
* lapack_complex_float_imag(z)
* lapack_complex_double_real(z)
* lapack_complex_double_imag(z)
*/
#elif defined(LAPACK_COMPLEX_STRUCTURE) || defined(_MSC_VER)
/* If user defines LAPACK_COMPLEX_STRUCTURE, then use a struct.
* Also use this for MSVC, as it has no C99 _Complex.
*/
typedef struct { float real, imag; } lapack_complex_float;
typedef struct { double real, imag; } lapack_complex_double;
#define lapack_complex_float_real(z) ((z).real)
#define lapack_complex_float_imag(z) ((z).imag)
#define lapack_complex_double_real(z) ((z).real)
#define lapack_complex_double_imag(z) ((z).imag)
#elif defined(LAPACK_COMPLEX_CPP) && defined(__cplusplus)
/* If user defines LAPACK_COMPLEX_CPP, then use C++ std::complex.
* This isn't compatible as a return type from extern C functions,
* so it may generate compiler warnings or errors. */
#include <complex>
typedef std::complex<float> lapack_complex_float;
typedef std::complex<double> lapack_complex_double;
#define lapack_complex_float_real(z) ((z).real())
#define lapack_complex_float_imag(z) ((z).imag())
#define lapack_complex_double_real(z) ((z).real())
#define lapack_complex_double_imag(z) ((z).imag())
#else
/* Otherwise, by default use C99 _Complex. */
#include <complex.h>
typedef float _Complex lapack_complex_float;
typedef double _Complex lapack_complex_double;
#define lapack_complex_float_real(z) (creal(z))
#define lapack_complex_float_imag(z) (cimag(z))
#define lapack_complex_double_real(z) (creal(z))
#define lapack_complex_double_imag(z) (cimag(z))
#endif
/* define so we can check later with ifdef */
#define lapack_complex_float lapack_complex_float
#define lapack_complex_double lapack_complex_double
#ifdef __cplusplus
extern "C" {
#endif /* __cplusplus */
lapack_complex_float lapack_make_complex_float( float re, float im );
lapack_complex_double lapack_make_complex_double( double re, double im );
#ifdef __cplusplus
}
#endif /* __cplusplus */
#ifndef LAPACK_malloc
#define LAPACK_malloc( size ) malloc( size )
#endif
#ifndef LAPACK_free
#define LAPACK_free( p ) free( p )
#endif
#endif /* LAPACK_CONFIG_H */
|
