tnt_math_utils.h

#ifndef MATH_UTILS_H
#define MATH_UTILS_H
 
 
/* needed for abs(), sqrt() below */
#include <cmath>
 
 
 
namespace TNT
{
 
using namespace std;
template <class Real>
Real hypot(const Real &a, const Real &b)
{
    
    if (a== 0)
        return abs(b);
    else
    {
        Real c = b/a;
        return abs(a) * sqrt(1 + c*c);
    }
}
} /* TNT namespace */
 
 
 
#endif
/* MATH_UTILS_H */