libclc/generic/lib/math/fract.inc

   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
  51 FRACT_DEF(local);
  52 FRACT_DEF(global);
  53
  54 #undef MIN_CONSTANT
  55 #undef ZERO