softfloat/f64_lt_quiet.c

   1
   2 #include <stdbool.h>
   3 #include <stdint.h>
   4 #include "platform.h"
   5 #include "internals.h"
   6 #include "specialize.h"
   7 #include "softfloat.h"
   8
   9 bool f64_lt_quiet( float64_t a, float64_t b )
  10 {
  11     union ui64_f64 uA;
  12     uint_fast64_t uiA;
  13     union ui64_f64 uB;
  14     uint_fast64_t uiB;
  15     bool signA, signB;
  16
  17     uA.f = a;
  18     uiA = uA.ui;
  19     uB.f = b;
  20     uiB = uB.ui;
  21     if (
  22            ( ( expF64UI( uiA ) == 0x7FF ) && fracF64UI( uiA ) )
  23         || ( ( expF64UI( uiB ) == 0x7FF ) && fracF64UI( uiB ) )
  24     ) {
  25         if (
  26             softfloat_isSigNaNF64UI( uiA ) || softfloat_isSigNaNF64UI( uiB )
  27         ) {
  28             softfloat_raiseFlags( softfloat_flag_invalid );
  29         }
  30         return false;
  31     }
  32     signA = signF64UI( uiA );
  33     signB = signF64UI( uiB );
  34     return
  35         ( signA != signB )
  36             ? signA && ( ( uiA | uiB ) & UINT64_C( 0x7FFFFFFFFFFFFFFF ) )
  37             : ( uiA != uiB ) && ( signA ^ ( uiA < uiB ) );
  38
  39 }
  40