1 /* Return value with minimum magnitude. 2 Copyright (C) 2016-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 #include <math.h> 20 21 FLOAT M_DECL_FUNC(__fminmag)22M_DECL_FUNC (__fminmag) (FLOAT x, FLOAT y) 23 { 24 FLOAT ax = M_FABS (x); 25 FLOAT ay = M_FABS (y); 26 if (isless (ax, ay)) 27 return x; 28 else if (isgreater (ax, ay)) 29 return y; 30 else if (ax == ay) 31 return x < y ? x : y; 32 else if (issignaling (x) || issignaling (y)) 33 return x + y; 34 else 35 return isnan (y) ? x : y; 36 } 37 38 declare_mgen_alias (__fminmag, fminmag); 39