sysdeps/loongarch/fpu/s_fmaximum_numf.c

   1 /* fmaximum_numf().  LoongArch version.
   2    Copyright (C) 2022-2023 Free Software Foundation, Inc.
   3
   4    This file is part of the GNU C Library.
   5
   6    The GNU C Library is free software; you can redistribute it and/or
   7    modify it under the terms of the GNU Lesser General Public
   8    License as published by the Free Software Foundation; either
   9    version 2.1 of the License, or (at your option) any later version.
  10
  11    The GNU C Library is distributed in the hope that it will be useful,
  12    but WITHOUT ANY WARRANTY; without even the implied warranty of
  13    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14    Lesser General Public License for more details.
  15
  16    You should have received a copy of the GNU Lesser General Public
  17    License along with the GNU C Library; if not, see
  18    <https://www.gnu.org/licenses/>.  */
  19
  20 #define NO_MATH_REDIRECT
  21 #include <math.h>
  22 #include <libm-alias-float.h>
  23 #include <fpu_control.h>
  24
  25 float
  26 __fmaximum_numf (float x, float y)
  27 {
  28   int x_cond;
  29   int y_cond;
  30   asm volatile ("fclass.s \t%0, %1" : "=f" (x_cond) : "f" (x));
  31   asm volatile ("fclass.s \t%0, %1" : "=f" (y_cond) : "f" (y));
  32
  33   if (__glibc_unlikely((x_cond & _FCLASS_NAN) && !(y_cond & _FCLASS_NAN)))
  34     {
  35       asm volatile ("fmax.s \t%0, %1, %2" : "=f" (x) : "f" (x), "f" (y));
  36       return y;
  37     }
  38   else if (__glibc_unlikely(!(x_cond & _FCLASS_NAN) && (y_cond & _FCLASS_NAN)))
  39     {
  40       asm volatile ("fmax.s \t%0, %1, %2" : "=f" (y) : "f" (x), "f" (y));
  41       return x;
  42     }
  43   else
  44     {
  45       asm volatile ("fmax.s \t%0, %1, %2" : "=f" (x) : "f" (x), "f" (y));
  46       return x;
  47     }
  48 }
  49 libm_alias_float (__fmaximum_num, fmaximum_num)