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)