Add quad floating point to float lib

lib/float/common.sail

@@ -32,8 +32,8 @@ $define _FLOAT_COMMON
 
 /* Floating point types definition */
 type fp_exception_flags = bits(5) /* Floating-point exception flags. */
-type fp_bits = { 'n, 'n in {16, 32, 64}. bits('n) } /* Floating-point in bits. */
-type fp_bits_x2 = { 'n, 'n in {16, 32, 64}. (bits('n), bits('n)) } /* Floating point x2 tuple */
+type fp_bits = { 'n, 'n in {16, 32, 64, 128}. bits('n) } /* Floating-point in bits. */
+type fp_bits_x2 = { 'n, 'n in {16, 32, 64, 128}. (bits('n), bits('n)) } /* Floating point x2 tuple */
 type fp_bool_and_flags = (bool, fp_exception_flags) /* Floating point bool and exception flags tuple */
 
 /* Floating point constants */
@@ -47,16 +47,20 @@ struct float_bits('n : Int) = {
  then 5
  else (if 'n == 32
  then 8
- else 11)),
+ else (if 'n == 64
+ then 11
+ else 15))),
  mantissa : bits(if 'n == 16
  then 10
  else (if 'n == 32
  then 23
- else 52)),
+ else (if 'n == 64
+ then 52
+ else 112))),
 }
 
 /* The val func implementations */
-val float_decompose : forall 'n, 'n in { 16, 32, 64 }. bits('n) -> float_bits('n)
+val float_decompose : forall 'n, 'n in { 16, 32, 64, 128 }. bits('n) -> float_bits('n)
 function float_decompose(op) = {
  match 'n {
  16 => struct {
@@ -73,6 +77,11 @@ function float_decompose(op) = {
  sign = op[63..63],
  exp = op[62..52],
  mantissa = op[51..0],
+ },
+ 128 => struct {
+ sign = op[127..127],
+ exp = op[126..112],
+ mantissa = op[111..0],
  }
  }
 }

test/float/eq_test.sail

@@ -75,6 +75,21 @@ function test_float_is_eq () -> unit = {
 
  assert(float_is_eq((0x7ff7000000000000, 0x0000234db0000000)) == (false, fp_eflag_invalid));
  assert(float_is_eq((0x7ff0000000000001, 0xfff0000003000001)) == (false, fp_eflag_invalid));
+
+ /* Double floating point */
+ assert(float_is_eq((0x00000000000000000000000000000001, 0x00000000000000000000000000000001)) == (true, fp_eflag_none));
+ assert(float_is_eq((0x0f00000000000000000000000000000f, 0x0f00000000000000000000000000000f)) == (true, fp_eflag_none));
+ assert(float_is_eq((0x80000000000000000000000000000000, 0x00000000000000000000000000000000)) == (true, fp_eflag_none));
+ assert(float_is_eq((0x7fff0000000000000000000000000000, 0x7fff0000000000000000000000000000)) == (true, fp_eflag_none));
+
+ assert(float_is_eq((0x00000000000000000000000000000001, 0x80000000000000000000000000000001)) == (false, fp_eflag_none));
+ assert(float_is_eq((0x0f00000000000000000000000000000f, 0x3f00000000000000000000000000000f)) == (false, fp_eflag_none));
+ assert(float_is_eq((0x7fff8000000000000000000000000000, 0x0000234db00000000000000000000000)) == (false, fp_eflag_none));
+ assert(float_is_eq((0x7fff8000000000000000000000000000, 0xffff8000000000000000000000000000)) == (false, fp_eflag_none));
+
+ assert(float_is_eq((0x7fff7000000000000000000000000000, 0x0000234db00000000000000000000000)) == (false, fp_eflag_invalid));
+ assert(float_is_eq((0x7fff0000000000000000000000000001, 0xfff00000030000000000000000000001)) == (false, fp_eflag_invalid));
+}
 }
 
 function main () -> unit = {

test/float/inf_test.sail

@@ -65,6 +65,17 @@ function test_float_is_inf () -> unit = {
  assert(float_is_inf(0xfff8000000000000) == false);
  assert(float_is_inf(0xfff0000000000001) == false);
  assert(float_is_inf(0xffc0000000000000) == false);
+
+ /* Quad floating point */
+ assert(float_is_inf(0x7fff0000000000000000000000000000));
+ assert(float_is_inf(0xffff0000000000000000000000000000));
+
+ assert(float_is_inf(0x7fff8000000000000000000000000000) == false);
+ assert(float_is_inf(0x7fff0000000000000000000000000001) == false);
+ assert(float_is_inf(0x7ffc0000000000000000000000000000) == false);
+ assert(float_is_inf(0xffff8000000000000000000000000000) == false);
+ assert(float_is_inf(0xffff0000000000000000000000000001) == false);
+ assert(float_is_inf(0xfffc0000000000000000000000000000) == false);
 }
 
 function main () -> unit = {

test/float/nan_test.sail

@@ -86,6 +86,24 @@ function test_float_is_nan () -> unit = {
  assert(float_is_nan(0xfff0000000000000) == false);
  assert(float_is_nan(0xffe0000000000000) == false);
  assert(float_is_nan(0xffc0000000000000) == false);
+
+ /* Quad floating point */
+ assert(float_is_nan(0x7fff8000000000000000000000000000));
+ assert(float_is_nan(0x7fff0000000000000000000000000001));
+ assert(float_is_nan(0x7fff8000000000000000000000000100));
+ assert(float_is_nan(0x7fffc000000000001000000000000000));
+
+ assert(float_is_nan(0xffff8000000000000000000000000000));
+ assert(float_is_nan(0xffff0000000000000000000000000001));
+ assert(float_is_nan(0xffff8000000000000000000000000100));
+ assert(float_is_nan(0xffffc000000000001000000000000000));
+
+ assert(float_is_nan(0x7fff0000000000000000000000000000) == false);
+ assert(float_is_nan(0x7ffe0000000000000000000000000000) == false);
+ assert(float_is_nan(0x7ffc0000000000000000000000000000) == false);
+ assert(float_is_nan(0xffff0000000000000000000000000000) == false);
+ assert(float_is_nan(0xfffe0000000000000000000000000000) == false);
+ assert(float_is_nan(0xfffc0000000000000000000000000000) == false);
 }
 
 function test_float_is_snan () -> unit = {

test/float/normal_test.sail

@@ -65,6 +65,17 @@ function test_float_is_normal () -> unit = {
  assert(float_is_normal(0x7ff0000000000000) == false);
  assert(float_is_normal(0x0008000000000000) == false);
  assert(float_is_normal(0x8008000000000000) == false);
+
+ /* Quad floating point */
+ assert(float_is_normal(0x7ffe0000000000000000000000000000));
+ assert(float_is_normal(0x7ffe0000000000000000000000000001));
+ assert(float_is_normal(0xfffc000000000000000000000000000f));
+ assert(float_is_normal(0x80030000000000000001000000000000));
+
+ assert(float_is_normal(0x7fff8000000000000000000000000000) == false);
+ assert(float_is_normal(0x7fff0000000000000000000000000000) == false);
+ assert(float_is_normal(0x00008000000000000000000000000000) == false);
+ assert(float_is_normal(0x80008000000000000000000000000000) == false);
 }
 
 function test_float_is_denormal () -> unit = {
@@ -100,6 +111,17 @@ function test_float_is_denormal () -> unit = {
  assert(float_is_denormal(0x7ff0000000000000) == false);
  assert(float_is_denormal(0xffc000000000000f) == false);
  assert(float_is_denormal(0x8030000000100000) == false);
+
+ /* Quad floating point */
+ assert(float_is_denormal(0x00008000000000000000000000000000));
+ assert(float_is_denormal(0x80008000000000000000000000000000));
+ assert(float_is_denormal(0x80000000000000000000000000000001));
+ assert(float_is_denormal(0x80008000000000000001000000000000));
+
+ assert(float_is_denormal(0x7fff8000000000000000000000000000) == false);
+ assert(float_is_denormal(0x7fff0000000000000000000000000000) == false);
+ assert(float_is_denormal(0xfffc000000000000000000000000000f) == false);
+ assert(float_is_denormal(0x80030000000000000001000000000000) == false);
 }
 
 function main () -> unit = {

test/float/sign_test.sail

@@ -62,6 +62,16 @@ function test_float_is_positive () -> unit = {
  assert(float_is_positive(0xf8f0000000023000) == false);
  assert(float_is_positive(0xfce0000000000000) == false);
  assert(float_is_positive(0xe8c0000000000000) == false);
+
+ /* Quad floating point */
+ assert(float_is_positive(0x7ff080000f0f00000000000000000000));
+ assert(float_is_positive(0x0000000000000fedc000000000000001));
+ assert(float_is_positive(0x7fffffffffffffffffffffffffffffff));
+ assert(float_is_positive(0x00000000003000000000010000000000));
+
+ assert(float_is_positive(0xf8ff0000000000000023000000000000) == false);
+ assert(float_is_positive(0xffce0000000000000000000000000000) == false);
+ assert(float_is_positive(0xe88c0000000000000000000000000000) == false);
 }
 
 function test_float_is_negative () -> unit = {
@@ -94,6 +104,16 @@ function test_float_is_negative () -> unit = {
  assert(float_is_negative(0x78f0000000023000) == false);
  assert(float_is_negative(0x7ce0000000000000) == false);
  assert(float_is_negative(0x68c0000000000000) == false);
+
+ /* Quad floating point */
+ assert(float_is_negative(0xfff0000000800000000f0f0000000000));
+ assert(float_is_negative(0x8000000000000fedc000000000000001));
+ assert(float_is_negative(0xffffffffffffffffffffffffffffffff));
+ assert(float_is_negative(0x80000000003000000000010000000000));
+
+ assert(float_is_negative(0x78ff0000000000000000023000000000) == false);
+ assert(float_is_negative(0x7fce0000000000000000000000000000) == false);
+ assert(float_is_negative(0x688c0000000000000000000000000000) == false);
 }
 
 function main () -> unit = {

test/float/zero_test.sail

@@ -59,6 +59,15 @@ function test_float_is_zero () -> unit = {
  assert(float_is_zero(0x7ff0000000000001) == false);
  assert(float_is_zero(0xfff8000000000000) == false);
  assert(float_is_zero(0xfff0000000000001) == false);
+
+ /* Quad floating point */
+ assert(float_is_zero(0x80000000000000000000000000000000));
+ assert(float_is_zero(0x00000000000000000000000000000000));
+
+ assert(float_is_zero(0x7fff8000000000000000000000000000) == false);
+ assert(float_is_zero(0x7fff0000000000000000000000000001) == false);
+ assert(float_is_zero(0xffff8000000000000000000000000000) == false);
+ assert(float_is_zero(0xffff0000000000000000000000000001) == false);
 }
 
 function main () -> unit = {