softfloat.hpp

Go to the documentation of this file.

#pragma once
#ifdef SYS_VM_SOFTFLOAT
#include <sysio/vm/exceptions.hpp>
 
#include <cstdint>
 
#include <softfloat.hpp>
 
namespace sysio { namespace vm {
inline constexpr uint32_t inv_float_eps = 0x4B000000;
inline constexpr uint64_t inv_double_eps = 0x4330000000000000;
 
inline bool is_nan( const float32_t f ) {
   return f32_is_nan( f );
}
inline bool is_nan( const float64_t f ) {
   return f64_is_nan( f );
}
 
inline float _sysio_f32_add( float a, float b ) {
   return ::from_softfloat32( ::f32_add( ::to_softfloat32(a), ::to_softfloat32(b) ) );
}
 
inline float _sysio_f32_sub( float a, float b ) {
   return ::from_softfloat32( ::f32_sub( ::to_softfloat32(a), ::to_softfloat32(b) ) );
}
 
inline float _sysio_f32_div( float a, float b ) {
   return ::from_softfloat32( ::f32_div( ::to_softfloat32(a), ::to_softfloat32(b) ) );
}
 
inline float _sysio_f32_mul( float a, float b ) {
   return ::from_softfloat32( ::f32_mul( ::to_softfloat32(a), ::to_softfloat32(b) ) );
}
 
inline float _sysio_f32_min( float af, float bf ) {
   float32_t a = to_softfloat32(af);
   float32_t b = to_softfloat32(bf);
   if (is_nan(a)) {
      return af;
   }
   if (is_nan(b)) {
      return bf;
   }
   if ( f32_sign_bit(a) != f32_sign_bit(b) ) {
      return f32_sign_bit(a) ? af : bf;
   }
   return ::f32_lt(a,b) ? af : bf;
}
 
inline float _sysio_f32_max( float af, float bf ) {
   float32_t a = to_softfloat32(af);
   float32_t b = to_softfloat32(bf);
   if (is_nan(a)) {
      return af;
   }
   if (is_nan(b)) {
      return bf;
   }
   if ( f32_sign_bit(a) != f32_sign_bit(b) ) {
      return f32_sign_bit(a) ? bf : af;
   }
   return ::f32_lt( a, b ) ? bf : af;
}
 
inline float _sysio_f32_copysign( float af, float bf ) {
   float32_t a = to_softfloat32(af);
   float32_t b = to_softfloat32(bf);
   a.v &= ~(1 << 31);             // clear the sign bit
   a.v = a.v | ((b.v >> 31) << 31); // add the sign of b
   return from_softfloat32(a);
}
 
inline float _sysio_f32_abs( float af ) {
   float32_t a = to_softfloat32(af);
   a.v &= ~(1 << 31);
   return from_softfloat32(a);
}
 
inline float _sysio_f32_neg( float af ) {
   float32_t a = to_softfloat32(af);
   uint32_t sign = a.v >> 31;
   a.v &= ~(1 << 31);
   a.v |= (!sign << 31);
   return from_softfloat32(a);
}
 
inline float _sysio_f32_sqrt( float a ) {
   float32_t ret = ::f32_sqrt( to_softfloat32(a) );
   return from_softfloat32(ret);
}
 
inline float _sysio_f32_ceil( float af ) {
   float32_t a = to_softfloat32(af);
   int e = (int)(a.v >> 23 & 0xFF) - 0X7F;
   uint32_t m;
   if (e >= 23)
      return af;
   if (e >= 0) {
      m = 0x007FFFFF >> e;
      if ((a.v & m) == 0)
         return af;
      if (a.v >> 31 == 0)
         a.v += m;
      a.v &= ~m;
   } else {
      if (a.v >> 31)
         a.v = 0x80000000; // return -0.0f
      else if (a.v << 1)
         a.v = 0x3F800000; // return 1.0f
   }
 
   return from_softfloat32(a);
}
 
inline float _sysio_f32_floor( float af ) {
   float32_t a = to_softfloat32(af);
   int e = (int)(a.v >> 23 & 0xFF) - 0X7F;
   uint32_t m;
   if (e >= 23)
      return af;
   if (e >= 0) {
      m = 0x007FFFFF >> e;
      if ((a.v & m) == 0)
         return af;
      if (a.v >> 31)
         a.v += m;
      a.v &= ~m;
   } else {
      if (a.v >> 31 == 0)
         a.v = 0;
      else if (a.v << 1)
         a.v = 0xBF800000; // return -1.0f
   }
   return from_softfloat32(a);
}
 
inline float _sysio_f32_trunc( float af ) {
   float32_t a = to_softfloat32(af);
   int e = (int)(a.v >> 23 & 0xff) - 0x7f + 9;
   uint32_t m;
   if (e >= 23 + 9)
      return af;
   if (e < 9)
      e = 1;
   m = -1U >> e;
   if ((a.v & m) == 0)
      return af;
   a.v &= ~m;
   return from_softfloat32(a);
}
 
inline float _sysio_f32_nearest( float af ) {
   float32_t a = to_softfloat32(af);
   int e = a.v>>23 & 0xff;
   int s = a.v>>31;
   float32_t y;
   if (e >= 0x7f+23)
      return af;
   if (s)
      y = ::f32_add( ::f32_sub( a, float32_t{inv_float_eps} ), float32_t{inv_float_eps} );
   else
      y = ::f32_sub( ::f32_add( a, float32_t{inv_float_eps} ), float32_t{inv_float_eps} );
   if (::f32_eq( y, {0} ) )
      return s ? -0.0f : 0.0f;
   return from_softfloat32(y);
}
 
inline bool _sysio_f32_eq( float a, float b ) {  return ::f32_eq( to_softfloat32(a), to_softfloat32(b) ); }
inline bool _sysio_f32_ne( float a, float b ) { return !::f32_eq( to_softfloat32(a), to_softfloat32(b) ); }
inline bool _sysio_f32_lt( float a, float b ) { return ::f32_lt( to_softfloat32(a), to_softfloat32(b) ); }
inline bool _sysio_f32_le( float a, float b ) { return ::f32_le( to_softfloat32(a), to_softfloat32(b) ); }
inline bool _sysio_f32_gt( float af, float bf ) {
   float32_t a = to_softfloat32(af);
   float32_t b = to_softfloat32(bf);
   if (is_nan(a))
      return false;
   if (is_nan(b))
      return false;
   return !::f32_le( a, b );
}
 
inline bool _sysio_f32_ge( float af, float bf ) {
   float32_t a = to_softfloat32(af);
   float32_t b = to_softfloat32(bf);
   if (is_nan(a))
      return false;
   if (is_nan(b))
      return false;
   return !::f32_lt( a, b );
}
 
inline double _sysio_f64_add( double a, double b ) {
   float64_t ret = ::f64_add( to_softfloat64(a), to_softfloat64(b) );
   return from_softfloat64(ret);
}
 
inline double _sysio_f64_sub( double a, double b ) {
   float64_t ret = ::f64_sub( to_softfloat64(a), to_softfloat64(b) );
   return from_softfloat64(ret);
}
 
inline double _sysio_f64_div( double a, double b ) {
   float64_t ret = ::f64_div( to_softfloat64(a), to_softfloat64(b) );
   return from_softfloat64(ret);
}
 
inline double _sysio_f64_mul( double a, double b ) {
   float64_t ret = ::f64_mul( to_softfloat64(a), to_softfloat64(b) );
   return from_softfloat64(ret);
}
 
inline double _sysio_f64_min( double af, double bf ) {
   float64_t a = to_softfloat64(af);
   float64_t b = to_softfloat64(bf);
   if (is_nan(a))
      return af;
   if (is_nan(b))
      return bf;
   if (f64_sign_bit(a) != f64_sign_bit(b))
      return f64_sign_bit(a) ? af : bf;
   return ::f64_lt( a, b ) ? af : bf;
}
 
inline double _sysio_f64_max( double af, double bf ) {
   float64_t a = to_softfloat64(af);
   float64_t b = to_softfloat64(bf);
   if (is_nan(a))
      return af;
   if (is_nan(b))
      return bf;
   if (f64_sign_bit(a) != f64_sign_bit(b))
      return f64_sign_bit(a) ? bf : af;
   return ::f64_lt( a, b ) ? bf : af;
}
 
inline double _sysio_f64_copysign( double af, double bf ) {
   float64_t a = to_softfloat64(af);
   float64_t b = to_softfloat64(bf);
   a.v &= ~(uint64_t(1) << 63);     // clear the sign bit
   a.v = a.v | ((b.v >> 63) << 63); // add the sign of b
   return from_softfloat64(a);
}
 
inline double _sysio_f64_abs( double af ) {
   float64_t a = to_softfloat64(af);
   a.v &= ~(uint64_t(1) << 63);
   return from_softfloat64(a);
}
 
inline double _sysio_f64_neg( double af ) {
   float64_t a = to_softfloat64(af);
   uint64_t sign = a.v >> 63;
   a.v &= ~(uint64_t(1) << 63);
   a.v |= (uint64_t(!sign) << 63);
   return from_softfloat64(a);
}
 
inline double _sysio_f64_sqrt( double a ) {
   float64_t ret = ::f64_sqrt( to_softfloat64(a) );
   return from_softfloat64(ret);
}
 
inline double _sysio_f64_ceil( double af ) {
   float64_t a = to_softfloat64( af );
   float64_t ret;
   int e = a.v >> 52 & 0x7ff;
   float64_t y;
   if (e >= 0x3ff+52 || ::f64_eq( a, { 0 } ))
      return af;
   /* y = int(x) - x, where int(x) is an integer neighbor of x */
   if (a.v >> 63)
      y = ::f64_sub( ::f64_add( ::f64_sub( a, float64_t{inv_double_eps} ), float64_t{inv_double_eps} ), a );
   else
      y = ::f64_sub( ::f64_sub( ::f64_add( a, float64_t{inv_double_eps} ), float64_t{inv_double_eps} ), a );
   /* special case because of non-nearest rounding modes */
   if (e <= 0x3ff-1) {
      return a.v >> 63 ? -0.0 : 1.0; //float64_t{0x8000000000000000} : float64_t{0xBE99999A3F800000}; //either -0.0 or 1
   }
   if (::f64_lt( y, to_softfloat64(0) )) {
      ret = ::f64_add( ::f64_add( a, y ), to_softfloat64(1) ); // 0xBE99999A3F800000 } ); // plus 1
      return from_softfloat64(ret);
   }
   ret = ::f64_add( a, y );
   return from_softfloat64(ret);
}
 
inline double _sysio_f64_floor( double af ) {
   float64_t a = to_softfloat64( af );
   float64_t ret;
   int e = a.v >> 52 & 0x7FF;
   float64_t y;
   if ( a.v == 0x8000000000000000) {
      return af;
   }
   if (e >= 0x3FF+52 || a.v == 0) {
      return af;
   }
   if (a.v >> 63)
      y = ::f64_sub( ::f64_add( ::f64_sub( a, float64_t{inv_double_eps} ), float64_t{inv_double_eps} ), a );
   else
      y = ::f64_sub( ::f64_sub( ::f64_add( a, float64_t{inv_double_eps} ), float64_t{inv_double_eps} ), a );
   if (e <= 0x3FF-1) {
      return a.v>>63 ? -1.0 : 0.0; //float64_t{0xBFF0000000000000} : float64_t{0}; // -1 or 0
   }
   if ( !::f64_le( y, float64_t{0} ) ) {
      ret = ::f64_sub( ::f64_add(a,y), to_softfloat64(1.0));
      return from_softfloat64(ret);
   }
   ret = ::f64_add( a, y );
   return from_softfloat64(ret);
}
 
inline double _sysio_f64_trunc( double af ) {
   float64_t a = to_softfloat64( af );
   int e = (int)(a.v >> 52 & 0x7ff) - 0x3ff + 12;
   uint64_t m;
   if (e >= 52 + 12)
      return af;
   if (e < 12)
      e = 1;
   m = -1ULL >> e;
   if ((a.v & m) == 0)
      return af;
   a.v &= ~m;
   return from_softfloat64(a);
}
 
inline double _sysio_f64_nearest( double af ) {
   float64_t a = to_softfloat64( af );
   int e = (a.v >> 52 & 0x7FF);
   int s = a.v >> 63;
   float64_t y;
   if ( e >= 0x3FF+52 )
      return af;
   if ( s )
      y = ::f64_add( ::f64_sub( a, float64_t{inv_double_eps} ), float64_t{inv_double_eps} );
   else
      y = ::f64_sub( ::f64_add( a, float64_t{inv_double_eps} ), float64_t{inv_double_eps} );
   if ( ::f64_eq( y, float64_t{0} ) )
      return s ? -0.0 : 0.0;
   return from_softfloat64(y);
}
 
inline bool _sysio_f64_eq( double a, double b ) { return ::f64_eq( to_softfloat64(a), to_softfloat64(b) ); }
inline bool _sysio_f64_ne( double a, double b ) { return !::f64_eq( to_softfloat64(a), to_softfloat64(b) ); }
inline bool _sysio_f64_lt( double a, double b ) { return ::f64_lt( to_softfloat64(a), to_softfloat64(b) ); }
inline bool _sysio_f64_le( double a, double b ) { return ::f64_le( to_softfloat64(a), to_softfloat64(b) ); }
inline bool _sysio_f64_gt( double af, double bf ) {
   float64_t a = to_softfloat64(af);
   float64_t b = to_softfloat64(bf);
   if (is_nan(a))
      return false;
   if (is_nan(b))
      return false;
   return !::f64_le( a, b );
}
 
inline bool _sysio_f64_ge( double af, double bf ) {
   float64_t a = to_softfloat64(af);
   float64_t b = to_softfloat64(bf);
   if (is_nan(a))
      return false;
   if (is_nan(b))
      return false;
   return !::f64_lt( a, b );
}
 
inline double _sysio_f32_promote( float a ) {
   return from_softfloat64(f32_to_f64( to_softfloat32(a)) );
}
 
inline float _sysio_f64_demote( double a ) {
   return from_softfloat32(f64_to_f32( to_softfloat64(a)) );
}
 
inline int32_t _sysio_f32_trunc_i32s( float af ) {
   float32_t a = to_softfloat32(af);
   SYS_VM_ASSERT(!(_sysio_f32_ge(af, 2147483648.0f) || _sysio_f32_lt(af, -2147483648.0f)), wasm_interpreter_exception, "Error, f32.convert_s/i32 overflow" );
 
   SYS_VM_ASSERT(!is_nan(a), wasm_interpreter_exception, "Error, f32.convert_s/i32 unrepresentable");
   return f32_to_i32( to_softfloat32(_sysio_f32_trunc( af )), 0, false );
}
 
inline int32_t _sysio_f64_trunc_i32s( double af ) {
   float64_t a = to_softfloat64(af);
   SYS_VM_ASSERT(!(_sysio_f64_ge(af, 2147483648.0) || _sysio_f64_lt(af, -2147483648.0)), wasm_interpreter_exception, "Error, f64.convert_s/i32 overflow");
   SYS_VM_ASSERT(!is_nan(a), wasm_interpreter_exception, "Error, f64.convert_s/i32 unrepresentable");
   return f64_to_i32( to_softfloat64(_sysio_f64_trunc( af )), 0, false );
}
 
inline uint32_t _sysio_f32_trunc_i32u( float af ) {
   float32_t a = to_softfloat32(af);
   SYS_VM_ASSERT(!(_sysio_f32_ge(af, 4294967296.0f) || _sysio_f32_le(af, -1.0f)),wasm_interpreter_exception, "Error, f32.convert_u/i32 overflow");
   SYS_VM_ASSERT(!is_nan(a), wasm_interpreter_exception, "Error, f32.convert_u/i32 unrepresentable");
   return f32_to_ui32( to_softfloat32(_sysio_f32_trunc( af )), 0, false );
}
 
inline uint32_t _sysio_f64_trunc_i32u( double af ) {
   float64_t a = to_softfloat64(af);
   SYS_VM_ASSERT(!(_sysio_f64_ge(af, 4294967296.0) || _sysio_f64_le(af, -1.0)), wasm_interpreter_exception, "Error, f64.convert_u/i32 overflow");
   SYS_VM_ASSERT(!is_nan(a), wasm_interpreter_exception, "Error, f64.convert_u/i32 unrepresentable");
   return f64_to_ui32( to_softfloat64(_sysio_f64_trunc( af )), 0, false );
}
 
inline int64_t _sysio_f32_trunc_i64s( float af ) {
   float32_t a = to_softfloat32(af);
      SYS_VM_ASSERT(!(_sysio_f32_ge(af, 9223372036854775808.0f) || _sysio_f32_lt(af, -9223372036854775808.0f)), wasm_interpreter_exception, "Error, f32.convert_s/i64 overflow");
   SYS_VM_ASSERT(!is_nan(a), wasm_interpreter_exception, "Error, f32.convert_s/i64 unrepresentable");
   return f32_to_i64( to_softfloat32(_sysio_f32_trunc( af )), 0, false );
}
 
inline int64_t _sysio_f64_trunc_i64s( double af ) {
   float64_t a = to_softfloat64(af);
   SYS_VM_ASSERT(!(_sysio_f64_ge(af, 9223372036854775808.0) || _sysio_f64_lt(af, -9223372036854775808.0)), wasm_interpreter_exception, "Error, f64.convert_s/i64 overflow");
   SYS_VM_ASSERT(!is_nan(a), wasm_interpreter_exception, "Error, f64.convert_s/i64 unrepresentable");
 
   return f64_to_i64( to_softfloat64(_sysio_f64_trunc( af )), 0, false );
}
 
inline uint64_t _sysio_f32_trunc_i64u( float af ) {
   float32_t a = to_softfloat32(af);
   SYS_VM_ASSERT(!(_sysio_f32_ge(af, 18446744073709551616.0f) || _sysio_f32_le(af, -1.0f)), wasm_interpreter_exception, "Error, f32.convert_u/i64 overflow");
   SYS_VM_ASSERT(!is_nan(a), wasm_interpreter_exception, "Error, f32.convert_u/i64 unrepresentable");
   return f32_to_ui64( to_softfloat32(_sysio_f32_trunc( af )), 0, false );
}
 
inline uint64_t _sysio_f64_trunc_i64u( double af ) {
   float64_t a = to_softfloat64(af);
   SYS_VM_ASSERT(!(_sysio_f64_ge(af, 18446744073709551616.0) || _sysio_f64_le(af, -1.0)), wasm_interpreter_exception, "Error, f64.convert_u/i64 overflow");
   SYS_VM_ASSERT(!is_nan(a), wasm_interpreter_exception, "Error, f64.convert_u/i64 unrepresentable");
   return f64_to_ui64( to_softfloat64(_sysio_f64_trunc( af )), 0, false );
}
 
inline float _sysio_i32_to_f32( int32_t a ) {
   return from_softfloat32(i32_to_f32( a ));
}
 
inline float _sysio_i64_to_f32( int64_t a ) {
   return from_softfloat32(i64_to_f32( a ));
}
 
inline float _sysio_ui32_to_f32( uint32_t a ) {
   return from_softfloat32(ui32_to_f32( a ));
}
 
inline float _sysio_ui64_to_f32( uint64_t a ) {
   return from_softfloat32(ui64_to_f32( a ));
}
 
inline double _sysio_i32_to_f64( int32_t a ) {
   return from_softfloat64(i32_to_f64( a ));
}
 
inline double _sysio_i64_to_f64( int64_t a ) {
   return from_softfloat64(i64_to_f64( a ));
}
 
inline double _sysio_ui32_to_f64( uint32_t a ) {
   return from_softfloat64(ui32_to_f64( a ));
}
 
inline double _sysio_ui64_to_f64( uint64_t a ) {
   return from_softfloat64(ui64_to_f64( a ));
}
}} //ns sysio::vm
#endif