1/* 2 * Copyright (c) 2015 Advanced Micro Devices, Inc. 3 * 4 * Permission is hereby granted, free of charge, to any person obtaining a copy 5 * of this software and associated documentation files (the "Software"), to deal 6 * in the Software without restriction, including without limitation the rights 7 * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell 8 * copies of the Software, and to permit persons to whom the Software is 9 * furnished to do so, subject to the following conditions: 10 * 11 * The above copyright notice and this permission notice shall be included in 12 * all copies or substantial portions of the Software. 13 * 14 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR 15 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, 16 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE 17 * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER 18 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, 19 * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN 20 * THE SOFTWARE. 21 */ 22 23#if __CLC_FPSIZE == 64 24#define MIN_CONSTANT 0x1.fffffffffffffp-1 25#define ZERO 0.0 26#elif __CLC_FPSIZE == 32 27#define MIN_CONSTANT 0x1.fffffep-1f 28#define ZERO 0.0f 29#elif __CLC_FPSIZE == 16 30#define MIN_CONSTANT 0x1.ffcp-1h 31#define ZERO 0.0h 32#endif 33 34_CLC_OVERLOAD _CLC_DEF __CLC_GENTYPE fract(__CLC_GENTYPE x, private __CLC_GENTYPE *iptr) { 35 *iptr = floor(x); 36 __CLC_GENTYPE r = fmin(x - *iptr, MIN_CONSTANT); 37 r = isinf(x) ? ZERO : r; 38 r = isnan(x) ? x : r; 39 return r; 40} 41 42 43#define FRACT_DEF(addrspace) \ 44 _CLC_OVERLOAD _CLC_DEF __CLC_GENTYPE fract(__CLC_GENTYPE x, addrspace __CLC_GENTYPE *iptr) { \ 45 __CLC_GENTYPE private_iptr; \ 46 __CLC_GENTYPE ret = fract(x, &private_iptr); \ 47 *iptr = private_iptr; \ 48 return ret; \ 49 } 50 51FRACT_DEF(local); 52FRACT_DEF(global); 53 54#undef MIN_CONSTANT 55#undef ZERO 56