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