tnt_math_utils.h

#ifndef MATH_UTILS_H
#define MATH_UTILS_H

/* needed for fabs, sqrt() below */
#include <cmath>



namespace TNT
{
template <class Real>
Real hypot(const Real &a, const Real &b)
{
        
        if (a== 0)
                return abs(b);
        else
        {
                Real c = b/a;
                return fabs(a) * sqrt(1 + c*c);
        }
}
} /* TNT namespace */



#endif
/* MATH_UTILS_H */