Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/externals/g4tools/include/tools/eqT

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ] 

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