July 2010 - PPL-devel - BUGSENG Mailman

[GIT] ppl/ppl(floating_point): Removed two probable unnecessary inclusions.
by Fabio Bossi 22 Jul '10

22 Jul '10

Module: ppl/ppl Branch: floating_point Commit: 598a0f5c035c8c0ebdb93d7a390b21c6500a8263 URL: http://www.cs.unipr.it/git/gitweb.cgi?p=ppl/ppl.git;a=commit;h=598a0f5c035c… Author: Fabio Bossi <bossi(a)cs.unipr.it> Date: Thu Jul 22 11:30:13 2010 +0200 Removed two probable unnecessary inclusions. --- src/Grid.defs.hh | 2 -- 1 files changed, 0 insertions(+), 2 deletions(-) diff --git a/src/Grid.defs.hh b/src/Grid.defs.hh index c51ffd0..bf4e617 100644 --- a/src/Grid.defs.hh +++ b/src/Grid.defs.hh @@ -40,9 +40,7 @@ site: http://www.cs.unipr.it/ppl/ . */ #include "Poly_Gen_Relation.defs.hh" #include "Grid_Certificate.types.hh" #include "Box.types.hh" -#include "Polyhedron.defs.hh" #include "Polyhedron.types.hh" -#include "Polyhedron.inlines.hh" #include "BD_Shape.types.hh" #include "Octagonal_Shape.types.hh" #include <vector>

1 0

[GIT] ppl/ppl(floating_point): Fixed a typo in the inclusion guard.
by Fabio Bossi 22 Jul '10

22 Jul '10

Module: ppl/ppl Branch: floating_point Commit: cc3822c311bf9c585c8a73c0bbdfed071de8277b URL: http://www.cs.unipr.it/git/gitweb.cgi?p=ppl/ppl.git;a=commit;h=cc3822c311bf… Author: Fabio Bossi <bossi(a)cs.unipr.it> Date: Thu Jul 22 10:29:16 2010 +0200 Fixed a typo in the inclusion guard. --- src/Float.templates.hh | 2 +- 1 files changed, 1 insertions(+), 1 deletions(-) diff --git a/src/Float.templates.hh b/src/Float.templates.hh index 4912b9a..139a153 100644 --- a/src/Float.templates.hh +++ b/src/Float.templates.hh @@ -21,7 +21,7 @@ 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 +#ifndef PPL_Float_Templates_hh #define PPL_Float_Templates_hh 1 #include <cmath>

1 0

[GIT] ppl/ppl(floating_point): Removed class absolute_error.
by Fabio Bossi 21 Jul '10

21 Jul '10

Module: ppl/ppl Branch: floating_point Commit: 0decead53d0c887298b1dc309e7003d2ce59698f URL: http://www.cs.unipr.it/git/gitweb.cgi?p=ppl/ppl.git;a=commit;h=0decead53d0c… Author: Fabio Bossi <bossi(a)cs.unipr.it> Date: Wed Jul 21 17:18:30 2010 +0200 Removed class absolute_error. Cache results locally in compute_absolute_error instead. --- src/Float.defs.hh | 60 +---------------------------------------------- src/Float.templates.hh | 47 ++++++++++++++++++++++++++++++++---- 2 files changed, 43 insertions(+), 64 deletions(-) diff --git a/src/Float.defs.hh b/src/Float.defs.hh index 43a0399..e68c3bd 100644 --- a/src/Float.defs.hh +++ b/src/Float.defs.hh @@ -361,64 +361,8 @@ public: #endif template <typename FP_Interval_Type> -FP_Interval_Type compute_absolute_error(Floating_Point_Format analyzed_format); - -template <typename FP_Interval_Type, Floating_Point_Format> -struct absolute_error; - -template <typename FP_Interval_Type> -struct absolute_error<FP_Interval_Type, IEEE754_HALF> { - static const FP_Interval_Type value; -}; - -template <typename FP_Interval_Type> -const FP_Interval_Type absolute_error<FP_Interval_Type, IEEE754_HALF>::value = - compute_absolute_error<FP_Interval_Type>(IEEE754_HALF); - -template <typename FP_Interval_Type> -struct absolute_error<FP_Interval_Type, IEEE754_SINGLE> { - static const FP_Interval_Type value; -}; - -template <typename FP_Interval_Type> -const FP_Interval_Type absolute_error<FP_Interval_Type, IEEE754_SINGLE>::value = - compute_absolute_error<FP_Interval_Type>(IEEE754_SINGLE); - -template <typename FP_Interval_Type> -struct absolute_error<FP_Interval_Type, IEEE754_DOUBLE> { - static const FP_Interval_Type value; -}; - -template <typename FP_Interval_Type> -const FP_Interval_Type absolute_error<FP_Interval_Type, IEEE754_DOUBLE>::value = - compute_absolute_error<FP_Interval_Type>(IEEE754_DOUBLE); - -template <typename FP_Interval_Type> -struct absolute_error<FP_Interval_Type, IBM_SINGLE> { - static const FP_Interval_Type value; -}; - -template <typename FP_Interval_Type> -const FP_Interval_Type absolute_error<FP_Interval_Type, IBM_SINGLE>::value = - compute_absolute_error<FP_Interval_Type>(IBM_SINGLE); - -template <typename FP_Interval_Type> -struct absolute_error<FP_Interval_Type, IEEE754_QUAD> { - static const FP_Interval_Type value; -}; - -template <typename FP_Interval_Type> -const FP_Interval_Type absolute_error<FP_Interval_Type, IEEE754_QUAD>::value = - compute_absolute_error<FP_Interval_Type>(IEEE754_QUAD); - -template <typename FP_Interval_Type> -struct absolute_error<FP_Interval_Type, INTEL_DOUBLE_EXTENDED> { - static const FP_Interval_Type value; -}; - -template <typename FP_Interval_Type> -const FP_Interval_Type absolute_error<FP_Interval_Type, INTEL_DOUBLE_EXTENDED>::value = - compute_absolute_error<FP_Interval_Type>(INTEL_DOUBLE_EXTENDED); +const FP_Interval_Type& compute_absolute_error( + Floating_Point_Format analyzed_format); } // namespace Parma_Polyhedra_Library diff --git a/src/Float.templates.hh b/src/Float.templates.hh index feaa8bb..4912b9a 100644 --- a/src/Float.templates.hh +++ b/src/Float.templates.hh @@ -29,41 +29,75 @@ site: http://www.cs.unipr.it/ppl/ . */ namespace Parma_Polyhedra_Library { template <typename FP_Interval_Type> -FP_Interval_Type compute_absolute_error( - const Floating_Point_Format analyzed_format) { +const FP_Interval_Type& compute_absolute_error( + const Floating_Point_Format analyzed_format) { typedef typename FP_Interval_Type::boundary_type analyzer_format; + static const FP_Interval_Type ZERO_INTERVAL = FP_Interval_Type(0); + // Cached results for each different analyzed format. + static FP_Interval_Type ieee754_half_result = ZERO_INTERVAL; + static FP_Interval_Type ieee754_single_result = ZERO_INTERVAL; + static FP_Interval_Type ieee754_double_result = ZERO_INTERVAL; + static FP_Interval_Type ibm_single_result = ZERO_INTERVAL; + static FP_Interval_Type ieee754_quad_result = ZERO_INTERVAL; + static FP_Interval_Type intel_double_extended_result = ZERO_INTERVAL; + + FP_Interval_Type* to_compute = NULL; // 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: + if (ieee754_half_result != ZERO_INTERVAL) + return ieee754_half_result; + + to_compute = &ieee754_half_result; 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: + if (ieee754_single_result != ZERO_INTERVAL) + return ieee754_single_result; + + to_compute = &ieee754_single_result; 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: + if (ieee754_double_result != ZERO_INTERVAL) + return ieee754_double_result; + + to_compute = &ieee754_double_result; 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: + if (ibm_single_result != ZERO_INTERVAL) + return ibm_single_result; + + to_compute = &ibm_single_result; 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: + if (ieee754_quad_result != ZERO_INTERVAL) + return ieee754_quad_result; + + to_compute = &ieee754_quad_result; 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: + if (intel_double_extended_result != ZERO_INTERVAL) + return intel_double_extended_result; + + to_compute = &intel_double_extended_result; 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; @@ -72,16 +106,17 @@ FP_Interval_Type compute_absolute_error( throw std::runtime_error("PPL internal error"); } + PPL_ASSERT(to_compute != NULL); + 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; + to_compute->build(i_constraint(GREATER_OR_EQUAL, -omega), + i_constraint(LESS_OR_EQUAL, omega)); + return *to_compute; } } // namespace Parma_Polyhedra_Library

1 0

[GIT] ppl/ppl(floating_point): Correctly handle unexpected formats.
by Fabio Bossi 21 Jul '10

21 Jul '10

Module: ppl/ppl Branch: floating_point Commit: d7bfff5fe6a73fd3496b5c6af2cd1864fab1e4b9 URL: http://www.cs.unipr.it/git/gitweb.cgi?p=ppl/ppl.git;a=commit;h=d7bfff5fe6a7… Author: Fabio Bossi <bossi(a)cs.unipr.it> Date: Wed Jul 21 15:48:54 2010 +0200 Correctly handle unexpected formats. --- src/Float.templates.hh | 2 ++ src/Linear_Form.templates.hh | 2 ++ 2 files changed, 4 insertions(+), 0 deletions(-) diff --git a/src/Float.templates.hh b/src/Float.templates.hh index 5cc53fc..feaa8bb 100644 --- a/src/Float.templates.hh +++ b/src/Float.templates.hh @@ -68,6 +68,8 @@ FP_Interval_Type compute_absolute_error( f_exponent_bias = float_intel_double_extended::EXPONENT_BIAS; f_mantissa_bits = float_intel_double_extended::MANTISSA_BITS; break; + default: + throw std::runtime_error("PPL internal error"); } analyzer_format omega = std::max( diff --git a/src/Linear_Form.templates.hh b/src/Linear_Form.templates.hh index 666207e..f676626 100644 --- a/src/Linear_Form.templates.hh +++ b/src/Linear_Form.templates.hh @@ -412,6 +412,8 @@ Linear_Form<C>::relative_error( f_base = float_intel_double_extended::BASE; f_mantissa_bits = float_intel_double_extended::MANTISSA_BITS; break; + default: + throw std::runtime_error("PPL internal error"); } C error_propagator;

1 0

[GIT] ppl/ppl(floating_point): Made class absolute_error a struct.
by Fabio Bossi 21 Jul '10

21 Jul '10

Module: ppl/ppl Branch: floating_point Commit: 5f7ad90bdec26fb1710939a6887bf0b89277f91f URL: http://www.cs.unipr.it/git/gitweb.cgi?p=ppl/ppl.git;a=commit;h=5f7ad90bdec2… Author: Fabio Bossi <bossi(a)cs.unipr.it> Date: Wed Jul 21 15:28:22 2010 +0200 Made class absolute_error a struct. --- src/Float.defs.hh | 14 +++++++------- 1 files changed, 7 insertions(+), 7 deletions(-) diff --git a/src/Float.defs.hh b/src/Float.defs.hh index 21270da..43a0399 100644 --- a/src/Float.defs.hh +++ b/src/Float.defs.hh @@ -364,10 +364,10 @@ template <typename FP_Interval_Type> FP_Interval_Type compute_absolute_error(Floating_Point_Format analyzed_format); template <typename FP_Interval_Type, Floating_Point_Format> -class absolute_error; +struct absolute_error; template <typename FP_Interval_Type> -class absolute_error<FP_Interval_Type, IEEE754_HALF> { +struct absolute_error<FP_Interval_Type, IEEE754_HALF> { static const FP_Interval_Type value; }; @@ -376,7 +376,7 @@ const FP_Interval_Type absolute_error<FP_Interval_Type, IEEE754_HALF>::value = compute_absolute_error<FP_Interval_Type>(IEEE754_HALF); template <typename FP_Interval_Type> -class absolute_error<FP_Interval_Type, IEEE754_SINGLE> { +struct absolute_error<FP_Interval_Type, IEEE754_SINGLE> { static const FP_Interval_Type value; }; @@ -385,7 +385,7 @@ const FP_Interval_Type absolute_error<FP_Interval_Type, IEEE754_SINGLE>::value = compute_absolute_error<FP_Interval_Type>(IEEE754_SINGLE); template <typename FP_Interval_Type> -class absolute_error<FP_Interval_Type, IEEE754_DOUBLE> { +struct absolute_error<FP_Interval_Type, IEEE754_DOUBLE> { static const FP_Interval_Type value; }; @@ -394,7 +394,7 @@ const FP_Interval_Type absolute_error<FP_Interval_Type, IEEE754_DOUBLE>::value = compute_absolute_error<FP_Interval_Type>(IEEE754_DOUBLE); template <typename FP_Interval_Type> -class absolute_error<FP_Interval_Type, IBM_SINGLE> { +struct absolute_error<FP_Interval_Type, IBM_SINGLE> { static const FP_Interval_Type value; }; @@ -403,7 +403,7 @@ const FP_Interval_Type absolute_error<FP_Interval_Type, IBM_SINGLE>::value = compute_absolute_error<FP_Interval_Type>(IBM_SINGLE); template <typename FP_Interval_Type> -class absolute_error<FP_Interval_Type, IEEE754_QUAD> { +struct absolute_error<FP_Interval_Type, IEEE754_QUAD> { static const FP_Interval_Type value; }; @@ -412,7 +412,7 @@ const FP_Interval_Type absolute_error<FP_Interval_Type, IEEE754_QUAD>::value = compute_absolute_error<FP_Interval_Type>(IEEE754_QUAD); template <typename FP_Interval_Type> -class absolute_error<FP_Interval_Type, INTEL_DOUBLE_EXTENDED> { +struct absolute_error<FP_Interval_Type, INTEL_DOUBLE_EXTENDED> { static const FP_Interval_Type value; };

1 0

[GIT] ppl/ppl(floating_point): Implemented Linear_Form::intervalize ( still not compiling,
by Fabio Bossi 21 Jul '10

21 Jul '10

Module: ppl/ppl Branch: floating_point Commit: 4c1245f7e89c21d15863001bbc4e5bd9a3cce348 URL: http://www.cs.unipr.it/git/gitweb.cgi?p=ppl/ppl.git;a=commit;h=4c1245f7e89c… Author: Fabio Bossi <bossi(a)cs.unipr.it> Date: Wed Jul 21 15:18:30 2010 +0200 Implemented Linear_Form::intervalize (still not compiling, probably due to some sort of cyclic dependency problem between Grid and Box). --- src/Linear_Form.defs.hh | 3 +++ src/Linear_Form.templates.hh | 17 +++++++++++++++++ 2 files changed, 20 insertions(+), 0 deletions(-) diff --git a/src/Linear_Form.defs.hh b/src/Linear_Form.defs.hh index c80b605..f7ed492 100644 --- a/src/Linear_Form.defs.hh +++ b/src/Linear_Form.defs.hh @@ -26,6 +26,7 @@ site: http://www.cs.unipr.it/ppl/ . */ #include "Linear_Form.types.hh" #include "Linear_Expression.types.hh" #include "Variable.defs.hh" +#include "Box.types.hh" #include <vector> namespace Parma_Polyhedra_Library { @@ -329,6 +330,8 @@ public: void relative_error(Floating_Point_Format analyzed_format, Linear_Form& result) const; + void intervalize(const Box<C>& store, C& result) const; + private: //! The generic coefficient equal to the singleton zero. static C zero; diff --git a/src/Linear_Form.templates.hh b/src/Linear_Form.templates.hh index 83e7c1f..666207e 100644 --- a/src/Linear_Form.templates.hh +++ b/src/Linear_Form.templates.hh @@ -24,6 +24,7 @@ site: http://www.cs.unipr.it/ppl/ . */ #define PPL_Linear_Form_templates_hh 1 #include "Linear_Expression.defs.hh" +#include "Box.defs.hh" #include <stdexcept> #include <iostream> @@ -445,6 +446,22 @@ Linear_Form<C>::relative_error( return; } +template <typename C> +void +Linear_Form<C>::intervalize(const Box<C>& store, C& result) const { + result = C(inhomogeneous_term()); + dimension_type dimension = space_dimension(); + assert(dimension <= store.space_dimension()); + for (dimension_type i = 0; i < dimension; ++i) { + C current_addend = coefficient(Variable(i)); + const C& curr_int = store.get_interval(Variable(i)); + current_addend *= curr_int; + result += current_addend; + } + + return; +} + /*! \relates Parma_Polyhedra_Library::Linear_Form */ template <typename C> std::ostream&

1 0

[GIT] ppl/ppl(floating_point): Fixed iostream inclusion.
by Fabio Bossi 21 Jul '10

21 Jul '10

Module: ppl/ppl Branch: floating_point Commit: 62fcf8c0d4beb9d94c555f6c0310a4030b5f2312 URL: http://www.cs.unipr.it/git/gitweb.cgi?p=ppl/ppl.git;a=commit;h=62fcf8c0d4be… Author: Fabio Bossi <bossi(a)cs.unipr.it> Date: Wed Jul 21 14:35:04 2010 +0200 Fixed iostream inclusion. --- src/Concrete_Expression.cc | 4 +--- 1 files changed, 1 insertions(+), 3 deletions(-) diff --git a/src/Concrete_Expression.cc b/src/Concrete_Expression.cc index b06e439..e2416fd 100644 --- a/src/Concrete_Expression.cc +++ b/src/Concrete_Expression.cc @@ -21,9 +21,7 @@ For the most up-to-date information see the Parma Polyhedra Library site: http://www.cs.unipr.it/ppl/ . */ #include <ppl-config.h> -#ifndef NDEBUG -#include <stdio> -#endif +#include <iostream> #include "Concrete_Expression.defs.hh"

1 0

[GIT] ppl/ppl(floating_point): Fixed error in the use of NDEBUG.
by Fabio Bossi 21 Jul '10

21 Jul '10

Module: ppl/ppl Branch: floating_point Commit: df4645ba455fac92547b0bd050d43b20edeadfcb URL: http://www.cs.unipr.it/git/gitweb.cgi?p=ppl/ppl.git;a=commit;h=df4645ba455f… Author: Fabio Bossi <bossi(a)cs.unipr.it> Date: Wed Jul 21 14:31:34 2010 +0200 Fixed error in the use of NDEBUG. --- src/Concrete_Expression.cc | 16 ++++++++-------- 1 files changed, 8 insertions(+), 8 deletions(-) diff --git a/src/Concrete_Expression.cc b/src/Concrete_Expression.cc index cc7b34d..b06e439 100644 --- a/src/Concrete_Expression.cc +++ b/src/Concrete_Expression.cc @@ -21,7 +21,7 @@ For the most up-to-date information see the Parma Polyhedra Library site: http://www.cs.unipr.it/ppl/ . */ #include <ppl-config.h> -#ifdef NDEBUG +#ifndef NDEBUG #include <stdio> #endif @@ -37,7 +37,7 @@ PPL::Concrete_Expression_Type::OK() const { && impl.bounded_integer_type_width != BITS_32 && impl.bounded_integer_type_width != BITS_64 && impl.bounded_integer_type_width != BITS_128) { -#ifdef NDEBUG +#ifndef NDEBUG std::cerr << "bounded integer type has illegal width " << impl.bounded_integer_type_width << std::endl; @@ -47,7 +47,7 @@ PPL::Concrete_Expression_Type::OK() const { if (impl.bounded_integer_type_representation != UNSIGNED && impl.bounded_integer_type_representation != SIGNED_2_COMPLEMENT) { -#ifdef NDEBUG +#ifndef NDEBUG std::cerr << "bounded integer type has illegal representation " << impl.bounded_integer_type_representation << std::endl; @@ -58,7 +58,7 @@ PPL::Concrete_Expression_Type::OK() const { if (impl.bounded_integer_type_overflow != OVERFLOW_WRAPS && impl.bounded_integer_type_overflow != OVERFLOW_UNDEFINED && impl.bounded_integer_type_overflow != OVERFLOW_IMPOSSIBLE) { -#ifdef NDEBUG +#ifndef NDEBUG std::cerr << "bounded integer type has illegal overflow " << impl.bounded_integer_type_overflow << std::endl; @@ -67,7 +67,7 @@ PPL::Concrete_Expression_Type::OK() const { } if (impl.floating_point_format != IEEE754_HALF) { -#ifdef NDEBUG +#ifndef NDEBUG std::cerr << "bounded integer type has illegal (unused) fp format " << impl.floating_point_format << std::endl; @@ -83,7 +83,7 @@ PPL::Concrete_Expression_Type::OK() const { && impl.floating_point_format != IBM_SINGLE && impl.floating_point_format != IEEE754_QUAD && impl.floating_point_format != INTEL_DOUBLE_EXTENDED) { -#ifdef NDEBUG +#ifndef NDEBUG std::cerr << "floating point type has illegal format " << impl.floating_point_format << std::endl; @@ -101,7 +101,7 @@ PPL::Concrete_Expression_Type::OK() const { } if (impl.bounded_integer_type_representation != SIGNED_2_COMPLEMENT) { -#ifdef NDEBUG +#ifndef NDEBUG std::cerr << "floating point type has illegal (unused) bi representation " << impl.bounded_integer_type_representation << std::endl; @@ -110,7 +110,7 @@ PPL::Concrete_Expression_Type::OK() const { } if (impl.bounded_integer_type_overflow != OVERFLOW_IMPOSSIBLE) { -#ifdef NDEBUG +#ifndef NDEBUG std::cerr << "floating point type has illegal (unused) bi overflow " << impl.bounded_integer_type_overflow << std::endl;

1 0

[GIT] ppl/ppl(floating_point): Implemented Linear_Form::overflows and
by Fabio Bossi 21 Jul '10

21 Jul '10

Module: ppl/ppl Branch: floating_point Commit: c85b888fd3f403f17df903db4016411d2c199196 URL: http://www.cs.unipr.it/git/gitweb.cgi?p=ppl/ppl.git;a=commit;h=c85b888fd3f4… Author: Fabio Bossi <bossi(a)cs.unipr.it> Date: Wed Jul 21 13:35:29 2010 +0200 Implemented Linear_Form::overflows and Linear_Form::relative_error. --- src/Linear_Form.defs.hh | 7 ++++ src/Linear_Form.inlines.hh | 16 +++++++++ src/Linear_Form.templates.hh | 70 ++++++++++++++++++++++++++++++++++++++++++ 3 files changed, 93 insertions(+), 0 deletions(-) diff --git a/src/Linear_Form.defs.hh b/src/Linear_Form.defs.hh index c58fbbc..c80b605 100644 --- a/src/Linear_Form.defs.hh +++ b/src/Linear_Form.defs.hh @@ -322,6 +322,13 @@ public: //! Swaps \p *this with \p y. void swap(Linear_Form& y); + // Floating point analysis related methods. + + bool overflows() const; + + void relative_error(Floating_Point_Format analyzed_format, + Linear_Form& result) const; + private: //! The generic coefficient equal to the singleton zero. static C zero; diff --git a/src/Linear_Form.inlines.hh b/src/Linear_Form.inlines.hh index 60f03de..43ca4e6 100644 --- a/src/Linear_Form.inlines.hh +++ b/src/Linear_Form.inlines.hh @@ -204,6 +204,22 @@ Linear_Form<C>::ascii_load(std::istream& s) { return false; } +// Floating point analysis related methods. +template <typename C> +inline bool +Linear_Form<C>::overflows() const { + if (!inhomogeneous_term().is_bounded()) + return true; + + dimension_type dimension = space_dimension(); + for (dimension_type i = 0; i < dimension; ++i) { + if (!coefficient(Variable(i)).is_bounded()) + return true; + } + + return false; +} + } // namespace Parma_Polyhedra_Library diff --git a/src/Linear_Form.templates.hh b/src/Linear_Form.templates.hh index 9c39fc3..83e7c1f 100644 --- a/src/Linear_Form.templates.hh +++ b/src/Linear_Form.templates.hh @@ -375,6 +375,76 @@ Linear_Form<C>::OK() const { return true; } +// Floating point analysis related methods. +template <typename C> +void +Linear_Form<C>::relative_error( + const Floating_Point_Format analyzed_format, + Linear_Form& result) const { + typedef typename C::boundary_type analyzer_format; + + // Get the necessary information on the analyzed's format. + unsigned int f_base; + unsigned int f_mantissa_bits; + switch (analyzed_format) { + case IEEE754_HALF: + f_base = float_ieee754_half::BASE; + f_mantissa_bits = float_ieee754_half::MANTISSA_BITS; + break; + case IEEE754_SINGLE: + f_base = float_ieee754_single::BASE; + f_mantissa_bits = float_ieee754_single::MANTISSA_BITS; + break; + case IEEE754_DOUBLE: + f_base = float_ieee754_double::BASE; + f_mantissa_bits = float_ieee754_double::MANTISSA_BITS; + break; + case IBM_SINGLE: + f_base = float_ibm_single::BASE; + f_mantissa_bits = float_ibm_single::MANTISSA_BITS; + break; + case IEEE754_QUAD: + f_base = float_ieee754_quad::BASE; + f_mantissa_bits = float_ieee754_quad::MANTISSA_BITS; + break; + case INTEL_DOUBLE_EXTENDED: + f_base = float_intel_double_extended::BASE; + f_mantissa_bits = float_intel_double_extended::MANTISSA_BITS; + break; + } + + C error_propagator; + analyzer_format lb = -pow(f_base, + -static_cast<analyzer_format>(f_mantissa_bits)); + error_propagator.build(i_constraint(GREATER_OR_EQUAL, lb), + i_constraint(LESS_OR_EQUAL, -lb)); + + // Handle the inhomogeneous term. + const C* current_term = &inhomogeneous_term(); + assert(current_term->is_bounded()); + + C current_multiplier(std::max(std::abs(current_term->lower()), + std::abs(current_term->upper()))); + Linear_Form current_result_term(current_multiplier); + current_result_term *= error_propagator; + result = Linear_Form(current_result_term); + + // Handle the other terms. + dimension_type dimension = space_dimension(); + for (dimension_type i = 0; i < dimension; ++i) { + current_term = &coefficient(Variable(i)); + assert(current_term->is_bounded()); + current_multiplier = C(std::max(std::abs(current_term->lower()), + std::abs(current_term->upper()))); + current_result_term = Linear_Form(Variable(i)); + current_result_term *= current_multiplier; + current_result_term *= error_propagator; + result += current_result_term; + } + + return; +} + /*! \relates Parma_Polyhedra_Library::Linear_Form */ template <typename C> std::ostream&

1 0

[GIT] ppl/ppl(floating_point): Do not const-qualify by-copy parameters in declarations.
by Fabio Bossi 21 Jul '10

21 Jul '10

Module: ppl/ppl Branch: floating_point Commit: 152569b0c53bb429e915bfa37b6bf1c7a0ea2fdb URL: http://www.cs.unipr.it/git/gitweb.cgi?p=ppl/ppl.git;a=commit;h=152569b0c53b… Author: Fabio Bossi <bossi(a)cs.unipr.it> Date: Wed Jul 21 12:36:33 2010 +0200 Do not const-qualify by-copy parameters in declarations. --- src/Float.defs.hh | 3 +-- 1 files changed, 1 insertions(+), 2 deletions(-) diff --git a/src/Float.defs.hh b/src/Float.defs.hh index d2f2678..21270da 100644 --- a/src/Float.defs.hh +++ b/src/Float.defs.hh @@ -361,8 +361,7 @@ public: #endif template <typename FP_Interval_Type> -FP_Interval_Type compute_absolute_error(const Floating_Point_Format - analyzed_format); +FP_Interval_Type compute_absolute_error(Floating_Point_Format analyzed_format); template <typename FP_Interval_Type, Floating_Point_Format> class absolute_error;

1 0