blob: ae31765f37ce0367fdbdf99a34dc3d92fff1b626 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
|
/** \file \ingroup smoke
*/
#ifndef MATH_UTILS_H
#define MATH_UTILS_H
/* needed for fabs, sqrt() below */
#include <cmath>
namespace TNT
{
/**
@returns hypotenuse of real (non-complex) scalars a and b by
avoiding underflow/overflow
using (a * sqrt( 1 + (b/a) * (b/a))), rather than
sqrt(a*a + b*b).
*/
template <class Real>
Real hypot(const Real &a, const Real &b)
{
if (a== 0)
return fabs(b);
else
{
Real c = b/a;
return fabs(a) * sqrt(1 + c*c);
}
}
} /* TNT namespace */
#endif
/* MATH_UTILS_H */
|