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