diff options
Diffstat (limited to 'math/s_fmaximum_mag_num_template.c')
-rw-r--r-- | math/s_fmaximum_mag_num_template.c | 35 |
1 files changed, 35 insertions, 0 deletions
diff --git a/math/s_fmaximum_mag_num_template.c b/math/s_fmaximum_mag_num_template.c new file mode 100644 index 0000000000..5c9cef21fd --- /dev/null +++ b/math/s_fmaximum_mag_num_template.c @@ -0,0 +1,35 @@ +/* Return number with maximum magnitude. + Copyright (C) 2016-2021 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + <https://www.gnu.org/licenses/>. */ + +#include <math.h> + +FLOAT +M_DECL_FUNC (__fmaximum_mag_num) (FLOAT x, FLOAT y) +{ + FLOAT ax = M_FABS (x); + FLOAT ay = M_FABS (y); + if (isgreater (ax, ay)) + return x; + else if (isless (ax, ay)) + return y; + else if (ax == ay) + return (M_COPYSIGN (1, x) >= M_COPYSIGN (1, y) ? x : y); + else + return isnan (y) ? (isnan (x) ? x + y : x) : y; +} +declare_mgen_alias (__fmaximum_mag_num, fmaximum_mag_num); |