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
|
/*========================== begin_copyright_notice ============================
Copyright (C) 2017-2021 Intel Corporation
SPDX-License-Identifier: MIT
============================= end_copyright_notice ===========================*/
#include "../include/BiF_Definitions.cl"
#include "../../Headers/spirv.h"
#include "../IMF/FP32/asinh_s_la.cl"
#if defined(cl_khr_fp64)
#include "../IMF/FP64/asinh_d_la.cl"
#endif // defined(cl_khr_fp64)
INLINE float __intel_asinh_f32( float x, bool doFast )
{
float result;
if(__FastRelaxedMath && (!__APIRS) && doFast)
{
// Implemented as log(x + sqrt(x*x + 1)).
// Conformance test checks for this "overflow" case, but
// I don't think we should have to handle it.
if( x > 1500.0f )
{
result = SPIRV_OCL_BUILTIN(log, _f32, )(x) + M_LN2_F;
}
else
{
result = x * x + 1.0f;
result = SPIRV_OCL_BUILTIN(sqrt, _f32, )( result );
result = x + result;
result = SPIRV_OCL_BUILTIN(log, _f32, )( result );
}
}
else
{
result = __ocl_svml_asinhf(x);
}
return result;
}
INLINE float SPIRV_OVERLOADABLE SPIRV_OCL_BUILTIN(asinh, _f32, )( float x )
{
return __intel_asinh_f32( x, true );
}
GENERATE_SPIRV_OCL_VECTOR_FUNCTIONS_1ARG_LOOP( asinh, float, float, f32 )
#if defined(cl_khr_fp64)
INLINE double SPIRV_OVERLOADABLE SPIRV_OCL_BUILTIN(asinh, _f64, )( double x )
{
return __ocl_svml_asinh(x);
}
GENERATE_SPIRV_OCL_VECTOR_FUNCTIONS_1ARGS( asinh, double, double, f64 )
#endif // defined(cl_khr_fp64)
#if defined(cl_khr_fp16)
INLINE half SPIRV_OVERLOADABLE SPIRV_OCL_BUILTIN(asinh, _f16, )( half x )
{
return SPIRV_OCL_BUILTIN(asinh, _f32, )((float)x);
}
GENERATE_SPIRV_OCL_VECTOR_FUNCTIONS_1ARG_LOOP( asinh, half, half, f16 )
#endif // defined(cl_khr_fp16)
|