1 files changed, 70 insertions, 0 deletions
diff --git a/extra/libofa/tnt_math_utils.h b/extra/libofa/tnt_math_utils.h
new file mode 100644
index 000000000..40de63f05
--- /dev/null
+++ b/extra/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 */