softfloat/f64_eq.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_eq( 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
  16     uA.f = a;
  17     uiA = uA.ui;
  18     uB.f = b;
  19     uiB = uB.ui;
  20     if (
  21            ( ( expF64UI( uiA ) == 0x7FF ) && fracF64UI( uiA ) )
  22         || ( ( expF64UI( uiB ) == 0x7FF ) && fracF64UI( uiB ) )
  23     ) {
  24         if (
  25             softfloat_isSigNaNF64UI( uiA ) || softfloat_isSigNaNF64UI( uiB )
  26         ) {
  27             softfloat_raiseFlags( softfloat_flag_invalid );
  28         }
  29         return false;
  30     }
  31     return
  32         ( uiA == uiB ) || ! ( ( uiA | uiB ) & UINT64_C( 0x7FFFFFFFFFFFFFFF ) );
  33
  34 }
  35