forked from apache/mxnet
-
Notifications
You must be signed in to change notification settings - Fork 0
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Revert "Fix and optimize handling of vectorized memory accesses (apac…
- Loading branch information
1 parent
9d44086
commit dc8eb60
Showing
19 changed files
with
444 additions
and
1,344 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,143 @@ | ||
/*! | ||
* Copyright (c) 2017 by Contributors | ||
* \file half2.h | ||
* \brief definition of vector float16, half2 type. | ||
* | ||
* \author Antti-Pekka Hynninen | ||
*/ | ||
#ifndef MSHADOW_HALF2_H_ | ||
#define MSHADOW_HALF2_H_ | ||
|
||
#if (defined(__CUDACC__) && __CUDA_ARCH__ >= 530 && MSHADOW_USE_CUDA && CUDA_VERSION >= 7050) | ||
#define MSHADOW_CUDA_HALF2 1 | ||
#include <cuda_fp16.h> | ||
#else | ||
#define MSHADOW_CUDA_HALF2 0 | ||
#endif | ||
|
||
#include<math.h> | ||
|
||
/*! \brief namespace for mshadow */ | ||
namespace mshadow { | ||
/* \brief name space for host/device portable half-precision floats */ | ||
namespace half { | ||
|
||
#define MSHADOW_HALF2_ASSIGNOP(AOP, OP) \ | ||
template<typename T> \ | ||
MSHADOW_XINLINE half2_t operator AOP (const T& a) { \ | ||
return *this = half2_t(*this OP a); /* NOLINT(*)*/ \ | ||
} \ | ||
|
||
class MSHADOW_ALIGNED(4) half2_t { | ||
public: | ||
#if MSHADOW_CUDA_HALF2 | ||
half2 half2_; | ||
#else | ||
half_t half_t2[2]; | ||
#endif | ||
|
||
MSHADOW_XINLINE half2_t() {} | ||
|
||
#if MSHADOW_CUDA_HALF2 | ||
MSHADOW_XINLINE explicit half2_t(half2 a) : half2_(a) {} | ||
#else | ||
MSHADOW_XINLINE explicit half2_t(half_t a, half_t b) { | ||
half_t2[0] = a; | ||
half_t2[1] = b; | ||
} | ||
#endif | ||
|
||
MSHADOW_XINLINE explicit half2_t(int a) { | ||
#if MSHADOW_CUDA_HALF2 | ||
half2_ = __half2half2(__int2half_rz(a)); | ||
#else | ||
half_t2[0] = (half_t)a; | ||
half_t2[1] = (half_t)a; | ||
#endif | ||
} | ||
|
||
MSHADOW_XINLINE half2_t operator+() { | ||
return *this; | ||
} | ||
|
||
MSHADOW_XINLINE half2_t operator-() { | ||
#if MSHADOW_CUDA_HALF2 | ||
return half2_t(__hneg2(half2_)); | ||
#else | ||
return half2_t(-half_t2[0], -half_t2[1]); | ||
#endif | ||
} | ||
|
||
MSHADOW_XINLINE half2_t operator=(const half2_t& a) { | ||
#if MSHADOW_CUDA_HALF2 | ||
half2_ = a.half2_; | ||
#else | ||
half_t2[0] = a.half_t2[0]; | ||
half_t2[1] = a.half_t2[1]; | ||
#endif | ||
return a; | ||
} | ||
|
||
MSHADOW_HALF2_ASSIGNOP(+=, +) | ||
MSHADOW_HALF2_ASSIGNOP(-=, -) | ||
MSHADOW_HALF2_ASSIGNOP(*=, *) | ||
MSHADOW_HALF2_ASSIGNOP(/=, /) | ||
}; | ||
|
||
/*! \brief overloaded + operator for half2_t */ | ||
MSHADOW_XINLINE half2_t operator+(half2_t a, half2_t b) { | ||
#if MSHADOW_CUDA_HALF2 | ||
return half2_t(__floats2half2_rn(__low2float(a.half2_) + __low2float(b.half2_), | ||
__high2float(a.half2_) + __high2float(b.half2_))); | ||
#else | ||
return half2_t(a.half_t2[0] + b.half_t2[0], a.half_t2[1] + b.half_t2[1]); | ||
#endif | ||
} | ||
/*! \brief overloaded - operator for half2_t */ | ||
MSHADOW_XINLINE half2_t operator-(half2_t a, half2_t b) { | ||
#if MSHADOW_CUDA_HALF2 | ||
return half2_t(__floats2half2_rn(__low2float(a.half2_) - __low2float(b.half2_), | ||
__high2float(a.half2_) - __high2float(b.half2_))); | ||
#else | ||
return half2_t(a.half_t2[0] - b.half_t2[0], a.half_t2[1] - b.half_t2[1]); | ||
#endif | ||
} | ||
/*! \brief overloaded * operator for half2_t */ | ||
MSHADOW_XINLINE half2_t operator*(half2_t a, half2_t b) { | ||
#if MSHADOW_CUDA_HALF2 | ||
return half2_t(__floats2half2_rn(__low2float(a.half2_) * __low2float(b.half2_), | ||
__high2float(a.half2_) * __high2float(b.half2_))); | ||
#else | ||
return half2_t(a.half_t2[0] * b.half_t2[0], a.half_t2[1] * b.half_t2[1]); | ||
#endif | ||
} | ||
/*! \brief overloaded / operator for half2_t */ | ||
MSHADOW_XINLINE half2_t operator/(half2_t a, half2_t b) { | ||
#if MSHADOW_CUDA_HALF2 | ||
return half2_t(__floats2half2_rn(__low2float(a.half2_) / __low2float(b.half2_), | ||
__high2float(a.half2_) / __high2float(b.half2_))); | ||
#else | ||
return half2_t(a.half_t2[0] / b.half_t2[0], a.half_t2[1] / b.half_t2[1]); | ||
#endif | ||
} | ||
/*! \brief overloaded % operator for half2_t */ | ||
MSHADOW_XINLINE half2_t operator%(half2_t a, half2_t b) { | ||
#if MSHADOW_CUDA_HALF2 | ||
return half2_t(__floats2half2_rn(::fmod(__low2float(a.half2_), __low2float(b.half2_)), | ||
::fmod(__high2float(a.half2_), __high2float(b.half2_)))); | ||
#else | ||
return half2_t(::fmod(a.half_t2[0], b.half_t2[0]), ::fmod(a.half_t2[1], b.half_t2[1])); | ||
#endif | ||
} | ||
/*! \brief overloaded == operator for half2_t */ | ||
MSHADOW_XINLINE bool operator==(half2_t a, half2_t b) { | ||
#if MSHADOW_CUDA_HALF2 | ||
return __hbeq2(a.half2_, b.half2_); | ||
#else | ||
return (a.half_t2[0] == b.half_t2[0] && a.half_t2[1] == b.half_t2[1]); | ||
#endif | ||
} | ||
|
||
} // namespace half | ||
} // namespace mshadow | ||
#endif // MSHADOW_HALF2_H_ |
Oops, something went wrong.