src/math/fmax.c

   1 /*
   2  * fmax: floating point maximum function.
   3  * Copyright (C) 2010 Nick Bowler.
   4  *
   5  * License LGPLv3+: GNU LGPL version 3 or later.  See COPYING.LIB for terms.
   6  * See <http://www.gnu.org/licenses/lgpl.html> if you did not receive a copy.
   7  * This is free software: you are free to change and redistribute it.
   8  * There is NO WARRANTY, to the extent permitted by law.
   9  */
  10
  11 #include <cstd/mathcommon.h>
  12
  13 TYPE MF(fmax)(TYPE x, TYPE y)
  14 {
  15         if (isgreater(x, y) || isnan(y))
  16                 return x;
  17         if (isgreater(y, x) || isnan(x))
  18                 return y;
  19         if (signbit(x))
  20                 return y;
  21         return x;
  22 }