"MarlinReco"  1.32.0
tnt_math_utils.h
1 #ifndef MATH_UTILS_H
2 #define MATH_UTILS_H
3 
4 /* needed for fabs, sqrt() below */
5 #include <cmath>
6 
7 
8 
9 namespace TNT
10 {
17 template <class Real>
18 Real hypot(const Real &a, const Real &b)
19 {
20 
21  if (a== 0)
22  return abs(b);
23  else
24  {
25  Real c = b/a;
26  return fabs(a) * sqrt(1 + c*c);
27  }
28 }
29 } /* TNT namespace */
30 
31 
32 
33 #endif
34 /* MATH_UTILS_H */