Narayan Kamath | c981c48 | 2012-11-02 10:59:05 +0000 | [diff] [blame] | 1 | /* |
| 2 | Copyright (c) 2011, Intel Corporation. All rights reserved. |
| 3 | |
| 4 | Redistribution and use in source and binary forms, with or without modification, |
| 5 | are permitted provided that the following conditions are met: |
| 6 | |
| 7 | * Redistributions of source code must retain the above copyright notice, this |
| 8 | list of conditions and the following disclaimer. |
| 9 | * Redistributions in binary form must reproduce the above copyright notice, |
| 10 | this list of conditions and the following disclaimer in the documentation |
| 11 | and/or other materials provided with the distribution. |
| 12 | * Neither the name of Intel Corporation nor the names of its contributors may |
| 13 | be used to endorse or promote products derived from this software without |
| 14 | specific prior written permission. |
| 15 | |
| 16 | THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND |
| 17 | ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED |
| 18 | WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE |
| 19 | DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR |
| 20 | ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES |
| 21 | (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; |
| 22 | LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON |
| 23 | ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT |
| 24 | (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS |
| 25 | SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. |
| 26 | |
| 27 | ******************************************************************************** |
| 28 | * Content : Eigen bindings to Intel(R) MKL |
| 29 | * Include file with common MKL declarations |
| 30 | ******************************************************************************** |
| 31 | */ |
| 32 | |
| 33 | #ifndef EIGEN_MKL_SUPPORT_H |
| 34 | #define EIGEN_MKL_SUPPORT_H |
| 35 | |
| 36 | #ifdef EIGEN_USE_MKL_ALL |
| 37 | #ifndef EIGEN_USE_BLAS |
| 38 | #define EIGEN_USE_BLAS |
| 39 | #endif |
| 40 | #ifndef EIGEN_USE_LAPACKE |
| 41 | #define EIGEN_USE_LAPACKE |
| 42 | #endif |
| 43 | #ifndef EIGEN_USE_MKL_VML |
| 44 | #define EIGEN_USE_MKL_VML |
| 45 | #endif |
| 46 | #endif |
| 47 | |
| 48 | #ifdef EIGEN_USE_LAPACKE_STRICT |
| 49 | #define EIGEN_USE_LAPACKE |
| 50 | #endif |
| 51 | |
| 52 | #if defined(EIGEN_USE_BLAS) || defined(EIGEN_USE_LAPACKE) || defined(EIGEN_USE_MKL_VML) |
| 53 | #define EIGEN_USE_MKL |
| 54 | #endif |
| 55 | |
| 56 | #if defined EIGEN_USE_MKL |
Carlos Hernandez | 7faaa9f | 2014-08-05 17:53:32 -0700 | [diff] [blame] | 57 | # include <mkl.h> |
| 58 | /*Check IMKL version for compatibility: < 10.3 is not usable with Eigen*/ |
| 59 | # ifndef INTEL_MKL_VERSION |
| 60 | # undef EIGEN_USE_MKL /* INTEL_MKL_VERSION is not even defined on older versions */ |
| 61 | # elif INTEL_MKL_VERSION < 100305 /* the intel-mkl-103-release-notes say this was when the lapacke.h interface was added*/ |
| 62 | # undef EIGEN_USE_MKL |
| 63 | # endif |
| 64 | # ifndef EIGEN_USE_MKL |
| 65 | /*If the MKL version is too old, undef everything*/ |
| 66 | # undef EIGEN_USE_MKL_ALL |
| 67 | # undef EIGEN_USE_BLAS |
| 68 | # undef EIGEN_USE_LAPACKE |
| 69 | # undef EIGEN_USE_MKL_VML |
| 70 | # undef EIGEN_USE_LAPACKE_STRICT |
| 71 | # undef EIGEN_USE_LAPACKE |
| 72 | # endif |
| 73 | #endif |
Narayan Kamath | c981c48 | 2012-11-02 10:59:05 +0000 | [diff] [blame] | 74 | |
Carlos Hernandez | 7faaa9f | 2014-08-05 17:53:32 -0700 | [diff] [blame] | 75 | #if defined EIGEN_USE_MKL |
Narayan Kamath | c981c48 | 2012-11-02 10:59:05 +0000 | [diff] [blame] | 76 | #include <mkl_lapacke.h> |
| 77 | #define EIGEN_MKL_VML_THRESHOLD 128 |
| 78 | |
| 79 | namespace Eigen { |
| 80 | |
| 81 | typedef std::complex<double> dcomplex; |
| 82 | typedef std::complex<float> scomplex; |
| 83 | |
| 84 | namespace internal { |
| 85 | |
| 86 | template<typename MKLType, typename EigenType> |
| 87 | static inline void assign_scalar_eig2mkl(MKLType& mklScalar, const EigenType& eigenScalar) { |
| 88 | mklScalar=eigenScalar; |
| 89 | } |
| 90 | |
| 91 | template<typename MKLType, typename EigenType> |
| 92 | static inline void assign_conj_scalar_eig2mkl(MKLType& mklScalar, const EigenType& eigenScalar) { |
| 93 | mklScalar=eigenScalar; |
| 94 | } |
| 95 | |
| 96 | template <> |
| 97 | inline void assign_scalar_eig2mkl<MKL_Complex16,dcomplex>(MKL_Complex16& mklScalar, const dcomplex& eigenScalar) { |
| 98 | mklScalar.real=eigenScalar.real(); |
| 99 | mklScalar.imag=eigenScalar.imag(); |
| 100 | } |
| 101 | |
| 102 | template <> |
| 103 | inline void assign_scalar_eig2mkl<MKL_Complex8,scomplex>(MKL_Complex8& mklScalar, const scomplex& eigenScalar) { |
| 104 | mklScalar.real=eigenScalar.real(); |
| 105 | mklScalar.imag=eigenScalar.imag(); |
| 106 | } |
| 107 | |
| 108 | template <> |
| 109 | inline void assign_conj_scalar_eig2mkl<MKL_Complex16,dcomplex>(MKL_Complex16& mklScalar, const dcomplex& eigenScalar) { |
| 110 | mklScalar.real=eigenScalar.real(); |
| 111 | mklScalar.imag=-eigenScalar.imag(); |
| 112 | } |
| 113 | |
| 114 | template <> |
| 115 | inline void assign_conj_scalar_eig2mkl<MKL_Complex8,scomplex>(MKL_Complex8& mklScalar, const scomplex& eigenScalar) { |
| 116 | mklScalar.real=eigenScalar.real(); |
| 117 | mklScalar.imag=-eigenScalar.imag(); |
| 118 | } |
| 119 | |
| 120 | } // end namespace internal |
| 121 | |
| 122 | } // end namespace Eigen |
| 123 | |
| 124 | #endif |
| 125 | |
| 126 | #endif // EIGEN_MKL_SUPPORT_H |