newlib/libm/machine/amdgcn/v64df_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/e_tgamma.c in Newlib. */
  28
  29 #include "amdgcnmach.h"
  30
  31 v64df v64df_exp_aux (v64df x, v64di __mask);
  32 v64df v64df_lgamma_r_aux (v64df x, v64si *signgamp, v64di __mask);
  33
  34 DEF_VD_MATH_FUNC (v64df, tgamma, v64df x)
  35 {
  36   v64si signgam_local;
  37   v64df y = v64df_exp_aux(v64df_lgamma_r_aux(x, &signgam_local, __mask), __mask);
  38   VECTOR_COND_MOVE(y, -y, signgam_local < 0);
  39         return y;
  40 }
  41
  42 DEF_VARIANTS (tgamma, df, df)