--- lib/JAMA/tnt_math_utils.h.orig 2006-06-17 01:46:22.000000000 +0300 +++ lib/JAMA/tnt_math_utils.h 2006-06-17 01:47:02.000000000 +0300 @@ -20,11 +20,20 @@ namespace TNT { /** + @returns the absolute value of a real (no-complex) scalar. +*/ +template +Real abs(const Real &a) +{ + return (a > 0 ? a : -a); +} +/** @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 Real hypot(const Real &a, const Real &b) { @@ -56,15 +65,6 @@ } */ -/** - @returns the absolute value of a real (no-complex) scalar. -*/ -template -Real abs(const Real &a) -{ - return (a > 0 ? a : -a); -} - } #endif /* MATH_UTILS_H */