Module: ppl/ppl Branch: floating_point Commit: 663dfab0dee076d37deb1f0afadb7dbcf4bebd42 URL: http://www.cs.unipr.it/git/gitweb.cgi?p=ppl/ppl.git;a=commit;h=663dfab0dee07...
Author: Fabio Bossi bossi@cs.unipr.it Date: Fri Sep 4 11:16:05 2009 +0200
Added a workaround for the fact that we don't have an explicit way to divide a Linear_Form by an Interval.
---
...Division_Floating_Point_Expression.templates.hh | 14 ++++++++++---- 1 files changed, 10 insertions(+), 4 deletions(-)
diff --git a/src/Division_Floating_Point_Expression.templates.hh b/src/Division_Floating_Point_Expression.templates.hh index 5fcfab9..ca464f0 100644 --- a/src/Division_Floating_Point_Expression.templates.hh +++ b/src/Division_Floating_Point_Expression.templates.hh @@ -43,10 +43,16 @@ typename Division_Floating_Point_Expression<FP_Interval_Type, FP_Format> FP_Interval_Type abs_error(-this->absolute_error); // FIXME: this may be incorrect for some policies. abs_error.join_assign(this->absolute_error); - FP_Linear_Form result = linearized_first_operand / intervalized_second_operand - + relative_error(linearized_first_operand) / intervalized_second_operand - + abs_error; - + /* + FIXME: since we currently lack an explicit way to divide a linear form + by a scalar, we temporarily multiply by 1/scalar. + */ + FP_Interval_Type reversed_intervalized_second_operand = + (FP_Interval_Type(boundary_type(1)) / intervalized_second_operand); + FP_Linear_Form result = linearized_first_operand * + reversed_intervalized_second_operand + + relative_error(linearized_first_operand) * + reversed_intervalized_second_operand + abs_error; return result; }