newlib/libm/machine/amdgcn/v64sf_lgamma.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
  27 /* Based on newlib/libm/math/wf_lgamma.c in Newlib.  */
  28
  29 #include "amdgcnmach.h"
  30
  31 v64si v64sf_finitef_aux (v64sf x, v64si __mask);
  32 v64sf v64sf_lgammaf_r_aux (v64sf x, v64si *signgamp, v64si __mask);
  33
  34 DEF_VS_MATH_FUNC (v64sf, lgammaf, v64sf x)
  35 {
  36   v64sf y = v64sf_lgammaf_r_aux(x, &(_REENT_V64SI_SIGNGAM(_V64_REENT)), __mask);
  37   if (ALL_ZEROES_P(v64sf_finitef_aux(y, __mask)) & !ALL_ZEROES_P(v64sf_finitef_aux(x, __mask))) {
  38     /* lgamma(finite) overflow */
  39     errno = ERANGE;
  40   }
  41   return y;
  42 }
  43
  44 DEF_VARIANTS (lgammaf, sf, sf)