newlib/libm/machine/amdgcn/v64df_ilogb.c

   1 /*
   2  * Copyright 2023 Siemens
   3  *
   4  * The authors hereby grant permission to use, copy, modify, distribute,
   5  * and license this software and its documentation for any purpose, provided
   6  * that existing copyright notices are retained in all copies and that this
   7  * notice is included verbatim in any distributions.  No written agreement,
   8  * license, or royalty fee is required for any of the authorized uses.
   9  * Modifications to this software may be copyrighted by their authors
  10  * and need not follow the licensing terms described here, provided that
  11  * the new terms are clearly indicated on the first page of each file where
  12  * they apply.
  13  */
  14
  15 /*
  16  * ====================================================
  17  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
  18  *
  19  * Developed at SunPro, a Sun Microsystems, Inc. business.
  20  * Permission to use, copy, modify, and distribute this
  21  * software is freely granted, provided that this notice
  22  * is preserved.
  23  * ====================================================
  24  */
  25
  26 /* Based on newlib/libm/common/s_ilogb.c in Newlib.  */
  27
  28 #include "amdgcnmach.h"
  29
  30 DEF_VD_MATH_PRED (v64si, ilogb, v64df x)
  31 {
  32   FUNCTION_INIT(v64si);
  33   v64si hx, lx, ix;
  34   EXTRACT_WORDS (hx, lx, x);
  35   hx &= 0x7fffffff;
  36   VECTOR_IF (hx < 0x00100000, cond)
  37     VECTOR_RETURN (VECTOR_INIT (-__INT_MAX__), cond & ((hx | lx) == 0));  // FP_ILOGB0
  38     VECTOR_IF2 (hx == 0, cond2, cond)
  39       ix = VECTOR_INIT (-1043);
  40       for (v64si i = lx;
  41             !ALL_ZEROES_P (cond2 & (i > 0));
  42             i <<= 1)
  43         VECTOR_COND_MOVE (ix, ix - 1, cond2 & (i > 0));
  44     VECTOR_ELSE2 (cond2, cond)
  45       ix = VECTOR_INIT (-1022);
  46       for (v64si i = (hx << 11);
  47             !ALL_ZEROES_P (cond2 & (i > 0));
  48             i <<= 1)
  49         VECTOR_COND_MOVE (ix, ix - 1, cond2 & (i > 0));
  50     VECTOR_ENDIF
  51     VECTOR_RETURN (ix, cond);
  52   VECTOR_ENDIF
  53   VECTOR_RETURN ((hx >> 20) - 1023, hx < 0x7ff00000);
  54   VECTOR_RETURN (VECTOR_INIT (__INT_MAX__), NO_COND);
  55
  56   FUNCTION_RETURN;
  57 }
  58
  59 DEF_VARIANTS (ilogb, si, df)