summaryrefslogtreecommitdiffstats
path: root/test/libofa/tnt_math_utils.h
diff options
context:
space:
mode:
authortnut <thierryn1 at hispeed dot ch>2010-02-26 23:22:14 +0100
committertnut <thierryn1 at hispeed dot ch>2010-02-26 23:22:14 +0100
commit1a6a3c4fe8ea50b6084546a4abd40e9f0e609b28 (patch)
tree871c615de837a33299aa620b78d92aaab9859b76 /test/libofa/tnt_math_utils.h
parenta7243fae5a8fa9a649453a49da27145810d5bc47 (diff)
downloadnutyx-extra-1a6a3c4fe8ea50b6084546a4abd40e9f0e609b28.tar.gz
nutyx-extra-1a6a3c4fe8ea50b6084546a4abd40e9f0e609b28.tar.bz2
nutyx-extra-1a6a3c4fe8ea50b6084546a4abd40e9f0e609b28.tar.xz
nutyx-extra-1a6a3c4fe8ea50b6084546a4abd40e9f0e609b28.zip
libofa dans test
Diffstat (limited to 'test/libofa/tnt_math_utils.h')
-rw-r--r--test/libofa/tnt_math_utils.h70
1 files changed, 70 insertions, 0 deletions
diff --git a/test/libofa/tnt_math_utils.h b/test/libofa/tnt_math_utils.h
new file mode 100644
index 000000000..40de63f05
--- /dev/null
+++ b/test/libofa/tnt_math_utils.h
@@ -0,0 +1,70 @@
+#ifndef MATH_UTILS_H
+#define MATH_UTILS_H
+
+#include <math.h>
+/* needed for sqrt() below */
+
+#ifdef PREANSI
+template <class _Tp>
+inline const _Tp& min(const _Tp& __a, const _Tp& __b) {
+ return __b < __a ? __b : __a;
+}
+
+template <class _Tp>
+inline const _Tp& max(const _Tp& __a, const _Tp& __b) {
+ return __a < __b ? __b : __a;
+}
+#endif
+
+
+namespace TNT
+{
+/**
+ @returns the absolute value of a real (no-complex) scalar.
+*/
+template <class Real>
+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 <class Real>
+Real hypot(const Real &a, const Real &b)
+{
+
+ if (a== 0)
+ return abs(b);
+ else
+ {
+ Real c = b/a;
+ return a * sqrt(1 + c*c);
+ }
+}
+
+/**
+ @returns the minimum of scalars a and b.
+template <class Scalar>
+Scalar min(const Scalar &a, const Scalar &b)
+{
+ return a < b ? a : b;
+}
+*/
+
+/**
+ @returns the maximum of scalars a and b.
+template <class Scalar>
+Scalar max(const Scalar &a, const Scalar &b)
+{
+ return a > b ? a : b;
+}
+*/
+
+}
+#endif
+/* MATH_UTILS_H */