aaron/ispc
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

first commit for {rsqrt,rcp}d knc support. going to test on other node now

Browse Source
This commit is contained in:
Evghenii
2014-02-05 13:43:07 +01:00
parent f225b558ec
commit eb01ffd4e6
3 changed files with 93 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

32
examples/intrinsics/knc.h
View File

@@ -1472,6 +1472,21 @@ static FORCEINLINE __vec16_f __rcp_varying_float(__vec16_f v) {
return _mm512_recip_ps(v);
#endif
}
static FORCEINLINE __vec16_d __rcp_varying_double(__vec16_d x) {
__vec16_i64 ex = __and(__cast_bits(__vec16_i64(), x), __smear_i64<__vec16_i64>(0x7fe0000000000000));
__vec16_d exp = __cast_bits(__vec16_d(), __sub(__smear_i64<__vec16_i64>(0x7fd0000000000000), ex));
__vec16_f xf = __cast_fptrunc(__vec16_f(), __mul(x, exp));
__vec16_f yf = __rcp_varying_float(xf);
__vec16_d y = __mul(__cast_fpext(__vec16_d(), yf), exp);
y = __add(y, __mul(y, __sub(__smear_double<__vec16_d>(2.0), __mul(x, y))));
y = __add(y, __mul(y, __sub(__smear_double<__vec16_d>(2.0), __mul(x, y))));
return y;
}
static FORCEINLINE double __rcp_uniform_double(double v)
{
return __extract_element(__rcp_varying_double(__smear_double<__vec16_d>(v)),0);
}
static FORCEINLINE __vec16_f __rsqrt_varying_float(__vec16_f v) {
#ifdef ISPC_FAST_MATH
@@ -1480,6 +1495,23 @@ static FORCEINLINE __vec16_f __rsqrt_varying_float(__vec16_f v) {
return _mm512_invsqrt_ps(v);
#endif
}
static FORCEINLINE __vec16_d __rsqrt_varying_double(__vec16_d x) {
__vec16_i64 ex = __and(__cast_bits(__vec16_i64(), x), __smear_i64<__vec16_i64>(0x7fe0000000000000));
__vec16_d exp = __cast_bits(__vec16_d(), __sub(__smear_i64<__vec16_i64>(0x7fd0000000000000), ex));
__vec16_d exph = __cast_bits(__vec16_d(), __sub(__smear_i64<__vec16_i64>(0x5fe0000000000000), __lshr(ex,1)));
__vec16_f xf = __cast_fptrunc(__vec16_f(), __mul(x, exp));
__vec16_f yf = __rsqrt_varying_float(xf);
__vec16_d y = __mul(__cast_fpext(__vec16_d(), yf), exph);
__vec16_d xh = __mul(x, __smear_double<__vec16_d>(0.5));
y = __add(y, __mul(y, __sub(__smear_double<__vec16_d>(0.5), __mul(xh, __mul(y,y)))));
y = __add(y, __mul(y, __sub(__smear_double<__vec16_d>(0.5), __mul(xh, __mul(y,y)))));
return y;
}
static FORCEINLINE double __rsqrt_uniform_double(double v)
{
return __extract_element(__rsqrt_varying_double(__smear_double<__vec16_d>(v)),0);
}
static FORCEINLINE __vec16_f __exp_varying_float(__vec16_f v) {
return _mm512_exp_ps(v);