| 1 |
--- lib/JAMA/tnt_math_utils.h.orig 2006-06-17 01:46:22.000000000 +0300
|
| 2 |
+++ lib/JAMA/tnt_math_utils.h 2006-06-17 01:47:02.000000000 +0300
|
| 3 |
@@ -20,11 +20,20 @@
|
| 4 |
namespace TNT
|
| 5 |
{
|
| 6 |
/**
|
| 7 |
+ @returns the absolute value of a real (no-complex) scalar.
|
| 8 |
+*/
|
| 9 |
+template <class Real>
|
| 10 |
+Real abs(const Real &a)
|
| 11 |
+{
|
| 12 |
+ return (a > 0 ? a : -a);
|
| 13 |
+}
|
| 14 |
+/**
|
| 15 |
@returns hypotenuse of real (non-complex) scalars a and b by
|
| 16 |
avoiding underflow/overflow
|
| 17 |
using (a * sqrt( 1 + (b/a) * (b/a))), rather than
|
| 18 |
sqrt(a*a + b*b).
|
| 19 |
*/
|
| 20 |
+
|
| 21 |
template <class Real>
|
| 22 |
Real hypot(const Real &a, const Real &b)
|
| 23 |
{
|
| 24 |
@@ -56,15 +65,6 @@
|
| 25 |
}
|
| 26 |
*/
|
| 27 |
|
| 28 |
-/**
|
| 29 |
- @returns the absolute value of a real (no-complex) scalar.
|
| 30 |
-*/
|
| 31 |
-template <class Real>
|
| 32 |
-Real abs(const Real &a)
|
| 33 |
-{
|
| 34 |
- return (a > 0 ? a : -a);
|
| 35 |
-}
|
| 36 |
-
|
| 37 |
}
|
| 38 |
#endif
|
| 39 |
/* MATH_UTILS_H */
|