[PPL-devel] [GIT] ppl/ppl(floating_point): Implemented compute_absolute_error.

[Fabio Bossi]

Module: ppl/ppl
Branch: floating_point
Commit: 95846dca5e7fe0afec95c2489192e4f3eab81f61
URL:    http://www.cs.unipr.it/git/gitweb.cgi?p=ppl/ppl.git;a=commit;h=95846dca5e7fe0afec95c2489192e4f3eab81f61

Author: Fabio Bossi <bossi at cs.unipr.it>
Date:   Wed Jul 21 10:54:28 2010 +0200

Implemented compute_absolute_error.

---

 src/Float.defs.hh      |   12 +++++++
 src/Float.templates.hh |   87 ++++++++++++++++++++++++++++++++++++++++++++++++
 src/Makefile.am        |    1 +
 3 files changed, 100 insertions(+), 0 deletions(-)

diff --git a/src/Float.defs.hh b/src/Float.defs.hh
index 1e3e413..fff7b68 100644
--- a/src/Float.defs.hh
+++ b/src/Float.defs.hh
@@ -23,6 +23,7 @@ site: http://www.cs.unipr.it/ppl/ . */
 #ifndef PPL_Float_defs_hh
 #define PPL_Float_defs_hh 1
 
+#include "globals.types.hh"
 #include "meta_programming.hh"
 #include "compiler.hh"
 #include <gmp.h>
@@ -359,8 +360,19 @@ public:
 };
 #endif
 
+template <typename FP_Interval_Type>
+FP_Interval_Type compute_absolute_error(const Floating_Point_Format
+                                        analyzed_format);
+
+// TODO.
+/*
+template <typename FP_Interval_Type>
+FP_Interval_Type absolute_error = compute_absolute_error();
+*/
+
 } // namespace Parma_Polyhedra_Library
 
 #include "Float.inlines.hh"
+#include "Float.templates.hh"
 
 #endif // !defined(PPL_Float_defs_hh)
diff --git a/src/Float.templates.hh b/src/Float.templates.hh
new file mode 100644
index 0000000..5cc53fc
--- /dev/null
+++ b/src/Float.templates.hh
@@ -0,0 +1,87 @@
+/* IEC 559 floating point format related functions:
+   non-inline template functions.
+   Copyright (C) 2001-2010 Roberto Bagnara <bagnara at cs.unipr.it>
+
+This file is part of the Parma Polyhedra Library (PPL).
+
+The PPL is free software; you can redistribute it and/or modify it
+under the terms of the GNU General Public License as published by the
+Free Software Foundation; either version 3 of the License, or (at your
+option) any later version.
+
+The PPL is distributed in the hope that it will be useful, but WITHOUT
+ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+for more details.
+
+You should have received a copy of the GNU General Public License
+along with this program; if not, write to the Free Software Foundation,
+Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02111-1307, USA.
+
+For the most up-to-date information see the Parma Polyhedra Library
+site: http://www.cs.unipr.it/ppl/ . */
+
+#ifndef PPL_Float_templates_hh
+#define PPL_Float_Templates_hh 1
+
+#include <cmath>
+
+namespace Parma_Polyhedra_Library {
+
+template <typename FP_Interval_Type>
+FP_Interval_Type compute_absolute_error(
+		 const Floating_Point_Format analyzed_format) {
+  typedef typename FP_Interval_Type::boundary_type analyzer_format;
+
+  // Get the necessary information on the analyzed's format.
+  unsigned int f_base;
+  int f_exponent_bias;
+  unsigned int f_mantissa_bits;
+  switch (analyzed_format) {
+    case IEEE754_HALF:
+      f_base = float_ieee754_half::BASE;
+      f_exponent_bias = float_ieee754_half::EXPONENT_BIAS;
+      f_mantissa_bits = float_ieee754_half::MANTISSA_BITS;
+      break;
+    case IEEE754_SINGLE:
+      f_base = float_ieee754_single::BASE;
+      f_exponent_bias = float_ieee754_single::EXPONENT_BIAS;
+      f_mantissa_bits = float_ieee754_single::MANTISSA_BITS;
+      break;
+    case IEEE754_DOUBLE:
+      f_base = float_ieee754_double::BASE;
+      f_exponent_bias = float_ieee754_double::EXPONENT_BIAS;
+      f_mantissa_bits = float_ieee754_double::MANTISSA_BITS;
+      break;
+    case IBM_SINGLE:
+      f_base = float_ibm_single::BASE;
+      f_exponent_bias = float_ibm_single::EXPONENT_BIAS;
+      f_mantissa_bits = float_ibm_single::MANTISSA_BITS;
+      break;
+    case IEEE754_QUAD:
+      f_base = float_ieee754_quad::BASE;
+      f_exponent_bias = float_ieee754_quad::EXPONENT_BIAS;
+      f_mantissa_bits = float_ieee754_quad::MANTISSA_BITS;
+      break;
+    case INTEL_DOUBLE_EXTENDED:
+      f_base = float_intel_double_extended::BASE;
+      f_exponent_bias = float_intel_double_extended::EXPONENT_BIAS;
+      f_mantissa_bits = float_intel_double_extended::MANTISSA_BITS;
+      break;
+  }
+
+  analyzer_format omega = std::max(
+  static_cast<analyzer_format>(pow(f_base,
+                               static_cast<analyzer_format>(1) -
+                               f_exponent_bias - f_mantissa_bits)),
+  std::numeric_limits<analyzer_format>::denorm_min());
+
+  FP_Interval_Type result;
+  result.build(i_constraint(GREATER_OR_EQUAL, -omega),
+               i_constraint(LESS_OR_EQUAL, omega));
+  return result;
+}
+
+} // namespace Parma_Polyhedra_Library
+
+#endif // !defined(PPL_Float_Templates_hh)
diff --git a/src/Makefile.am b/src/Makefile.am
index 62efcad..e4dbcb7 100644
--- a/src/Makefile.am
+++ b/src/Makefile.am
@@ -182,6 +182,7 @@ Octagonal_Shape.types.hh \
 fpu.types.hh \
 Float.defs.hh \
 Float.inlines.hh \
+Float.templates.hh \
 C_Integer.hh \
 fpu.defs.hh \
 fpu-c99.inlines.hh \