MathFunctions.h 2.0 KB
Newer Older
Z
zhangjinchao01 已提交
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
/* Copyright (c) 2016 Baidu, Inc. All Rights Reserve.

Licensed under the Apache License, Version 2.0 (the "License");
you may not use this file except in compliance with the License.
You may obtain a copy of the License at

    http://www.apache.org/licenses/LICENSE-2.0

Unless required by applicable law or agreed to in writing, software
distributed under the License is distributed on an "AS IS" BASIS,
WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
See the License for the specific language governing permissions and
limitations under the License. */

#ifndef MATHFUNCTIONS_H_
#define MATHFUNCTIONS_H_

#ifdef PADDLE_USE_MKL
#include <mkl.h>
L
lzhao4ever 已提交
20
#include <mkl_lapacke.h>
Z
zhangjinchao01 已提交
21 22 23 24
#else
extern "C" {
#include <cblas.h>
}
L
lzhao4ever 已提交
25 26 27 28 29 30 31
#ifdef PADDLE_USE_ATLAS
extern "C" {
#include <clapack.h>
}
#else
#include <lapacke.h>
#endif
Z
zhangjinchao01 已提交
32 33
#endif

L
liaogang 已提交
34 35
#include <cmath>

Z
zhangjinchao01 已提交
36 37 38 39 40 41 42 43 44
namespace paddle {

template<class T>
void gemm(const CBLAS_TRANSPOSE transA, const CBLAS_TRANSPOSE transB,
          const int M, const int N, const int K,
          const T alpha, const T* A, const int lda,
          const T* B, const int ldb,
          const T beta, T* C, const int ldc);

L
lzhao4ever 已提交
45 46 47 48 49 50 51 52
template<class T>
int getrf(const CBLAS_ORDER Order, const int M, const int N,
          T *A, const int lda, int *ipiv);

template<class T>
int getri(const CBLAS_ORDER Order, const int N, T *A,
          const int lda, const int *ipiv);

Z
zhangjinchao01 已提交
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
template<class T>
void axpy(const int n, const T alpha, const T* x, T* y);

template<class T>
T dotProduct(const int n, const T* x, const T* y);

template<class T>
void vExp(const int n, const T* a, T* r);

template<class T>
void vPow(const int n, const T* a, const T b, T* r);

template<class T>
void vLog(const int n, const T* a, T* r);

template<class T>
void vAdd(const int n, const T* a, const T* b, T* r);

template<class T>
void vInvSqrt(const int n, const T* a, T* r);

template<class T>
void vLog1p(const int n, const T* a, T* r);

template<class T>
void vTanh(const int n, const T* a, T* r);

}  // namespace paddle

#endif  // MATHFUNCTIONS_H_