newlib/libm/machine/amdgcn/v64sf_tgamma.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/ef_tgamma.c in Newlib. */
  28
  29 #include "amdgcnmach.h"
  30
  31 v64sf v64sf_expf_aux (v64sf x, v64si __mask);
  32 v64sf v64sf_lgammaf_r_aux (v64sf x, v64si *signgamp, v64si __mask);
  33
  34 DEF_VS_MATH_FUNC (v64sf, tgammaf, v64sf x)
  35 {
  36   v64si signgam_local;
  37   v64sf y = v64sf_expf_aux(v64sf_lgammaf_r_aux(x, &signgam_local, __mask), __mask);
  38   VECTOR_COND_MOVE(y, -y, signgam_local < 0);
  39         return y;
  40 }
  41
  42 DEF_VARIANTS (tgammaf, sf, sf)