diff options
| author | tnut <thierryn1 at hispeed dot ch> | 2010-02-26 23:22:14 +0100 |
|---|---|---|
| committer | tnut <thierryn1 at hispeed dot ch> | 2010-02-26 23:22:14 +0100 |
| commit | 1a6a3c4fe8ea50b6084546a4abd40e9f0e609b28 (patch) | |
| tree | 871c615de837a33299aa620b78d92aaab9859b76 /test/libofa/tnt_math_utils.h | |
| parent | a7243fae5a8fa9a649453a49da27145810d5bc47 (diff) | |
| download | nutyx-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.h | 70 |
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 */
|