w_lgamma_r_template.c 1.7 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748
  1. /* Wrapper to set errno for lgamma_r.
  2. Copyright (C) 2017-2019 Free Software Foundation, Inc.
  3. This file is part of the GNU C Library.
  4. The GNU C Library is free software; you can redistribute it and/or
  5. modify it under the terms of the GNU Lesser General Public
  6. License as published by the Free Software Foundation; either
  7. version 2.1 of the License, or (at your option) any later version.
  8. The GNU C Library is distributed in the hope that it will be useful,
  9. but WITHOUT ANY WARRANTY; without even the implied warranty of
  10. MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
  11. Lesser General Public License for more details.
  12. You should have received a copy of the GNU Lesser General Public
  13. License along with the GNU C Library; if not, see
  14. <http://www.gnu.org/licenses/>. */
  15. /* Only build wrappers from the templates for the types that define the macro
  16. below. This macro is set in math-type-macros-<type>.h in sysdeps/generic
  17. for each floating-point type. */
  18. #if __USE_WRAPPER_TEMPLATE
  19. # include <errno.h>
  20. # include <fenv.h>
  21. # include <math.h>
  22. # include <math_private.h>
  23. #define M_DECL_FUNC_R_X(x) x ## _r
  24. #define M_DECL_FUNC_R_S(x) M_DECL_FUNC_R_X (x)
  25. #define M_DECL_FUNC_R(x) M_DECL_FUNC_R_S (M_SUF (x))
  26. #define M_CALL_FUNC_R_X(x) x ## _r
  27. #define M_CALL_FUNC_R_S(x) M_CALL_FUNC_R_X (x)
  28. #define M_CALL_FUNC_R(x) M_CALL_FUNC_R_S (M_SUF (x))
  29. FLOAT
  30. M_DECL_FUNC_R (__lgamma) (FLOAT x, int *signgamp)
  31. {
  32. FLOAT y = M_CALL_FUNC_R (__ieee754_lgamma) (x, signgamp);
  33. if (__glibc_unlikely (!isfinite (y)) && isfinite (x))
  34. /* Pole error: lgamma_r(integer x<0). Or overflow. */
  35. __set_errno (ERANGE);
  36. return y;
  37. }
  38. declare_mgen_alias_r (__lgamma, lgamma)
  39. #endif /* __USE_WRAPPER_TEMPLATE. */