| Geant4 Cross Reference |
1 // Copyright (C) 2010, Guy Barrand. All rights 1 // Copyright (C) 2010, Guy Barrand. All rights reserved.
2 // See the file tools.license for terms. 2 // See the file tools.license for terms.
3 3
4 #ifndef tools_eqT 4 #ifndef tools_eqT
5 #define tools_eqT 5 #define tools_eqT
6 6
7 namespace tools { 7 namespace tools {
8 8
9 template <class NUMBER,class PREC> 9 template <class NUMBER,class PREC>
10 inline bool numbers_are_equal(const NUMBER& a_ 10 inline bool numbers_are_equal(const NUMBER& a_left,const NUMBER& a_right,const PREC& a_prec,PREC(*a_fabs)(const NUMBER&)) {
11 NUMBER diff = a_left - a_right; 11 NUMBER diff = a_left - a_right;
12 if(a_fabs(diff)>=a_prec) return false; 12 if(a_fabs(diff)>=a_prec) return false;
13 return true; 13 return true;
14 } 14 }
15 15
16 template <class NUMBER,class PREC> 16 template <class NUMBER,class PREC>
17 inline bool is_zero(const NUMBER& a_left,const 17 inline bool is_zero(const NUMBER& a_left,const PREC& a_prec,PREC(*a_fabs)(const NUMBER&)) {
18 if(a_fabs(a_left)>=a_prec) return false; 18 if(a_fabs(a_left)>=a_prec) return false;
19 return true; 19 return true;
20 } 20 }
21 21
22 template <class VEC,class PREC> 22 template <class VEC,class PREC>
23 inline bool vectors_are_equal(const VEC& a_1,c 23 inline bool vectors_are_equal(const VEC& a_1,const VEC& a_2,const PREC& a_prec,PREC(*a_fabs)(const PREC&)) {
24 if(a_1.size()!=a_2.size()) return false; 24 if(a_1.size()!=a_2.size()) return false;
25 typedef typename VEC::size_type sz_t; 25 typedef typename VEC::size_type sz_t;
26 sz_t sz = a_1.size(); 26 sz_t sz = a_1.size();
27 //bool status = true; 27 //bool status = true;
28 for(sz_t index=0;index<sz;index++) { 28 for(sz_t index=0;index<sz;index++) {
29 if(!numbers_are_equal(a_1[index],a_2[index 29 if(!numbers_are_equal(a_1[index],a_2[index],a_prec,a_fabs))
30 //{ ::printf("debug : vectors_are_equals : 30 //{ ::printf("debug : vectors_are_equals : %lu : %g %g\n",index,a_1[index],a_2[index]);
31 return false; 31 return false;
32 //status = false; 32 //status = false;
33 //} 33 //}
34 } 34 }
35 return true; 35 return true;
36 //return status; 36 //return status;
37 } 37 }
38 38
39 template <class VECVEC,class PREC> 39 template <class VECVEC,class PREC>
40 inline bool vecvecs_are_equal(const VECVEC& a_ 40 inline bool vecvecs_are_equal(const VECVEC& a_1,const VECVEC& a_2,const PREC& a_prec,PREC(*a_fabs)(const PREC&)) {
41 if(a_1.size()!=a_2.size()) return false; 41 if(a_1.size()!=a_2.size()) return false;
42 typedef typename VECVEC::size_type sz_t; 42 typedef typename VECVEC::size_type sz_t;
43 sz_t sz = a_1.size(); 43 sz_t sz = a_1.size();
44 for(sz_t index=0;index<sz;index++) { 44 for(sz_t index=0;index<sz;index++) {
45 if(!vectors_are_equal(a_1[index],a_2[index 45 if(!vectors_are_equal(a_1[index],a_2[index],a_prec,a_fabs)) return false;
46 } 46 }
47 return true; 47 return true;
48 } 48 }
49 49
50 ////////////////////////////////////////////// 50 ///////////////////////////////////////////////////////////////////////////////////////////////////////////////
51 /// PREC(*a_fabs)(PREC) : //////////////////// 51 /// PREC(*a_fabs)(PREC) : /////////////////////////////////////////////////////////////////////////////////////
52 ////////////////////////////////////////////// 52 ///////////////////////////////////////////////////////////////////////////////////////////////////////////////
53 53
54 // for histo equals functions. 54 // for histo equals functions.
55 55
56 template <class NUMBER,class PREC> 56 template <class NUMBER,class PREC>
57 inline bool numbers_are_equal(const NUMBER& a_ 57 inline bool numbers_are_equal(const NUMBER& a_left,const NUMBER& a_right,const PREC& a_prec,PREC(*a_fabs)(NUMBER)) {
58 NUMBER diff = a_left - a_right; 58 NUMBER diff = a_left - a_right;
59 if(a_fabs(diff)>=a_prec) return false; 59 if(a_fabs(diff)>=a_prec) return false;
60 return true; 60 return true;
61 } 61 }
62 62
63 template <class VEC,class PREC> 63 template <class VEC,class PREC>
64 inline bool vectors_are_equal(const VEC& a_1,c 64 inline bool vectors_are_equal(const VEC& a_1,const VEC& a_2,const PREC& a_prec,PREC(*a_fabs)(PREC)) {
65 if(a_1.size()!=a_2.size()) return false; 65 if(a_1.size()!=a_2.size()) return false;
66 typedef typename VEC::size_type sz_t; 66 typedef typename VEC::size_type sz_t;
67 sz_t sz = a_1.size(); 67 sz_t sz = a_1.size();
68 //bool status = true; 68 //bool status = true;
69 for(sz_t index=0;index<sz;index++) { 69 for(sz_t index=0;index<sz;index++) {
70 if(!numbers_are_equal(a_1[index],a_2[index 70 if(!numbers_are_equal(a_1[index],a_2[index],a_prec,a_fabs))
71 //{ ::printf("debug : vectors_are_equals : 71 //{ ::printf("debug : vectors_are_equals : %lu : %g %g\n",index,a_1[index],a_2[index]);
72 return false; 72 return false;
73 //status = false; 73 //status = false;
74 //} 74 //}
75 } 75 }
76 return true; 76 return true;
77 //return status; 77 //return status;
78 } 78 }
79 79
80 template <class VECVEC,class PREC> 80 template <class VECVEC,class PREC>
81 inline bool vecvecs_are_equal(const VECVEC& a_ 81 inline bool vecvecs_are_equal(const VECVEC& a_1,const VECVEC& a_2,const PREC& a_prec,PREC(*a_fabs)(PREC)) {
82 if(a_1.size()!=a_2.size()) return false; 82 if(a_1.size()!=a_2.size()) return false;
83 typedef typename VECVEC::size_type sz_t; 83 typedef typename VECVEC::size_type sz_t;
84 sz_t sz = a_1.size(); 84 sz_t sz = a_1.size();
85 for(sz_t index=0;index<sz;index++) { 85 for(sz_t index=0;index<sz;index++) {
86 if(!vectors_are_equal(a_1[index],a_2[index 86 if(!vectors_are_equal(a_1[index],a_2[index],a_prec,a_fabs)) return false;
87 } 87 }
88 return true; 88 return true;
89 } 89 }
90 90
91 ////////////////////////////////////////////// 91 ///////////////////////////////////////////////////////////////////////////////////////////////////////////////
92 ////////////////////////////////////////////// 92 ///////////////////////////////////////////////////////////////////////////////////////////////////////////////
93 ////////////////////////////////////////////// 93 ///////////////////////////////////////////////////////////////////////////////////////////////////////////////
94 94
95 } 95 }
96 96
97 #endif 97 #endif