1 /* Difference_Floating_Point_Expression class implementation:
2    non-inline template functions.
3    Copyright (C) 2001-2010 Roberto Bagnara <bagnara@cs.unipr.it>
4    Copyright (C) 2010-2016 BUGSENG srl (http://bugseng.com)
5 
6 This file is part of the Parma Polyhedra Library (PPL).
7 
8 The PPL is free software; you can redistribute it and/or modify it
9 under the terms of the GNU General Public License as published by the
10 Free Software Foundation; either version 3 of the License, or (at your
11 option) any later version.
12 
13 The PPL is distributed in the hope that it will be useful, but WITHOUT
14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
15 FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
16 for more details.
17 
18 You should have received a copy of the GNU General Public License
19 along with this program; if not, write to the Free Software Foundation,
20 Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02111-1307, USA.
21 
22 For the most up-to-date information see the Parma Polyhedra Library
23 site: http://bugseng.com/products/ppl/ . */
24 
25 #ifndef PPL_Difference_Floating_Point_Expression_templates_hh
26 #define PPL_Difference_Floating_Point_Expression_templates_hh 1
27 
28 namespace Parma_Polyhedra_Library {
29 
30 template <typename FP_Interval_Type, typename FP_Format>
31 bool Difference_Floating_Point_Expression<FP_Interval_Type, FP_Format>
linearize(const FP_Interval_Abstract_Store & int_store,const FP_Linear_Form_Abstract_Store & lf_store,FP_Linear_Form & result) const32 ::linearize(const FP_Interval_Abstract_Store& int_store,
33             const FP_Linear_Form_Abstract_Store& lf_store,
34             FP_Linear_Form& result) const {
35   if (!first_operand->linearize(int_store, lf_store, result)) {
36     return false;
37   }
38   FP_Linear_Form rel_error;
39   relative_error(result, rel_error);
40   result += rel_error;
41   FP_Linear_Form linearized_second_operand;
42   if (!second_operand->linearize(int_store, lf_store,
43                       linearized_second_operand)) {
44     return false;
45   }
46   result -= linearized_second_operand;
47   relative_error(linearized_second_operand, rel_error);
48   result += rel_error;
49   result += this->absolute_error;
50   return !this->overflows(result);
51 }
52 
53 } // namespace Parma_Polyhedra_Library
54 
55 #endif // !defined(PPL_Difference_Floating_Point_Expression_templates_hh)
56