1 /* Wrapper to set errno for j1 and y1. 2 Copyright (C) 2017-2022 Free Software Foundation, Inc. 3 This file is part of the GNU C Library. 4 5 The GNU C Library is free software; you can redistribute it and/or 6 modify it under the terms of the GNU Lesser General Public 7 License as published by the Free Software Foundation; either 8 version 2.1 of the License, or (at your option) any later version. 9 10 The GNU C Library is distributed in the hope that it will be useful, 11 but WITHOUT ANY WARRANTY; without even the implied warranty of 12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 13 Lesser General Public License for more details. 14 15 You should have received a copy of the GNU Lesser General Public 16 License along with the GNU C Library; if not, see 17 <https://www.gnu.org/licenses/>. */ 18 19 /* Only build wrappers from the templates for the types that define the macro 20 below. This macro is set in math-type-macros-<type>.h in sysdeps/generic 21 for each floating-point type. */ 22 #if __USE_WRAPPER_TEMPLATE 23 24 # include <errno.h> 25 # include <fenv.h> 26 # include <math.h> 27 # include <math_private.h> 28 29 FLOAT M_DECL_FUNC(__j1)30M_DECL_FUNC (__j1) (FLOAT x) 31 { 32 return M_SUF (__ieee754_j1) (x); 33 } declare_mgen_alias(__j1,j1)34declare_mgen_alias (__j1, j1) 35 36 FLOAT 37 M_DECL_FUNC (__y1) (FLOAT x) 38 { 39 if (__glibc_unlikely (islessequal (x, M_LIT (0.0)))) 40 { 41 if (x < 0) 42 /* Domain error: y1(x<0). */ 43 __set_errno (EDOM); 44 else if (x == 0) 45 /* Pole error: y1(0). */ 46 __set_errno (ERANGE); 47 } 48 return M_SUF (__ieee754_y1) (x); 49 } 50 declare_mgen_alias (__y1, y1) 51 52 #endif /* __USE_WRAPPER_TEMPLATE. */ 53