[Mesa-dev] [PATCH v3 09/44] util: add fp64 -> fp32 conversion support for RTNE and RTZ rounding modes

Samuel Iglesias Gonsálvez siglesias at igalia.com
Wed Feb 6 10:44:38 UTC 2019


Signed-off-by: Samuel Iglesias Gonsálvez <siglesias at igalia.com>
---
 src/util/Makefile.sources |   2 +
 src/util/double.c         | 197 ++++++++++++++++++++++++++++++++++++++
 src/util/double.h         |  46 +++++++++
 src/util/meson.build      |   2 +
 4 files changed, 247 insertions(+)
 create mode 100644 src/util/double.c
 create mode 100644 src/util/double.h

diff --git a/src/util/Makefile.sources b/src/util/Makefile.sources
index f09b89b3be5..1b998bf26c4 100644
--- a/src/util/Makefile.sources
+++ b/src/util/Makefile.sources
@@ -11,6 +11,8 @@ MESA_UTIL_FILES := \
 	debug.h \
 	disk_cache.c \
 	disk_cache.h \
+	double.c \
+	double.h \
 	fast_idiv_by_const.c \
 	fast_idiv_by_const.h \
 	format_r11g11b10f.h \
diff --git a/src/util/double.c b/src/util/double.c
new file mode 100644
index 00000000000..61f23e3b57b
--- /dev/null
+++ b/src/util/double.c
@@ -0,0 +1,197 @@
+/*
+ * Mesa 3-D graphics library
+ *
+ * Copyright (C) 2018 Intel Corporation
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining a
+ * copy of this software and associated documentation files (the "Software"),
+ * to deal in the Software without restriction, including without limitation
+ * the rights to use, copy, modify, merge, publish, distribute, sublicense,
+ * and/or sell copies of the Software, and to permit persons to whom the
+ * Software is furnished to do so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included
+ * in all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
+ * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.  IN NO EVENT SHALL
+ * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR
+ * OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,
+ * ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
+ * OTHER DEALINGS IN THE SOFTWARE.
+ */
+
+#include <stdio.h>
+#include <assert.h>
+#include <math.h>
+#include <string.h>
+
+#include "rounding.h"
+#include "double.h"
+
+
+
+
+typedef union { double f; int64_t i; uint64_t u; } fi_type;
+
+/**
+ * Convert a 8-byte double to a 4-byte float.
+ *
+ * Not all float64 values can be represented exactly as a float32 value. We
+ * round such intermediate float64 values to the nearest float32. When the
+ * float64 lies exactly between two float32 values, we round to the one with
+ * an even mantissa.
+ */
+
+float
+_mesa_double_to_float(double val)
+{
+   const fi_type fi = {val};
+   const int64_t flt_m = fi.i & 0x0fffffffffffff;
+   const int64_t flt_e = (fi.i >> 52) & 0x7ff;
+   const int64_t flt_s = (fi.i >> 63) & 0x1;
+   int s, e, m = 0;
+   float result;
+
+   /* sign bit */
+   s = flt_s;
+
+   /* handle special cases */
+   if ((flt_e == 0) && (flt_m == 0)) {
+      /* zero */
+      /* m = 0; - already set */
+      e = 0;
+   }
+   else if ((flt_e == 0) && (flt_m != 0)) {
+      /* denorm -- denorm float64 maps to 0 */
+      /* m = 0; - already set */
+      e = 0;
+   }
+   else if ((flt_e == 0x7ff) && (flt_m == 0)) {
+      /* infinity */
+      /* m = 0; - already set */
+      e = 255;
+   }
+   else if ((flt_e == 0x7ff) && (flt_m != 0)) {
+      /* NaN */
+      m = 1;
+      e = 255;
+   }
+   else {
+      /* regular number */
+      const int new_exp = flt_e - 1023;
+      if (new_exp < -126) {
+         /* The float64 lies in the range (0.0, min_normal32) and is rounded
+          * to a nearby float32 value. The result will be either zero, subnormal,
+          * or normal.
+          */
+         e = 0;
+         m = _mesa_lroundeven(((double)((uint64_t)1 << 54)) * fabs(fi.f));
+      }
+      else if (new_exp > 127) {
+         /* map this value to infinity */
+         /* m = 0; - already set */
+         e = 255;
+      }
+      else {
+         /* The float64 lies in the range
+          *   [min_normal32, max_normal32 + max_step32)
+          * and is rounded to a nearby float32 value. The result will be
+          * either normal or infinite.
+          */
+         e = new_exp + 127;
+         m = _mesa_lroundeven((double)flt_m / (double) (1 << 29));
+      }
+   }
+
+   assert(0 <= m && m <= (1 << 23));
+   if (m == (1 << 23)) {
+      /* The float64 was rounded upwards into the range of the next exponent,
+       * so bump the exponent. This correctly handles the case where f64
+       * should be rounded up to float32 infinity.
+       */
+      ++e;
+      m = 0;
+   }
+
+   unsigned result_int = (s << 31) | (e << 23) | m;
+   memcpy(&result, &result_int, sizeof(float));
+   return result;
+}
+
+float
+_mesa_double_to_float_rtz(double val)
+{
+   const fi_type fi = {val};
+   const int64_t flt_m = fi.i & 0x0fffffffffffff;
+   const int64_t flt_e = (fi.i >> 52) & 0x7ff;
+   const int64_t flt_s = (fi.i >> 63) & 0x1;
+   int s, e, m = 0;
+   float result;
+
+   /* sign bit */
+   s = flt_s;
+
+   /* handle special cases */
+   if ((flt_e == 0) && (flt_m == 0)) {
+      /* zero */
+      /* m = 0; - already set */
+      e = 0;
+   }
+   else if ((flt_e == 0) && (flt_m != 0)) {
+      /* denorm -- denorm float64 maps to 0 */
+      /* m = 0; - already set */
+      e = 0;
+   }
+   else if ((flt_e == 0xff) && (flt_m == 0)) {
+      /* infinity */
+      /* m = 0; - already set */
+      e = 255;
+   }
+   else if ((flt_e == 0xff) && (flt_m != 0)) {
+      /* NaN */
+      m = 1;
+      e = 255;
+   }
+   else {
+      /* regular number */
+      const int new_exp = flt_e - 1023;
+      if (new_exp < -126) {
+         /* The float64 lies in the range (0.0, min_normal32) and is rounded
+          * to a nearby float32 value. The result will be either zero, subnormal,
+          * or normal.
+          */
+         e = 0;
+         m = trunc((double)((uint64_t)1 << 54) * fabs(fi.f));
+      }
+      else if (new_exp > 127) {
+         /* map this value to infinity */
+         /* m = 0; - already set */
+         e = 255;
+      }
+      else {
+         /* The float64 lies in the range
+          *   [min_normal32, max_normal32 + max_step32)
+          * and is rounded to a nearby float32 value. The result will be
+          * either normal or infinite.
+          */
+         e = new_exp + 127;
+         m = trunc((double)flt_m / (double) (1 << 29));
+      }
+   }
+
+   assert(0 <= m && m <= (1 << 23));
+   if (m == (1 << 23)) {
+      /* The float64 was rounded upwards into the range of the next exponent,
+       * so bump the exponent. This correctly handles the case where f64
+       * should be rounded up to float32 infinity.
+       */
+      ++e;
+      m = 0;
+   }
+
+   unsigned result_int = (s << 31) | (e << 23) | m;
+   memcpy(&result, &result_int, sizeof(float));
+   return result;
+}
diff --git a/src/util/double.h b/src/util/double.h
new file mode 100644
index 00000000000..0314c8ea0b5
--- /dev/null
+++ b/src/util/double.h
@@ -0,0 +1,46 @@
+/*
+ * Mesa 3-D graphics library
+ *
+ * Copyright (C) 2018 Intel Corporation
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining a
+ * copy of this software and associated documentation files (the "Software"),
+ * to deal in the Software without restriction, including without limitation
+ * the rights to use, copy, modify, merge, publish, distribute, sublicense,
+ * and/or sell copies of the Software, and to permit persons to whom the
+ * Software is furnished to do so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included
+ * in all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
+ * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.  IN NO EVENT SHALL
+ * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR
+ * OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,
+ * ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
+ * OTHER DEALINGS IN THE SOFTWARE.
+ */
+
+#ifndef _DOUBLE_H_
+#define _DOUBLE_H_
+
+
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+float _mesa_double_to_float(double val);
+float _mesa_double_to_float_rtz(double val);
+
+static inline float
+_mesa_double_to_float_rtne(double val)
+{
+   return _mesa_double_to_float(val);
+}
+
+#ifdef __cplusplus
+} /* extern C */
+#endif
+
+#endif /* _DOUBLE_H_ */
diff --git a/src/util/meson.build b/src/util/meson.build
index 397c2228129..23371ef2742 100644
--- a/src/util/meson.build
+++ b/src/util/meson.build
@@ -35,6 +35,8 @@ files_mesa_util = files(
   'debug.h',
   'disk_cache.c',
   'disk_cache.h',
+  'double.c',
+  'double.h',
   'fast_idiv_by_const.c',
   'fast_idiv_by_const.h',
   'format_r11g11b10f.h',
-- 
2.19.1



More information about the mesa-dev mailing list